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

Dispersion-free optical coherence depth sensing with a spatial frequency comb generated by an angular spectrum modulator

Open Access Open Access

Abstract

As an alternative to the conventional optical frequency comb technique, a spatial frequency comb technique is proposed for dispersionfree optical coherence depth sensing. Instead of generating an optical frequency comb over a wide range of time spectrum, we generate a spatial frequency comb by modulating the incident angle of a monochromatic plane wave with a spatial light modulator (SLM). The use of monochromatic light combined with the SLM enables dispersion-free depth sensing that is free from mechanical moving components.

©2006 Optical Society of America

1. Introduction

Optical coherence tomography and profilometry using an optical frequency comb (OFC) has the marked advantage of fast and direct measurements being free from mechanical moving components as well as being free from complicated signal processing [1–3]. By generating the OFC with stepwise frequency modulation of the light from a tunable laser and synchronous phase modulation faster than the integration time of the image sensor, Hotate et al. have successfully demonstrated the synthesis of a temporal coherence function that has high coherence peaks at arbitrarily specified depth locations [4, 5], and more recently Choi et al. have proposed improved OFC-based interferometry [6]. These methods presented a new concept of coherence synthesis in which a temporal coherence function was controlled by time modulation of the source spectrum. However, because OFC is composed of multiple spectral components whose spectral range must be extended for higher depth resolution, the technique suffers from spectral absorption and/or index dispersion problems, particularly when the object and/or the propagation medium have inhomogeneous spectral response as in the case of biological samples submerged in a liquid medium. This is a fundamental problem inherent to many other techniques making use of a light source with a broadband spectrum, among which are white-light interferometry [7–9] or low-coherence interferometry [10–12] and spectral interferometry [13–18].

A hint for a solution for this dispersion problem can be obtained from the analogy between space and time that can be found in many principles of optical metrology [19–21]. As an alternative to the use of OFC, we propose the use of a spatial frequency comb (SFC), in which the angular spectrum of quasi-monochromatic light is tailored to have a comb shape in the spatial frequency domain with a spatial-frequency-tunable source made of a spatial light modulator (SLM). The proposed technique enables spatial coherence depth sensing that is completely free from dispersion problems and mechanical moving components. Experimental results are presented that demonstrate the validity of the proposed principle.

2. Principle

The technique for generating SFC may be regarded as a natural extension of the idea of the angular spectrum scanning technique recently proposed by these authors [22] and also the tilt scanning technique proposed by Ruiz and Huntley [23]. For convenience of explanation, let us first briefly review the principle of the angular spectrum control for two-beam interferometry with a Michelson interferometer illustrated in Fig. 1. A point source S is placed on the focal plane of a lens L1 whose optical axis is normal to the surface of a reference mirror MR. One of the collimated rays (shown in red) exiting from lens L1 is reflected by beam splitter BS and reaches an observation point A on the surface of the object Obj; the point A is imaged onto a point à on an image sensor by lens L2. Another ray (shown in blue) comes to the same point à on the image sensor after being reflected at point B on the surface of a reference mirror. In Fig. 1, point A’ is the mirror image of point A with respect to the virtual reference mirror M’R. The propagation vector k (which will be referred to as the k-vector for short) of the collimated beam and the height vector h are in the direction of the vectors BA and AA , respectively. The phase difference Δφ between these two rays is given by

Δφ=k·2h=2khh=2hkcosθ,

where kh =k cos θ is the height component of the vector k, and θ is the angle of incidence to the reference surface defined by the angle between the vector k and the height vector h.

It is common practice in optical tomography and profilometry using OFC to adjust the k-vector k to be parallel to the height vector h so as to maximize the fringe sensitivity such that Δφ=-2hk with θ=0. The OFC is formed with equally-spaced multiple line spectra corresponding to wavenumbers k=nΔω/c, with n and Δω being an integer and the modefrequency separation, respectively. This characteristic of the OFC is illustrated in the k-vector space of Ewald sphere shown in Fig. 2.

 figure: Fig. 1.

Fig. 1. Two-beam interferometry with oblique illumination.

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. Ewald sphere representation of optical frequency comb (a) and spatial frequency comb (b).

Download Full Size | PDF

As shown in Fig. 2, the object is illuminated with the beam whose k-vector k(0) is parallel to the height vector h such that θ=0. Under this illumination, the OFC corresponds to a set of collinear k-vectors aligned parallel to the vector k(0) with their arrow tips equally spaced at an interval Δkω/c on the line through the center of the Ewald sphere. These radially distributed k-vectors inside the Ewald sphere cause the dispersion problems as they correspond to multiple optical frequencies. If we take a closer look at Eq. (1), we note an alternative solution in which we change the angle θ while keeping the optical frequency constant. In the k-space shown in Fig. 2, this operation corresponds to changing the cone angle θ of the k(θ) vector while keeping the radius of the k-sphere unchanged. The projected height component -kh of the k(θ) vector plays the role of the k(0) vector in the OFC. For example, if one can change θ over 0~30 degrees for the wavelength of 633nm, one can in principle realize the dispersion-free measurement with the performance comparable to the OFC with the wavelength range as wide as 98nm. We use a set of collimated monochromatic beams with different k-vectors angles whose longitudinal spatial frequency components kh are equally spaced in the direction of h, and refer to this technique as the spatial frequency comb (SFC) technique to differentiate it from the conventional optical frequency comb (OFC) technique. To generate the equally spaced SFC spectra equivalent to those of the OFC, the lateral component of the k-vector, k (θ), which is the projection of k(θ) onto the plane normal to the k(0) vector, should have a concentric Fresnel zone-plate structure as shown in Fig. 2(b). Just as the OFC generates a temporal coherence comb function with periodic coherence peaks with a time interval Δτ=πω, the SFC generates a longitudinal spatial coherence comb function with periodic coherence peaks with an axial distance interval Δh=πkh (where the double optical path in the reflection interferometer has been taken into account). The Fourier-transform relation of Wiener-Khinchin theorem between the temporal coherence function and the optical frequency spectrum is now replaced by the Fourier-transform relation of the generalized van Cittert-Zernike theorem between the longitudinal spatial coherence function and the longitudinal component of the spatial frequency spectrum, similarly to the McChutchen theorem [24]. To show this, let the fields at point A and point A’ in Fig. 1 be uA and uA′ , respectively. They are created as a superposition of the angular spectra of collimated beams from an extended quasimonochromatic source placed on the focal plane of the collimator lens:

uA=A(k)exp(ik·rA)dk
uA=A(k)exp(ik·rA)dk,

where rA and rA′ are position vectors pointing at point A and point A’ from an arbitrarily chosen origin, and the integration is taken over the surface of the Ewald sphere assuming a quasi-monochromatic source. The mutual intensity between the two points A and A’ are given by

JAA(h)=A(k)A*(k)exp[i(k·rAk·rA)]dkdk
=A(k)2exp[ik·(rA-rA)]dk=2πkA(kh)2exp(i2khh)dkh,

where we have assumed that the each elementary point source constituting the extended source is perfectly incoherent such that 〈A(k)A*(k′)〉=|A(k)|2 Δ(k-k′) with 〈〉 being an ensemble average, and have used the relation dk=2πksin θkdθ=2πkdkh about the ring surface element on the Ewald sphere. The complex degree of coherence is therefore given by

μAA(h)=JAA(h)JAA(0)JAA(0)=A(kh)2exp(i2khh)dkhA(kh)2dkh

For the SFC with the maximum angular spatial frequency k hmax=k and the comb spatial frequency interval Δkh

A(kh)2=n=0N1δ[kh(khmaxnΔkh)],

the complex degree of coherence becomes

μAA(h)=sin(NΔkhh)Nsin(Δkhh)exp{i[2k(N1)Δkh]h}..

The modulus of the complex degree of coherence has a comb shape with multiple coherence peaks separated by πkh as in the Hotate scheme, and the depth resolution will be given by Δh=2π/(NΔkh ), which is inversely proportional to the total comb frequency bandwidth.

The practical implementation of the SFC is a straightforward extension of the technique we developed for angular spectrum scanning [22]. Because tuning angular spectrum is equivalent to controlling the angle of the illumination beam, the SFC can be generated by a set of thin ring sources, with appropriate radii, placed coaxially on the focal plane of a collimator lens L1, which are created, for example, with a spatial light modulator (SLM). Suppose the optical axis of the collimator lens L1 is normal to the reference mirror surface MR and is parallel to the height vector by virtue of reflection at beam splitter BS. Then the h-component of the angular spectrum kh can be written as

kh=kcosθ=kfr2+f2k(1r22f2),

where r is the radius of the ring source, and in arriving at the last expression, the paraxial assumption has been made that fr. One can generate a SFC with the maximum angular spatial frequency k hmax=k and the comb spatial frequency interval Δkh by adjusting the radius of the n-th ring source as

rn=f2Δkhk×n.

Note that these ring sources correspond to the concentric Fresnel zone-plate structure shown in Fig. 2(b). In common with the angular spectrum scanning technique [22], the use of the ring source has the advantages in the amount of usable light and also in the robustness to the shading problem that occurs when a high and/or deep object is illuminated at a large incidence angle. The shadow-free illumination from the circular source solves the shadowing problem intrinsic to the SFC technique. To avoid the spurious interference fringe noises arising from the interference between the beams from different points on the ring source, the ring source should be a spatially incoherent (and yet temporally coherent) source. Such a source can be realized by placing a rotating ground glass on the source plane to destroy the spatial coherence (but preserve the temporal coherence).

 figure: Fig. 3.

Fig. 3. Optical geometry with tilt misalignment in reference mirror.

Download Full Size | PDF

Unlike the OFC made of aligned collinear k-vectors, the SFC has a freedom in the choice of the direction of k-vectors. This can be used for the compensation of misalignment of the optical system. Let us consider a misaligned optical geometry as shown in Fig. 3, in which reference mirror MR has a tilt with the normal vectors n, and n′ are not aligned to the optical axis. For this geometry, k·(rA-rA′) in Eq. (4) becomes k·2h′=2kh′h′, with h′ being a new vector 2h=AA shown Fig. 3. Since the k-vectors are now projected onto the new direction of the height vector h′, the projected spatial frequency spectrum |A(kh′ )|2 no longer has the spectral form of an equally spaced comb function. As the result of projection of the k-vectors of the beams from the ring sources, |A(kh′ )|2 can even have a continuous broadband spectrum, as shown in green color in Fig. 4(a), and gives rise to low spatial coherence for all points except for the point h′=0. In this way, the SFC generated by ring sources is extremely sensitive to a tilt alignment error of the reference mirror. This sensitivity to tilt may first look a disadvantage but our experience in experiments has revealed that this is in fact a great advantage in monitoring and adjusting the optical alignment. This characteristic of a ring source has also been used successfully for reducing the effects of coherent artifacts in interferometry by Kuechel [25]. Furthermore, this influence of mirror tilt can easily be compensated for by rotating the Ewald sphere to align its axis to the new height-vector h′, as shown in Fig. 4(b). To generate the equally spaced comb spectra in the direction of the h’- spatial frequency component, the lateral component of the k-vector k (θ), which is the projection of k(θ) onto the plane normal to the k(0) vector, should now have a decentered ellipse structure as shown in Fig. 4(b); the center of the shifted ellipse for the n-th SFC spectrum is given by (k-nΔkh′ )n, with n being the lateral component of the unit normal vector of the virtual reference mirror MR’. It should be noted that by combining the tilt of the reference mirror with the adequately tailored spatially incoherent elliptical sources generated by SLM, we can control the direction of the depth sensing. This is one of the unique characteristics of the SFC technique as compared with the optical frequency comb technique.

 figure: Fig. 4.

Fig. 4. Effect of the tilt of reference mirror and the depth sensing in h′ direction. (a) Misaligned geometry in which the tilt causes the projected concentric ring sources to function as a continuous broadband source. (b) Tilt-compensated geometry in which decentered ring sources can generate a spatial frequency comb.

Download Full Size | PDF

3. Relation to other techniques

Although the proposed SFC technique is inspired by the OFC technique and also by the analogy between space and time, it naturally has close relationship to several techniques that have been proposed for spatial coherence control. In this section, we clarify the relation of our technique to some of the relevant techniques.

Let us consider a special case where the SFC spectrum is replaced by a non-negative sinusoidal function of the same period Δkh and the same bandwidth Kh =(N-1) Δkh such that

A(kh)2=1+cos(2πkhΔkh)(kKhkhk).

We call this as a sinusoidal spatial frequency comb (SSFC). One can easily see from Eq. (5) that this will give a longitudinal coherence function

μAA(h)={sinc[Khh]+12sinc[Kh(h+πΔkh)]exp{i[2kKh](πΔkh)}
+12sinc[Kh(hπΔkh)]exp{i[2kKh](πΔkh)}}exp[i(2kKh)h]

which has three coherence peaks at h=0 and hπkh . When projected onto the lateral plane normal to the h vector as in Fig. 2(b), the SSFC becomes a zone plate in terms of the lateral component of the k-vector k=k2kh2 ;

A(k)2=1+cos{2π[k(k22k)]Δkh}.

It has now become clear that the SSFC technique corresponds to the technique for longitudinal coherence control proposed by Rosen and Takeda [26]. However, there are several important differences in their characteristics between the previously proposed SSFC and the newly proposed SFC. The SSFC requires the precise gray level control of SLM, and has only a pair of side coherence peaks, with their maximum modulus value being reduced down to a half, as seen in Eq. (11). The SFC, on the other hand, needs only simple binary control of SLM, and has equally spaced multiple high coherence peaks, with their modulus value being kept unity; this gives the possibility of the extending the dynamic range of depth sensing by making use of the higher order coherence peaks. Recently Gokhler and Rosen [27] proposed an alternative technique to generate multiple coherence peaks. They divided the ring sources into several sectors and assigned SSFCs of Eq. (10) with different spatial frequency intervals Δkh to these sectors, so that each sector generates coherence peaks at different locations. The advantage of the Gokhler-Rosen scheme is that the multiple coherence peaks need not be equally spaced as in the case of the SFC technique. However, since their technique adheres to the sinusoidal comb, the modulus value of the coherence peaks is halved as compared with the central peak, as seen in Eq. (11), and requires gray level control of SLM for its ideal performance. Furthermore, because of the sectored structure of the multiple zone plates, directional anisotropy is introduced to the sensitivity to tilt and the immunity to shading.

4. Experiments

4.1. Experimental setup

The schematic diagram of the depth-sensing system based on the SFC is illustrated in Fig. 5. Light from a He-Ne laser (10mW) is expanded and collimated by a beam expander EX and lens L1. A spatial light modulator (DMD Discovery TM 1100 Starter Kit, Texas Instruments) based on DMD (Digital Micromirror Device) is used to project the light source pattern generated by a personal computer PC. A diaphragm D1 is placed at the entry side of DMD to adjust the proper area of the illumination. The spatially-modulated light reflected from DMD enters a confocal lens pair consisting of lens L2, L3 and a diaphragm D2, which functions as a spatial filter to remove the effect of discrete pixels of DMD and to adjust the illumination level. A SFC displayed on DMD is relayed by the confocal lens pair and imaged onto a rotating ground glass GG placed on the front focal plane of lens L4, to generate a spatially incoherent source. The light from the source created on this rotating ground glass is introduced into a Michelson interferometer composed of a beam splitter BS, a reference mirror MR, and an object GB made of a pair of gauge brocks. Lens L5 images the interference fringe pattern on the object surface onto CCD. A quarter-wave plate QWP is placed directly behind the laser to generate a circularly polarized light for uniformity of reflection at DMD. The object is a pair of gauge blocks with heights h1=1.600mm and h2=2.000mm, which are placed on an optical plate. Lens L5 is focused on the surfaces of the gauge blocks. The magnification of the imaging lens L5 is adjusted to ~1×. The focal length of lens L2, L3 and L4 are f2 =120mm, f3 =250mm and f4 =150mm, respectively. Before starting the measurement, we first performed alignment by placing the center of the ring sources on the axis of lens L4 so as to ensure that the angular spectrum from each point source on the ring gives the same spatial comb-frequency component kh =kcosθ. The detail of the alignment procedure has been described in our previous paper [22], and, as mentioned in Sec. 2, we made use of the fact that the ring source is highly sensitive to misalignment.

 figure: Fig. 5.

Fig. 5. Experimental setup.

Download Full Size | PDF

4.2. Spatial frequency comb of fireworks

For convenience of explanation, we have so far assumed that the SFC is implemented in the form of concentric thin ring sources displayed on SLM. However, if the thin ring sources have the same thickness and the same irradiance, the projected spatial frequency spectra |A(kh′ )|2 cannot create such an ideal SFC with uniform comb heights as modeled in Eq. (6), but it will generate a tapered comb whose comb height varies in proportion to the peripheral length of the ring sources. One could solve this problem by adjusting the irradiance of the ring sources by introducing the gray level control of SLM, but this would compromise the advantage of SFC over SSFC. One could also control the thickness of the ring sources, but inner rings would become too thick to function as the ideal sources for comb spectra. To get around this problem, we propose a new implementation scheme which we call a spatial frequency comb of fireworks (SFCF). As its name stands for, the source structure looks like a firework as shown in Fig. 6(a). In stead of a continuous ring structure, we adopted an array structure of discrete point sources of the same irradiance on a polar mesh. Specifically, we placed 18 discrete point sources at equal azimuthal angular intervals on the periphery of the ring. In Fig. 6(b) shows the radial distribution of the irradian ce on the source plane, and (c) shows the the corresponding comb spectrum in the spatial frequency domain. We excluded the central point source at r 0 and its spectrum k 0, because 18 point sources on the ring degenerate into a single point source at the center, and the comb spectrum for this point is reduced down to 1/18.

 figure: Fig. 6:

Fig. 6: (a) Source structure for a spatial frequency comb of firework (SFCF). (b) Comb spectrum in the source plane. (c) Comb spectrum in spatial frequency domain.

Download Full Size | PDF

The full bandwidth of the SFCF can be calculated from the predefined experimental parameters. We generated the SFCF on a M×M pixel matrix of the DMD, where M=600, and the pixel size of the DMD is ΔxΔy=13.7µm×13.7µ. Considering the magnification of the confocal lens system composed of L2 (f 2=120mm) and L3 (f 3=250mm), and the field angle of the Fourier transform lens L4 (f 4=150mm), the full spatial frequency bandwidth of the SFCF can be calculated as k(f 3/f 2)2 (MΔx/2)2/(2f42). For the SFCF generated by a set of discrete point sources on the polar mesh composed of N rings, the spatial comb frequency interval is given by

Δkh=k(f3f2)2(MΔx)28f42N=16.19N[radians/mm],

which is equivalent to an OFC with the optical frequency interval of

Δν=cΔkh2π=kc(f3f2)2(MΔx)216πf42N=0.773N[THz],

where c is light velocity. The corresponding coherence comb interval becomes

Δh=πΔkh=c2Δν=4λ(f2f4f3MΔx)2N=0.194N0.2N[mm].

In the following proof-of-principle experiments, we use the number of rings N in the firework source as a parameter specifying SFCF.

4.3. Characteristic of longitudinal coherence comb function

To demonstrate the characteristic of the longitudinal coherence comb function predicted from Eq. (7), we observed the fringe contrast on one of the gauge block surfaces for the varying height h by scanning the reference mirror. We carried out experiments for two different SFCFs with the ring parameters N=8 and N=16. The results are shown in Fig. 7. As predicted from Eq. (15), high contrast fringes were observed when the height becomes integer multiple of Δh, which means h=1.6mm for the SFCF with N=8, and h=3.2mm for the SFCF with N=16.

 figure: Fig. 7.

Fig. 7. Variation of fringe contrast with the object height relative to the reference mirror. Two different parameters N=8 and N=16 were adopted for the SFCF.

Download Full Size | PDF

4.4. Depth ssensing by spatial comb frequency scanning

Next we fixed the offset height of the gauge blocks and scanned the spatial comb frequency parameter N of SFCF from 1 to 32. Several examples of fringe patters selected from the total 32 steps are shown in Fig. 8. During 32-step scanning, high-contrast fringes appeared at NL=2, 4, 8, and 16 on the left surface, and at NR=7 and 14 on the right surface of the gauge block set. From Eq. (15), the heights of the left and right surfaces, hL and hR , relative to the reference mirrors are given by the integer multiple of the coherence comb intervals for the largest SFCF parameters NL =16 and NR =14, such that hL =3.2mL [mm] and hR =2.8mR [mm], with mL and mR being integer numbers. Our a priori knowledge that the heights do not exceed 6.4mm and their difference hL -hR is less than 0.5mm permits us to uniquely determine the integers as mL =1 and mR =1. We therefore have hL =3.2mm and hR =2.8mm, and their height difference 0.4mm is in agreement to the nominal height-difference (400 µm) of the block gauge set.

 figure: Fig. 8.

Fig. 8. Examples of fringe patterns selected from the total 32 steps of SFCF scanning.

Download Full Size | PDF

4.5. Relation of SFC to angular spectrum scanning

While the angular spectrum scanning technique in our previous paper was based on the sequential display of a single angular spectrum, the SFC and SFCF techniques proposed in this paper are based on a set of multiple spatial angular spectra displayed simultaneously. In order for the SFC or the SFCF to produce a high coherence peak for the specified height, mutually incoherent point sources constituting SFC or SFCF should individually produce inphase fringe signals that are superposed constructively on intensity basis at the specified position of observation, much in the same manner as a mode-locked laser produces high-peak pulses at a specified time interval as a result of in-phase superposition of multiple waves. We call this situation as SFC or SFCF being tuned to the height. On the other hand, when SFC or SFCF is detuned from the specified height, the fringes will no longer be phase-locked, and will vary with the position of the point source. Then the fringes are averaged out as the result of superposition.

 figure: Fig. 9.

Fig. 9. (a). Fringe pattern for the illumination with a SFCF (N=8) tuned to left surface. (b) Fringe pattern for the illumination with the first ring source of the SFCF (N=8). (c) Fringe pattern for the illumination with the third ring source of the SFCF (N=8). Note the phase shift of the fringes on the right surface as compared to (b). (d) The fringe intensity on the SFCF-tuned left surface does not oscillate with the shift of the ring source. (f) The fringe intensity on the SFCFdetuned right surface varies with the shift of the ring source. (g) The sinusoidal variation of the fringe intensity on the left surface with the shift of the ring source for the SFCF (N=32) tuned to 6.4mm height. The red circles indicate the intensity at the ring position for the SFCF (N=8) tuned to the left surface.

Download Full Size | PDF

To confirm this observation by experiment, we adjusted the height of the left surface of the gauge block to hL =3.2mm with respect to the reference mirror, and generated a SFCF with eight rings N=8, which is tuned to the left surface by its secondary coherence peak as shown in Fig. 7, but which is detuned to the right surface whose height is hR =2.8mm. As seen in Fig. 9(a), the contrast of fringes on the left surface is high, but the fringe contrast on the right surface is much lower. Instead of illuminating all the rings simultaneously, we illuminated only the point sources on a selected single ring at a time. The selected ring is shifted sequentially from the inner ring to the outer ring, and the fringes on the left and right surfaces are observed at each step of the sequential ring shift. In Fig. 9(b) and 9(c) show fringe patterns observed by the illumination with the first and third ring sources of the SFCF (N=8), respectively. Note the phase shift of the fringes on the right surface between (b) and (c). As shown in Fig. 9(d), the fringe intensity on the left surface does not change significantly with the shift of the ring source because the SFCF (N=8) is tuned on this surface; the slight decrease of the intensity is due to the vignetting of the light beam by the limited acceptance angle of the interferometer. On the other hand, the fringe intensity on the SFCF-detuned right surface varies significantly with the shift of the ring source, as seen in Fig. 9(e), serving as the evidence for the predicted large phase change for the SFCF-detuned height. Next we increased the comb parameter of N of the SFCF to 32, which is tuned to the distance 6.4mm by its first coherence peak, with 1/4 comb frequency interval. Fig. 9(g) shows the sinusoidal variation of the fringe intensity on the left surface with the shift of the ring source by the (four times finer) comb frequency interval for the SFCF (N=32); this operation is similar to the angular spectrum scanning technique in our previous paper. The red circles indicate the intensity at the ring position for the SFCF (N=8) tuned to the left surface. Note that the fringes are all in phase at the ring position of the tuned SFCF (N=8), which is in agreement to our theoretical observation.

4.6. Comparison between SFCF and SSFC

Next we filled the interval of the spatial frequency comb spectra with one period of sinusoidal function to generate SSFC with a Fresnel-zone-plate (FZP) light source distribution. We first illuminated the interferometer with a SFCF by varying N=4, 8, 16, and then switched DMD from a binary to a grey mode of operation using an ALP controlling board (ViALUX, GmbH) and sequentially generated a SSFC of the corresponding parameters N=4, 8, 16. The result is shown in Fig. 10. Whereas the SFCF produces the multiple coherence peaks (the primary coherence peak for N=8, and the secondary coherence peak for N=4), the SSFC produces only one coherence peak (the primary coherence peak for N=8), as predicted from Eq. (11).

 figure: Fig. 10.

Fig. 10. Comparison between SFCF and SSFC.

Download Full Size | PDF

5. Conclusion

Inspired by the optical frequency comb (OFC) technique and the analogy between the roles played by space and time in many optical systems, we have proposed a new technique of depth sensing based on the spatial frequency comb (SFC) and clarified its relation to other relevant techniques. We have presented experimental results that demonstrate the validity of the principle of the proposed SFC technique, which can serve an alternative to the conventional OFC technique. The use of monochromatic light combined with the SLM enables dispersion-free absolute interferometry that is free from mechanical moving components. The use of SFCF with a firework-like structure was proposed, which can improve the uniformity of the comb heights and solve the shadowing problem.

Acknowledgments

Part of this work was supported by Grant-in-Aid of JSPS B (2) No. 18360034, and The 21st Century Center of Excellence (COE) Program on “Innovation of Coherent Optical Science” granted to The University of Electro-Communications, from Japanese Government. Zhihui Duan gratefully acknowledges the scholarship given from Ministry of Education, Culture, Sports, Science and Technology of Japan.

1. K. Hotate and T. Okugawa, “Selective extraction of a two-dimensional optical image by synthesis of the coherence function,” Opt. Lett. 17, 1529–1531 (1992). [CrossRef]   [PubMed]  

2. K. Hotate and O. Kamatani, “Optical coherence domain reflectometry by synthesis of coherence function,” J. Lightwave Technol. 11, 1701–1709 (1993). [CrossRef]  

3. S.-J. Lee, B. Widiyatmoko, M. Kourogi, and M. Ohtsu, “Ultrahigh scanning speed optical coherence tomography using optical frequency comb generators,” Jpn. J. Appl. Phys. Part II , 40, L878–L880 (2001). [CrossRef]  

4. K. Hotate and T. Okugawa, “Optical information processing by synthesis of the coherence function,” J. Lightwave Technol. 12, 1247–1255 (1994). [CrossRef]  

5. Z. He and K. Hotate, “Synthesized optical coherence tomography for imaging of scattering objects by use of a stepwise frequency-modulated tunable laser diode,” Opt. Lett. 24, 1502–1504 (1999). [CrossRef]  

6. S. Choi, M. Yamamoto, D. Moteki, T. Shioda,Y. Tanaka, and T. Kurokawa, “Frequency-comb-based interferometer for profilometry and tomography,” Opt. Lett. 31, 1976–1978 (2006). [CrossRef]   [PubMed]  

7. P. A. Flournoy, R. W. McClure, and G. Wyntjes, “White-light interferometric thickness gauge,” Appl. Opt. , 11, 1907–1915 (1972). [CrossRef]   [PubMed]  

8. N. Tanaka, M. Takeda, and K. Matsumoto, “Interferometrically measuring the physical properties of test object,” United States Patent 4,072,422 (filed October 20 1976, 1978).

9. M. Davidson, K. Kaufman, I. Mazor, and F. Cohen, “An application of interference microscopy to integrated circuit inspection and metrology,” in Integrated Circuit Metrology, Inspection, and Process Control, K. M. Monahan, ed., Proc. SPIE775, 233–247 (1987).

10. B. S. Lee and T. C. Strand, “Profilometry with a coherence scanning microscope,” Appl. Opt. 29, 3784–3788 (1990). [CrossRef]   [PubMed]  

11. T. Dresel, G. Hausler, and H. Venzke, “Three-dimensional sensing of rough surfaces by coherence radar,” Appl. Opt. 31, 919–925 (1992). [CrossRef]   [PubMed]  

12. D. Huang, E. A. Swanson, C. P. Lin, J. S. Schuman, W. G. Stinson, W. Chang, M. R. Hee, T. Flotte, K. Gregory, C. A. Puliafito, and J. G. Fujimoto, “Optical coherence tomography,” Science 254, 1178–1181 (1991). [CrossRef]   [PubMed]  

13. J. Schwider and L. Zhou, “Dispersive interferometric profilometer,” Opt. Lett. 19, 995–997 (1994). [CrossRef]   [PubMed]  

14. M. Takeda and H. Yamamoto, “Fourier-transform speckle profilometry: Three-dimensional shape measurements of diffuse objects with large height steps and/or spatially isolated surfaces,” Appl. Opt. 33, 7829–7837 (1994). [CrossRef]   [PubMed]  

15. S. Kuwamura and I. Yamaguchi, “Wavelength scanning profilometry for real-time surface shape measurement,” Appl. Opt. 36, 4473–4482 (1997). [CrossRef]   [PubMed]  

16. H. J. Tiziani, B. Franze, and P. Haible, “Wavelength-shift speckle interferometry for absolute profilometry using mode-hope free external cavity diode laser,” J. Mod. Opt. 44, 1485–1496 (1997). [CrossRef]  

17. M. Kinoshita, M. Takeda, H. Yago, Y. Watanabe, and T. Kurokawa, “Optical frequency-domain microprofilometry with a frequency-tunable liquid-crystal Fabry-Perot etalon device,” Appl. Opt. 38, 7063–7068 (1999). [CrossRef]  

18. D. S. Mehta, M. Sugai, H. Hinosugi, S. Saito, M. Takeda, T. Kurokawa, H. Takahashi, M. Ando, M. Shishido, and T. Yoshizawa, “Simultaneous three-dimensional step-height measurement and high-resolution tomographic imaging with a spectral interferometric microscope,” Appl. Opt. 41, 3874–3885 (2002). [CrossRef]   [PubMed]  

19. M. Takeda, “The philosophy of fringes: Analogies and dualities in fringe generation and analysis,” in Fringe ’97 Automatic Processing of Fringe Patterns, W. Jueptner and W. Osten, eds., Akademie Verlag Series in Optical Metrology, (Akademie Verlag, Berlin, 1997), pp.17–26.

20. M. Takeda, “Space-time analogy in absolute optical profilometry using frequency-scan techniques: A tutorial review,” in Perspectives in Engineering Optics, K. Singh and P. K. Rastogi, eds, (Anita Publications, New Delhi, 2003), pp.192–198.

21. J. Ch. Vienot, R. Ferriere, and J. P. Goedgebur, “Conjugation of space and time variables in optics,” in Optics in Four Dimensions, M. A. Machado and L. M. Narducci, eds., No. 65 of AIP Conference Proceedings, (American Institute of Physics, New York, 1981), pp.49–62.

22. Z. Duan, Y. Miyamoto, and M. Takeda, “Dispersion-free absolute interferometry based on angular spectrum scanning,” Opt. Express 14, 655–663 (2006). [CrossRef]   [PubMed]  

23. P. D. Ruiz and J. M. Huntley, “Depth-resolved displacement measurement using tilt scanning speckle interferometry,” in Fringe 2005 Automatic Processing of Fringe Patterns, W. Osten, ed. (Springer, Berlin, 2005), pp. 238–241.

24. C. W. McCutchen, “Generalized source and the van Cittert-Zernike theorem: A study of spatial coherence required for interferometry,” J. Opt. Soc. Am. 56, 727–733 (1966). [CrossRef]  

25. M. Kuechel, “Apparatus and method(s) for reducing the effects of coherent artifacts in an interferometry,” US Patent 6804011 B2, (2004) or US Patent Application 2003/0030819 A1, (2003).

26. J. Rosen and M. Takeda, “Longitudinal spatial coherence applied for surface profilometry,” Appl. Opt. 37, 4107–4111 (2000). [CrossRef]  

27. M. Gokhler and J. Rosen, “Synthesis of a multiple-peak spatial degree of coherence for imaging through absorbing media,” Appl. Opt. 44, 2921–2927 (2005). [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 (10)

Fig. 1.
Fig. 1. Two-beam interferometry with oblique illumination.
Fig. 2.
Fig. 2. Ewald sphere representation of optical frequency comb (a) and spatial frequency comb (b).
Fig. 3.
Fig. 3. Optical geometry with tilt misalignment in reference mirror.
Fig. 4.
Fig. 4. Effect of the tilt of reference mirror and the depth sensing in h ′ direction. (a) Misaligned geometry in which the tilt causes the projected concentric ring sources to function as a continuous broadband source. (b) Tilt-compensated geometry in which decentered ring sources can generate a spatial frequency comb.
Fig. 5.
Fig. 5. Experimental setup.
Fig. 6:
Fig. 6: (a) Source structure for a spatial frequency comb of firework (SFCF). (b) Comb spectrum in the source plane. (c) Comb spectrum in spatial frequency domain.
Fig. 7.
Fig. 7. Variation of fringe contrast with the object height relative to the reference mirror. Two different parameters N=8 and N=16 were adopted for the SFCF.
Fig. 8.
Fig. 8. Examples of fringe patterns selected from the total 32 steps of SFCF scanning.
Fig. 9.
Fig. 9. (a). Fringe pattern for the illumination with a SFCF (N=8) tuned to left surface. (b) Fringe pattern for the illumination with the first ring source of the SFCF (N=8). (c) Fringe pattern for the illumination with the third ring source of the SFCF (N=8). Note the phase shift of the fringes on the right surface as compared to (b). (d) The fringe intensity on the SFCF-tuned left surface does not oscillate with the shift of the ring source. (f) The fringe intensity on the SFCFdetuned right surface varies with the shift of the ring source. (g) The sinusoidal variation of the fringe intensity on the left surface with the shift of the ring source for the SFCF (N=32) tuned to 6.4mm height. The red circles indicate the intensity at the ring position for the SFCF (N=8) tuned to the left surface.
Fig. 10.
Fig. 10. Comparison between SFCF and SSFC.

Equations (17)

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

Δ φ = k · 2 h = 2 k h h = 2 h k cos θ ,
u A = A ( k ) exp ( i k · r A ) d k
u A = A ( k ) exp ( i k · r A ) d k ,
J AA ( h ) = A ( k ) A * ( k ) exp [ i ( k · r A k · r A ) ] d k d k
= A ( k ) 2 exp [ i k · ( r A - r A ) ] d k = 2 π k A ( k h ) 2 exp ( i 2 k h h ) d k h ,
μ AA ( h ) = J AA ( h ) J AA ( 0 ) J A A ( 0 ) = A ( k h ) 2 exp ( i 2 k h h ) d k h A ( k h ) 2 d k h
A ( k h ) 2 = n = 0 N 1 δ [ k h ( k h max n Δ k h ) ] ,
μ AA ( h ) = sin ( N Δ k h h ) N sin ( Δ k h h ) exp { i [ 2 k ( N 1 ) Δ k h ] h } . .
k h = k cos θ = k f r 2 + f 2 k ( 1 r 2 2 f 2 ) ,
r n = f 2 Δ k h k × n .
A ( k h ) 2 = 1 + cos ( 2 π k h Δ k h ) ( k K h k h k ) .
μ AA ( h ) = { sinc [ K h h ] + 1 2 sinc [ K h ( h + π Δ k h ) ] exp { i [ 2 k K h ] ( π Δ k h ) }
+ 1 2 sinc [ K h ( h π Δ k h ) ] exp { i [ 2 k K h ] ( π Δ k h ) } } exp [ i ( 2 k K h ) h ]
A ( k ) 2 = 1 + cos { 2 π [ k ( k 2 2 k ) ] Δ k h } .
Δ k h = k ( f 3 f 2 ) 2 ( M Δ x ) 2 8 f 4 2 N = 16.19 N [ radians/ mm ] ,
Δ ν = c Δ k h 2 π = k c ( f 3 f 2 ) 2 ( M Δ x ) 2 16 π f 4 2 N = 0.773 N [ THz ] ,
Δ h = π Δ k h = c 2 Δ ν = 4 λ ( f 2 f 4 f 3 M Δ x ) 2 N = 0.194 N 0.2 N [ mm ] .
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.