High Precision Signal Processing Algorithm for White Light Interferometry

A new signal processing algorithm for absolute temperature measurement using white light interferometry has been proposed and investigated theoretically. The proposed algorithm determines the phase delay of an interferometer with very high precision (≪ one fringe) by identifying the zero order fringe peak of cross-correlation of two fringe scans of white light interferometer. The algorithm features cross-correlation of interferometer fringe scans, hypothesis testing and fine tuning. The hypothesis test looks for a zero order fringe peak candidate about which the cross-correlation is symmetric minimizing the uncertainty of mis-identification. Fine tuning provides the proposed algorithm with high precision sub-sample resolution phase delay estimation capability. The shot noise limited performance of the proposed algorithm has been analyzed using computer simulations. Root-mean-square (RMS) phase error of the estimated zero order fringe peak has been calculated for the changes of three different parameters (SNR, fringe scan sample rate, coherence length of light source). Computer simulations showed that the proposed signal processing algorithm identified the zero order fringe peak with a miss rate of 3 × 10-4 at 31 dB SNR and the extrapolated miss rate at 35 dB was 3 × 10-8. Also, at 35 dB SNR, RMS phase error less than 10-3 fringe was obtained. The proposed signal processing algorithm uses a software approach that is potentially inexpensive, simple and fast.


Introduction
Although fiber optic interferometric sensors offer the possibility of performing measurements with very high sensitivity and resolution [1], they suffer from problems such as complex signal processing techniques, quadrature point stabilization, and uncertainty as to whether an increase or decrease in the value of measurands has occurred. In order to fully utilize the capability of fiber optic sensors, a different sensing principle termed as "white light interferometer" or "white light interferometry (WLI)" was developed [2]. From the beginning of this research, an all fiber white light interferometry (AFWLI) absolute temperature measurement system as shown in Figure 1 was selected as the application area of the proposed signal processing algorithm. White light interferometry departs from the conventional interferometry in that it uses a broadband light source. SLD in Figure 1 represents a superluminescent diode (SLD) used as a broadband light source and PD1, PD2 denote photodetectors 1 and 2, respectively. In the AFWLI two fiber Fabry-Perot interferometers (FFPI, sensing FFPI and reference FFPI in Figure 1) and a Mach-Zehnder Interferometer as a processing interferometer (scanning interferometer in Figure 1, hereafter termed as MZI), which has two piezoelectric transducers (hereafter termed as PZT) in its two arms, are connected in tandem. The sensing and the reference interferometer output signals from PD1 (Equation 1 for PD1) and PD2 (Equation 2 for PD2) are given respectively by: where L C is the coherence length of light source, Φ P is the OPD (Optical Path-length Difference) between two arms of Mach-Zehnder processing interferometer and Φ S , Φ R are the round trip phase shifts for the sensing and reference FFPIs, respectively. In Equation 1 and Equation 2 it is assumed that light source has a Gaussian power spectrum. A constant d.c. voltage V DC (100~150 volt, Figure 1) is applied to the PZT1 in one arm to coarsely match the OPD of the MZI to that of the sensing FFPI. And an alternating ramp voltage V SAW (Figure 1) is applied to PZT2 in the other arm to scan the processing interferometer so that OPD Φ P between two arms of MZI can be varied over a certain range. This AFWLI for temperature measurement produces two fringe scans, one from the sensing FFPI and another one from the reference FFPI, as shown Figure 2. In AFWLI, the sensing FFPI is exposed to the temperature T S to be measured and the reference FFPI is protected from environmental disturbances but exposed to the known reference temperature T R . Now, assume that the known temperature of the sensing FFPI and the reference FFPI are T S and T R , respectively. When the phase Φ P of MZI is scanned and exactly matched to that of the sensing FFPI (the reference FFPI), a zero order fringe peak of the sensing FFPI (the reference FFPI) is produced at certain Φ P = Φ P,S (Φ P = Φ P,R ), as shown in Figure 2. In Figure 2, Φ P,S (Φ P,R ) denotes the phase of the processing interferometer producing the zero order fringe peak of the sensing FFPI (the reference FFPI). If we can identify the phase difference Φ d = Φ P,S -Φ P,R (this is possible by the proposed signal processing algorithm which will be shown later in this article) then Φ d is mapped to the temperature T S and absolute temperature measurement is possible. This problem has been known as "Time Delay Estimation (TDE) [3] (or Phase Delay Estimation)". In this article, a new signal processing algorithm to estimate the phase delay Φ d of AFWLI is proposed. This article consists of five sections. Section 2 describes the previous related works for time delay estimation methods. In Section 3, a signal processing algorithm to measure the phase difference Φ d of AFWLI is proposed. In Section 4 the performance of the proposed signal processing algorithm has been demonstrated using computer simulations. Section 5 shows the comparison to the previous literature and Section 6 is the conclusion of this article.

Previous works
Two major classes of signal processing algorithms for WLI are the hardware approach and the software approach. Both approaches have a more or less "tracking zero order fringe peak" feature. Gerges proposed a hardware approach which locks the zero order fringe position of interferometer by a feedback loop [4]. An improvement of the sensitivity up to 1/240 fringe was claimed. To the author's best knowledge this method, while still dependent on the incremental characteristic of laser interferometry and not fully taking advantage of WLI's potential to identify the interference fringe, demonstrated the feasibility of locking and tracking the fringe peak for absolute measurement for the first time.
There are many software algorithms to estimate phase delay Φ d [3,5]. Among many detection methods, the cross-correlation method dominates the field of phase delay estimation in practice due to its easier implementation [6]. Many other TDE methods are based on this algorithm. The crosscorrelation method cross-correlates the two fringe signals i S (n) and i R (n), which are sampled versions of I S (Φ P ) and I R (Φ P ) respectively, into i (n) and considers the sample number argument n=n d that corresponds to the maximum peak in cross-correlation i (n) as the estimated time delay [7].
While WLI has the potential to identify the interference fringe order from the output pattern of an interferometer [8], it is difficult to distinguish the zero order fringe peak from its adjacent first order fringe peaks when noise is present in the interferometer output, as shown in Figure 3. From Equations 1 and 2, the amplitude difference ∆I between the zero order fringe peak and the first order fringe peaks is represented as: Clearly, if a system has a noise level which is equal or greater than ∆I the zero order fringe peak cannot be identified directly, simply through inspection of its amplitude. If the normalized zero order fringe peak value is defined as unit signal, a minimum signal-to-noise ratio SNR min required to identify the zero order fringe peak [9] is given by Representative values of the coherence length of different light sources like white light lamp, light emitting diode (LED), superluminescent diode (SLD), are about 10, 20, and 40 in terms of optical fringes. The SNR min required to identify the zero order fringe peak by amplitude difference ∆I is given from Equation 3 as 28 dB, 40 dB, and 52 dB respectively [9]. One disadvantage related to crosscorrelation is the broadening of envelope (or peak) (hereafter termed as envelope) from L C to C L 2 due to the cross-correlation of two almost identical Gaussian envelope signals [17]. This results in higher SNR min in Equation 3 requiring 6 dB more than before cross correlation. Then, with this broadening, the attainable resolution is often not better than one fringe and a precise peak location may be somewhat questionable. This difficulty has inhibited the application of fiber optic sensors using WLI, for example, absolute OPD determination [9].  Coherence length L C  8  Additionally, once the zero order fringe peak is identified, then for a more accurate sub-sample resolution time delay estimation we will have to use interpolation which is possible by either quadratic interpolation in time domain or frequency domain zero-padding [10,11]. But quadratic interpolation uses three cross-correlation values centered at the estimated peak of cross-correlation. This method has a shortcoming of bias caused by time domain sample rate [12,13] and the difference between true peak shape and the fitted quadratic polynomial. In frequency domain zero-padding, the number of zeros are padded in the middle of Fourier Transform of cross-correlation i(n). Notice that zero-padding in frequency domain increases the discrete frequency by a certain factor which eventually results in time-domain interpolation in cross-correlation signal i(n). Zero-padding in frequency domain is a useful tool to improve the peak location accuracy, but it increases the computational complexity [14] and storage requirement associated with inverse Fast Fourier Transform (FFT) operations [15].
Notwithstanding the above mentioned shortcomings, cross-correlation is a still useful tool for time delay estimation as shown that cross-correlation with no pre-filtering is an optimal maximum likelihood estimator to estimate the time delay between two similar signals if the noise processes w R (n), w S (n) of signal i R (n), i S (n) are white noises and at least one of signals has high signal-to-noise ratio greater than 10 dB [16].
The outcome of the time delay estimation depends on the combined performance of coarse estimation (zero order fringe peak identification) and sub-sample resolution estimation of time delay. In this article, a new signal processing algorithm which can accurately identify the zero order fringe peak of cross-correlation i (n) of two fringe scan output signals of a WLI is proposed. This algorithm still uses a cross-correlation technique taking advantage of simple implementation. But this algorithm combines the hypothesis test as a coarse estimation to reduce the possibility of mis-identification of peak with fine tuning algorithm as a sub-sample resolution peak estimation to overcome the shortcomings of quadratic interpolation or frequency domain zero-padding.
The proposed signal processing algorithm uses a software approach, which is potentially inexpensive, simple and fast. And, this proposed signal processing algorithm has a low peak misidentification rate of 3 x 10 -4 at 31 dB SNR and has a high precision fine tuning capability down to 5 x 10 -4 fringe as will be shown from the computer simulation results.

Proposed Signal Processing Algorithm
The proposed signal processing algorithm consists of five steps applied to sampled signal of WLI fringe scans. They are: 1) Normalization and cross-correlation 2) Peak and zero crossing detection 3) Matched filtering 4) Hypothesis test 5) Fine Tuning Each procedure is explained briefly below.

Normalization and cross-correlation
As a preliminary procedure, the output of photodetector signals I S (Φ P ) and I R (Φ P ) are sampled and normalized respectively into fringe scan i S (n) and i R (n): and: where n S is zero order fringe peak sample number for sensing FFPI, n R is zero order fringe peak sample number for reference FFPI, f S is the sample rate in samples per fringe (or samples/fringe). A normalization procedure is carried out by removing the d.c. component (constant "1" in Equation 1 and Equation 2) of each fringe scan. w S (n) and w R (n) are the white noise related to the fringe scan i S (n) and i R (n) respectively with zero mean and variance 2 w  . The phase delay between the sensing FFPI and the reference FFPI is defined in terms of samples as: After normalization, i S (n) and i R (n) are cross-correlated into i(n) and normalized again. Crosscorrelation i(n) can be expressed in mathematical form with its zero order fringe peak p 0 at n=n d as shown in Equation 7.
In Equation 7 the effective coherence length L C,eff is given as is a white noise of ) (n i . Exponential term in Equation 7 is termed as "the envelope" of cross-correlation i(n). This cross-correlation improves SNR at zero order fringe peak up to 14 dB [17]. One disadvantage related to cross-correlation is the broadening of peak from L C to C L 2 due to the cross-correlation of two (almost) identical Gaussian envelope signals [17]. For example, this results in higher SNR min in Equation 6 requiring 6 [dB] more than before cross-correlation when L C =26λ, but this will be compensated by 14 dB SNR improvement at the zero order fringe peak [17]. After cross-correlation, the task of the signal processing algorithm becomes to find a zero order fringe peak p 0 =n d (global maximum) of cross-correlation i (n) correctly.

Peak and Zero-crossing Detection
At this stage, all the peaks p i , all the minima q i and all the zero crossings z i of the cross-correlation i(n) are detected and registered. Every peak is labeled as p i where its value p 0 =n d is zero order fringe peak position, p -1 is negative first order fringe peak position, p 1 is positive first order fringe peak position in terms of sample number and so on. Also every zero crossing between peak p i and minimum q i is detected and labeled as z j where j=1,2,3…in terms of sample number. Then zero crossing position z j with sub-sample resolution is calculated [9] by using interpolation formula as shown in Equation 8, Equation 9 and Figure 4.
Linear interpolation is the reasonable method because the cosine function crosses the zero essentially as a straight line as shown in Figure 4.   . .
Then the zero crossing period b is calculated by fitting all the z j into the form of: by applying least square fitting to the distance between all pairs of zero (z j and z j+1 ) as shown Figure 5 where N Z is the total number of zero crossings in i(n). In Figure 5 the slope b of the y (n) is the estimated number of samples inside half fringe (or half the sample rate) and the estimated sample rate is given by: where [2b] is the closest integer number of 2b.

Matched Filtering
A matched filter is the optimum filter to maximize the SNR of signal out of the matched filter in the presence of additive stochastic noise at input signal [18]. Consider the case where the input signal is s(t) and n(t) is the white noise with zero mean, variance 2 N  related to s(t). In this case, matched filter theory states that the maximum SNR at the output will occur when the matched filter has an impulse response h(t)=s(t 0 -t) that is equal to the time-reversed version of the signal waveform s(t) to which matched filter is matched and time-delayed by certain t 0 [second]. When input signal s(t) is timelimited signal existing only when 0< t <T then the matched filter is defined by [18]: Note that the output s 0 (t) of the matched filter, which is the convolution of s(t) and h(t)=s(T-t), is given by: . But s(-t) and h(t) are identical two waveforms with time distance T separated and s 0 (t)=R s(-t)h(t) (t) becomes R ss (t-T) when R ss (t) denotes the autocorrelation R ss (t) of input signal s(t). Then s 0 (t) has its maximum value when t=T and this maximum value is same as the R ss (0) of autocorrelation R ss (t) of input signal s(t).
Thus if we process a signal-plus-noise with a matched filter, the largest peak outputs due to the signal will correspond to R ss (0) and the SNR [18] of the largest peak output of matched filter is given by: If we define the SNR of input signal s(t) as the ratio between the largest value  a of input signal s(t) and σ N , then SNR improvement due to the matched filtering is given as: At this stage, following the concept of matched filter, only one fringe of cross-correlation signal i (n) (i.e. signal between i(p i ) and i(p i+1 ) of i (n)) is considered as the signal to be matched and one period of cosine function, Then, for a given one fringe of cross-correlation signal i(n), signal out of the matched filter is maximized when the time delay between two signals is S fˆ and maximum value is calculated as in for given J i 's (i=0, ±1, ±2,…). The envelope function of i(n) is slowly changing in the vicinity of the zero order fringe peak of i (n) and Equation 15 can be rewritten as: where h i in Equation 16 is the average value of exponential envelope function between p i and p i+1 as shown in Equation 16b.
For the case of a zero order fringe peak i(p 0 ) of i(n), average value h 0 of envelope between p 0 and p 1 is approximately "1" because i (p 0 ) =1 and i (p 0 ) ≈ i (p 1 ). Then, Equation 13 becomes: and SNR of i(p 0 ) is improved by the factor of: In the above Equation 16a any J i value is the weighted integration of one fringe between two peaks of i(n) and also can be considered as maximum value out of matched filtering which has a improved SNR over i(p i ) as shown Equation 13. So, J i values will be used in the hypothesis test instead of peak value i(p i ) of i(n) as will be shown later (see Section 3.4.). One useful property of J i is J 0 = J -1 , J 1 = J -2 … (J i = J -i-1 ) due to the even symmetric property of i (n) about n=p 0 = n d .

Hypothesis Test
In hypothesis test, signal processing algorithm chooses the nine biggest peaks of i (n) as zero order fringe peak candidates. Ideally, p j 's of j=0, ±1, ±2, ±3, ±4 are selected. Hypothesis test presumes that each candidate peak is the zero order fringe peak and calculates the parameter g(p j ): where j of notation d(j,i) corresponds to the candidate peak p j on the hypothesis test and i is the distance from the zero order fringe peak candidate in terms of fringe. Then g(p 0 ) is expressed as: Ideally all the values of d(j=0,i) for the zero order fringe peak p 0 is zero (Figure 7) due to the symmetric property of J i = J -i-1 and the zero order fringe peak candidate producing│g(p j )│=0 is announced as the estimated zero order fringe peak  p 0 . But, practically the zero order fringe peak candidate producing minimum │g(p j )│ is announced as the estimated zero order fringe peak  p 0 . Note that ideally zero order fringe peak p 0 happens at n= n d , first positive order fringe peak p 1 at n= n d + f S and so on.  i Monotonically increasing or decreasing d (j,i) Figure 8 shows the vicinity of zero order fringe peak of cross-correlation i(n). From Figure 8 it is clear that the discrete sample zero order fringe peak position p d of the cross-correlation i(n) identified by the hypothesis test is not necessarily same as the "true zero order fringe peak position", p t and the goal of fine tuning algorithm is to find the distance between p d and p t . First, fine tuning algorithm assumes that the cross-correlation i (n) is represented as:

Fine Tuning
Note that p t in Equation 20 is not necessarily a integer sample number. Also, fine tuning algorithm generates test cross-correlation, i test (n) defined by: In Equation 21 it is assumed that the discrete sample zero order fringe peak position p d happens ideally at n= n d and also that the approximate value of the coherence length of the light source, C L is known priori to us. This is possible by counting the number of fringes within full width at which the intensity of interferogram using a particular light source decreases down to e -1 of its maximum value. And value of S fˆ is given from as shown in Section 3.2.
Then, the distance between p d and p t can be found by calculating J(M) for varying integer number M: where i test (n-MΔn) is a computer-generated test cross-correlation with its zero order fringe peak at p d +MΔn and Δn is the desired fine tuning resolution in sample.
In Equation 23 N sub is the number of sub-divisions in one fringe (the desired fine tuning resolution in fringe). Then J(M) becomes: Interpretation of Equation 24 is that fine tuning algorithm generates its own test cross-correlation, i test (n) with the true zero order fringe peak positioned at n d +MΔn and calculate J(M) for varying M. Then J(M) has a peak value at the certain M f which makes M f Δn ≈ p t -n d because the best similarity between i (n) and i test (n) is attained when n d +M f Δn ≈ p t . Then, fine tuning algorithm announces (n d +M f Δn) as the estimated true zero order fringe peak location t p .
Calculating J(M) for the range of -∞ ≤ M ≤ + ∞ takes many calculations and is time-consuming. The true zero order fringe peak p t normally is within several samples of the discrete sample zero order fringe peak p d . Especially, when SNR of the fringe scans is high enough, the true zero order fringe peak p t is ideally within the half sample of the discrete zero order fringe peak p d . So, for the above case  Additionally it must be emphasized that the principle of fine tuning in Equation (24) and Equation (25) can be used for matched filtering in order to obtain the best similarity between the signal inside the one fringe of cross-correlation (signal between any peak i(p i ) and i(p i+1 )) and the matched filter i M (n) and to maximize the SNR improvement at peaks by matching the phase between the discrete sample i-th order fringe peak i(p i ) and the very first value of the matched filer i M (n).

Computer Simulations
The proposed signal processing algorithm was verified using computer simulations. To see the shotnoise limited performance of the proposed signal processing algorithm, the normalized AFWLI fringe scans, i S (n) and i R (n) were computer-generated using Equation 4 and Equation 5 respectively and shot noise was added to the AFWLI fringe scans, i S (n) and i R (n) instead of white noise. In the computer simulations the zero order fringe peak position p 0,S (and p 0,R ) for i S (n) (and i R (n)) were randomly selected as real number. i S (n) and i R (n) were cross-correlated into i(n). Then, p t is calculated as (p 0,Sp 0,R ) and zero order fringe peak p d = n d of i(n) becomes the integer part of p t . The coherence length of i S (n) and i R (n) were chosen as L C =26λ to simulate the coherence length of the commercially available SLD like OKI OE350S from Oki semiconductor. For all the computer simulations presented in this article parameters are fixed as follows, unless the specified parameter is varied for a certain range and circumflex notation  on the top of the parameter denotes the estimated value of that particular parameter calculated by the proposed signal processing algorithm.

1. Simulation: Miss rate simulation
In the first simulation miss rate (misidentification rate) of the proposed signal processing algorithm was tested at different shot noise levels. The SNR was varied from 1 dB to 28 dB with 1 dB separation and a set of 10000 simulations was carried out at different SNR. When the position difference between the estimated zero order fringe peak d p and the computer-generated zero order fringe peak p d is bigger than half fringe (8 samples), the zero order fringe peak is considered to have been misidentified. Table 2 shows the miss rate of the proposed signal processing algorithm and miss rate is defined as the ratio of number of miss to 10000. In Figure 9, miss rate was plotted along with Bit Error Rate (BER) of binary fiber optic communication system [19]. In Figure 9 (or Table 2) we have a rule of thumb that every dB improvement in SNR (over the range of 26 ~ 31 dB) produces approximately one order of magnitude improvement in error rate. This kind of behavior is also the case for binary fiber optic communication system (two orders of magnitude improvement in error rate for the binary fiber optic communication system).
To extrapolate the miss rate beyond the range of 10 -4 on the BER curve, data points in abscissa in Figure 9 were left-shifted by ~14.2 dB (by trial and error) until data points of miss rate between 16 dB ~ 31 dB were visually fitted on the BER curve. Then four more data points were extrapolated beyond the data point of miss rate at 31 dB SNR on the BER curve as shown in Figure 10. It is predicted that miss rate will be 3 x 10 -8 at SNR of 35 dB.  Figure 10. Extrapolation of miss rate on BER curve.

Simulation: Root mean square (RMS) error of the zero order fringe peak identification
After the zero order fringe peak was identified, fine tuning was calculated for resolution enhancement. Phase error Φ error,i between computer generated zero order fringe peak p t (or p 0 ) and fine-tuned zero order fringe peak t p (or 0 p ) was calculated. Phase error Φ error,i was averaged over 30 simulations at a given SNR and this average gives out the root mean square (RMS) error of the zero order fringe peak identification.   Figure 11 shows the change of RMS error along with the SNR in the range of 10 dB to 40 dB and Figure 12 shows the change of RMS error along with the SNR in the range of 30 dB to 40 dB. It is shown that the minimum SNR required to achieve RMS error less than 10 -3 [fringe] (which is the fine tuning step size) must be greater than 35 dB SNR.

Simulation: High Precision Performance of Fine Tuning
In this simulation, the SNR was fixed at 30 dB and RMS error was compared between two cases. One is the case where the estimated zero order fringe peak  p 0 was fine-tuned for resolution enhancement and the other case was where the estimated zero order fringe peak  p 0 was not fine-tuned.  Figure 13, without fine tuning, RMS error of  p 0 is approximately 0.3/f S fringe and RMS error of  p 0 totally depends on the sample rate. But this RMS error was reduced down to ~0.0015 fringe when  p 0 was fine tuned. Note that RMS error was not sensitive to the sample rate f S over the range between 10~20 samples/fringe. This feature predicts that signal processing algorithm can use lower sample rate f S for the faster signal processing, but still maintain lower RMS error by using fine tuning, which is comparable to the error of higher sample rate f S . In the proposed algorithm, the crosscorrelation between i S (n) and i R (n) is most time-consuming part of signal processing algorithm, but the sample rate f S can be lowered to speed up the signal processing without losing higher sample rate RMS error. One strong benefit of the proposed algorithm is that once zero order fringe is identified correctly then the resolution of the fine tuning can be adjusted as low as desired simply by changing the value of N sub and high precision phase delay estimation can be obtained.

Simulation: Effects of b ( S fˆ)
In this simulation, estimated effective coherence length eff C L , was set as 36 fringes, SNR=30 dB and estimation error of the zero crossing period b was varied from -10 % to 10 % when b=10 (f S =20 samples/fringe) to see the effect of estimation error of the zero crossing period b on the performance of the proposed signal processing algorithm. The estimation error of the zero crossing period b is defined by: where b is the half the sample rate f S . In Figure 14, it can be shown that RMS error of fine tuning was not sensitive to the estimation error of b within the range of 6  % estimation error of b and especially, RMS error of 0.001 [fringe] was obtained within the range of 2% estimation error of b . But, on the range of 7  % ~ 10  % estimation error b RMS error increased dramatically to 0.25 fringe. This is presumably due to the fact that i(n) and test cross-correlation i test (n) (Equation 24 or Equation 25) started to be out of phase and RMS error increased fast.

Simulation: Effects of C L
This simulation is to show the effect of estimated coherence length C L on RMS error assuming that b b  , L C =26λ(L C,eff ≈ 36λ) and SNR was fixed as 30 dB. Estimated coherence length C L was varied from 21λ to 30λ . In other words, estimated effective coherence length eff C L , was varied from 30λ to 42λ and RMS error was observed. Figure 15 shows that RMS error turned out to be not sensitive to the estimation error in coherence length . This is due to the property of cross-correlation. Crosscorrelation is maximized when p t and p d +MΔn are in phase as long as both cross-correlation i (n) and test cross-correlation i test (n) are symmetric.

Comparison to the literature
In this section a comparison to the previous literature data regarding the resolution of zero order fringe peak detection is given. Interestingly enough, reference [20] calculated the theoretical limit of scanning white light interferometry signal evaluation algorithm. In this reference the theoretical limit of resolution of the fringe order detection was given as: where S i =S(z i ) consists of a fixed number of intensity values, typically taken at equidistant positions z i , I(z) is the ideal input signal (i.e. fringe scan) to be subject to intensity noise N(z), z c of I(z) is the higher SNR range (30~40 [dB]). This is not a surprising result because in higher SNR fine tuning algorithm will enhance the fringe order detection resolution down to the theoretical limit once zero order fringe peak is identified correctly. But, over low SNR range, fine tuning algorithm is not effective in enhancing the fringe order detection resolution because probably zero order fringe peak is misidentified and fine tuning is searched within the half sample of the misidentified zero order fringe peak.  Increasing the fine tuning range will help to locate the zero order fringe peak correctly and lower the RMS error down to the theoretical limit, although this is time consuming. Optimum fine tuning range of M over low SNR range to find the correct zero order fringe peak will be the subject of further research.

Conclusions
A new signal processing algorithm for white light interferometry has been proposed. The goal of signal processing algorithm was to find the time delay (phase shift) between two fringe scans which makes it possible to measure the absolute optical path length of a sensing interferometer. This new signal processing algorithm can be used for absolute temperature measure measurement by mapping the zero order fringe peak position of cross-correlation i(n) to the time delay between two fringe scan. Cross-correlation between two fringe scans utilizes all the photons in fringe scans effectively and the uncertainty of the zero order fringe peak mis-identification was reduced. Monte-Carlo simulations showed that the proposed signal processing algorithm identified the zero order fringe peak with a miss rate of 3x10 -4 at 31 dB shot noise and the extrapolated miss rate at 35 dB shot noise was 3x10 -8 . Also resolution of less than 10 -3 fringe was obtained at 35 dB shot noise ( Figure 12). The fine tuning of signal processing algorithm requires some prior knowledge on the coherence length L C of the SLD and sample rate f S . But the performance of proposed signal processing algorithm turned out to be not sensitive to the estimation error in C L and S fˆ (within the range from -6% to 6% estimation error of b ). Also the proposed signal processing reached the theoretical limit of fringe order detection over the higher SNR range (30~40 [dB]). The proposed signal processing algorithm uses a software approach which is potentially inexpensive, simple and fast. As a whole, the proposed signal processing algorithm has proven to be a high precision signal processing algorithm for AFWLI phase (time) delay estimation.
In this appendix, SNR improvement at the peak of cross-correlation i(t) is derived. First, i S (Φ) and i R (Φ) are given respectively as: Here L P , L S , L R are the path length difference of the Mach-Zehnder processing interferometer, sensing FFPI (twice the Fabry-Perot cavity length), reference FFPI and n is the effective refractive index of core of the optical fiber (not to be confused with sample number n in i(n)). Notation    , and are reserved to denote the cross-correlation, convolution and Fourier transform relationship.
When L P is scanned L P is a function of time t and also Φ P is a function of time t. Then: where f (or w 0 ) is a modulation frequency (or angular velocity) of fringe scan. Following the same token, Φ S , Φ R are given by: -5) and the relation between t S and t R is given by: where t d is the time delay between sensing FFPI fringe scan and reference FFPI fringe scan. Using the following substitutions: and considering the noise in i S (Φ P ) and i R (Φ P ), time domain representation of the Equation A-1 and Equation A-2 can be rewritten as: In Equation A-10 and A-11 we used a substitution a = (w 0 /πN f ) 2 and it is assumed that Φ R =0 without loss of generality and that w S (t), w R (t) are white noise with zero mean and variance 2 i R (t) respectively and independent each other. Also, w 0 was set to "1" for normalization in order to make i S (t) = e -1 at t = (N f /2) x 2π.
The first part of Appendix A shows that the increase of coherence length from L C to C L 2 due to cross-correlation. We know that cross-correlation and convolution are identical operations for the case where either i S (t) or i R (t) is an even symmetric function. Then, the cross-correlation i(t) of i S (t) and i R (t) is given as follows:  And the coherence length of the cross-correlation(effective coherence length) has increased by a factor of 2 . In other words, The second part of appendix A shows that 14 dB SNR improvement is obtained at zero order fringe peak after cross-correlation. Peak value (zero-order fringe peak) of cross-correlation i(t) happens when t = t d and its value i S * i R (t=t d ) is calculated by substituting t = t d in Equation Then SNR at the peak of cross-correlation i (t) is given as the ratio between magnitude of peak signal to the standard deviation of noise as shown in Equation A-25.