Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Multi-frequency and multiple phase-shift sinusoidal fringe projection for 3D profilometry

Open Access Open Access

Abstract

In this paper, we report on a laser fringe projection set-up, which can generate fringe patterns with multiple frequencies and phase shifts. Stationary fringe patterns with sinusoidal intensity distributions are produced by the interference of two laser beams, which are frequency modulated by a pair of acousto-optic modulators (AOMs). The AOMs are driven by two RF signals with the same frequency but a phase delay between them. By changing the RF frequency and the phase delay, the fringe spatial frequency and phase shift can be electronically controlled, which allows high-speed switching from one frequency or phase to another thus makes a dynamic 3D profiling possible.

©2005 Optical Society of America

1. Introduction

Non-contact and full-field optical surface profilometry is becoming an important measurement tool in both scientific research and industrial applications including automated manufacturing, solid modeling, computer graphics, machine vision systems, etc. Different methods have been investigated and developed, among them those based on fringe projection and triangulation have attracted a lot of attention and have been widely applied in practical applications due to their superior performance (such as full-field and 3D fast measurement) as well as simple implementation [15]. In the fringe projection based profilometry, a set of fringes with a preferably sinusoidal intensity distribution is projected onto the surface to be measured and a CCD camera is positioned at a location with a viewing angle to record the distorted fringes. From the recorded images, one can calculate the phase map of the distorted fringes by using different methods, such as Fourier transform or phase shifting methods. The 3D profile of the measured object can be obtained from the phase map by using a transform determined by the geometry of the optical set-up. However, since the phase calculation involves the arctangent operation, the phase map calculated from the detected intensity distributions through either Fourier transform or phase shifting methods has principal values ranging from -π to π, and has discontinuities with 2π phase jumps. Therefore, a so-called phase unwrapping process is necessary to remove the phase jumps and to obtain continuous phase distribution. The phase unwrapping process is normally carried out using the so-called spatial unwrapping algorithm, which compares neighboring pixels and adjusts their relative phases by addition of integral multiples of 2π until the difference lies within the range of (-π, π). However, the spatial unwrapping algorithm suffers from the drawback that large phase errors can propagate across the image. Moreover it could be impossible to unwrap the phases correctly when the surfaces of interest contain discontinuities or the recorded images contain speckle noise [6]. Intensive studies have been devoted to develop reliable and robust phase unwrapping methods and algorithms, and a review on this topic can be found in [7].

Among the proposed phase unwrapping methods, the one based on multi-frequency fringe projection can be applied to the objects with discontinuities [8, 9]. In this method, multiple images are recorded at fringe patterns with different periods (from long to short) or frequencies (from low to high) in time series, and the phase at each pixel is unwrapped along the time axis. The longest fringe period is selected so that the phase across the whole object is within (-π, π) range. The phase map obtained from a lower fringe frequency will be used to guide the phase unwrapping process at a higher fringe frequency. In order to realize this phase unwrapping process, a fringe projection set-up which can generate fringe patterns with variable and controllable frequencies is desirable. On the other hand, in the phase shifting profilometry (PSP), fringe patterns with sinusoidal intensity distribution and multiple phase shifts are also necessary.

Multi-frequency and multiple phase-shift fringe patterns can be generated by using mechanical devices, such as a translation stage [10]. Due to hysteresis, backlash and wear, inaccurate shift such as unequal or incomplete shift often results, leading to undesirable errors. Therefore, it is still a challenging task to develop methods which can generate sinusoidal fringe patterns with controllable frequencies and phase shifts. In a paper published in 2000, a fringe projector using an acousto-optic modulator (AOM) has been demonstrated [11]. In this method, an AOM driven by two RF frequencies generates two sets of gratings inside the acousto-optic crystal and diffracts an incident laser beam into two first-order beams which are separated by a angle determined by the frequency difference of the two RF driving signals. The two first-order beams are brought together and expended by a lens. As the two diffracted beams are of different frequency shifts, the fringe pattern generated by the interference of the two first-order beams are moving with a velocity determined by the frequency difference. To make the fringes look stationary, a laser diode which amplitude is modulated by a signal synchronized with the RF driving signals is employed as the laser source. By controlling the frequency difference between the two RF driving signals, the period of the generated fringe pattern can be changed. Phase shifting is realized by controlling the time when the laser diode is switched on.

In this paper, we propose a fringe projection set-up also based on acousto-optic modulation. We demonstrate that by using two AOMs, a stationary fringe pattern with variable periods and phase shifts can be generated. In the following section, the principle of the proposed fringe projector is first introduced. Section 3 presents a theoretical analysis of the process of the fringe formation, followed by a description of the experiments carried out and the experimental results. Conclusions are given in Section 5.

 figure: Fig. 1.

Fig. 1. Schematic diagram of the proposed fringe projector. BS - beam splitter; AOM - acousto-optic modulator; WP- wedge prism, PH - pinholes; BC - beam combiner; MO - microscope objective.

Download Full Size | PDF

2. Experimental set-up and principle of operation

A schematic diagram of the projector proposed in this study is shown in Fig. 1. A beam from a CW laser is first split by a beam-splitter into two parallel beams with equal intensity. The two parallel beams then enter two acousto-optic modulators (AOMs) respectively. The AOMs are driven by two RF signals with identical frequency but different initial phases. The AOMs are arranged in such a way that diffractions generated by both AOMs occur in the horizontal plane and the diffracted beams are symmetrical about the optical axis of the system. The first-order diffracted beams from the AOMs are combined by the second beam-splitter (as a beam-combiner) and focused by a lens. Since the two beams are of the same frequency shift, a set of stationary interference fringes are produced in the intersecting area of the two beams as the result of laser interference. The interference fringes are projected onto the surface to be measured by a microscope objective. The AOMs are normally designed to work in the so-called Bragg regime where only the first-diffraction order is generated. In practice, however, the zero- and higher-order diffracted beams always exist. In order to filter out the beams other than the first orders, two pinholes can be used. Normally the AOMs need to work at a frequency higher than 40 MHz to achieve a high diffraction efficiency. In this case, the diffracted angles generated at such a frequency could be larger than the angles needed for the lowest fringe frequency. To adjust the angle between the first-order diffracted beams before combining them, two wedge prisms with a small angle are used.

When the frequency of the RF driving signals is varied, the diffraction angles of the first-order beams will change correspondingly. As a consequence, the period (or spatial frequency) of the fringe pattern generated by the interference will be changed. The phase shifts can be realized by introducing phase delay between the two RF driving signals, which can be done by using the well-developed phase shifting or time delay electronic circuitry.

As shown in Fig. 1, the optical layout of the projector is symmetrical about its optical axis. Therefore, the influence of the environmental conditions on the fringe patterns will be minimized. As demonstrated below, the intensity distribution of the fringe pattern generated by the interference of laser beams is sinusoidal.

3. Analysis

Suppose that the beam from the laser is a fundamental mode Gaussian beam. Mathematically the electric field of a Gaussian beam can be described by [12]

E(x,y,z,t)=E0wow(z)exp{x2+y2w2(z)}exp{i[kzωttan1(zz0)+kx2+y22R(z)]},

where, E 0 is the amplitude of the electric field of the beam; k is the propagation constant; w 0 is the radius of the beam waist; z is the distance along the beam propagation direction; ω is the frequency of the laser beam. Other parameters are defined as follows:

z0=πw02λ,
w(z)=w0[1+(zz0)2]12,
R(z)=z[1+(z0+z)2].

When a Gaussian beam is well collimated and the optical path is short (as in the case considered in this study), it can be simplified to a plane wave with Guassian intensity distribution, ie

E(x,y,z,t)=E0exp{x2+y2w02}exp{i[kzωt]}.

We assume that the acoustic waves inside the crystals of the AOMs are

A1=A10exp[i(K1X1Ωt+Φ1)],
A2=A20exp[i(K2X2Ωt+Φ2)].

where, A 10 and A 20 are the amplitudes of the acoustic waves inside the crystals; K 1 and K 2 are the propagation constants of the acoustic waves along the X 1 - and X 2 -directions respectively; Ω is the frequency of the acoustic waves and is equal to the frequency of the RF driving signals; Φ1 and Φ2 are the initial phases of the RF driving signals.

The electric fields of the first-order diffracted beams before leaving the acousto-optic crystals are

E1(x1,y1,z1,t)=E10exp{x12+y12w02}exp{i[k1z1(ω0+Ω)t+Φ1]},
E2(x2,y2,z2,t)=E20exp{x22+y22w02}exp{i[k2z2(ω0+Ω)t+Φ2]},

where, k 1and k 2 are propagation constants of the first-order diffracted beams produced by AOM1 and AOM2 respectively.

The propagation directions of the diffracted beams are adjusted by the wedge prisms. After the beam combiner, the beams are brought together by the focusing lens. At the focal point of the lens, two beams intersect with an angle, θ(Ω), which is determined by the focal length of the focusing lens, F, and the driving frequency of the AOMs, Ω. The microscope objective converts the input beams into spherical waves. In the space after the microscope objective, interference between two point sources located at (d(Ω)/2, 0, 0) and (- d(Ω)/2, 0, 0) happens. d(Ω) is determined by the focal lengths of the focusing lens and the microscope objective, and the AOM driving frequency, Ω.

The electric field at a point (x, y, z) in the fringe area is

E(x,y,z,t)=E10w0w(z)exp{r2w2(z)}exp{i[kz+r22R1(ω0+Ω)t+Φ1]}
+E20w0w(z)exp{r2w2(z)}exp{i[kz+r22R2(ω0+Ω)t+Φ2]},

where,

r 2=x 2+y 2,

R 1=[(x-d(Ω)/2)2+y 2+z 2]1/2 ;

R 2=[(x+d(Ω)/2)2+y 2+z2]1/2,

and w0 is the beam waist radius after the focusing lens and the microscope objective;

w(z)=w0[1+(λzπw02)2]12..

The interference intensity at (x, y, z) is given by

I(x,y,z)=[E(x,y,z,t)·E(x,y,z,t)*].

As well known and implied by the above equations, the interference patterns produced by the two point sources are planes parallel to the y-axis and the fringe spacing in a plane perpendicular to the z-axis linearly increases with the z-distance. The fringe spacing is also a function of the distance between the two point sources, d(Ω), meaning that the fringe spacing can be varied by changing the AOM driving frequency, Ω. At a point (x, y, z), the phase of the fringe pattern is determined by the phase difference between the two driving signals, Φ21. Therefore, one can introduce a phase shift to the fringe pattern by simply controlling the phase delay between the AOM driving signals.

By using the equations given above, one can calculate the fringe patterns produced by the projector. Shown in Fig. 2 is an example of the calculated fringe pattern projected on a convex surface with a plane as reference.

 figure: Fig. 2.

Fig. 2. Simulated fringe pattern projected on a convex surface with a plane as reference.

Download Full Size | PDF

4. Experiments

An experimental set-up was designed and constructed to test the fringe projector as shown in Fig. 1 and to demonstrate the controllability of the fringe spacing and the phase of the fringe pattern generated. In the experiments, a laser diode with an output power of 10 mW at 670 nm was used. Two identical integrated prisms were used as the beam-splitter and beam-combiner. The beam spacing is 25 mm. The crystals used in the AOMs are PbMoO4. The AOMs have a central frequency of 65 MHz and a bandwidth of 20 MHz. The focal length of the focusing lens is 10 mm, and the power of the microscope objective is 20X.

 figure: Fig. 3.

Fig. 3. Schematic diagram of the electronics developed for driving the AOMs, and for controlling the fringe spacing and the phase shift.

Download Full Size | PDF

Figure 3 shows a schematic diagram of the electronics developed for driving the AOMs, and for controlling the fringe spacing and the phase shift. A sinusoidal signal was generated by a voltage-controlled oscillator (VCO), and split into two channels. One of them was amplified to the level required by the AOM, and another was fed into a phase delay circuit where the input signal was phase shifted. Four fixed phase delays were generated with the current design. Using a multi-channel switching IC, one of the signals with different phases was selected by a digital control. The phase-shifted signal was amplified before feeding to the second AOM. The control of the VCO frequency was realized by either a DC voltage or a D/A converter.

 figure: Fig. 4.

Fig. 4. Image recorded by a CCD camera when the fringe pattern was projected on a dome-shaped object with a reference plane.

Download Full Size | PDF

Shown in Fig. 4 is an image recorded by a CCD camera when the fringe pattern was projected on a dome shape with a reference plane. The viewing angle was 30 degrees, corresponding to the simulated fringe pattern as given in Fig. 2.

In order to test the phase shift function of the projector, the CCD camera with imaging lens removed was directly positioned in the fringe area. Depicted in Fig. 5 is the fringes recorded at two states with a nominal phase shift of 180 degrees. In the experiments, the phase delays between the two RF driving signals to the AOMs were monitored and measured by using a digital oscilloscope.

 figure: Fig. 5.

Fig. 5. Fringes recorded with a nominal phase shift of 180 degrees.

Download Full Size | PDF

The developed projector was tested in a 3D measurement of a rectangle block. Shown in Fig. 6 are the deformed fringe patterns with three different fringe periods. The reconstructed 3D surface is given in Fig. 7.

 figure: Fig. 6.

Fig. 6. Recorded fringe patterns with three different periods projected on a rectangle block

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Reconstructed 3D surface plot of the rectangle block

Download Full Size | PDF

5. Conclusion

We have proposed and demonstrated a laser fringe projection method, which can generate fringe patterns with multiple frequencies and phase shifts. Stationary fringe patterns with sinusoidal intensity distributions can be produced by the interference of two laser beams, which are frequency shifted by using two AOMs. By controlling the frequency shifts, phase delays and the pointing directions of the two laser beams, the fringe spacing and the phase of the fringe pattern can be controlled electronically at a fast switching speed. Thus, precision mechanical motions can be replaced with precision time delays of electronic signals. No moving part is necessary.

The proposed AOM projector is an enabling technique to realize dynamic 3D profile measurement. Because the interference happens in whole intersection area, the AOM projector does not have the defocusing problem suffered by other grating projectors. In another word, it has unlimited depth of focus. The AOM projector can be applied to both large-scale and micro-scale objects.

The AOM projector will be further tested in our ongoing projects. By constructing a 3D profile measurement system, we will evaluate the overall performance of the system including the accuracy and resolution. This work is currently under way, and the results will be reported in a separate paper.

Acknowledgements

The work was supported by Natural Science Foundation of China (NSFC) under the grant No. 60275012 and No. 60472107, and partly supported by Natural Science Foundation of Gaungdong province (No. 031804) and Science & Technology Bureau of Shenzhen (No. 200341). The financial support from the Australian Research Council (ARC) is also acknowledged. The authors would like to thank the reviewers of the manuscript for their suggestions for revising the manuscript of the paper.

References and links

1. X. Su and W. Chen, “Fourier transform profilometry: a review,” Opt. Lasers Eng. 35, 263–284 (2001). [CrossRef]  

2. V. Srinivasan, H. C. Liu, and M. Halioua, “Automated phase-measuring profilometry of 3-D diffuse objects,” Appl. Opt. 23, 3105–3108 (1984). [CrossRef]   [PubMed]  

3. X. Su, W. S. Zhou, V. von Bally, and D. Vukicevic, “Automated phase-measuring profilometry using defocused projection of a Ronchi grating,” Opt. Commun. 94, 561–73 (1992). [CrossRef]  

4. M. Takeda, H. Ina, and S. Koboyashi, “Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry,” J. Opt. Soc. Am. 72, 156–160 (1982). [CrossRef]  

5. M. Takeda and H. Motoh, “Fourier transform profilometry for the automatic measurement of 3-D object shapes,” Appl. Opt. 22, 3977–3982 (1983). [CrossRef]   [PubMed]  

6. C. A. Hobson, J. T. Atkinson, and F. Lilley, “The application of digital filtering to phase recovery when surface contouring using fringe projection techniques,” Opt. Lasers Eng. 27, 355–368 (1997). [CrossRef]  

7. X. Su and W. Chen, Opt. “Reliability-guided phase unwrapping algorithm: a review,” Opt. Lasers Eng. 42, 245–261 (2004). [CrossRef]  

8. J. M. Huntley and H. O. Saldner, “Temporal phase unwrapping algorithm for automated interferogram analysis,” Appl. Opt. 32, 3047–3052 (1993). [CrossRef]   [PubMed]  

9. H. O. Saldner and J. M. Huntley, “Temporal phase unwrapping: application to surface profiling of discontinuous objects,” Appl. Opt. 36, 2770–2775 (1997). [CrossRef]   [PubMed]  

10. G. C. Jin, N.K. Bao, and P.S. Chung, “Applications of a novel phase-shift method using a computer-controlled polarization mechanism,” Opt. Eng. 33, 2733–2737 (1994). [CrossRef]  

11. M. S. Mermeistein, D. L. Feldkhum, and L. G. Shirley, “Video-rate surface profiling with acousto-optic accordion fringe interferometry,” Opt. Eng. 39, 106–113 (2000). [CrossRef]  

12. H. Kogelnik and T. Li, “Laser beams and resonators,” Appl. Opt. 5, 1550–1567 (1966). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1.
Fig. 1. Schematic diagram of the proposed fringe projector. BS - beam splitter; AOM - acousto-optic modulator; WP- wedge prism, PH - pinholes; BC - beam combiner; MO - microscope objective.
Fig. 2.
Fig. 2. Simulated fringe pattern projected on a convex surface with a plane as reference.
Fig. 3.
Fig. 3. Schematic diagram of the electronics developed for driving the AOMs, and for controlling the fringe spacing and the phase shift.
Fig. 4.
Fig. 4. Image recorded by a CCD camera when the fringe pattern was projected on a dome-shaped object with a reference plane.
Fig. 5.
Fig. 5. Fringes recorded with a nominal phase shift of 180 degrees.
Fig. 6.
Fig. 6. Recorded fringe patterns with three different periods projected on a rectangle block
Fig. 7.
Fig. 7. Reconstructed 3D surface plot of the rectangle block

Equations (12)

Equations on this page are rendered with MathJax. Learn more.

E ( x , y , z , t ) = E 0 w o w ( z ) exp { x 2 + y 2 w 2 ( z ) } exp { i [ k z ω t tan 1 ( z z 0 ) + k x 2 + y 2 2 R ( z ) ] } ,
z 0 = π w 0 2 λ ,
w ( z ) = w 0 [ 1 + ( z z 0 ) 2 ] 1 2 ,
R ( z ) = z [ 1 + ( z 0 + z ) 2 ] .
E ( x , y , z , t ) = E 0 exp { x 2 + y 2 w 0 2 } exp { i [ k z ω t ] } .
A 1 = A 10 exp [ i ( K 1 X 1 Ω t + Φ 1 ) ] ,
A 2 = A 20 exp [ i ( K 2 X 2 Ω t + Φ 2 ) ] .
E 1 ( x 1 , y 1 , z 1 , t ) = E 10 exp { x 1 2 + y 1 2 w 0 2 } exp { i [ k 1 z 1 ( ω 0 + Ω ) t + Φ 1 ] } ,
E 2 ( x 2 , y 2 , z 2 , t ) = E 20 exp { x 2 2 + y 2 2 w 0 2 } exp { i [ k 2 z 2 ( ω 0 + Ω ) t + Φ 2 ] } ,
E ( x , y , z , t ) = E 10 w 0 w ( z ) exp { r 2 w 2 ( z ) } exp { i [ k z + r 2 2 R 1 ( ω 0 + Ω ) t + Φ 1 ] }
+ E 20 w 0 w ( z ) exp { r 2 w 2 ( z ) } exp { i [ k z + r 2 2 R 2 ( ω 0 + Ω ) t + Φ 2 ] } ,
I ( x , y , z ) = [ E ( x , y , z , t ) · E ( x , y , z , t ) * ] .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.