Method for estimating refraction of an eye of an individual, the method comprising the following steps: a) Providing a set of N acquired eccentric photorefraction images of the eye, where N is an integer higher than or equal to one; b) Initializing values for a set of parameters including at least low order aberration such as sphere; c) Providing a set of N simulated eccentric photorefraction images using a simulation model according to the values of the set of parameters; d) Determining an estimator of a difference between the set of N acquired eccentric photorefraction images and the set of N simulated eccentric photorefraction images; e) Performing an optimization using an optimization algorithm so as to optimize said estimator by adjusting values for the set of parameters and iterating steps c) and d), and f) Deducing an estimation of at least one refraction parameter of the eye from the adjusted values for the set of parameters, wherein step a) includes the following steps g) to i), and wherein each step c) includes the following step j): g) placing a system comprising an image capturing device having an aperture and a plurality of M light sources in front of the eye, where M is an integer higher than or equal to two, the plurality of M light sources being arranged eccentrically around the aperture of the image capturing device at determined positions along at least two directions transverse to an optical axis of the image capturing device, each light source of the M light sources being adapted and configured to illuminate the eye with a light pulse; h) illuminating the eye using each light source of the plurality of M light sources; i) recording the set of N acquired eccentric photorefraction images of the eye using the image capturing device, where N is lower than or equal to M, each image of the set of N acquired eccentric photorefraction images representing the eye illuminated by the light pulse of at least one light source of the plurality of M light sources; j) generating the set of N simulated eccentric photorefraction images using the simulation model according to the values of the set of parameters, wherein the steps b), c), d), e), f), j) are carried out by a processor. a) Providing a set of N acquired eccentric photorefraction images of the eye, where N is an integer higher than or equal to one; b) Initializing values for a set of parameters including at least low order aberration such as sphere; c) Providing a set of N simulated eccentric photorefraction images using a simulation model according to the values of the set of parameters; d) Determining an estimator of a difference between the set of N acquired eccentric photorefraction images and the set of N simulated eccentric photorefraction images; e) Performing an optimization using an optimization algorithm so as to optimize said estimator by adjusting values for the set of parameters and iterating steps c) and d), and f) Deducing an estimation of at least one refraction parameter of the eye from the adjusted values for the set of parameters, g) placing a system comprising an image capturing device having an aperture and a plurality of M light sources in front of the eye, where M is an integer higher than or equal to two, the plurality of M light sources being arranged eccentrically around the aperture of the image capturing device at determined positions along at least two directions transverse to an optical axis of the image capturing device, each light source of the M light sources being adapted and configured to illuminate the eye with a light pulse; h) illuminating the eye using each light source of the plurality of M light sources; i) recording the set of N acquired eccentric photorefraction images of the eye using the image capturing device, where N is lower than or equal to M, each image of the set of N acquired eccentric photorefraction images representing the eye illuminated by the light pulse of at least one light source of the plurality of M light sources; j) generating the set of N simulated eccentric photorefraction images using the simulation model according to the values of the set of parameters, wherein the steps b), c), d), e), f), j) are carried out by a processor.