Abstract:
An error correction for effective-wavelength variations is implemented by adjusting the filter parameters of the quadrature demodulation algorithm of a high definition vertical scanning process using a phase step that accounts for phase-step changes associated with variations in the effective wavelength irradiating the sample when the surface is curved. The nominal phase step is replaced in the filter with an actual phase step size that includes a phase parameter generated for each pixel by calibration or modeling of a specific type of surface. This substitution eliminates all errors produced by surface-dependent variations in the effective wavelength of the irradiating light.

Description:
RELATED APPLICATIONS 
     This is a continuation-in-part application of U.S. application Ser. No. 12/082,846, filed Apr. 15, 2008. 
    
    
     BACKGROUND OF THE INVENTION 
     1. Field of the Invention 
     This invention relates in general to vertical-scanning interferometry (VSI) for surface characterization. In particular, it relates to an error-correction method for the high-definition vertical-scan interferometric (HDVSI) procedure described in Ser. No. 11/473,447, herein incorporated by reference. 
     2. Description of the Related Art 
     The well known techniques generally classified as phase-shifting interferometry (PSI) and vertical scanning interferometry (VSI) make it possible to measure the profile of most samples. However, they do not allow low-noise measurement of samples that combine smooth surfaces with large profile gradients and discontinuities. Measuring the profile of such samples requires the large scanning range of VSI (with vertical resolution in the order of 1/100 of a wavelength and scanning range in the order of millimeters), while characterizing a smooth surface texture requires the vertical resolution normally afforded by PSI (in the order of 1/1000 of a wavelength or less). 
     This problem has been addressed by the development of enhanced VSI algorithms (named EVSI in the art) that combine both PSI and VSI. In particular, Ser. No. 11/473,447 disclosed a new approach (named HDVSI, from high-definition VSI) based on performing a coarse calculation of surface profile with a conventional VSI method such as center of mass (COM), quadrature center of mass (QCOM), or zero crossing detection (ZCD), and concurrently performing a phase calculation carried out using the same irradiance frame data acquired for VSI. However, the phase calculation does not utilize the conventional n-frame phase shifting approach of the prior art. Instead, it utilizes a quadrature-demodulation (QD) algorithm applied to the irradiance data contained in the VSI correlogram. 
     As a result of this concurrent procedure, the phase calculation is independent of the position of any particular interferometric fringe and, therefore, it is more accurate and its results are more certain than those produced by the combined VSI/PSI methods of the prior art. Once both calculations are accomplished, the phase data are incorporated into the coarse profile data through a unique “unwrapping” method that yields a final surface map with sub-nanometer resolution within a large z-height range. 
     The HDVSI algorithm decouples the calculation of phase from the calculation of fringe-intensity peak. This results from the fact that QD manipulation of the correlogram data does not require a prior determination of the position of the coherence intensity peak. Therefore, the calculated phase map of the surface does not inherit the errors produced by the calculation of the coarse map. Such errors are removed from the coarse map, prior to combining it with the phase map, by consistently rounding all VSI measurements to an integer multiple of 2π. In other words, the VSI coarse map is rounded to an integer multiple of λ/2. (Typically, the accuracy of a VSI coarse map is much better than λ/2.) As a result, the VSI errors are completely removed from the VSI coarse map through the rounding process, an achievement that is not possible with prior-art approaches. Thus, the HDVSI algorithm allows profiling to sub-nanometer accuracy with conventional equipment, the only necessary changes being in the computational components required to round the VSI data, to implement the QD algorithm, and to combine the phase data with the corrected VSI data in real time, as data are acquired during the scan. 
     Ideally, the scanner steps and the effective wavelength of the light source are constant and can be determined through precise calibration procedures. However, in reality the scanner position suffers from system errors, such as scanner nonlinearity and random deviations due to mechanical vibrations. Similarly, the effective wavelength of the irradiance produced by light source changes with the slope of the sample surface and also as a result of scanner nonlinearities. Therefore, the scanner sampling step and the effective wavelength are not constant and cause errors that distort the measured sample-surface profile. 
     A techniques for real-time correction of scanner-nonlinearities is described in a separate disclosure. The present invention is directed at a novel approach for the correction of errors caused by changes in the system&#39;s effective wavelength. 
     BRIEF SUMMARY OF THE INVENTION 
     The invention is based on combining the HDVSI procedure disclosed in Ser. No. 11/473,447 with a different quadrature-demodulation parameter found suppress phase errors due to effective-wavelength variations in the light irradiating curved sample surfaces. In HDVSI, the fringe phase and the peak of the coherent envelope of the correlogram are calculated independently and in parallel. The fringe phase is calculated using a quadrature demodulation technique wherein the nominal scanner step size and the effective wavelength are filter parameters. 
     According to the present invention, the error correction for effective-wavelength variations is implemented by adjusting the filter parameters of the quadrature demodulation algorithm using a phase step that accounts for phase-step changes associated with variations in the effective wavelength irradiating the sample when the surface is curved. The nominal phase step is replaced in the filter with an actual phase step size that includes a phase parameter generated for each pixel by calibration or modeling of a specific type of surface. As detailed below, this substitution eliminates all errors produced by surface-dependent variations in the effective wavelength of the irradiating light. 
     Various other advantages of the invention will become clear from its description in the specification that follows and from the novel features particularly pointed out in the appended claims. Therefore, to the accomplishment of the objectives described above, this invention consists of the features hereinafter illustrated in the drawings, fully described in the detailed description of the preferred embodiment, and particularly pointed out in the claims. However, such drawings and description disclose but a few of the various ways in which the invention may be practiced. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  is an illustration of the typical fringe noise produced by variations in the effective wavelength of an interferometric profiler when an HDVSI algorithm is used to measure a spherical sample surface. 
         FIG. 2  shows the effect of the error-correction method of the invention when applied to the interferometric data of  FIG. 1 . 
         FIG. 3  is a diagram representing the application the filter parameter of the invention to the HDVSI quadrature-demodulation algorithm. 
         FIG. 4  is a schematic representation of a conventional interferometric profiler suitable to practice the invention. 
     
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     The heart of the invention lies in the use of a demodulation parameter that accounts for phase errors caused by variations in the effective wavelength of light in the calculation of phase in the HDVSI algorithm described in Ser. No. 11/473,447. In particular, the phase change produced at each pixel by the effective-wavelength variation is determined and added to the nominal phase step used as a parameter in the demodulation algorithm. This substitution has been found to substantially correct the phase errors introduced by wavelength variations during the data-acquisition scan of the object. 
     The invention for convenience is described throughout in terms of VSI (i.e., the process of determining the coherence peak of an interferogram produced by a mechanical scan), but the term is intended to refer as well and include any coherence-peak sensing method that may be applicable to the procedure. Therefore, the broader notion of producing a coarse map of the sample surface by finding the coherence peak of interference fringes, however produced, is intended to be part of and covered by this invention. Thus, the invention is intended to include any method that produces interference fringes by changing the optical path difference between a test and a reference beam, such as, without limitation, mechanical scanning and polarization shifting methods. 
     An analysis of the error produced by variations in the effective wavelength of the irradiance illuminating the sample shows that, when constant filter parameters (k 0  and ΔΦ 0 , defined below) are used and when wavelength-variation errors (τ 2 ) are present, the phase obtained from the HDVSI demodulation filter is a combination of the true sample phase (Φ) and a phase error (δ Φ ) related to the wavelength-variation error. In particular, as the VSI scan is implemented through a plurality of scanning steps (ΔΦ or Δz, with reference to phase or distance, respectively) of predetermined nominal size (referred to herein interchangeably as nominal scan steps or phase steps), each detector pixel produces a signal with an amplitude that passes through and drops off very rapidly from its maximum value [I max , corresponding to zero optical path difference (OPD)]. The interferometric signal can be expressed as
 
 I ( n )= G ( z   n   −z   0 )cos(2 nk   0   Δz −Φ)  (1)
 
where G(z n −z 0 ) is the Gaussian envelope of the correlogram (i.e., the amplitude of irradiance at each scanner position); z 0  and z n  are the positions of the scanning objective at zero OPD and at the n th  scanning step, respectively; k 0  is the wave number (i.e., k 0 =2π/λ 0 , λ 0  being the nominal effective wavelength of the illumination); and Δz, and Φ are the n th  scan step and the phase of the correlogram, respectively.
 
     As taught in Ser. No. 11/473,447, the interferometric signals I(n) are manipulated with the HDVSI algorithm to produce a low-noise sub-nanometer resolution surface map. The HDVSI algorithm comprises several steps. After undergoing pre-processing filtration, if needed, the input to the algorithm (i.e., the interferometric irradiance signals produced by the VSI scan) is processed simultaneously in two ways that allow for the independent calculation of a coarse height z 0  of each surface point and of the phase Φ of the correlogram corresponding to that point. In the aggregate, the values of z 0  and Φ for all surface points imaged onto the light detector yield a coarse surface map, z VSI (x,y), and a high-definition phase map, Φ(x,y), of the test surface, respectively. As one skilled in the art would readily understand, the phase map Φ(x,y) can be equivalently expressed as a fine surface (or height) map, z HD (x,y), because the phase difference between two surface points is directly proportional to the difference in their heights. Therefore, for the purposes of this disclosure, the maps Φ(x,y) and z HD (x,y) and the corresponding terms “phase map,” “surface map,” “height map” and “map” may be used interchangeably. As a separate step of the HDVSI algorithm, the two maps are judiciously combined to generate the final high-definition surface map, z HDVSI (x,y). 
     The phase Φ of the correlogram is determined by applying a quadrature-demodulation approach. The signal I(n) is passed through a QD filter that performs several operations. First, the irradiance signal is differentiated to eliminate the DC component in the original correlogram. Then, it is split into “in-phase” and “quadrature” signal components, J and Q, respectively, wherein
 
 Q   n   =G ( z   n   −z   0 )cos( n 2 kΔz   n −Φ)sin( n 2 k   0   Δz   0 ) and  (2a)=
 
 J   n   =G ( z   n   −z   0 )cos( n 2 kΔz   n −Φ)cos( n 2 k   0   Δz   0 )  (2b)
 
where k=2π/λ, λ being the actual effective wavelength of the irradiance; Δz n  is the actual size of the n th  scan step; and Δz 0  is the parameter of the filter, chosen to match the nominal step size of the scan [typically equal to mλ 0 /8, m=0, 1, 2 . . . ].
 
     At each pixel x,y of the sample surface, the deviations of the wavelength from the nominal effective value λ 0  may be represented by a generalized wavelength error τ(x,y) as follows,
 
λ( x,y )=λ 0 +τ( x,y )  (3)
 
Assuming an actual scan step equal to the nominal scan step of size λ 0 /8 (i.e., 90°) and substituting for Δz n , Δz 0  and k, Equations 2a, 2b, and 3 yield the following:
 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           Q 
                           n 
                         
                         = 
                         
                           
                             G 
                             ⁡ 
                             
                               ( 
                               n 
                               ) 
                             
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   n 
                                   ⁢ 
                                   
                                     
                                       4 
                                       ⁢ 
                                       π 
                                     
                                     λ 
                                   
                                   ⁢ 
                                   
                                     
                                       λ 
                                       0 
                                     
                                     8 
                                   
                                 
                                 - 
                                 Φ 
                               
                               ) 
                             
                           
                           ⁢ 
                           
                             sin 
                             ⁡ 
                             
                               ( 
                               
                                 n 
                                 ⁢ 
                                 
                                   π 
                                   2 
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           = 
                           
                             
                               G 
                               ⁡ 
                               
                                 ( 
                                 n 
                                 ) 
                               
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   
                                     n 
                                     ⁢ 
                                     
                                       π 
                                       2 
                                     
                                     ⁢ 
                                     
                                       
                                         λ 
                                         0 
                                       
                                       
                                         
                                           λ 
                                           0 
                                         
                                         - 
                                         
                                           τ 
                                           xy 
                                         
                                       
                                     
                                   
                                   - 
                                   Φ 
                                 
                                 ) 
                               
                             
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 
                                   n 
                                   ⁢ 
                                   
                                     π 
                                     2 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                         , 
                         and 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     4 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
                           J   n     =       G   ⁡     (   n   )       ⁢     cos   ⁡     (       n   ⁢       4   ⁢   π     λ     ⁢       λ   0     8       -   Φ     )       ⁢     cos   ⁡     (     n   ⁢     π   2       )                     =       G   ⁡     (   n   )       ⁢     cos   ⁡     (       n   ⁢     π   2     ⁢       λ   0         λ   0     -     τ   xy           -   Φ     )       ⁢       cos   ⁡     (     n   ⁢     π   2       )       .                     (     4   ⁢   b     )               
When τ xy /λ 0  is much smaller than 1,
 
                     Q   n     =       G   ⁡     (   n   )       ⁢     cos   ⁡     (       n   ⁢     π   2     ⁢     (     1   +     δ   xy       )       -   Φ     )       ⁢     sin   ⁡     (     n   ⁢     π   2       )       ⁢           ⁢   and             (     5   ⁢   a     )               
Defining δ xy  as τ xy /λ 0 , the phase error due to effective-wavelength variations, Equations 4a and 4b may be written as
 
     
       
         
           
             
               
                 
                   
                     Q 
                     n 
                   
                   = 
                   
                     
                       G 
                       ⁡ 
                       
                         ( 
                         n 
                         ) 
                       
                     
                     ⁢ 
                     
                       cos 
                       ⁡ 
                       
                         ( 
                         
                           
                             n 
                             ⁢ 
                             
                               π 
                               2 
                             
                             ⁢ 
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   δ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                           - 
                           Φ 
                         
                         ) 
                       
                     
                     ⁢ 
                     
                       sin 
                       ⁡ 
                       
                         ( 
                         
                           n 
                           ⁢ 
                           
                             π 
                             2 
                           
                         
                         ) 
                       
                     
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     and 
                   
                 
               
               
                 
                   ( 
                   
                     5 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
     
       
         
           
             
               
                 
                   
                     J 
                     n 
                   
                   = 
                   
                     
                       G 
                       ⁡ 
                       
                         ( 
                         n 
                         ) 
                       
                     
                     ⁢ 
                     
                       cos 
                       ⁡ 
                       
                         ( 
                         
                           
                             n 
                             ⁢ 
                             
                               π 
                               2 
                             
                             ⁢ 
                             
                               ( 
                               
                                 1 
                                 + 
                                 
                                   δ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                           - 
                           Φ 
                         
                         ) 
                       
                     
                     ⁢ 
                     
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           
                             n 
                             ⁢ 
                             
                               π 
                               2 
                             
                           
                           ) 
                         
                       
                       . 
                     
                   
                 
               
               
                 
                   ( 
                   
                     5 
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     The quadrature demodulation procedure involves calculating the ratio 
                     (         ∑     n   =   1     N     ⁢           ⁢     Q   n           ∑     n   =   1     N     ⁢           ⁢     J   n         )     .           (   6   )               
Accordingly, from Equation 5a,
 
                       ∑     n   =   1     N     ⁢           ⁢     Q   n       =       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     cos   ⁡     (       n   ⁢     π   2     ⁢     (     1   +     δ   xy       )       -   Φ     )       ⁢       sin   ⁡     (     n   ⁢     π   2       )       .                 (     7   ⁢   a     )               
Equation 2b similarly yields:
 
     
       
         
           
             
               
                 
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         1 
                       
                       N 
                     
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       J 
                       n 
                     
                   
                   = 
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         1 
                       
                       N 
                     
                     ⁢ 
                     
                       
                         G 
                         ⁡ 
                         
                           ( 
                           n 
                           ) 
                         
                       
                       ⁢ 
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           
                             
                               n 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                               ⁢ 
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     δ 
                                     xy 
                                   
                                 
                                 ) 
                               
                             
                             - 
                             Φ 
                           
                           ) 
                         
                       
                       ⁢ 
                       
                         
                           sin 
                           ⁡ 
                           
                             ( 
                             
                               n 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                             
                             ) 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     7 
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     Applying to Equation 7a the general trigonometric formulas
 
sin( a+b )=sin( a )cos( b )+cos( a )sin( b ) and  (8a)
 
sin( a−b )=sin( a )cos( b )−cos( a )sin( b ),  (8b)
 
which, combined, produce the relation
 
sin( a+b )−sin( a−b )=2 cos( a )sin( b ),  (8c)
 
Equation 7a becomes
 
                     2   ⁢       ∑     n   =   1     N     ⁢           ⁢     Q   n         =         ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (       n   ⁢           ⁢   π     +     n   ⁢     π   2     ⁢     δ   xy       -   Φ     )           -       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (       n   ⁢     π   2     ⁢     δ   xy       -   Φ     )                     (   9   )               
The first summation term is approximately zero when N&gt;&gt;1; therefore,
 
                           2   ⁢       ∑     n   =   1     N     ⁢           ⁢     Q   n         =     -       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (       n   ⁢     π   2     ⁢     δ   xy       -   Φ     )                         =         -     cos   ⁡     (   Φ   )         ⁢       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (     n   ⁢     π   2     ⁢     δ   xy       )             +       sin   ⁡     (   Φ   )       ⁢       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     cos   ⁡     (     n   ⁢     π   2     ⁢     δ   xy       )                             (     10   ⁢   a     )               
which follows from Equation 7b. Similarly,
 
     
       
         
           
             
               
                 
                   
                     2 
                     ⁢ 
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       
                         J 
                         n 
                       
                     
                   
                   = 
                   
                     
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           Φ 
                           ) 
                         
                       
                       ⁢ 
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                         ⁢ 
                         
                           
                             G 
                             ⁡ 
                             
                               ( 
                               n 
                               ) 
                             
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 n 
                                 ⁢ 
                                 
                                   π 
                                   2 
                                 
                                 ⁢ 
                                 
                                   δ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                     + 
                     
                       
                         sin 
                         ⁡ 
                         
                           ( 
                           Φ 
                           ) 
                         
                       
                       ⁢ 
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                         ⁢ 
                         
                           
                             G 
                             ⁡ 
                             
                               ( 
                               n 
                               ) 
                             
                           
                           ⁢ 
                           
                             sin 
                             ⁡ 
                             
                               ( 
                               
                                 n 
                                 ⁢ 
                                 
                                   π 
                                   2 
                                 
                                 ⁢ 
                                 
                                   δ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     10 
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     Equations 10a and 10b may be further simplified by setting 
     
       
         
           
             
               
                 
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         1 
                       
                       N 
                     
                     ⁢ 
                     
                       
                         G 
                         ⁡ 
                         
                           ( 
                           n 
                           ) 
                         
                       
                       ⁢ 
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           
                             n 
                             ⁢ 
                             
                               π 
                               2 
                             
                             ⁢ 
                             
                               δ 
                               xy 
                             
                           
                           ) 
                         
                       
                     
                   
                   = 
                   
                     a 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     and 
                   
                 
               
               
                 
                   ( 
                   
                     11 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
                         ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (     n   ⁢     π   2     ⁢     δ   xy       )           =   b     ⁢                   (     11   ⁢   b     )               
which yields
 
     
       
         
           
             
               
                 
                   
                     
                       2 
                       ⁢ 
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           Q 
                           n 
                         
                       
                     
                     = 
                     
                       
                         
                           - 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               Φ 
                               ) 
                             
                           
                         
                         ⁢ 
                         b 
                       
                       + 
                       
                         
                           sin 
                           ⁡ 
                           
                             ( 
                             Φ 
                             ) 
                           
                         
                         ⁢ 
                         a 
                       
                     
                   
                   , 
                   
                       
                   
                   ⁢ 
                   and 
                 
               
               
                 
                   ( 
                   
                     12 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
                     2   ⁢       ∑     n   =   1     N     ⁢           ⁢     J   n         =         cos   ⁡     (   Φ   )       ⁢   a     +       sin   ⁡     (   Φ   )       ⁢     b   .                 (     12   ⁢   b     )               
Defining the phase error in terms of a new variable ζ xy  at each x,y pixel, such that
 
     
       
         
           
             
               
                 
                   
                     cos 
                     ⁡ 
                     
                       ( 
                       
                         ζ 
                         xy 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       a 
                       
                         
                           
                             a 
                             2 
                           
                           + 
                           
                             b 
                             2 
                           
                         
                       
                     
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     and 
                   
                 
               
               
                 
                   ( 
                   
                     13 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
                       sin   ⁡     (     ζ   xy     )       =     b         a   2     +     b   2             ,           (     13   ⁢   b     )               
and further dividing Equations 12a and 12b by √{square root over (a 2 +b 2 )} and substituting Equations 13a and 13b produces the following relations in terms of the error ζ xy :
 
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           2 
                           ⁢ 
                           
                             
                               ∑ 
                               
                                 n 
                                 = 
                                 1 
                               
                               N 
                             
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             
                               
                                 Q 
                                 n 
                               
                               / 
                               
                                 
                                   
                                     a 
                                     2 
                                   
                                   + 
                                   
                                     b 
                                     2 
                                   
                                 
                               
                             
                           
                         
                         = 
                         
                           
                             
                               - 
                               
                                 cos 
                                 ⁡ 
                                 
                                   ( 
                                   Φ 
                                   ) 
                                 
                               
                             
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 
                                   ζ 
                                   xy 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 Φ 
                                 ) 
                               
                             
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 
                                   ζ 
                                   xy 
                                 
                                 ) 
                               
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         
                           = 
                           
                             sin 
                             ⁡ 
                             
                               ( 
                               
                                 Φ 
                                 - 
                                 
                                   ζ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                         
                         , 
                         
                             
                         
                         ⁢ 
                         and 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     14 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
     
       
         
           
             
               
                 
                   
                     
                       
                         
                           2 
                           ⁢ 
                           
                             
                               ∑ 
                               
                                 n 
                                 = 
                                 1 
                               
                               N 
                             
                             ⁢ 
                             
                               
                                 J 
                                 n 
                               
                               / 
                               
                                 
                                   
                                     a 
                                     2 
                                   
                                   + 
                                   
                                     b 
                                     2 
                                   
                                 
                               
                             
                           
                         
                         = 
                         
                           
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 Φ 
                                 ) 
                               
                             
                             ⁢ 
                             
                               cos 
                               ⁡ 
                               
                                 ( 
                                 
                                   ζ 
                                   xy 
                                 
                                 ) 
                               
                             
                           
                           + 
                           
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 Φ 
                                 ) 
                               
                             
                             ⁢ 
                             
                               sin 
                               ⁡ 
                               
                                 ( 
                                 
                                   ζ 
                                   xy 
                                 
                                 ) 
                               
                             
                           
                         
                       
                     
                   
                   
                     
                       
                         = 
                         
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 Φ 
                                 - 
                                 
                                   ζ 
                                   xy 
                                 
                               
                               ) 
                             
                           
                           . 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     14 
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     Finally, dividing Equations 14a and 14b yields 
                 (         ∑     n   =   1     N     ⁢     Q   n           ∑     n   =   1     N     ⁢     J   n         )     =       sin   ⁡     (     Φ   -     ζ   xy       )         cos   ⁡     (     Φ   -     ζ   xy       )           ,         
or
 
                     Φ   -     ζ   xy       =         tan     -   1       ⁡     (         ∑     n   =   1     N     ⁢     Q   n           ∑     n   =   1     N     ⁢     J   n         )       .             (   15   )               
Form Equation 15 it is apparent that, when constant filter parameters (k 0 , Δz n ) are used and wavelength-related errors τ xy  are present during the scan, the phase obtained from the demodulation filter is the combination of the true sample phase Φ and a phase error ζ xy  related to the wavelength-variation error.
 
     Thus, in practice, when the effective wavelength changes (λ=λ 0 −τ) across a sloped surface, it distorts the relative phase between surface pixels. If the distorted phases are converted to z heights using the nominal effective wavelength (λ=λ 0 ), a distorted z surface map is obtained in the form of 2π jump error, as illustrated in  FIG. 1  for a spherical surface. According to the invention, if the demodulation filters 
             cos   ⁡     (     n   ⁢     π   2       )           
and
 
             sin   ⁢           ⁢     (     n   ⁢     π   2       )           
of Equations 7a and 7b above are replaced by
 
             sin   ⁡     (     n   ⁢     π   2     ⁢     (     1   +     δ   xy       )       )           
and
 
               cos   ⁡     (     n   ⁢     π   2     ⁢     (     1   +     δ   xy       )       )       ,         
then the quadrature Equations 7a and 7b become:
 
     
       
         
           
             
               
                 
                   
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                       ⁢ 
                       
                         Q 
                         n 
                       
                     
                     = 
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                       ⁢ 
                       
                         
                           G 
                           ⁡ 
                           
                             ( 
                             n 
                             ) 
                           
                         
                         ⁢ 
                         
                           cos 
                           ⁡ 
                           
                             ( 
                             
                               
                                 n 
                                 ⁢ 
                                 
                                   π 
                                   2 
                                 
                                 ⁢ 
                                 
                                   ( 
                                   
                                     1 
                                     + 
                                     
                                       ζ 
                                       xy 
                                     
                                   
                                   ) 
                                 
                               
                               - 
                               Φ 
                             
                             ) 
                           
                         
                         ⁢ 
                         
                           sin 
                           ⁡ 
                           
                             ( 
                             
                               n 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                               ⁢ 
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     δ 
                                     xy 
                                   
                                 
                                 ) 
                               
                             
                             ) 
                           
                         
                       
                     
                   
                   , 
                   
                       
                   
                   ⁢ 
                   and 
                 
               
               
                 
                   ( 
                   
                     16 
                     ⁢ 
                     a 
                   
                   ) 
                 
               
             
           
         
       
     
     
       
         
           
             
               
                 
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         1 
                       
                       N 
                     
                     ⁢ 
                     
                       J 
                       n 
                     
                   
                   = 
                   
                     
                       ∑ 
                       
                         n 
                         = 
                         1 
                       
                       N 
                     
                     ⁢ 
                     
                       
                         G 
                         ⁡ 
                         
                           ( 
                           n 
                           ) 
                         
                       
                       ⁢ 
                       
                         cos 
                         ⁡ 
                         
                           ( 
                           
                             
                               n 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                               ⁢ 
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     ζ 
                                     xy 
                                   
                                 
                                 ) 
                               
                             
                             - 
                             Φ 
                           
                           ) 
                         
                       
                       ⁢ 
                       
                         
                           cos 
                           ⁡ 
                           
                             ( 
                             
                               n 
                               ⁢ 
                               
                                 π 
                                 2 
                               
                               ⁢ 
                               
                                 ( 
                                 
                                   1 
                                   + 
                                   
                                     δ 
                                     xy 
                                   
                                 
                                 ) 
                               
                             
                             ) 
                           
                         
                         . 
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     16 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     Applying again the trigonometric formula of Equation 8c to Equation 16a yields: 
                     2   ⁢       ∑     n   =   1     N     ⁢     Q   n         =         ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (       n   ⁢           ⁢     π   ⁡     (     1   +     δ   xy       )         -   Φ     )           +       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (   Φ   )                     (   17   )               
The first summation term is approximately zero when N&gt;&gt;1; therefore,
 
                     2   ⁢       ∑     n   =   1     N     ⁢     Q   n         ≈       ∑     n   =   1     N     ⁢       G   ⁡     (   n   )       ⁢     sin   ⁡     (   Φ   )                   (     18   ⁢   a     )               
which follows from Equation 7b. Similarly,
 
     
       
         
           
             
               
                 
                   
                     2 
                     ⁢ 
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                       ⁢ 
                       
                         J 
                         n 
                       
                     
                   
                   = 
                   
                     
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                         ⁢ 
                         
                           
                             G 
                             ⁡ 
                             
                               ( 
                               n 
                               ) 
                             
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               
                                 
                                   n 
                                   ⁢ 
                                   
                                       
                                   
                                   ⁢ 
                                   
                                     π 
                                     ⁡ 
                                     
                                       ( 
                                       
                                         1 
                                         + 
                                         
                                           δ 
                                           xy 
                                         
                                       
                                       ) 
                                     
                                   
                                 
                                 - 
                                 Φ 
                               
                               ) 
                             
                           
                         
                       
                       + 
                       
                         
                           ∑ 
                           
                             n 
                             = 
                             1 
                           
                           N 
                         
                         ⁢ 
                         
                           
                             G 
                             ⁡ 
                             
                               ( 
                               n 
                               ) 
                             
                           
                           ⁢ 
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               Φ 
                               ) 
                             
                           
                         
                       
                     
                     ≈ 
                     
                       
                         ∑ 
                         
                           n 
                           = 
                           1 
                         
                         N 
                       
                       ⁢ 
                       
                         
                           G 
                           ⁡ 
                           
                             ( 
                             n 
                             ) 
                           
                         
                         ⁢ 
                         
                           
                             cos 
                             ⁡ 
                             
                               ( 
                               Φ 
                               ) 
                             
                           
                           . 
                         
                       
                     
                   
                 
               
               
                 
                   ( 
                   
                     18 
                     ⁢ 
                     b 
                   
                   ) 
                 
               
             
           
         
       
     
     From Equations 18a and 18b it follows that 
                       tan     -   1       ⁡     (         ∑     n   =   1     N     ⁢     Q   n           ∑     n   =   1     N     ⁢     J   n         )       =     Φ   .             (   19   )               
Equation 19 shows that, if the quadrature demodulation filter of the HDVSI algorithm uses a parameter equal
 
             (     n   ⁢     π   2     ⁢     (     1   +     δ   xy       )       )         
instead of the standard parameter
 
               (     n   ⁢     π   2       )     ,         
the phase error δ xy  produced by effective-wavelength variations is suppressed.
 
     Thus, in order to correct errors caused by effective-wavelength deviation according to the invention, an effective-wavelength deviation map, τ xy , is needed (from which a corresponding phase-error map, δ xy , is advantageously derived). Such a map can be obtained empirically by a calibration process in which the spectral amplitude of the correlogram produced at each pixel is calculated using the quadrature-demodulation method. From the spectral amplitude, the effective wavelength, and thus the wavelength variation, or error, relative to the nominal wavelength at each pixel can be extracted. The extracted effective wavelength variation error can then used in the QD filter to correct the effective wavelength at each pixel. 
     An alternative approach to produce the δ xy  map is to use a parameterized model fitted empirically through calibration. For example, given the fact that variations in the effective wavelength are known to result mainly from curvatures in the sample surface, a model that reflects the slope of the surface can be used. For example, it can be assumed that δ xy  is equal to some function of the surface slope [i.e., δ xy =F(slope)]. One example that proved to be effective in the error correction procedure of the invention for a spherical surface is the function
 
δ xy =αr β ,  (20)
 
where r is the distance between pixel (x,y) and a center pixel (x 0 ,y 0 ). This center pixel (x 0 ,y 0 ) is usually a surface location where the surface slope equal to zero. Accordingly,
 
 r =[( x−x   0 ) 2 +( y−y   0 ) 2 ] 1/2 .  (21)
 
The coefficient α and the exponent β are empirical parameters that can be determined through some kind of calibration process. For example, for a given spherical surface shape, the α and β parameters can be determined by finding the values that minimize the surface 2π error caused by the effective-wavelength variation.
 
     Once a δ xy  map is obtained, it is used according to the invention in the demodulation filter, as described above, to correct the image distortion caused by wavelength variations.  FIG. 2  shows the corrected spherical surface image produced by applying correction procedure to the data of  FIG. 1 . The error correction algorithm is illustrated schematically in  FIG. 3 . 
     Note that the wavelength error correction of the invention may be applied to the phase-demodulation module concurrently with the scanner-nonlinearity error correction described in a separate disclosure for real-time HDVSI data processing. 
     The invention may be implemented with conventional interferometric equipment modified to perform the correction procedure described above. With reference to  FIG. 4 , wherein like parts are designated with the same numerals and symbols used throughout the drawings, such an interferometer  10  may comprise a light source  12  (such as white-light) directing a beam of light through an aperture  14  and field stop  16  toward a beam splitter  18  that reflects the light in the direction of a test surface S. The light reflected by the beam splitter  18  passes through a microscope objective  20  focused on the test surface S. The objective incorporates an interferometer  22 , such as Mirau, comprising a beam splitter  24  and a reference mirror adapted for relative movement with respect to the test surface, so that two light beams are generated for producing interference fringes as a result of the optical path difference between the reference mirror and the test surface S. The beams reflected from the reference mirror  26  and the test surface S pass back up through the optics of the microscope objective  20  and through the beam splitter  18  to a solid-state detector array  28  in a camera in coaxial alignment with the objective  20 . 
     Typically, the detector array  28  consists of individual CCD cells or other sensing apparatus adapted to produce a two-dimensional array of digitized intensity data corresponding to light signals received at each sensor cell as a result of the interference of the coherent light beams reflected from individual x-y coordinates or pixels in the surface S and from corresponding coordinates in the reference mirror  26 . Appropriate electronic hardware (not shown) is provided to transmit the digitized intensity data generated by the detector to a microprocessor for processing. The microscope objective  20 , as well as the interferometer typically incorporated within it, is adapted for vertical movement to focus the image of the test surface on the detector array  28 . 
     According to the invention described in U.S. Pat. No. 6,987,570, an additional reference-signal interferometer may be built into the conventional measurement interferometer of the instrument and be dedicated to monitor OPD changes through a separate reference-signal channel. Such a reference signal is necessary to correct errors produced by scanner nonlinearities, as described in a separate disclosure. As illustrated in  FIG. 4 , a laser  62  may be used to provide enough light for a good-quality reference signal and a beamsplitter  64  is added in the optical measurement path to direct the laser beam toward the main splitter  18 . A portion of the light is passed through to a reference mirror  66  for the reference signal, while the balance of the light is reflected toward a mirror or other reflective element  68  (such as a corner cube or a reflecting grating) placed on the moving part of the scanner  38  of the measurement interferometer. 
     While the present invention has been shown and described herein in what is believed to be the most practical and preferred embodiments, it is recognized that departures can be made therefrom within the scope of the invention. For example, neither the VSI-scan nor the phase parameter of the QD-filter used for calculation of the phase of the correlogram has to equal π/2, but it should be chosen to equal the value that is most appropriate for conducting the VSI scan. Also, the invention has been described in terms of a quadrature demodulation algorithm used to calculate phase in real time. However, it is recognized that quadrature demodulation is a particular case of wavelet filtering and that, as such, other wavelet filters may be used as well. 
     Therefore, the invention is not to be limited to the details disclosed herein but is to be accorded the full scope of the claims so as to embrace any and all equivalent processes and products.