Determination of two-dimensional characteristics of the anterior surface of the eye is becoming increasingly important in modern optometry and ophthalmology practice. In particular, accurate estimation of the pupil size and centration is crucial in customized refractive surgery, corneal transplantation, and advanced contact lens fitting. The pupil parameters change under different lighting conditions so they often need to be related to some fixed reference such as the limbus outline. However, current commercial pupillometers do not estimate limbus position. We present a novel algorithm for automatic extraction of pupil parameters from digital images that takes the relative limbus information into account. The algorithm utilizes several customized image processing techniques that form a robust procedure which performs well for a wide range of clinical images. We apply the developed algorithm to images obtained by a standard digital camera, and specialized ophthalmic instruments such as a wavefront sensor and a high-speed imaging system.