Code for measurement of pupil radius and ellipticity