Traceable, highly accurate calibration of wavefront sensors such as Shack–Hartmann sensors is a challenging task and an active field of research. We have developed a measurement system for the traceable calibration of wavefront sensors employing spherical wavefronts and a point light source in combination with a three-axis linear stage. We obtain an absolute calibration of sensor errors, including the reference spot position and wavefront gradient-dependent errors for each microlens individually. We discuss the main error influences and present an initial measurement uncertainty budget for the calibration. The calibration can be performed with an expanded measurement uncertainty (95% coverage interval) of better than 5 μrad for the wavefront gradient deviation. |
1.IntroductionShack–Hartmann wavefront sensors (SHSs)1 are widely used in different fields that require accurate wavefront sensing. These fields include ophthalmology,2 adaptive optics for astronomy,3 laser beam characterization,4 and quality control in optics fabrication. To obtain a highly accurate measurement, the wavefront sensors must be calibrated. Sensor errors occur due to residual misalignment between the microlens array (MLA) and the image sensor5 or imperfections such as aberrations of the MLA. The calibration should include the reference spot position for each microlens, which creates a constant offset error for each measurement and the wavefront-dependent errors of the sensor. A common approach for calibration uses a collimator to generate a plane reference wavefront. This approach relies heavily on the quality of the employed optics, thus limiting the calibration accuracy. Another approach to obtaining a high-quality spherical reference wavefront is to use a single-mode fiber with a known distance to the sensor. The main problem, in this case, is how to determine the exact absolute distance between the sensor and the fiber tip: an erroneous value of this distance would result in a power error due to the spherical reference wavefront employed. This can be overcome by performing several measurements at different distances between the fiber tip and the sensor and by measuring only the relative movements.6–8 Additional sensor parameters such as the distance between the MLA and the image sensor can be obtained using this approach. However, microlens-specific errors that are dependent on the incident wavefront slope cannot be characterized. The setup proposed in this paper can be used to measure all these errors with high accuracy. This is possible by additionally shifting the light source perpendicular to the optical axis and by consecutively recording wavefront gradients in the respective measurement positions. 1.1.Spot Displacement Error in Shack–Hartman Wavefront SensorsThe SHS measures the incoming wavefront gradient using an MLA. For small angles, the wavefront gradient in the and directions is well approximated by the local wavefront tilt angles and . The indices employed are allocated by separating the wavefront under test into patches of index , and assuming an MLA with the lens index , . For ease of notation, the direction-indicating index of the angle is employed only when an explicit notation of the considered direction is required, since the underlying theory is completely analog for the and directions. The same is true for the spot displacement . An image sensor is placed in the focal plane of the MLA. When the incident wavefront passes through the MLA, each lens produces a spot on the image sensor. The wavefront-dependent focal spot displacement relative to its reference position and its associated microlens and wavefront-dependent deviation indicate the wavefront angle over the size of the respective microlens. Figure 1 shows the basic principle of the SHS. The determination of the local wavefront angle as a function of the spot displacement is given in Eq. (1) The shape of the whole wavefront can be reconstructed from the measured spot displacements using numerical reconstruction algorithms. Commonly these algorithms retrieve the wavefront by least-square fitting.9 For an ideal SHS with no displacement errors and no microlens aberrations, the reference positions of the focal spots would be on a predefined regular grid. Different errors that can occur during sensor manufacturing lead to a systematic deviation. A nonideal substrate of the microlenses or image sensors can result in uneven spacing between pixels or lenses. Furthermore, a curvature of the substrates or a thickness variation will cause the distance between the image sensor and the lens array to vary, where indicates the distance variation at the respective microlens. Additional errors can occur during the alignment between the image sensor and the MLA;5 these errors also contribute to and . The image recording and signal evaluation of the focal spot using the image sensor introduce additional errors. These deviations are dependent on the local wavefront gradient and the microlens parameters. Quantization effects caused by the image sensor10 or variations between pixels in terms of their gain can substantially influence the accuracy of the spot localization and cause a microlens-specific spot displacement deviation. All these errors are represented by in Eq. (1). 2.Setup of the Measurement SystemThe setup uses a fiber-coupled light source. The fiber tip is mounted on a three-axis linear stage and the wavefront sensor under test. The light source is an LED in combination with a multimode fiber with a core diameter of . Alternatively, a laser coupled into a single-mode fiber can be used. Figure 2 shows a schematic drawing of the measurement setup. 2.1.Measurement ProcedureFor the calibration, the fiber tip is moved stepwise in a plane perpendicular to the optical axis at a fixed radial distance . The step width equals the microlens pitch . At each step, each spot displacement on the image sensor is measured and recorded. Since the step width of the fiber displacement equals the microlens pitch , the same sections of the wavefront are measured with different microlenses. The curvature (and thus, the local gradient and angle of the incoming spherical wavefront) depends on the distance of the light source to the sensor . Therefore, if the fiber end is shifted at a fixed distance parallel to the orientation of the MLA by increments of the microlens pitch , the wavefront sections associated with the starting position will be consecutively shifted to each individual microlens. When the measurement is complete, each lens has measured each wavefront section considered. This results in a linear equation system of equations with , for each of the wavefront patches and , and for each of the microlenses in accordance with Eq. (1). The double index indicates the measurement of the ’th wavefront patch using the ’th microlens. Because the same ’th wavefront section is measured times by the respective microlenses, an overdetermined subsystem of equations for a fixed index set results.This procedure is repeated with a different distance . The distance between and equals the movement of the traceably calibrated axis and is referred to as . 2.2.Calibration of the Image Sensor Distance VariationsTo determine the unknown distance variations between the microlenses and the image sensor, measurements of the wavefront using the same ’th microlens at two different radii and are employed. Thus, each microlens measures two spherical wavefronts with different unknown radii but with a known axial displacement between them. A spherical wavefront of radius produces a linear slope on the measured spot displacement6 proportional to the position of the center of the respective microlens (Fig. 3). The local wavefront angle for each lens position is given by Eq. (3) and the according spot displacement by Eq. (4). The index counts the patch difference in the direction from the central patch From Eqs. (3) and (4), it follows as The value in Eq. (5) represents the measured spot displacement for the ’th microlens that records the ’th wavefront patch. For an ideal system, the measurement value would equal the linear increase when scanning the wavefront patches in the or direction with the ’th microlens. However, in the real system, the spot displacement is distorted. To reduce the influence of these deviations, the information of all wavefront patches measured with the same microlens is employed to yield an estimate of the slope . The spot displacement in direction is measured for the two different radii and . The measurement of the undisturbed spot displacement is affected by the microlens-specific error ; this error consists of a systematic error due to the shift of the reference spot position , which is calibrated in Sec. 2.3, and stochastic errors.As outlined in Sec. 2.1, the wavefront sections are all measured with the same ’th microlens, which introduces a constant offset error . Thus, the linear slope of the spot displacement associated with a spherical reference wave impinging on an ideal wavefront sensor illustrated in Fig. 3 may be approximated by linear regression of the spot displacement measurements. Since is represented by the slope of the fitted line, the constant offset does not influence the result and the stochastic errors are reduced via the averaging effect of the approach. The estimate of the slope in the ’th microlens employing the complete set of the measured ’th wavefront patches results from the following linear regression: In Eqs. (6)–(9), denotes the mean value of the measured spot displacement and refers to the estimated slope of the spot displacement.For the two measurements at radius and , we obtain The distance between and can be written as , which is measured using the glass scales of the axis. Thus, Eq. (11) can be written as Equations (10) and (12) can be solved for the radial distance . This distance is the same for all microlenses, thus enabling averaging over all microlenses to further reduce measurement noise. This value for the radial distance can be used to determine the distance from Eq. (10) The step size in Eq. (14) and the distance change in Eq. (13) are applied via the calibrated linear axis and characterized using the uncertainty of the axis. By considering Fig. 1, we can immediately conclude that represents a measure of the sensitivity of the spot displacement to changes in the wavefront gradient.In this step, we assume the pixel pitch from the image sensor to be ideal. An erroneous pixel pitch that is constant in the image sensor region considered for the ’th microlens would cause a linear error in . However, it would also affect the measurement of the focal spot in the same way. From Eq. (1), we can observe that these elements cancel each other out and do not affect the measurement of the wavefront angle. It follows that the determined value also corrects an erroneous pixel pitch, although in this case, it does not represent the real distance between the microlens and the image sensor. 2.3.Reference Spot Position CalibrationAs described in Sec. 2.1, of each microlens is measured with the radial distance and the associated local wavefront angle of each wavefront patch. The systematic deviation is identified in advance using the procedure introduced in Sec. 2.1 and included as a compensation value for the values in this section. Since the undisturbed local wavefront angle and the associated spot displacement depend on the ’th wavefront patch considered, they do not change if measured with different microlenses. To determine the constant systematic, microlens-specific spot displacement error , the same ’th wavefront section is measured with each ’th microlens. The term represents the spot displacements measured in the experiment. Using pairwise comparison of the spot displacements measured for the same wavefront section and by averaging the differences, the microlens-specific spot displacement coefficient is determined by employing the following equation: The resulting calibration parameter describes the systematic spot displacement error of the respective single microlens. After application of for the correction of the spot displacement of each of the ’th microlenses, all microlenses show the same unknown offset deviation from the absolute reference spot position. This unknown offset results in a constant tilt error in the wavefront reconstruction. However, the wavefront form is not affected. By referencing the microlenses against each other, the calibration of the reference spot position is independent of the quality of the incident wavefront; it is not necessary to measure its form (and thus, the absolute position of the fiber tip on the optical axis defining the radial distance ) with high accuracy. In general, for the calibration, the wavefront emitted by the fiber tip must be stable and the relative movement of the axis must be measured accurately. Furthermore, it is assumed that the microlens has been manufactured with an ideal pitch . An error in the pitch would result in different parts of the wavefront being compared and would therefore create an error in that would depend on the incident wavefront shape. Significant errors in the microlens pitch would become apparent through differences in the respective calibration maps when measurements with different wavefront curvatures are conducted. The result of Eq. (15) for the MLA is confirmed when the calculation with the second set of measurements is repeated at the position. The above steps may be taken to preclude the usage of MLAs with an erroneous microlens pitch. A feasible result for may also be obtained employing only a subset of the microlenses for comparison. This setup can also be used to retrieve the wavefront angles of the incident wavefront. As shown in Eq. (1), the measured spot displacement is dependent on the two variables and . The constant systematic spot displacement error and the distance variation are determined above. To obtain the wavefront gradient, we average all measurements associated with the respective ’th wavefront section, taking the correction values and into account By averaging over all microlenses, stochastic errors in can be substantially reduced.As with the calibration of the reference spot position, only the relative tilt and form of the wavefront can be reconstructed. The light source for this setup can be easily replaced. This enables us to analyze the effects of different wavelengths and the degree of coherence of the light source on the measurement. We have published some of these effects in Ref. 11 using a different calibration setup. The traceability chain is shown in Fig. 4. The method relies on the axes being correctly positioned to ensure that the microlenses are illuminated with the identical wavefront section. A displacement-measuring interferometer is used to validate the linear scales of the axes. The tilt errors of the axes are obtained via measurements with an autocollimator. The wavefront measurement is self-referencing because the wavefront is measured during the calibration as outlined in Eq. (15) and only the differences between the lenses affect the results. The laboratory temperature influences the measurement in different ways. The sensor characteristics ( and ) are particularly temperature dependent due to thermal expansion and must therefore be measured as well. 3.Measurement ResultsIn this section, we present the measurement results of a Shack–Hartmann sensor obtained with the setup described above. As a light source, we used a green LED (). The distance between the sensor and the light source is varied between and . Figure 5 shows the effect of the systematic spot displacement error on the local wavefront angle measurement as given in Eq. (1), , for each microlens. Furthermore, we obtained the image sensor distance variation . This is shown in Fig. 6 for each microlens. The microlens-to-image-sensor distance of this sensor is on a low scale and ranges between 17.7 and 17.74 mm with a standard deviation of . Our setup measures the same part of a wavefront with different lenses. After the correction of the systematic sensor errors using the results above, the difference between the measurements of the same wavefront part should be minimized. This can be used to verify our calibration system. These results are shown in Fig. 7. The color bar gives the rms deviation between the lenses for the spot displacement measurement. The mean value of this error is about . Without the correction of the reference spot position, the mean rms value is at , which shows the capability of our calibration. The accuracy is limited by the sensor noise of the image sensor. Additional errors are generated by our calibration setup, including positioning errors of the axes and errors due to the stability of the laboratory environment. Another error is due to the pixelization effect of the image sensor, which influences the spot measurement; this creates a dependency of the wavefront measured on the position of the focal spot on a pixel (see Ref. 10). These pixelization effects become apparent by comparing the wavefront measured through all microlenses with the measurement of a single microlens. In this way, we obtain detailed information on the tilt error as a function of the incoming wavefront gradient. Figure 8 shows the corresponding data for a single microlens. To enhance the visibility of the systematic parts, the noise is reduced by smoothing the data with a low-pass filter. The and directions denote the incoming wavefront gradient and the color bar shows the measurement error for this microlens. The error has a periodic structure equal to the periodicity of the image sensor pixels. A shift of one pixel on the sensor equals a wavefront gradient of 0.72 mrad. This effect is caused by pixel effects of the sensor. This systematic error is about peak-to-valley (PV) and therefore on a low scale. Additionally, we obtain the gradient of the incident wavefront. Since we used a point-light source for illumination, the wavefront has a good spherical shape, which results in a linear slope in the map of the wavefront gradient. Figures 9(a) and 9(b) show the measured wavefront gradient. In Figs. 9(c) and 9(d), a best-fit plane is subtracted to show the nonspherical part of the wavefront. The remaining differences from an ideal sphere have a low PV value in the range of and are again caused by the pixelization effects of the sensor. 4.Uncertainty AnalysisThe measurement uncertainty is influenced by various properties of the system and the specific sensor. Overall, we determine the wavefront angle error with an expanded uncertainty (). For the uncertainty of the calibration, the main input quantity is the exact determination of the position of the light source. The positioning accuracy depends on the axis used. We estimated the positioning accuracy to be . This general positioning error consists of the positioning errors of the single axis and errors due to a constant tilt between the axes and due to movement-dependent tilt effects. At a distance of about 1.7 m between the sensor and the light source, this general positioning error contributes a standard uncertainty of to . Further input quantities include the adjustment of the setup, the temperature stability of the setup and random errors due to noise on the image sensor vibrations and air turbulence. The repeatability of the calibration is determined by comparing several measurements of the same wavefront sensor, which leads to a standard deviation of . This is employed as an estimate of the standard uncertainty. The repeatability includes the random errors caused by electric, thermal, pixel-related, and mechanical noise. These random deviations are estimated from the results shown in Fig. 7 and yield a standard deviation for a single measurement in the interval of for the respective microlenses. Their impact is significantly reduced via averaging over the 2601 wavefront patches, which results in a standard deviation of the mean value below . The alignment between the three-movement axes and the sensor axes also contributes to the measurement uncertainty. The standard uncertainty of the adjustment errors is estimated to be . Temperature effects are analyzed using a cooled wavefront sensor with varying temperatures. The temperature variations during the measurement are 0.2 K; the resulting contribution to the uncertainty is below determined by comparison measurements at different temperatures. Table 1 lists these input quantities together with the corresponding uncertainties. Errors due to the tilt of the fiber tip during measurement of the wavefront are neglectable due to the low deviation of the wavefront from an ideal sphere and the accuracy of the axes. Table 1Measurement uncertainty budget.
This uncertainty can be further improved by improving the measurement of the light source position or by increasing the distance between the sensor and the light source. However, this is limited by the increased sensor noise caused by the reduced light intensity. 5.Conclusion and DiscussionIn this paper, we present a setup for a highly accurate, traceable characterization of SHS using spherical wavefronts and linear stages. This calibration setup allows measurements of the microlens-specific reference spot position deviation , the microlens-specific distance variation to the image sensor , and further wavefront-dependent effects . An SHS was calibrated using this setup with an expanded uncertainty of () for the angular reference spot position deviation of each microlens. The main input quantities for this uncertainty budget are taken from the positioning of the light source. The setup could potentially be further enhanced by better measurement of the light source position. The setup presented could easily be expanded to include other wavefront sensors such as shearing-based sensors. AcknowledgmentsThe funding of the project by the TransMeT program of the German Federal Ministry for Economic Affairs is gratefully acknowledged. ReferencesB. Platt et al.,
“History and principles of Shack–Hartmann wavefront sensing,”
J. Refract. Surg., 17
(5), S573
–S577
(2001). https://doi.org/10.3928/1081-597X-20010901-13 JRSUEY 0883-0444 Google Scholar
J. Liang et al.,
“Objective measurement of the wave aberrations of the human eye with the use of a Hartmann–Shack wave-front sensor,”
J. Opt. Soc. Am. A, 11 1949
–1957
(1994). https://doi.org/10.1364/JOSAA.11.001949 JOAOD6 0740-3232 Google Scholar
A. Glindemann et al.,
“Adaptive optics on large telescopes,”
Exp. Astron., 10 5
–47
(2000). https://doi.org/10.1023/A:1008116831367 EXASER 0922-6435 Google Scholar
B. Schäfer, M. Lübbecke and K. Mann,
“Hartmann–Shack wave front measurements for real time determination of laser beam propagation parameters,”
Rev. Sci. Instrum., 77 053103
(2006). https://doi.org/10.1063/1.2198795 RSINAK 0034-6748 Google Scholar
J. Pfund et al.,
“Misalignment effects of the Shack–Hartmann sensor,”
Appl. Opt., 37 22
–27
(1998). https://doi.org/10.1364/AO.37.000022 APOPAI 0003-6935 Google Scholar
A. Chenishov et al.,
“Calibration of a Shack–Hartmann sensor for absolute measurements of wavefronts,”
Appl. Opt., 44 6419
–6425
(2005). https://doi.org/10.1364/AO.44.006419 APOPAI 0003-6935 Google Scholar
A. Nikitin et al.,
“Absolute calibration of a Shack–Hartmann wavefront sensor for measurements of wavefronts,”
Proc. SPIE, 10925 109250K
(2019). https://doi.org/10.1117/12.2510047 PSISDG 0277-786X Google Scholar
J. Yang et al.,
“Absolute calibration of Hartmann–Shack wavefront sensor by spherical wavefronts,”
Opt. Commun., 283
(6), 910
–916
(2010). https://doi.org/10.1016/j.optcom.2009.11.022 OPCOB8 0030-4018 Google Scholar
W. Southwell,
“Wave-front estimation from wave-front slope measurements,”
J. Opt. Soc. Am., 70
(8), 998
–1006
(1980). https://doi.org/10.1364/JOSA.70.000998 JOSAAH 0030-3941 Google Scholar
B. Neal et al.,
“Shack–Hartmann wavefront sensor precision and accuracy,”
Proc. SPIE, 4779 148
–160
(2002). https://doi.org/10.1117/12.450850 PSISDG 0277-786X Google Scholar
J. Bautsch et al.,
“DGaO Jahrestagung,”
(2018) https://www.dgao-proceedings.de/download/119/119_p34.pdf Google Scholar
|