Fast algorithm for streaming wavefront
View Patent ↗The invention is generally directed to the field of image processing, and more particularly to a method and an apparatus for determining a wavefront of an object, in particular a human eye. The invention discloses a method and an apparatus for real-time wavefront sensing of an optical system utilizing two different algorithms for detecting centroids of a centroid image as provided by a Hartmann-Shack wavefront sensor. A first algorithm detects an initial position of all centroids and a second algorithm detects incremental changes of all centroids detected by said first algorithm.
1. A method for real-time wavefront sensing of an optical system comprising the following steps: (a) acquiring a first centroid image and determining the positions of the centroids in said first centroid image using a first centroid determining algorithm; (b) acquiring a next centroid image and determining the positions of the centroids in the next centroid image using a second centroid determining algorithm, wherein said second centroid determining algorithm determines the position of each centroid with reference to the position of the respective centroid in the first centroid image, wherein the second centroid determining algorithm is an algorithm faster than the first centroid determining algorithm.
2. The method of claim 1 , wherein the second centroid determining algorithm determines the position of a centroid in the next centroid image by searching in an area adjacent to the position of the respective centroid in the first centroid image.
3. The method of claim 2 , wherein said area adjacent to the position of a centroid in the first centroid image is defined as the direct neighbourhood being determined by a deviation dx and dy.
4. The method of claim 3 , wherein the deviation dx, dy determining the neighbourhood of a position of a centroid in the first centroid image amounts to 0.1% to 10% of the distance X of two neighbouring centroids in the x-direction and the distance Y of two neighbouring centroids in the y-direction, respectively.
5. The method of any of claims 1 to 4 , further comprising the steps:
(c) acquiring a further centroid image and determining the positions of the centroids in the further centroid image using the second centroid determining algorithm, wherein the position of a centroid in the further centroid image is determined by reference to the position of the respective centroid in a previously acquired centroid image.
6. The method of claim 5 , further comprising the step of comparing the positions of the centroids in a succeeding centroid image with the positions of the respective centroids in the first centroid image and if the deviation of the positions of the respective centroids is larger than a pre-given tolerance or if some centroids are lacking using the next succeeding centroid image.
7. The method of claim 6 , further comprising the step of determining whether within a pre-given time any further succeeding centroid image deviates more than the pre-given tolerance, and evaluating the next centroid image with the first centroid determining algorithm if all further succeeding centroid images deviate more than the pre-given tolerance.
8. The method of claim 7 , further comprising the step of determining changes of the wavefront ΔW(x, y) on the basis of the changes of the centroid positions for each following centroid image.
9. The method of claim 8 , wherein said first centroid determining algorithm operates on a 1 Hz-2 Hz base.
10. The method of claim 9 , wherein said second centroid determining algorithm operates on a 10 Hz-15 Hz base.
11. The method of claim 6 , further comprising the step of determining changes of the wavefront ΔW(x, y) on the basis of the changes of the centroid positions for each following centroid image.
12. The method of claim 11 , wherein said first centroid determining algorithm operates on a 1 Hz-2 Hz base.
13. The method of claim 12 , wherein said second centroid determining algorithm operates on a 10 Hz-15 Hz base.
14. The method of claim 5 , further comprising the step of determining changes of the wavefront ΔW(x, y) on the basis of the changes of the centroid positions for each following centroid image.
15. The method of claim 14 , wherein said first centroid determining algorithm operates on a 1 Hz-2 Hz base.
16. The method of claim 15 , wherein said second centroid determining algorithm operates on a 10 Hz-15 Hz base.
17. Apparatus for real-time wavefront sensing of an optical system comprising:
(a) acquiring means for acquiring centroid images;
(b) first determining means for determining the positions of the centroids in a first centroid image using a first centroid determining algorithm;
(c) second determining means for determining the positions of the centroids in a next centroid image using a second centroid determining algorithm, wherein said second centroid determining algorithm determines the position of each centroid with reference to the position of the respective centroid in the first centroid image, wherein the second centroid determining algorithm is an algorithm faster than the first centroid determining algorithm.
18. The apparatus of claim 17 , wherein the second determining means is adapted to determine the position of a centroid in a further centroid image by reference to the position of the respective centroid in a previously acquired centroid image.
19. The apparatus of claim 17 or 18 , wherein said first determining means operates on a 1 Hz-2 Hz base.
20. The apparatus of claim 19 , wherein said second determining means operates on a 10 Hz-15 Hz base.
21. The apparatus of claim 18 , wherein said second determining means operates on a 10 Hz-15 Hz base.
22. The apparatus of claim 17 wherein said second determining means operates on a 10 Hz-15 Hz base.