An adjustable, high sensitivity, wide dynamic range two channel wave-front sensor based on moiré deflectometry has been constructed for measuring distortions of light wave-front transmitted through the atmosphere. In this approach, a slightly divergent laser beam is passed through the turbulent ground level atmosphere and then a beam-splitter divides it into two beams. The beams pass through a pair of moiré deflectometers which are installed parallel and close together. From deviations in the moiré fringes we calculate the two orthogonal components of angle of arrival at each location across the wave-front. The deviations have been deduced in successive frames which allows evolution of the wave-front shape to be determined. The dynamic range and sensitivity of detection are adjustable by merely changing the separation of the gratings and the angle between the rulings of the gratings in both of channels. The spatial resolution of the method is also adjustable by means of bright, dark, and virtual traces for given moiré fringes without paying a toll in the measurement precision.
© 2010 Optical Society of America
Various wave-front sensing techniques have been developed for use in a variety of applications ranging from measuring the wave aberrations of human eyes  to adaptive optics in astronomy . The most commonly used wave-front sensors are the Shack-Hartmann [3,4], curvature sensing , shearing interferometry , phase retrieval methods  and Pyramid wave-front sensor . The Shack-Hartmann sensor is also the most commonly used technique for measurement of turbulence-induced phase distortions for various applications in atmospheric studies and adaptive optics. But, the dynamic range of the Shack-Hartmann sensor is limited by the optical parameters of its microlenses, namely, the spacing and the focal length of the microlens array.
Although not so well known, Talbot effect and moiré deflectometry have been employed for wave-front sensing in different fields and in several different ways [9–13]. For example, a technique for aspheric surface metrology by means of dual-path moiré deflectometry has been suggested by . The method used by them is based on measurements of the distribution of deflection angles for a line near the center of the moiré deflectometer. A limitation of the set-up is that the method is applicable only for stationary cases. To the best of our knowledge, moiré deflectometry has not been used for sensing time-varying distortions of wave-fronts, although by means of moiré technique we have estimated atmospheric turbulence parameters in various schemes . Based on the above experience, we recently suggested a two channel wave-front sensor based on moiré deflectometry . However, it was recognized that the implementation of the method is not very straightforward. In this paper we present a new, reliable, and quite simple method for wave-front sensing. Also, using an algorithm developed for moiré fringe processing, the spatial resolution of the method is improved. In this paper, the theory of the approach by means of which angle of arrival (AA) fluctuations across the wave-front can be determined is presented.
A schematic diagram of the experimental setup is shown in Fig. 1. A slightly divergent laser beam passes through a turbulent ground level atmosphere and enters the wave-front sensor. Then a beam splitter splits it into two beams, a mirror reflects the second beam into a direction parallel to the first beam’s propagation direction, and the two beams pass through a pair of moiré deflectometers. The moiré deflectometers are installed quite close to each other. Directions of the gratings’ rulings are almost parallel in each moiré deflectometer but are perpendicular in the two beams. Moiré patterns are formed at a plane where the second gratings of the moiré deflectometers and a diffuser are installed. The moiré patterns from both the beams are projected on a CCD camera. Successive moiré patterns are recorded by the CCD camera and transferred to a computer, to allow temporal fluctuations of the wave-front phase to be measured highly accurately. Displacements of the moiré fringes in the recorded patterns correspond to the fluctuations of two orthogonal components of the AA across the wave-front. The fluctuations have been deduced in successive frames, and the evolution of the wave-front shape is determined.
2. Theoretical framework
In a usual configuration of moiré deflectometer consisting of two Ronchi gratings, the set up is sensitive to the incoming wave-front changes only in one direction normal to its rulings. In this work, we have used two moiré deflectometers in two channels. We choose a coordinate system so that in the first channel the grating rulings are in the x-direction and those in the second channel are in the y-direction, as shown in Fig. 1. Theoretical considerations are similar for both channels. We present the theory for determining the AA fluctuations for a single channel first. Later the information from both channels are combined.
We assume that the wave-front illuminates two gratings G1 and G2 of equal periods d, separated by Zk along the optical axis, where Zk is the kth Talbot distance of the first grating. Due to wave-front distortions, the self-image of the first grating fluctuates on the second grating. The changes in the moiré patterns are related to grating separation in the experimental setup Zk, the relative rotation between the grating rulings in the moiré deflectometer, and the y component of the AA fluctuations of the incident wave-front αy(x, y). According to moiré formulation, for displacement δy(x,y) of the self-image in a direction normal to the grating rulings, moiré fringes are shifted by δxm(x, y) given by
According to Eq. (2), by increasing the distance between the gratings, increasing the moiré fringe spacing, or decreasing the period of the gratings, the measurement precision of AA fluctuations can be improved. In this work, we have used , and Zk=22.4 cm, and by using , the minimum measurable AA fluctuations, is determined to be 5.9 × 10−6 rad or 1.2 arc sec. The minimum measurable value of is influenced by the number of CCD pixels which were covered by the moiré fringe spacing. On the other hand, if we stipulate that the maximum measurable displacement of a trace is equal to half the moiré fringe spacing, one can measure AA fluctuations in the range −30 to +30 arc sec without having any adjustment in the set-up.
Now, we can generalize the formulation for the second channel. The component αx(x, y) of the AA fluctuation in the direction perpendicular to the rulings of gratings G3 and G4 (parallel to the corresponding moiré traces) on the G3 grating is given by:
Since αx(x, y) and αy(x, y) are equal to the incident wave-front gradients in x and y-directions, respectively, the incident wave-front gradients at a point (x,y) are determined by
It should be mentioned that, the presented formulation is valid for a collimated incoming laser beam. In principal, for a divergent incoming laser beam some modifications are needed . In this work, we have used L=360 m, , and Zk=22.4 cm, and by using Eqs. (2) and (4) of , the moiré fringe spacing, is determined to be 10.6 cm when the relative rotation between the grating rulings, θ, is equal to zero. This is several times larger than the size of the gratings that we have used. As a result the beam divergence in the incoming laser beam is not important and the above simple formulation is applicable.
3. Moiré fringes processing
In the previous work , we used only the displacements of the bright moiré fringes for evaluation of the wave-front shape, as it is the most commonly employed approach. According to Eq. (4), by increasing the moiré fringe spacing the precision of AA fluctuations measurements can be improved. But, then the number of moiré fringes in the field of view as well as the spatial resolution of the method is decreased. In this paper we have used an improved algorithm for processing the moiré fringes to overcome this limitation. Here we have introduced a new concept that we call virtual traces in the moiré patterns. In this approach, the traces of bright moiré fringes, dark moiré fringes, and of points with intensities equal to the mean intensity of the adjacent bright and dark traces (first order virtual traces) were determined. One can potentially produce a large number of virtual traces between two adjacent bright and dark traces by using their intensities and locations. According to Eqs. (4) and (5) the displacements of the traces in the recorded patterns with respect to their reference positions correspond to the fluctuations of AA components which in turn represent the local wave-front gradients. With our approach proposed here, the precision of AA fluctuations can be improved by increasing the moiré fringe spacing (by merely changing the separation of the gratings and/or the angle between the rulings of the gratings in both of moiré deflectometers), while at the same time the desired spatial resolution can be achieved by using a number of virtual traces. Finally one can use the concept of virtual traces in all of the moiré based methods for atmospheric turbulence measurements  to improve their spatial resolution.
4. Wave-front reconstruction
Another important part of a wave-front sensor is the software to convert the 2-D wave-front gradient data into 2-D wave-front phase data. In order to perform the wave-front reconstruction from the measured moiré patterns, we consider the displacements of the bright, dark, and the first order virtual traces with respect to their reference positions, which represents an estimate of the local x-gradients or y-gradients of the wave-front phase. The reference positions of the traces are determined from a long-exposure frame. In practice, by considering two sets of vertical and horizontal moiré traces of a frame in a x–y coordinate system (the vertical and horizontal traces are overlapped in the x–y coordinate system), the intersection points of the vertical and horizontal bright, dark, and the first order virtual traces are determined. x-gradients and y-gradients of the wave-front are deduced from the displacements of the intersection points in successive frames. From this 2-D gradient field we performed an estimate of the wave-front.
5. Experimental results
In the experiment, a He-Ne laser beam propagates through a turbulent surface layer of the atmosphere and passes through a square aperture with the dimension 20 mm×20 mm, and then enters the wave-front sensor (Fig. 1). The source and the experimental set-up are installed at a height of 55 cm, over an asphalted area, at a distance 360 m from each other. Within the wave-front sensor, the laser beam splits into two beams by an amplitude beam splitter. The second beam is reflected by a mirror, M, to a direction parallel to the first beam propagation direction. The beams strike the first gratings G1 and G3, respectively. The gratings G1 and G3 are installed on a single holder and the second gratings G2 and G4 are installed on another holder. The holders are installed on an optical rail of changeable length to choose the desired Talbot’s distance. The grating holders are held on rotary mounts which can rotate about the optical axis for adjustment of the angle between the grating rulings. In the experiment G2 and G4 are installed at a distance of 22.4 cm from G1 and G3 respectively. All the gratings are identical and have a ruling period of . Moiré patterns are formed on a plane where the G2 and G4 and a diffuser, D, are installed. Using a projection lens, PL, both of the moiré patterns are projected on a CCD camera. After alignment of the setup, the beam intensity was reduced by a neutral density filter to a level below the saturation level of the CCD. The moiré fringes were recorded at a sampling rate of 60 frames/s with the CCD camera (model DMK 21AU), and fed to the PC. In the described experiment, the digitized frames consist of 480×640 pixels. Moiré pattern in each channel consists of about 276×276 pixels and dm was covered by 40 pixels in both of channels. The measurement precision of AA fluctuations was about 1.5 arc sec in the x and y directions. Several sets of experimental data were recorded and digitized. Each set of data was collected in 5sec and contained 300 frames. The experiment was performed at IASBS. A typical recorded frame is shown in Fig. 2. Each frame consists of two sets of orthogonal moiré patterns. The low frequency illumination of the pattern of the fringes were recovered by using a spatial fast Fourier transform.
Once all 300 data frames of the data set were processed, the traces of bright and dark moiré fringes, the first order virtual traces, and the mean position of the traces (the reference positions of the traces) were determined to within one pixel accuracy. In Fig. 3 (a) and (d) typical moiré fringes in the horizontal and vertical directions are shown. In Fig. 3 (b) and (e) the corresponding low-frequency illumination of the patterns are shown. Derived traces of bright and dark moiré fringes and the first order virtual traces are shown in 3 (c) and (f). The background movie ( Media 1) of Fig. 3, contains all of the recorded patterns and the corresponding bright, dark, and first order virtual traces.
For all the bright, dark and the first order virtual traces of the moiré fringes in Fig. 2, the intersection points are determined. Displacement vectors of the intersection points with respect to their mean positions are calculated and are shown in Fig. 4. A larger shift is visualized by a larger length of the arrow. From the magnitude of the displacements in x and y directions and using Eqs. (4) and (5) the wave-front gradients in the corresponding directions were obtained. The background movie ( Media 2) of Fig. 4, contains the displacement vectors of all of the intersection points for 300 successive frames. In the described experiment, in order to perform the wave-front reconstruction from the measured wave-front slopes, we have used Hudgin’s discretization [18, 19]. The reconstructed wave-front surface plot, corresponding to distortions generated by atmospheric turbulence are shown in Fig. 5. The background movie ( Media 3) of Fig. 5, shows the evolution of wave-front in 300 successive frames. Let us now consider the temporal and spatial evolution of the wave-front in the background movie ( Media 3) of Fig. 5. A good fit to all of the reconstructed wave-front surfaces over the area 20 × 20 mm2 is almost a plane wave. As a result, the major atmospheric distortions over the selected area are tip and tilt. Thus, it is easy to say that Fried’s parameter is larger than the size of the beam size entering the wave-front sensor. Also, considering the temporal evolution of the reconstructed wave-front, it is clear that the atmospheric distortions are changing faster than our frame sampling rate. By some changes in the set-up, for example increasing the size of the gratings and the size of the beam splitter to cover a large area of the wave-front and using a high speed CCD camera, one can investigate more precisely the effect of the atmospheric turbulence. By employing the mentioned modification to the wave-front sensor, the atmospheric turbulence parameters can be determined more precisely. This work is in progress and will be published in the future.
In this paper we have presented an adjustable, high-sensitivity, wide dynamic range two channel wave-front sensor based on moiré deflectometry. We have also developed software for analyzing the moiré fringes by introducing the concept of virtual traces. In the proposed wave-front sensor the dynamic range and sensitivity of detection are adjustable by merely changing the separation of the gratings and the angle between the rulings of the gratings in both the moiré deflectometers. The spatial resolution of the method is also adjustable by means of bright, dark, and virtual traces for a given set of moiré fringes without paying a toll in the measurement precision. The implementation of the technique is straightforward. Compared to our previously proposed two channel wave-front sensor, this method is reliable and quite simple. The measurement is relatively insensitive to the alignment of the beam into the sensor. This sensor has many practical applications ranging from wave aberrations of human eyes to adaptive optics in astronomy. For low light applications as one would normally expect in astronomy, high sensitive detectors are required. Also, to reduce the energy loss, moiré deflectometers can be performed with phase gratings.
S. Rasouli would like to acknowledge Professor Yousef Sobouti for his valuable support for accomplishing the research described in this paper. Also S. Rasouli gratefully acknowledge partial support by the Institute for Advanced Studies in Basic Sciences (IASBS) Research Council under grant No. G2009IASBS107.
References and links
2. F. Roddier, Adaptive optics in astronomy (Cambridge university press, Cambridge, United Kingdom) (1999). [CrossRef]
3. R. V. Shack and B. C. Platt, “Production and use of a lenticular Hartmann screen,” J. Opt. Soc. Am. 61, 656 (1971).
4. B. C. Platt and R. Shack, “History and principles of Shack-Hartmann wavefront sensing,” J. Refract. Surg. 17, S573–S577 (2001). [PubMed]
6. G. W. R. Leibbrandt, G. Harbers, and P. J. Kunst, “Wavefront analysis with high accuracy by use of a double-grating lateral shearing interferometer,” Appl. Opt. 35, 6151–6161 (1996). [CrossRef] [PubMed]
7. R. A. Gonsalves, “Phase retrieval and diversity in adaptive optics,” Opt. Eng. 21, 829–832 (1982).
8. R. Ragazzoni and J. Farinato, “Sensitivity of a pyramidic wave front sensor in closed loop adaptive optics,” Astron. Astrophys. 350, L23–L26 (1999).
9. R. Legarda-Saenz, “Robust wavefront estimation using multiple directional derivatives in moiré deflectometry,” Opt. Lasers Eng. 45, 915–921 (2007). [CrossRef]
10. J. A. Quiroga, D. Crespo, and E. Bernabeu, “Fourier transform method for automatic processing of moiré deflectograms,” Opt. Eng. 38, 974–982 (1999). [CrossRef]
11. N. H. Salama, D. Patrignani, L. D. Pasquale, and E. E. Sicre, “Wavefront sensor using the Talbot effect,” Opt. Laser Technol. 31, 269–272 (1999). [CrossRef]
12. Ch. Siegel, F. Loewenthal, and L. E. Balmer, “A wavefront sensor based on the fractional Talbot effect,” Opt. Commun. 194, 265–275 (2001). [CrossRef]
13. R. Sekine, T. Shibuya, K. Ukai, S. Komatsu, M. Hattori, T. Mihashi, N. Nakazawa, and Y. Hirohara, “Measurement of wavefront aberration of human eye using Talbot image of two-dimensional grating,” Opt. Rev. 13, 207–211 (2006). [CrossRef]
14. M. Rottenkolber and H. Podbielska, “Measuring ophthalmologic surfaces by means of moiré deflectometry,” Opt. Eng. 35, 1124–1133 (1996). [CrossRef]
16. S. Rasouli, Anamparambu N. Ramaprakash, H. K. Das, C. V. Rajarshi, Y. Rajabi, and M. Dashti, “Two channel wavefront sensor arrangement employing moiré deflectometry,” in Optics in Atmospheric Propagation and Adaptive Systems XII , A. Kohnle, K. Stein, and J. D. Gonglewski, ed., Proc. SPIE7476, 74760K-1–9 (2009).
17. S. Rasouli and M. T. Tavassoly, “Application of the moiré deflectometry on divergent laser beam to the measurement of the angle of arrival fluctuations and the refractive index structure constant in the turbulent atmosphere,” Opt. Lett. 33, 980–982 (2008). [CrossRef] [PubMed]
18. J. Herrmann, “Least-squares wave front errors of minimum norm,” J. Opt. Soc. Am. 70, 28–35 (1980). [CrossRef]
19. B. R. Hunt, “Matrix formulation of the reconstruction of phase values from phase differences,” J. Opt. Soc. Am. 69, 393–399 (1979). [CrossRef]