Fringe-Order Determination Method in White-Light Phase-Shifting Interferometry for the Compensation of the Phase Delay and the Suppression of Excessive Phase Unwrapping

Journal of the Optical Society of Korea.
2013.
Oct,
17(5):
415-422

- Received : July 16, 2013
- Accepted : September 24, 2013
- Published : October 25, 2013

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

White-light phase-shifting interferometry (WLPSI) is widely recognized as a standard method to measure shapes with high resolution over a long distance. In practical applications, WLPSI, however, is associated with some degree of ambiguity of its phase, which occurs due to a phase delay, which is the offset between the phase of the fringes and the fringe envelope peak position. In this paper, a new algorithm is proposed for the determination of a fringe order suitable for samples in which the phase delay mainly occurs due to noise, diffraction and a steep angle. The concepts of the decouple factor and the connectivity are introduced and a method for calculating the decouple factor and the connectivity is developed. With the phase-unwrapping procedure which considers these values, it is demonstrated that our algorithm determines the correct fringe order. To verify the performance of the algorithm, a simulation was performed with the virtual step height under noise. Some specimens such as step height standard and a column spacer with a steep angle are also measured with a Mirau interference microscope, after which the algorithm is shown to be effective and robust.
where 𝝓
_{wlpsi}
is the phase from WLPSI, 𝝓
_{wsi}
is the phase of the fringe envelope peak position, and 𝝓
_{psi}
is the phase of the fringe relative to the position. When the phase delay, which is the offset between the phase of the fringes and the fringe envelope peak position is greater than
π
, 𝝓
_{psi}
cannot be calculated correctly due to the ambiguity of the phase. We can then only calculate
between -
π
and
π
as shown in
Fig. 1.
This means that 𝝓
_{psi}
and
can be related as follows;
Phase definition of the interference signal for a low-coherence source. The upper line is the envelope of interference signal. The black circle point is the peak of the envelope which is WSI the phase. The lower line is the interference signal. The blank circle points denote the captured signal when we use a 4 bucket or 5 bucket algorithm for calculating the phase. The left rectangle point is the phase, which is calculated with the incorrect fringe-order determination. The right rectangle point indicates the correct phase.
When we assume the fringe order
n
to be zero, the equation for WLPSI, Eq. (1), can be written as:
where
is the phase from WLPSI when the fringe order
n
is assumed to be zero.
From Eq. (3), a reasonable proposition can be made such that the fringe order
can be assigned as the same fringe order of
, when the both sides of Eq. (3) are phase-unwrapped. We can express this condition as follows;
For neighboring points,
where
m
is the integer of the fringe-order difference between two neighboring points.
Equation (4), however, cannot be met at all points due to practical measurement conditions such as noise, diffraction and s steep angle. Therefore, 𝝓
_{d}
, a decouple factor, can be introduced to address this issue:
where 𝝓
_{wsi}
and
are compensated as 𝝓
_{wsi,d}
and
, respectively by the decouple factor 𝝓
_{d}
. Then Eq. (3) can be rewritten using Eq. (5).
If we can calculate 𝝓
_{d}
which makes 𝝓
_{wsi,d}
and
satisfy Eq. (4) at all points, the unwrapping method can become more robust to determine the fringe order.
2.2. Calculation of the Decouple Factor
In this section, we propose a method to calculate the decouple factor in practical applications. From our practical application, 𝝓
_{wsi,d}
can be considered such that the difference in 𝝓
_{wsi,d}
between the neighboring points falls within the range of -π and π. In addition, the fringe-order difference
m
can be limited between -2 and 2 such that the phase delay due to undesirable conditions such as noise and diffraction are limited between -5π and 5π .
Then Eq. (4) can be rewritten considering Eq. (5) and Eq. (6).
After Eq. (5) is applied to Eq. (7),
Let △𝝓
_{d}
(
i
) be the difference in the decouple factor between the neighboring points,
Apply Eq. (9) into Eq. (8),
From Eq. (10),△𝝓
_{d}
(
i
) can be ranged as follows;
For a valid △𝝓
_{d}
(
i
) in practical application, the range of △𝝓
_{d}
(
i
) can be obtained between the maximum of the two left-hand limits and the minimum of the two right-hand limits of Eq. (11),
In order to calculate the decouple factor 𝝓
_{d}
(
i
) , an iteration procedure is used. The initial
can be assigned as the median value of the limits of Eq. (12). That is,
will be used instead of △𝝓
_{d}
(
i
) in the subsequent calculations.
For an easy explanation, we consider the pixel and edge figure shown in
Fig. 2.
The pixel figure indicates the value defined at the pixel, in this case
,
, and 𝝓
_{d}
. The edge figure indicates the difference value between the pixel values, in this case Δ𝝓
_{d}
,
A
, and
B
in Eq. (9) and (11).
Definition of a pixel and an edge. We define the edge to explain the connectivity, which is defined between neighboring pixels.
We extend the previous equations from a one dimensional case to a two dimensional case as shown in
Fig. 3.
In a two dimensional case, one point has four edges on the upper, lower, left, and right sides. The difference values,
A
(
i
),
B
(
i
), Eq. (13) and △𝝓
_{d}
(
i
) can be rewritten as follows;
Decouple factor in a two dimensional case. This figure shows how the definition of the decouple factor changes in a two dimensional case. At the pixel, the decouple factor is defined. At the edge, the difference of the decouple factor is defined.
We already know the values of 𝝓
_{wsi}
and
at all pixels. Using Eq. (14), (15), and (18), we can calculate the values of
(
i, j
) at all edges.
Next, we explain how to calculate the decouple factor 𝝓
_{d}
from the difference of the decouple factor
. In
Fig. 3
, the values of △𝝓
_{d}
in the four gray edges are written as follows;
Each equation has one 𝝓
_{d}
(
i, j
). Let 𝝓
_{d}
(
i, j
) denote the average of four 𝝓
_{d}
(
i, j
) in Eq. (20). This leads to Eq. (21).
After this strategy is applied once more to the second term in Eq. (21), Eq. (21) is rewritten as follows;
The third term in Eq. (22) has a constant
This effect is almost 1% of the total effect. This term can be ignored.
_{d}
(
i, j
) to satisfy Eq. (12). Second, this prohibits the unwrapping of the phase at the disconnected edges.
The connectivity can be given at all edges as follows;
The connectivity initially is assigned as one at all edges, and then changed to a zero under the next three conditions
where
M_{x}
is the horizontal image size and
M_{y}
is the vertical image size. When we consider connectivity, the equation for the decouple factor, Eq. (22), is changed, as follows;
_{d}
(
i
,
j
), is calculated by Eq. (24) and the connectivity at the first pixel in
Fig. 5.
If the range conditions, Eq. (16) and (17), of the difference in 𝝓
_{d}
(
i
,
j
) are not satisfied or if the connectivity of the edges is not disconnected, the connectivity at the edge where the difference in
is greatest among the four edges is changed to “disconnected”. Then, 𝝓
_{d}
(
i
,
j
) is re-calculated with the updated connectivity. The reason to check only the range condition, the Eq. (16) and (17), at the left edge and the upper edge is that checking at the right edge and the bottom edge will be processed at the next pixels. For example, when the decouple factor and the connectivity at pixel 5 was checked as shown in
Fig. 5
, we don’t check at the right edge and the bottom edge. But when checked at pixel 8 and pixel 9, these edges will be checked. This means that the decouple factor and the connectivity has been checked at all edges of pixel 5. This calculation is repeated until 𝝓
_{d}
(
i
,
j
) satisfies the condition equations, Eq. (16) and (17). This sequence is performed following the order shown in
Fig. 5.
The points which the edge connectivity changes from “connected” to “disconnected” demonstrate a major difference in
between the neighboring points. Upon the completion this process, we obtain the decouple factor 𝝓
_{d}
(
i
,
j
) at all points and the connectivity
C_{H}
(
i
,
j
),
C_{V}
(
i
,
j
) at all edges.
Flow chart to check and modify the connectivity to satisfy Equation 16 and Equation 17. Upon the completion of this flow chart, we obtain one new decouple factor and the four new connectivity at point (i, j).
The order of the flow chart in Figure 4 at two dimensional points. After completing the flow chart at all points, we obtain the new decouple factor map and the two new connectivity maps of the vertical map and the horizontal map.
, which is corrected with the decouple factor in place of the original phase 𝝓
_{psi}
. The phase-unwrapping step is not performed at the disconnected edges identified by the connectivity. During the phase unwrapping procedure, stable phase unwrapping is possible because we use the decouple factor to compensate for the phase delay and the connectivity to suppress excessive phase unwrapping.
kπ
to
kπ
was simulated as follows;
where
rand
(
i
,
j
) denotes a random signal with range of -1 to 1 and
k
is the noise level. This noise can be considered for a simulated phase delay. The analysis proceeded with and without our fringe-order determination method. We discuss the fringe-order determination and if the fringeorder determination is incorrect at some pixel, the height error is exactly a multiple of the half of mean wavelength at that pixel. The virtual step height will be perfectly recovered if we obtain the correct fringe order. So the accuracy of noise simulation was not an important point in this noise simulation. We pointed out how many points fail to determinate the fringe order under noise. The results are shown in
Fig. 6
, showing the ratio of the points where unwrapping fails on the total points at a noise level ranging from
k
= 0 to
k
= 3.
We measured the simulation-surface induced random phase delay by Equation 25 with the conventional WLPSI and the new algorithm. The lower line shows the acceptable performance of the new algorithm.
When the noise level
k
is less than 1, the results are identical between the WLPSI and our algorithm because there is no phase ambiguity. When the noise level range is from 1.0 to 2.0, our algorithm gives a better result, as shown in
Fig. 7.
When the noise level,
k
, exceeds 2
π
, the fringe order is not determined well enough under our algorithm due to the limited condition of the phase unwrapping in which the phase difference between successive points cannot exceed 2
π
on a continuous surface. When the noise level is over 2 in an actual measurement, it can be estimated that the object was measured under too much vibration or at a steep angle.
The noise simulated surface measurement at k=1.5: (a) WLPSI, (b) the proposed algorithm. Fig. 7(a) shows an instance of incorrect fringe-order determination. The noise’s height is half of the mean wavelength. In Fig. 7(b), the new algorithm corrects the incorrect fringe-order determination completely.
When the measurement height is identical to or lower than the wavelength of the light source, the measurement result can be distorted by diffraction at the sharp edge of the step height, as noted in earlier research by, for instance,
[5]
. We measured the step height standard, VLSI SHS-4500, at a height of 470.3 nm, which is quite close to the wavelength of light, with a 50X Mirau interference microscope. We intend to compare the result after we intentionally create the worst case with a large phase delay. In that case, the apparatus has a high-magnification lens and uses a Mirau-type lens. As shown in
Fig. 8(a)
, we observed that the measured result was distorted at the edge of the step height standard without our fringe-order determination method. However as shown in
Fig. 8(b)
, the measurement result was not distorted.
Profiles of VLSI SHS-4500 at a height of 470.3 nm as measured with (a) WLPSI and (b) the proposed algorism using a 50X Mirau interference microscope. At a sharp edge of the step sample, the new algorithm can correct the incorrect fringe-order determination.
As shown in
Fig. 9
, another measurement test was performed with a Mirau interference microscope. We measured the LCD column spacer sample as shown in
Fig. 9
, which has a gradient on its surface with a phase delay of
π
. In particular, on the left side of the sample, the phase delay is greater than that on the right side due to the optic misalignment of the microscope. When applying phase unwrapping without our fringe-order determination method, the measurement result is distorted at the boundary of the sample, with a phase delay of more than
π
, as shown in
Fig. 9(a)
. However, when applying with our fringe-order determination method, it was possible stably to unwrap a phase delay of more than
π
; in addition the measurement result is not distorted, leading to correct measurement result, as shown in
Fig. 9(b)
.
3D shapes and profiles of column spacer with a steep angle as measured by a) WLPSI and b) the proposed algorism using a 50X Mirau interference microscope. This result shows that the new algorithm corrects the incorrect fringe-order determination on the left side of the column spacer.
The proposed algorithm is summarized as shown in
Fig. 10.
This algorithm adds the decouple factor and the connectivity to the conventional WLPSI process which has the phaseunwrapping process. So our algorithm can be used in the conventional WLPSI which uses any phase method like the 4 bucket, 5 bucket and 7 bucket algorithm. The calculating time according to this algorithm is small. Some initial value calculation time and the flow chart calculation time of
Fig. 4
are needed additionally. The flow chart calculation of
Fig. 4
was processed one time for each pixel. So the number of total calculations is the same as the total number of pixels.
Summarization of the proposed algorithm. The proposed algorithm does not modify the conventional WLPSI algorithm.

White light interferometry
;
Fringe order determination
;
Phase unwrapping
;
Decouple factor
;
Connectivity

I. INTRODUCTION

White-light phase-shifting interferometry (WLPSI) proposed by
[2]
and
[3
,
4]
has high resolution over long distances. WLPSI, however, cannot overcome the limitation of phaseambiguity entirely, because it can also detect the fringe order incorrectly due to the problem of phase ambiguity. In order to address this issue, a specific strategy has been required to determine the correct fringe order in relation to WLPSI.
Previous research related to determining of fringe order was done by
[3
,
4]
,
[5
,
6]
and
[7]
. This work, however, yielded solutions only for limited samples.
[8]
suggested a method to determine the fringe order by phase unwrapping, though it is not sufficient in practical applications given undesirable conditions such as noise, diffraction and steep angles.
In this paper, a new algorithm that corrects errors in fringe-order determinations is proposed. Prior to the phase unwrapping process, parameters such as the decouple factor and the connectivity are introduced and calculated. The decouple factor can compensate for the effect of a phase delay. As the connectivity contains information on the edges at which phase unwrapping is prohibited, it is possible to suppress excessive phase unwrapping.
To verify the performance of the algorithm, a simulation was performed with the virtual step height under noise. Some specimens such as step height standard and a column spacer with a steep angle are also measured with a Mirau interference microscope, after which the algorithm is shown to be effective and robust.
II. METHODS

- 2.1. Introduction of the Decouple Factor

The WLPSI equation using both the phase of white-light scanning interferometry (WSI) and the phase of phase shifting interferometry (PSI) is
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.3. Introduction of the Concept of Connectivity

We introduce the concept of connectivity in order to consider the disconnected edges, serving two purposes. First, the connectivity can change the decouple factor 𝝓
PPT Slide

Lager Image

- 1. At the boundary of the measurement area: that is,i=0,i=Mx;j=0,j=My.
- At the edges of an point where the light intensity signal cannot be obtained.
- 3. If the difference inbetween the neighboring points exceed over 5π , as mentioned in relation to Eq. (7).

PPT Slide

Lager Image

- 2.4. Iterative Correction of the Decouple Factor and the Connectivity

We check and correct the decouple factor and the connectivity using the steps in
Fig. 4
and
Fig. 5
, where
Fig. 4
shows a flow chart for checking the decouple factor and changing the connectivity at each pixel and
Fig. 5
shows when using this flow chart. First 𝝓
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 2.5. Phase Unwrap

The remaining final step involves simply unwrapping the phase. In this paper, we use
[10]
a phase unwrapping method, which is modified to consider the decouple factor and the connectivity. A modified version of Herráez’s method uses a phase of
PPT Slide

Lager Image

III. EXPERIMENTS

To examine the robustness to noise of our algorithm, a simulation was performed using the following procedure. The virtual step height with white noise in a range of -
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

IV. CONCLUSION & DISCUSSION

In this paper, a new algorithm is proposed for the fringe-order determination in the white light phase-shifting interferometry to compensate for the phase delay and to suppress excessive phase unwrapping to overcome the phase ambiguity problem.
- (1)A decouple factor is proposed to compensate for the phase delay. When the phase is compensated for by the decouple factor prior to unwrapping, the WLPSI result can be successfully unwrapped, using the proposed Eq. (24) for calculating the decouple factor at all points.
- (2)The concept of connectivity is proposed to suppress excessive phase unwrapping. The connectivity is defined in the vertical and the horizontal directions, and proper initial connectivity values are devised beforehand as explained in Section 2.3. This process allows the determination of the edges at which phase unwrapping is prohibited.
- (3)The decouple factor and the connectivity are successfully calculated and corrected according to the proposed flow chart shown inFig. 4.
- (4)We performed a simulation under noise to verify the proposed method. The result showed that our method can reduce the errors induced by diffraction at a low step height close to the wavelength of light. We also measured a sample with a steep angle, with the result demonstrating that our method can successfully perform phase unwrapping for a phase delay of more thanπ.
- (5)The proposed algorithm is found to be very robust against noise, diffraction, and the angle of the surface through measurement tests using white-light phaseshifting interferometry.

PPT Slide

Lager Image

Kino G. S.
,
Chim S. S. C.
1990
“Mirau correlation microscope,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.29.003775
29
3775 -
3783

Larkin K. G.
1996
“Efficient nonlinear algorithm for envelope detection in white light interferometry,”
J. Opt. Soc. Am.
A 13
832 -
843

Sandoz P.
1996
“An algorithm for profilometry by white-light phase-shifting interferometry,”
J. Mod. Opt.
43
1545 -
1554

Sandoz P.
,
Devillers R.
,
Plata A.
1997
“Unambiguous profilometry by fringe-order identification in white-light phaseshifting interferometry,”
J. Mod. Opt.
http://dx.doi.org/10.1080/09500349708232918
44
519 -
534

Harasaki A.
,
Wyant J. C.
2000
“Fringe modulation skewing effect in white-light vertical scanning interferometry,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.39.002101
39
2101 -
2106

Harasaki A.
,
Schmit J.
,
Wyant J. C.
2000
“Improved verticalscanning interferometry,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.39.002107
39
2107 -
2115

Schmit J.
,
Olszak A.
2002
“High-precision shape measurement by white-light interferometry with real-time scanner error correction,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.41.005943
41
5943 -
5950

de Groot P.
,
de Lega X. C.
,
Kramer J.
,
Turzhitsky M.
2002
“Determination of fringe order in white-light interference microscopy,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.41.004571
41
4571 -
4578

Itoh K.
1982
“Analysis of the phase unwrapping problem,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.21.002470
21
2470 -
2470

Herráez M. A.
,
Burton D. R.
,
Lalor M. J.
,
Gdeisat M. A.
2002
“Fast two-dimensional phase-unwrapping algorithm based on sorting by reliability following a noncontinuous path,”
Appl. Opt.
http://dx.doi.org/10.1364/AO.41.007437
41
7437 -
7444

Kim J.-H.
,
Yoon S.-W.
,
Lee J.-H.
,
Ahn W.-J.
,
Pahk H.-J.
2008
“New algorithm of white-light phase-shifting interferometry pursing higher repeatability by using numerical phase error correction schemes of pre-processor, main processor, and post-processor,”
Optics and Lasers in Engineering
http://dx.doi.org/10.1016/j.optlaseng.2007.08.008
46
140 -
148

Citing 'Fringe-Order Determination Method in White-Light Phase-Shifting Interferometry for the Compensation of the Phase Delay and the Suppression of Excessive Phase Unwrapping
'

@article{ E1OSAB_2013_v17n5_415}
,title={Fringe-Order Determination Method in White-Light Phase-Shifting Interferometry for the Compensation of the Phase Delay and the Suppression of Excessive Phase Unwrapping}
,volume={5}
, url={http://dx.doi.org/10.3807/JOSK.2013.17.5.415}, DOI={10.3807/JOSK.2013.17.5.415}
, number= {5}
, journal={Journal of the Optical Society of Korea}
, publisher={Optical Society of Korea}
, author={Kim, SeongRyong
and
Kim, ungHwan
and
Pahk, HeuiJae}
, year={2013}
, month={Oct}