Next Article in Journal
Rethinking the Design of Low-Cost Point-of-Care Diagnostic Devices
Previous Article in Journal
Microfluidical Microwave Reactor for Synthesis of Gold Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Topography Measurement of Large-Range Microstructures through Advanced Fourier-Transform Method and Phase Stitching in Scanning Broadband Light Interferometry

1
State Key Laboratory of Optical Technologies for Microfabrication, Institute of Optics and Electronics, Chinese Academy of Sciences, Chengdu 610209, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Micromachines 2017, 8(11), 319; https://doi.org/10.3390/mi8110319
Submission received: 13 September 2017 / Revised: 11 October 2017 / Accepted: 23 October 2017 / Published: 26 October 2017

Abstract

:
Scanning broadband light interferometry (SBLI) has been widely utilized in surface metrology due to its non-contact and high-accuracy method. In SBLI, phase evaluation through Fourier Transform (FT) is a prevalent and efficient technique, where the topography measurement can often be achieved through one interferogram. Nevertheless, the accuracy of the FT method would be significantly influenced by intensity modulation depth: “the lower the modulation of the pixel, the higher the error probability of its phase assignment”. If the structure has a large enough range along the z-axis, several areas in an individual interferogram would be weakly modulated due to the limited depth of focus (DOF). In this paper, we propose an advanced FT-based method when it comes to large-height structures. Spatial modulation depth is first calculated for each interferogram independently. After that, a binary control mask is reasonably constructed to identify the pixels that are valid for phase unwrapping. Then, a phase stitching method along the z-axis is carried out to conduct the large-height topography measurement within a giving field of view. The theoretical principle, simulation, and experimental validation are elaborated to demonstrate that the method can achieve an improved robustness for the reconstruction of large-range microstructures, the advantages of which include the elimination of stepping errors, the suppression of light fluctuations, and the freedom of a limited DOF.

1. Introduction

Three-dimensional measurement of microstructure is an area with growing needs and interest since it plays an important role in addressing quality issues and observing the performance of micro products. Given its high-precision, non-contact, and full-field method, the “topography measurement,” which is based on scanning broadband light interferometry (SBLI), has been widely applied in measuring and calibrating microstructures ranging from the nanometer to the micrometer scale [1,2,3,4]. Nonetheless, the accuracy of most previously recovering algorithms, such as the centroid algorithm, the peak location method, and the envelope detection technique, are at times negatively influenced by moving errors and intensity fluctuations [5,6,7,8,9]. Typically, phase evaluation is a general and effective method for precise metrology in SBLI. In phase measuring profilometry, there are mainly two techniques to retrieve the phase data including phase-shifting interferometry (PSI) and Fourier transform (FT) in the spatial frequency domain. For PSI, the phase mapping from N frames to obtain a wrapped phase is pixel-to-pixel analysis. Nevertheless, the PSI-based method would be significantly limited in a complex environment, since the accuracy strongly depends on stepping precision and light stability [10,11,12,13,14,15].
Contrarily, the Fourier transform is always a reliable technique to retrieve the phase distribution without moving errors and light fluctuations, for the reason that every interferogram is analyzed independently and filtering is adopted in the spatial frequency domain [16,17,18,19]. Additionally, as was pointed out in [20], “to make reliable phase measurements, the incident intensity must modulate sufficiently at each detector point to yield an accurate phase.” In the FT method, the accuracy of phase measurement directly depends on the modulation depth, “the lower the modulation depth of the pixel, the higher the error probability (the possibility to cause errors) of its phase assignment” [21,22,23]. Actually, when the measurement sample has a large enough range along the z-axis, several parts in the interferogram would be weakly modulated due to the limited depth of focus (DOF) of the imaging system. Particularly, devices in many science domains, including biomedical engineering, chemistry, and integrated circuits (ICs), are at times large-height structures [24,25,26,27]. In that case, when the height range is beyond the DOF, only a part of areas can be recognized at a definite scanning position, causing difficulty to reshape the topography within a giving field of view through one interferogram.
In this paper, we propose an advanced FT-based method and phase stitching technique when it comes to large-height structures. Spatial modulation depth is first calculated in the {xy} plane (imaging plane) for each charge-coupled device (CCD) image independently. After that, a binary control mask, which is constructed based on modulation threshold, is developed to identify the pixels that are valid for phase unwrapping in each interferogram. The threshold limit of the intensity modulation is appropriately set, and all pixels that do not modulate sufficiently are excluded. Then, a phase stitching method along the z-axis is proposed to complete the large-range measurement. The principle has been elaborated in detail, which contains modulation extraction, identification for valid pixels, the phase stitching process, and reconstruction algorithms. Simulation on a slope plane is organized to explain the procedure and verify the effectiveness. Additionally, the experimental results on a micro lens with a height of nearly 3000 nm show that the proposed method can precisely reshape the surface topography, the advantages of which include freedom of stepping errors, the suppression of light fluctuations, and a significant improvement in the measurement range.

2. Principle

Figure 1 illustrates a Mirau correlation microscope-based system of SBLI for dimensional metrology of micro samples. In this system, a broadband light is utilized to illuminate the test and reference surface. Light transmits through the interferometry objective, which contains a beam splitter and a several-millimeter reference mirror, where the light is divided into two beams by the beam splitter (the transmittance ratio is about 50%). One is reflected on the reference mirror and the other is on the measurement structure. The superposition beams are relayed to the CCD camera and produce an interferogram of the sample topography that can be captured by the individual CCD pixels. This technique is achieved by vertically scanning the sample in a high-precision piezo-electric transducer (PZT) stage, where the interferograms can be rapidly detected at different scanning positions.
In broadband light interferometry, when intensities from the measurement arm and reference arm are equal, the interference intensity field can be described as
I ( z ) = I 0 { 1 + exp [ ( z h z 0 l c ) 2 ] cos ( 4 π z h z 0 λ 0 ) }
where I0 is the background irradiance, z presents the scanning distance, h describes the sample height, z0 denotes the length of reference arm, lc illustrates the coherence length of the light source, and λ0 is the central wavelength [5,28]. When an interferogram is obtained at a definite scanning position zn, the irradiance distribution can be expressed as
I ( x , y ) = A ( x , y ) + B ( x , y ) cos ( 4 π z n h ( x , y ) z 0 λ 0 ) , with   B ( x , y ) = A ( x , y ) exp [ ( z n h ( x , y ) z 0 l c ) 2 ]
where A(x, y) presents the background intensity in a definite position, B(x, y) can be deemed as the fringe contrast.
In this method, a measured sample is slightly tilted to develop a small angle between the light axis and sample surface to produce a carrier frequency, which gives
I ( x , y ) = A ( x , y ) + B ( x , y ) cos ( 4 π z n ( h ( x , y ) k x ) z 0 λ 0 )
where k x is the initiatively tilted slope along the x axis, while k is the slope factor. Then, the equation can be rewritten as
I ( x , y ) = A ( x , y ) + B ( x , y ) cos ( 4 π k x λ 0 + 4 π z n z 0 h ( x , y ) λ 0 )
To simplify the expression, the intensity distribution can be illustrated as
I ( x , y ) = A ( x , y ) + B ( x , y ) cos ( 2 π f 0 x + ϕ ( x , y ) )
where ϕ ( x , y ) = 4 π z n z 0 h ( x , y ) λ 0 , f 0 = 2 k λ 0 .
For extracting the sideband signal in the frequency domain, Equation (5) can also be rewritten in the following form according to Euler’s formula,
I ( x , y ) = A ( x , y ) + C ( x , y ) e i 2 π f 0 x + C ( x , y ) e i 2 π f 0 x ,   with   C ( x , y ) = 1 2 B ( x , y ) e i ϕ ( x , y )
where presents the complex conjugate. Then, we apply two-dimensional Fourier transform of Equation (6) to obtain the frequency spectra, as expressed in
G   ( f x , f y ) = G A ( f x , f y ) + D 1 ( f x , f y ) G 1 ( f x , f y ) + D 1 ( f x , f y ) G 1 ( f x , f y ) ) , with   D 1 ( f x , f y ) = F [ C ( x , y ) ]
where G(fx, fy) denotes the Fourier spectra, GA(fx, fy) is the zero frequency, G−1(fx, fy) and G1(fx, fy) are the sideband signals, is the convolution symbol, and F represents the 2D Fourier transform. We extract either of the two spectra on the carrier, say D−1(fx, fy), and apply inverse Fourier transform to obtain C(x, y) Then, the spatial modulation depth and phase distribution in the {xy} plane can be calculated as
M ( x , y ) = | B ( x , y ) | = 2 | C ( x , y ) |
ϕ ( x , y ) = Angle [ C ( x , y ) ]
where M(x, y) is expressed as the spatial modulation depth, ‘Angle’ is the function to obtain the wrapped phase values within (−π, π). After that, unwrapped phase algorithms are utilized to obtain the continuous phase values, through which the surface topography can be effectively reshaped. To guarantee the accuracy of the measurement, the pixels with lower modulation should be removed from consideration in phase unwrapping. Commonly, when the measurement sample has a large enough range along the z-axis, the intensity modulation depth will decrease significantly due to the limited DOF. Therefore, an estimation of modulation depth is essential to satisfactorily retrieve the phase assignment.
Initially, two-dimensional Fourier transform and signal processing of individual interferogram are conducted to calculate both intensity modulation depth and discrete phase. Then, an effective threshold limit is reasonably determined for distinguishing the pixels that are valid and invalid for phase unwrapping, where a binary mask is logically constructed (the valid pixels are set to be one, and invalid ones are generated as zero). If the threshold is set too high, excluded areas will be enlarged. In that case, although the obtained phase will have higher reliability, it will take more time and cause fewer intersections between adjacent frames. On the contrary, if the threshold is established too low, the phase unwrapping procedure will lack reliability.
In this technique, we utilize the Diamond phase-unwrapping algorithm to calculate continuous phase, as described in [29]. The algorithm firstly identifies a seed point, which will spread to four points nearby. After that, the four points will serve as the second group of seed points. These seed points will spread again to four points nearby the second group in turn and pass through all of the effective points by a diamond path. Seed points will eventually bring about the phase unwrapping in the whole image. In the method, the sample is scanned by the PZT scanner along the z-axis, where a sequence of interferograms will be obtained at different scanning positions. When the step pace is set appropriately, a number of pixels will be valid in both adjacent interferograms. In that case, the unwrapped phase of the valid areas in two frames can be utilized for phase stitching process along the z-axis.

3. Simulation

To explain the proposed procedure specifically, the simulation has been conducted on a large-range slope. The height difference of the slope sample ranges greatly along the z-axis, as shown in Figure 2a. In the process, the central wavelength of simulated broadband light source (enveloped by Gaussian curve, full width at half maximum (FWHM) is 200 nm) is 580 nm. Additionally, a sequence of interferograms will be captured in the scanning direction, from which we select two adjacent frames as the examples. In the interference frames, it can be noted that some areas are of a high fringe-pattern contrast, while others are lowly modulated, as illustrated in Figure 2b,c. Then, 2D Fourier transform and signal processing methods are applied to obtain the spatial modulation depth in the {xy} plane for each frame, as shown in Figure 2d,e. After that, a binary mask will be generated referring to the valid pixels, whose modulation depth is above the threshold limit (the quadrangle denotes the threshold). Consequently, absolute phase values will be calculated only for the valid pixels in each interferogram, as presented in Figure 2f,g.
Figure 3a shows the modulation distribution of two adjacent interferograms, illustrating the trend of fringe-pattern contrast. Additionally, with the purpose of making it clear to explain the stitching method, a cross section indicated by a rectangle is presented in Figure 3b, where there are two modulation curves, expressed as I and II. The green section denotes the intersection of both curves, which means that these areas are commonly valid in both interferograms. Note that the phase values of intersection pixels are calculated in both adjacent frames, but the values are not same because of the different scanning positions, as presented in Figure 3c. In the stitching step, the mean values of the intersection areas are first calculated. After that, the difference can be obtained as described in
D V = Mean ( Ι ) Mean ( Ι Ι )
in which Mean(I) and Mean(II) illustrate the mean values of the intersection in Frames I and II, respectively. Lastly, in order to eliminate the phase difference between the two adjacent frames, the DV is added to the phase values of Frame II and the stitched continuous phase is obtained. In that case, the phase distribution of the valid areas in both two interferograms will be calculated, as shown in Figure 3d. It should also be noted that additional frames are required to completely recover the large-height surface topography.

4. Experimental Results

An experiment on a large-range micro-dome silica structure (micro lens) is carried out to demonstrate the advanced FT technique and the phase stitching method in SBLI. The experiment setup is developed referring to Figure 1. The system mainly contains a broadband light source (LED illuminant enveloped by Gaussian function whose central wavelength is 560 nm and whose bandwidth is 160 nm), a CCD camera that contains 768 × 576 pixels with a resolution of 6.25 × 6.5 μm, a Mirau interferoscope (50×, NA = 0.55, DOF = 900 nm) produced by the Nikon Company (Tokyo, Japan), a PZT scanning stage (with a resolution of 1.0 nm and a scanning range of 20.0 μm, produced by PI Ceramic Co., Lederhose, Germany), and self-designed imaging components.
In the experiment, nine interferograms are used to reshape the topography of the micro lens. Figure 4a–d shows four of the interferograms, which are captured with a step pace of 200 nm. Then, the spatial modulation depth distributions in the {xy} plane are obtained for each interferogram independently and the modulation threshold limit is reasonably set according to different materials (with different reflection index) of samples. In order to balance the measurement accuracy and time, the threshold will be set higher if the reflection ability is strong enough; otherwise, the threshold will be adjusted lower. In the latter case, we will try different limit values for improved experimental results. Once the threshold limit is confirmed, the absolute phase values of the valid pixels will be calculated via the FT method, as illustrated in Figure 4e–h. It should be noted that the quadrilateral in the figure is just an example of the indication line, and the boundaries of “satisfactory areas” could actually be confirmed by any other curve.
Generally, only some of the pixels in the micro lens are valid in one interferogram due to the limited DOF. Therefore, a stitching process along the z-axis is required to calculate the absolute phase of the given field of view. Based on the identification of the modulation threshold, the valid pixels would be successfully recognized, as shown in Figure 5a–d. The red color illustrates the valid pixels, of which the continuous phase will be calculated, while the blank shows the invalid ones. Consequently, the continuous phase can be stitched frame by frame like a growing tree, as shown in Figure 5e–h.
In the procedure, when all the interferograms are analyzed in sequence, the absolute phase distribution of the sample surface will be obtained, as illustrated in Figure 6a. Since the technique uses mean values of the intersection areas to carry out the stitching method, the phase of individual pixels will be slightly deviated. Here, the stitching errors are defined as
E r r o r s = Phase ( Ι Ι ) Phase ( Ι )
where Phase(I) are the phase values of the intersection pixels in frame I, and Phase(II) are the phase distributions after adding the DV to Frame II of the intersection areas. Therefore, there will be stitching errors for each independent pixel, which are the main causes of measurement errors in the proposed method. It can be seen from Figure 6b that the stitching phase errors are within the range of −0.008–0.008 rad. Accordingly, when it comes to the topography information, Equation (12) will be applied in the Mirau interference system.
h ( x , y ) = ϕ ( x , y ) 2 π * λ 0 2
From which, the stitching errors of this technique are obtained to be ±0.36 nm.
In addition, a 3D plane-fitting algorithm is applied to obtain the tilting slope in the experiment according to
F ( x , y , z ) = A x + B y + C z + D
where F(x, y, z) indicates the tilting slope, while the coefficients A, B, C, and D can be calculated via the least square method. When the redundantly tilting slope has been removed, the surface topography of the large-height micro-dome structure can be obtained, as presented in Figure 6c,d. In the experiment, we also conducted four repetitions under the same experiment environment, where the results of the positions (x = 300 pixels, y = 200 pixels and x = 400 pixels, y = 100 pixels) were obtained, as shown in Table 1. Additionally, in order to demonstrate the accuracy of this method, the results measured by a commercial 3D profiler (New View 8000, Zygo Corporation, Middlefield, CT, USA) are tabulated in Table 1. The micro-dome lens, which is produced by photolithography method, is aspherical and asymmetrical in the surface, so it is feasible to match the topography in SBLI with the results of the 3D profiler through the aspherical shape. Additionally, since the CCD camera used in SBLI is the same as the one applied in the 3D profiler, we can compare the data pixel to pixel to verify the height information.
As a result, the measurement repeatability can be estimated to be within the range of 0.1%. Meanwhile, the mean errors between the proposed method and the 3D profiler are experimentally obtained to be less than 1.5 nm. We also compare the height information of various other pixels to demonstrate that the stitching errors are almost uniform at all lateral positions, for the reason that the stitching deviations are nearly unvarying for the giving surface. It can also be noted that the height difference of the sample is nearly 3000 nm, which is significantly larger than the DOF of the experiment system. Additionally, the light fluctuations and environmental disturbances can be suppressed in some extent, because the filtering process is adopted in the spatial frequency domain and the background intensity is removed for each interferogram.

5. Discussion and Conclusions

Here, the spatial modulation depth and discrete phase are first obtained via the FT method, and the modulation threshold is then reasonably set to identify the pixels that are valid for phase unwrapping. After that, the intersection areas between two adjacent interferograms are utilized to conduct the phase stitching process. In this technique, the background noises and light fluctuations are substantially suppressed due to the signal process on each interferogram, and the moving errors are effectively eliminated because of the independent analysis of each interferogram. Importantly, the height range can be improved to multiple DOFs, which can be applied in various samples with heights from the nanometer level to the micrometer level.
In conclusion, this paper has analyzed the theoretical principle, simulation and experimental validation of the phase stitching method for dimensional metrology of microstructures. The technique could achieve an improved robustness for the measurement of large-height microstructures, the advantages of which include the elimination of stepping errors, the suppression of light fluctuations, and the freedom of a limited DOF.

Acknowledgments

The authors thank Michael Stephen for his helpful comments to improve the quality of the paper and appreciate the financial support from the Natural Science Foundation of China (Grant No. 61675206 and 61605212) and Youth Innovation Promotion Association, Chinese Academy of Sciences.

Author Contributions

Yi Zhou designed the experiments and wrote the paper; Yan Tang performed the experiments; Yong Yang and Song Hu helped to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schmit, J.; Olszak, A. High-precision shape measurement by white-light interferometry with real-time scanner error correction. Appl. Opt. 2002, 41, 5943–5950. [Google Scholar] [CrossRef] [PubMed]
  2. De Groot, P.; De Lega, X.C. Angle-resolved three-dimensional analysis of surface films by coherence scanning interferometry. Opt. Lett. 2007, 32, 1638–1640. [Google Scholar] [CrossRef] [PubMed]
  3. Mandal, R.; Coupland, J.; Leach, R.; Mansfield, D. Coherence scanning interferometry: Measurement and correction of three-dimensional transfer and point-spread characteristics. Appl. Opt. 2014, 53, 1554–1563. [Google Scholar] [CrossRef] [PubMed]
  4. Cheymol, G.; Gusarov, A.; Gaillot, S.; Destouches, C.; Caron, N. Dimensional measurements under high radiation with optical fibre sensors based on white light interferometry—Report on irradiation tests. In Proceedings of the 2013 3rd International Conference on Advancements in Nuclear Instrumentation Measurement Methods and their Applications (ANIMMA), Marseille, France, 23–27 June 2013; pp. 1–7. [Google Scholar]
  5. De Groot, P. Principles of interference microscopy for the measurement of surface topography. Adv. Opt. Photonics 2015, 7, 1–65. [Google Scholar] [CrossRef]
  6. Chen, S.; Palmer, A.W.; Grattan, K.T.V.; Meggitt, B.T. Fringe order identification in optical fibre white-light interferometry using centroid algorithm method. Electron. Lett. 1992, 28, 553–555. [Google Scholar] [CrossRef]
  7. Kitagawa, K. Simultaneous measurement of film surface topography and thickness variation using white-light interferometry. Proc. SPIE Int. Soc. Opt. Eng. 2006, 6375, 637507–637510. [Google Scholar]
  8. Lehmann, P. Systematic effects in coherence peak and phase evaluation of signals obtained with a vertical scanning white-light mirau interferometer. In Proceedings of the SPIE Photonics Europe Conference on Optical Micro- and Nanometrology in Microsystems Technology, Strasbourg, France, 3–7 April 2006; p. 618811. [Google Scholar]
  9. Ullmann, V.; Emam, S.; Manske, E. White-light interferometers with polarizing optics for length measurements with an applicable zero-point detection. Measur. Sci. Technol. 2015, 26, 084010. [Google Scholar] [CrossRef]
  10. Schmit, J.; Olszak, A.G.; Creath, K.; Schmit, J. Challenges in white-light phase-shifting interferometry. In Proceedings of the International Optical Science and Technology Symposium on Interferometry XI: Techniques and Analysis, Seattle, WA, USA, 7–11 July 2002; pp. 118–127. [Google Scholar]
  11. Li, W.; Wang, L.X.; Zheng, J.Y.; Li, M.; Zhu, N.H. Photonic generation of ultrawideband signals with large carrier frequency tunability based on an optical carrier phase-shifting method. IEEE Photonics J. 2013, 5, 5502007. [Google Scholar] [CrossRef]
  12. Shen, M.H.; Hwang, C.H.; Wang, W.C. Using higher steps phase-shifting algorithms and linear least-squares fitting in white-light scanning interferometry. Opt. Lasers Eng. 2015, 66, 165–173. [Google Scholar] [CrossRef]
  13. Terasawa, T.; Wada, A.; Ishii, Y. Multiple phase-shifting white-light interferometry for quantitative detection. In Proceedings of the Photonics Asia on Advanced Materials and Devices for Sensing and Imaging III, Beijing, China, 11–15 November 2007; p. 682928. [Google Scholar]
  14. Arieli, Y.; Epshtein, S.; Yakubov, I.; Weitzman, Y.; Locketz, G.; Harris, A. Surface measurements by white light spatial-phase-shift imaging interferometry. Opt. Express 2014, 22, 15632–15638. [Google Scholar] [CrossRef] [PubMed]
  15. Sung, Y.; Choi, W.; Fang-Yen, C.; Badizadegan, K.; Dasari, R.R.; Feld, M.S. Optical diffraction tomography for high resolution live cell imaging. Opt. Express 2009, 17, 266–277. [Google Scholar] [CrossRef] [PubMed]
  16. Takeda, M.; Ina, H.; Kobayashi, S. Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry. J. Opt. Soc. Am. 1982, 72, 156–160. [Google Scholar] [CrossRef]
  17. Qian, K. Two-dimensional windowed fourier transform for fringe pattern analysis: Principles, applications and implementations. Opt. Lasers Eng. 2007, 45, 304–317. [Google Scholar]
  18. Gao, W.; Huyen, N.T.; Loi, H.S.; Kemao, Q. Real-time 2d parallel windowed fourier transform for fringe pattern analysis using graphics processing unit. Opt. Express 2009, 17, 23147–23152. [Google Scholar] [CrossRef] [PubMed]
  19. Massig, J.H.; Heppner, J. Fringe-pattern analysis with high accuracy by use of the fourier-transform method: Theory and experimental tests. Appl. Opt. 2001, 40, 2081–2088. [Google Scholar] [CrossRef] [PubMed]
  20. Yue, H.M.; Su, X.Y.; Liu, Y.Z. Fourier transform profilometry based on composite structured light pattern. Opt. Laser Technol. 2007, 39, 1170–1175. [Google Scholar] [CrossRef]
  21. Su, X.Y.; Bally, G.V.; Vukicevic, D. Phase-stepping grating profilometry: Utilization of intensity modulation analysis in complex objects evaluation. Opt. Commun. 1993, 98, 141–150. [Google Scholar] [CrossRef]
  22. Wada, N.; Boku, S.; Yoda, T. High Resolution Optical Waveform Shaping Device Having Phase Shift Compensation Associated with Optical Intensity Modulation. U.S. Patent 8,571,418 B2, 29 October 2013. [Google Scholar]
  23. Zhou, Y.; Tang, Y.; Deng, Q.; Liu, J.; Wang, J.; Zhao, L. Dimensional metrology of smooth micro structures utilizing the spatial modulation of white-light interference fringes. Opt. Laser Technol. 2017, 93, 187–193. [Google Scholar] [CrossRef]
  24. Wippermann, F.; Bräuer, A. Fabrication of microstructures with continuous surface profiles and very large sag heights by laser lithography. In Proceedings of the SPIE MOEMS-MEM Conference on Micromachining and Microfabrication Process Technology XVII, San Francisco, CA, USA, 21–26 January 2012; pp. 824807–824811. [Google Scholar]
  25. Silvestri, S.; Schena, E. Micromachined flow sensors in biomedical applications. Micromachines 2012, 3, 225–243. [Google Scholar] [CrossRef]
  26. Khoo, H.S.; Lin, C.; Huang, S.H.; Tseng, F.G. Self-assembly in micro- and nanofluidic devices: A review of recent efforts. Micromachines 2011, 2, 17–48. [Google Scholar] [CrossRef]
  27. Kim, M.S.; Allegre, L.; Sunarjo, J.; Noell, W.; Voelkel, R. Metrology techniques for refractive microlenses and microlens array manufacturing. In Proceedings of the Microoptics Conference, Fukuoka, Japan, 25–28 October 2015; pp. 1–2. [Google Scholar]
  28. Born, M.; Wolf, E.; Hecht, E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light; Pergamon Press: Oxford, UK, 1980; pp. 30–38. [Google Scholar]
  29. Wan, W.B.; Jun-Hong, S.U.; Yang, L.H.; Jun-Qi, X.U. Phase unwrapping algorism for image processing of interferogram. J. Appl. Opt. 2011, 1, 018. [Google Scholar]
Figure 1. Schematic of Mirau correlation microscope-based scanning broadband light interferometry (SBLI) system.
Figure 1. Schematic of Mirau correlation microscope-based scanning broadband light interferometry (SBLI) system.
Micromachines 08 00319 g001
Figure 2. Process for continuous phase calculation through the modulation identification. (a) The simulated large-height slope. (b,c) Fringe patterns at two different scanning positions. (d,e) The corresponding spatial modulation of the obtained interferograms, with quadrangle indicating the threshold limit. (f,g) The continuous phase values only for valid pixels.
Figure 2. Process for continuous phase calculation through the modulation identification. (a) The simulated large-height slope. (b,c) Fringe patterns at two different scanning positions. (d,e) The corresponding spatial modulation of the obtained interferograms, with quadrangle indicating the threshold limit. (f,g) The continuous phase values only for valid pixels.
Micromachines 08 00319 g002
Figure 3. An example showing the steps of the phase stitching method. (a) The spatial modulation of two adjacent interferograms. (b) Cross section at the marked position in (a) and the identification of valid pixels. (c) The phase stitching process through intersection areas, with a red line indicating the different values of the same areas. (d) Absolute phase distribution through the stitching process.
Figure 3. An example showing the steps of the phase stitching method. (a) The spatial modulation of two adjacent interferograms. (b) Cross section at the marked position in (a) and the identification of valid pixels. (c) The phase stitching process through intersection areas, with a red line indicating the different values of the same areas. (d) Absolute phase distribution through the stitching process.
Micromachines 08 00319 g003
Figure 4. Identification based on spatial modulation. (ad) The captured interference frames at different scanning positions. (eh) Spatial modulation depth distributions with the red line indicating the satisfactory areas (whose modulation is higher than 0.3).
Figure 4. Identification based on spatial modulation. (ad) The captured interference frames at different scanning positions. (eh) Spatial modulation depth distributions with the red line indicating the satisfactory areas (whose modulation is higher than 0.3).
Micromachines 08 00319 g004
Figure 5. Phase stitching through valid areas. (ad) The valid pixels are indicated by the red color, while the invalid ones are blank. (eh) The stitched continuous phase obtained frame by frame.
Figure 5. Phase stitching through valid areas. (ad) The valid pixels are indicated by the red color, while the invalid ones are blank. (eh) The stitched continuous phase obtained frame by frame.
Micromachines 08 00319 g005
Figure 6. Experimental results and stitching errors. (a) The stitched phase values through the nine interferograms. (b) Stitching errors. (c) The reshaped topography of the micro-dome structure when redundant slope has been removed. (d) Cross-section profile at the marked position (y = 200) in (c).
Figure 6. Experimental results and stitching errors. (a) The stitched phase values through the nine interferograms. (b) Stitching errors. (c) The reshaped topography of the micro-dome structure when redundant slope has been removed. (d) Cross-section profile at the marked position (y = 200) in (c).
Micromachines 08 00319 g006
Table 1. Experimental results using the proposed method and the 3D profiler.
Table 1. Experimental results using the proposed method and the 3D profiler.
Repetitions(x = 300, y = 200)(x = 400, y = 100)
Height/nmHeight/nm
12639.22126.6
22638.12125.5
32637.92126.9
42639.92128.5
Average2638.82126.9
3D Profiler2640.22128.1

Share and Cite

MDPI and ACS Style

Zhou, Y.; Tang, Y.; Yang, Y.; Hu, S. Topography Measurement of Large-Range Microstructures through Advanced Fourier-Transform Method and Phase Stitching in Scanning Broadband Light Interferometry. Micromachines 2017, 8, 319. https://doi.org/10.3390/mi8110319

AMA Style

Zhou Y, Tang Y, Yang Y, Hu S. Topography Measurement of Large-Range Microstructures through Advanced Fourier-Transform Method and Phase Stitching in Scanning Broadband Light Interferometry. Micromachines. 2017; 8(11):319. https://doi.org/10.3390/mi8110319

Chicago/Turabian Style

Zhou, Yi, Yan Tang, Yong Yang, and Song Hu. 2017. "Topography Measurement of Large-Range Microstructures through Advanced Fourier-Transform Method and Phase Stitching in Scanning Broadband Light Interferometry" Micromachines 8, no. 11: 319. https://doi.org/10.3390/mi8110319

APA Style

Zhou, Y., Tang, Y., Yang, Y., & Hu, S. (2017). Topography Measurement of Large-Range Microstructures through Advanced Fourier-Transform Method and Phase Stitching in Scanning Broadband Light Interferometry. Micromachines, 8(11), 319. https://doi.org/10.3390/mi8110319

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop