Abstract
Stitching interferometry is performed by collecting interferometric data from overlapped sub-apertures and stitching these data together to provide a full surface map. The propagation of the systematic error in the measured subset data is one of the main error sources in stitching interferometry for accurate reconstruction of the surface topography. In this work, we propose, using the redundancy of the captured subset data, two types of two-dimensional (2D) self-calibration stitching algorithms to overcome this issue by in situ estimating the repeatable high-order additive systematic errors, especially for the application of measuring X-ray mirrors. The first algorithm, called CS short for “Calibrate, and then Stitch”, calibrates the high-order terms of the reference by minimizing the de-tilted discrepancies of the overlapped subsets and then stitches the reference-subtracted subsets. The second algorithm, called SC short for “Stitch, and then Calibrate”, stitches a temporarily result and then calibrates the reference from the de-tilted discrepancies of the measured subsets and the temporarily stitched result. In the implementation of 2D scans in $x$- and $y$-directions, step randomization is introduced to generate nonuniformly spaced subsets which can diminish the periodic stitching errors commonly observed in evenly spaced subsets. The regularization on low-order terms enables a highly flexible option to add the curvature and twist acquired by another system. Both numerical simulations and experiments are carried out to verify the proposed method. All the results indicate that 2D high-order repeatable additive systematic errors can be retrieved from the 2D redundant overlapped data in stitching interferometry.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The surface topographic measurement based on interferometry is a major tool for optical surfaces inspection [1,2]. As the accuracy and precision requirement is getting higher while the measuring range and size being larger, the “standard” interferometer Field Of View (FOV) becomes insufficient for some industrial applications, e.g. to measure large or long and curved mirrors. To access a bigger measuring size with a limited interferometer FOV, the Sub-Aperture (SA) stitching interferometry [3–5] was developed as a low cost solution to flexibly extend the measuring range and size, while preserving the native lateral resolution. Stitching interferometry is one of the high-precision optical surface measurement techniques and it is applied in many fields [3,4,6–10].
Stitching interferometry acquires many highly-overlapped SAs. Using the high redundancy of the overlapped data in the acquisition process, it is possible to retrieve additional information [11–13]. In some applications, the Surface Under Test (SUT) has the same quality as the reference surface of the interferometer [8,9]. Therefore, it is a practical and meaningful task to carefully design the stitching strategy from data acquisition to stitching algorithm. By doing this, we can fully use the redundant measurement data to get an accurate estimation of the systematic error. Polack et. al. use Legendre polynomial interpolation of the reference to solve the 1D stitching problem in the case where the stitching steps are not an integer number of pixels [14]. Nicolas et. al. provided a solution to reconstruct a one-dimensional (1D) profile of the surface and extract the 1D additive systematic error in the stitching process [15]. These 1D solutions [14,15] based on Legendre polynomials or cubic B-splines offer the capability of using the redundant stitching data in a practical manner and paths a way for two-dimensional (2D) cases. The 2D self-calibration stitching method for the general optics has been studied [16–20], which show the effectiveness in reducing the reference errors.
In this work, we present two types of 2D stitching algorithms capable of estimating high-order additive systematic errors of the interferometer with linear least squares method, especially for measuring X-ray mirrors. The scanning steps during the data acquisition are specially designed for reference retrieval to avoid unwanted periodic errors. We validate the proposed method with simulation using a known reference as a benchmark. In parallel, experiments are also carried out to demonstrate the feasibility of the proposed method to retrieve reference during the stitching measurement process. The retrieved 2D reference map is compared with the calibrated reference data. The capability and limitations of the proposed method are discussed for possible future improvements.
2. Principle and method
2.1 Mathematical model of 2D data stitching
In order to clearly describe the stitching problem, a mathematical model of the 2D stitching is established. As shown in Fig. 1, the measurement of $n$th subset among the total $N$ subsets $\mathbf {m}_{n}(x,y)$ is taken by an interferometer at location $(x_{n},y_{n})$ during the stitching acquisition process.
The $n$th subset can be described by the expression
2.2 Two self-calibration stitching algorithms
Two types of self-calibration stitching algorithms are presented. In the first self-calibration stitching algorithm, the reference $\mathbf {r}$ is first calibrated and then stitch the reference-removed subsets $\mathbf {m}-\mathbf {r}$ into a whole piece of $\mathbf {z}$. We name this algorithm CS which is short for “Calibrate, and then Stitch.” The second self-calibration stitching algorithm, called SC short for “Stitch, and then Calibrate”, stitches the $i$th updated reference-removed subsets $\mathbf {m}-\mathbf {r}_i$ to get $\mathbf {z}_i$ first, and then estimate the reference $\mathbf {r}_{i+1}$ for the next iteration. It is interesting to notice that the iterative operation has also been studied and applied for the reference estimation with three-flat test [21,22].
Both of CS and SC algorithms can be applied to any kind of underlying stitching strategy. Since we solve the stitching problem based on matrix operation in this work, the proposed self-calibration stitching algorithms (CS and SC) will be mainly addressed to cooperate with the pixel-relation-based stitching strategy and the subset-relation-based stitching strategy [23]. In pixel-relation-based stitching strategy, all relations between the measured data $\mathbf {m}$, $\mathbf {t}$, $\mathbf {z}$, and $\mathbf {r}$ (if involved) are established at each pixel, while the subset-relation-based stitching strategy establishes those dependant equations with each subset.
When the CS algorithm cooperating with the pixel-relation-based stitching strategy named CS-P as illustrated in Fig. 2(a), we can consider that the SUT shape is cancelled out in the overlapped height difference $\mathbf {d}$. Similar to the 1D case [13], its relations with $\mathbf {t}$ and $\mathbf {r}$ can be established by matching the corresponding pixels as
After that, the stitching result $\mathbf {z}$ can be obtained by simply merging the reference-subtracted data $\mathbf {m}-\mathbf {r}$ with the known $\mathbf {t}$.
If the subset-relation-based stitching strategy is used with the CS algorithm, named CS-S, the initial reference $\mathbf {r}_0$ can be set as zeros. From the de-tilted height difference $\mathbf {\hat {d}}_i$ in the overlaps with $\mathbf {m}-\mathbf {r}_i$ in the $i$th iteration, only the reference is estimated through a compensatory update $\mathbf {r}_{i+1} = \mathbf {r}_i+\mathbf {\hat {r}}_i$, where the reference compensation amount $\mathbf {\hat {r}}_i$ is estimated from the de-tilted discrepancies $\mathbf {\hat {d}}_i$ by
By solving this linear least squares problem, $\mathbf {\hat {r}}_i$ is calculated usingThe essential idea of the SC algorithm in Fig. 2(c) is based on the following consideration. If the estimated reference $\mathbf {r}$ is correct, the discrepancy $\mathbf {d^s}$ between the reference-removed subsets $\mathbf {m}-\mathbf {r}$ and the stitched height $\mathbf {z}$ in corresponding SA should only be a tilted plane with random noise. The superscribe “$\mathbf {s}$” denotes the operation in SA. The discrepancy $\mathbf {d^s}$ can be de-tilted and symbolized as $\mathbf {\hat {d^s}}$. We can make iterative estimations of the reference $\mathbf {r}$ by minimizing this de-tilted discrepancies $\mathbf {\hat {d^s}}$ with a least squares method. As illustrated in Fig. 2(c), the reference $\mathbf {r}_0$ is initialized with zeros. For the $i$th iteration, $\mathbf {z}_i$ is stitched from the reference-subtracted subsets $\mathbf {m}-\mathbf {r}_i$. The reference compensation amount $\mathbf {\hat {r}}_i$ is then estimated from the de-tilted discrepancies $\mathbf {\hat {d^s}}_i$ between each $\mathbf {m}_{n}, n\in [1, N]$ and the stitched height $\mathbf {z}_i$ with the corresponding pixels inside the $n$th SA as
where $\mathbf {D^s}$ is the relation matrix. In this matrix based algorithm, some matrix operations are carefully taken when composing $\mathbf {D^s}$ to match pixels in all subsets with the ones of the stitched map in the corresponding SAs. The size of $\mathbf {D^s}$ can be fixed as $M\times M$ by collapsing matrices describing different SAs into one. As the size of $\mathbf {D^s}$ becomes independent of $N$, this operation will be useful when stitching a large number of SAs. By solving this linear least squares problem, $\mathbf {\hat {r}}_i$ is calculated usingAs a result, the interferometric measurements can be best explained by the stitched result $\mathbf {z}$ and the estimated reference $\mathbf {r}$ using these two types of self-calibration stitching algorithms. In the next section, we will address some possible ambiguity issues in this stitching model and self-calibration estimation followed by our corresponding strategies to resolve them.
3. Ambiguities and regularization in data acquisition and processing
As shown in Eq. (1), the measured $m_{n}(x,y)$ is a combination of three unknowns (SUT shape $\mathbf {z}$, reference $\mathbf {r}$, and geometric parameters $\mathbf {t}$). Ambiguities may happen under certain conditions which makes the stitching problem become ill-posed. To resolve the issues, some regularization needs to be applied in the stitching strategy.
3.1 Regularization for periodic errors in two directions
In many stitching applications, the $x$- and $y$-translation step sizes $\Delta x$ and $\Delta y$ are selected as a uniform value in one direction during the stitching process. However, using these uniform stepped stitching, the stitched height map will suffer with periodic errors in both SUT stitching and reference reconstruction, if we want to reconstruct the reference from the redundant data.
Let’s consider there is situation in which $z(x,y)$ and $r(x,y)$ are replaced by another pair of SUT shape and reference shape $z^*(x,y)$ and $r^*(x,y)$ given by the following expression
where $p(x)$ and $q(y)$ are periodic functions with periods $T_x$ and $T_y$ in $x$ and $y$ as If the step sizes are uniform, $x_{n}=n\cdot \Delta x$ and $y_{n}=n\cdot \Delta y$. Then choosing $T_x=\Delta x/n_x$ with $n_x\in \mathbb {N}$ and $T_y=\Delta y/n_y$ with $n_y\in \mathbb {N}$, we have $x_{n}=nn_x\cdot T_x$ and $y_{n}=nn_y\cdot T_y$. Since $nn_x \in \mathbb {N}$ and $nn_y \in \mathbb {N}$, we get $p(x-x_{n})=p(x-nn_x\cdot T_x)=p(x)$, and $q(y-y_{n})=q(y-nn_y\cdot T_y)=q(y)$. According to the equations above, we haveThis periodic error issue is due to the uniform stitching steps, in order to overcome this problem, we suggest to make randomized nonuniform steps for both $x$- and $y$-directions during the data acquisition process to make this ill-posed problem solvable. As illustrated in Fig. 3, the periodic repeated pattern in Fig. 3(a) is “broken” by the nonuniform steps by randomization in Fig. 3(b).
In our implementation, the stitching step sizes are always integer pixels, and we first calculate the preferred step sizes by the user-preferred overlapping ratios and the desired stitching length in $x$- and $y$-directions to cover the entire SUT. We can make sure at least one randomization solution is available by keeping the step number larger than the absolute value of the residual pixel number. Row-by-row or column-by-column vectors of potential step size are added with random variations. Because the random variations are generated in a certain range determined from the previous steps, the step sizes in each vector are not absolutely random. Consequently, we randomly sample the step sizes from one calculated step vector to make the resultant steps completely random.
With this regularization implemented in data acquisition, the ambiguous periodic errors created on SUT and reference by the self-calibration stitching algorithms can be significantly reduced.
3.2 Regularization for curvature and twist ambiguities
Despite the periodic errors when using uniform stitching steps, the curvature and twist are also ambiguous in the mathematical model described in Eq. (1). We consider three second order terms ($x$-curvature, $y$-curvature and twist) adding to $z(x,y)$ and $r(x,y)$ to compose the new shapes of SUT $z^*(x,y)$ and reference $r^*(x,y)$ as
where the coefficients $C_{11}$, $C_{12}$, and $C_{22} \in \mathbb {R}$ can be arbitrary values. The linear terms $a_{n} x+b_{n} y+c_{n}$ in Eq. (1) are replaced by $a_{n}^* x+b_{n}^* y+c_{n}^*$ whose coefficients are The measured data acquired by the interferometer will be given byUnlike the periodic errors in Section 1 coupling with SUT and reference which occurs under uniform step condition only, the curvature and twist ambiguity always exists in the stitching model and couples not only with the curvature and twist on SUT and reference but also with the mechanical motions. Since this ambiguity exists without any extra hypothesis, it implies that the information from the stitching data acquisition is not enough to solve this ambiguity problem. A possible regularization can be a better knowledge of the reference curvature and twist or the mechanical motion. For example, the real curvature and twist of the reference can be utilized in the stitching algorithm, or the motion angles and straightness can be traced to reduce the uncertainties of the curvature and twist. In this work, our stitching algorithm is regularized to not estimate the piston, $x$-tilt, $y$-tilt, $x$-curvature, twist, and $y$-curvature for the reference shape.
The estimation of $\mathbf {r}$ in Eq. (3) is implemented by using the extended matrices $\mathbf {G_{ex}}$ and $\mathbf {D_{ex}}$, as well as the extended vector $\mathbf {d_{ex}}$ as
Similarly, the estimation of $\mathbf {\hat {r}}_i$ in Eq. (5) becomes
The estimation of $\mathbf {\hat {r}}_i$ in Eq. (7) can be regularized as
The regularization in Eqs. (19)–(21) enables a highly flexible option to add the curvature and twist acquired by another system. The high-order additive systematic errors are estimated with the redundant data.
4. Simulation-based numerical case study
A stitching measurement of a $40$mm $\times$ $10$mm sample is simulated to verify the proposed iterative self-calibration algorithms. The height distribution of the SUT in Fig. 4(a) is generated as
The SA of the interferometer is $64\times 64$ pixels with $0.12$ mm lateral resolution. The user-defined overlapping ratio is set as $80\%$ for both directions. As discussed in Section 1, randomized step sizes are performed to resolve the periodic errors. The range for the randomization is $\pm 3$ pixels. In total, $3\times 23$ subsets of interferometric height maps are simulated with the reference $\mathbf {r}$ shown in Fig. 4(b), $1$ $\mu$rad RMS angular adjustment in geometric parameters $\mathbf {t}$, and $0.5$ nm RMS additive normally distributed random noises $n_a$. The CS-P, CS-S, and SC algorithms are performed to demonstrate the performance of the proposed 2D self-calibration stitching.
The non-iterative CS-P result is shown in Fig. 5. The stitching result in Fig. 5(a) and estimated reference in Fig. 5(b) are very close to the simulated SUT and reference in Fig. 4. The up to the 2nd order removed stitching error in Fig. 5(c) is about 0.26 nm RMS, and corresponding reference estimation error is 0.07nm RMS as shown in Fig. 5(d).
The results using CS-S algorithm are shown in Fig. 6. The first estimated reference $\mathbf {r}_1$ in Fig. 6(b) is similar to the simulated truth with only 0.11 nm RMS in Fig. 6(c). However, low-frequency variations are shown on the stitching error in Fig. 6(d). With the compensation $\mathbf {r}_2 = \mathbf {r}_1 + \mathbf {\hat {r}}_1$ in the first iteration, the reference estimation error is 0.07 nm RMS in Fig. 6(g) and the stitching error is only 0.21 nm RMS in Fig. 6(h). Since the RMS value of the next compensation amount $\mathbf {\hat {r}}_2$ in Fig. 6(i) is less than the preset threshold $thr=0.01$ nm, the iteration terminates with the final stitching error shown in Fig. 6(l). In fact, the stitching calculations in iterations shown in Figs. 6(d) and 6(h) are not necessary in CS-S algorithm, and only the reference is estimated through iterative compensation.
By applying the SC algorithm, the reference is estimated during the iterative stitching process in Fig. 7. As shown in Figs. 7(a), 7(e), and 7(i), the RMS value of the up to the 2nd order removed stitching error is getting lower, as the reference $\mathbf {r}$ is updated during the stitching process. The stitching error with repeating patterns shown in Fig. 7(a) is very typical if the reference is not well calibrated. The reference compensation amount $\mathbf {\hat {r}}_0$ in Fig. 7(b) is estimated from the stitched result $\mathbf {z}_0$ in this step and the raw SA measures. The updated reference $\mathbf {r}_1 = \mathbf {r}_0 +\mathbf {\hat {r}}_0$ in Fig. 7(c) has a reference estimation errors (terms up to the 2nd order removed) only at the sub-nanometer level with a period correlated to the average step size. In the 1st iteration, the stitching errors are reduced from $3.14$ nm RMS down to $0.29$ nm RMS, which also indicates that the previous reference estimation are effective. The remaining errors are small but still evident with a periodic shape shown in Fig. 7(e). The estimated $\mathbf {\hat {r}}_1$ in Fig. 7(f) is very similar to the previous reference estimation error in Fig. 7(d). After updating the reference $\mathbf {r}_2 = \mathbf {r}_1 +\mathbf {\hat {r}}_1$, the high-order discrepancy from the true value is almost random noise as shown in Fig. 7(h). A small random noise like stitching error(terms up to the 2nd order removed) can be expected in the next iteration as demonstrated in Fig. 7(i) with $0.2$ nm RMS. The stitching algorithm converges after the 2nd iteration, as the reference compensation amount $\mathbf {\hat {r}}_2$ in Fig. 7(j) is less than the preset RMS threshold $thr=0.01$ nm.
The simulation demonstrates that the proposed method is effective in estimating the reference $\mathbf {r}$ and stitching the SUT shape $\mathbf {z}$. Next, we will show a real stitching interferometry measurement to verify the feasibility of the proposed method.
5. Experimental verification
Some experiments are implemented on the interferometric stitching platform developed at NSLS-II [23] to test and verify the proposed method. The setup of the stitching interferometer (see Fig. 8) is composed of a Fizeau interferometer and several translation and rotation stages to adjust the relative position and angle between the SUT and the interferometer reference surface. As marked in Fig. 8, we only need the $x$- and $y$-motions, $\theta _x$ and $\theta _y$ rotations to perform 2D scans and to null the fringes at each subset location in this experiment.
The SUT is a $190$-mm-long flat silicon mirror. The central $180$ mm $\times$ $10$ mm area is the region of interest. We use a $256\times 64$ pixels rectangular mask window as the interferometer SA, and the pixel lateral resolution is $0.12$ mm. The user-defined overlap ratios are set as $80\%$ for both $x$- and $y$-directions. The software algorithm automatically calculates the random step variation to cover the whole region of interest. It ensures that the random step size variations are within the preset limits ($[-2, +2]$ pixels) and the sizes are in integer pixel to avoid the sub-pixel matching issue. As a result, the actual overlap ratio is about $81\%$ in $x$-direction and $85\%$ in $y$-direction. In total, $3\times 26$ subsets are captured to fully cover the region of interest. We take 64-averaged measurement at each SA. The repeatability of a 64-averaged measurement in our environment is about $0.2$ nm RMS. The tip-tilt threshold is set as $1$ $\mu$rad to null fringes before capturing each subset. The RMS threshold for reference compensation $thr =0.01$ nm.
In order to make a quantitative evaluation of the proposed self-calibration, a well-polished $240$mm$\times 40$mm flat X-ray mirror with about $0.6$ nm RMS is used to calibrate the reference by averaging several uncorrelated regions. This calibrated reference (Fig. 9(a)) is used as the benchmark to judge the estimated reference during iterations. The stitched SUT shape in Fig. 9(b) is used as the benchmark to evaluate the stitching errors.
5.1 Self-calibration stitching results
We present self-calibration stitching results with CS-P and SC algorithms. The CS-P estimated reference is shown in Fig. 10(a) with a $0.11$ nm RMS high-order reference estimation error in Fig. 10(b). With the reference $\mathbf {r}$, the stitched height $\mathbf {z}$ is shown in Fig. 10(c).
The iterative stitching and reference estimation results with the SC algorithm are shown in Fig. 11. The calculation starts with $\mathbf {r}_0=0$. The stitched $\mathbf {z}_0$ shows repeating error patterns due to the lack of knowledge on the reference. It is obvious that their period is about the averaged step size.
The reference compensation amount $\mathbf {\hat {r}}_0$ is estimated from the difference between $\mathbf {m}$ and $\mathbf {z}_0$. The main topography of the reference $\mathbf {r}_1=\mathbf {r}_0+\mathbf {\hat {r}}_0$ is updated via this estimation. However, comparing to the calibrated reference, there are periodic errors shown in Fig. 11(d), which results in obvious waviness errors on the stitched height map in the next iteration shown in Fig. 11(e). In SC iteration 1, The reference compensation amount $\mathbf {\hat {r}}_1$ is estimated as shown in Fig. 11(f). After this compensation, the discrepancy between the updated reference $\mathbf {r}_2 =\mathbf {r}_1+\mathbf {\hat {r}}_1$ and the calibrated reference benchmark is very small (only a $0.11$ nm RMS) with no obvious periodic error shown in Fig. 11(h). The periodic patterns on the stitched $\mathbf {z}_2$ is reduced (see Fig. 11(i)) in the SC iteration 2. Finally, the calculated $\mathbf {\hat {r}}_2$ in Fig. 11(j) is less than the preset RMS threshold $thr=0.01$ nm. The iteration ends up with the final stitched SUT height $\mathbf {z}_2$ and the estimated reference $\mathbf {r}_2$. The final results of the SC algorithm in Fig. 11 is very similar to the result of the CS-P algorithm in Fig. 10. The following studies are based on results of the SC algorithm.
5.2 Repeatability study
To study the repeatability of the proposed self-calibration stitching method, 10 repeating scans are performed. These 10 SC algorithm stitched SUT maps are shown in Fig. 12(a). To highlight the details of the height variation on the SUT, the color range is only within $[-3\sigma , +3\sigma ]$ where $\sigma$ is the RMS value of the stitched height values.
Illustrated with the same color map, Fig. 12(b) shows their discrepancies from the average map. All the discrepancies give an overall repeatability value less than $0.1$ nm RMS.
The 10 estimated reference maps are displayed in Fig. 13(a) with only $0.015$ nm RMS overall repeatability in their discrepancies shown in Fig. 13(b), which illustrates excellent repeatability of the reference estimation.
5.3 Comparison with the stitching result using the calibrated reference
The stitched SUT height with the proposed SC algorithm (Fig. 11(i)) is evaluated by comparing with the stitching result using the calibrated reference as benchmark (Fig. 9). The stitching result without subtracting the systematic error (Fig. 11(a)) is also evaluated as a contrast. The comparison results are shown in Fig. 14.
The obvious periodic error in Fig. 14(a) is one of the drawbacks to use a stitching mechanism with systematic errors. Our proposed self-calibration method provide a way to in situ estimate the repeatable high-order additive systematic error during the analysis of the measurement data. As shown in Fig. 14(b), the artificial waviness diminishes in our stitched results without any pre-knowledge of the systematic error.
5.4 Self-consistency study using different SA windows
We perform another measurement of the same SUT using a new SA with the same size ($256\times 64$ pixels) but in a different location on our 4-inch pupil interferometer. The previous SA is labelled as SA$_1$ and the new one as SA$_2$.
Comparing to SA$_1$, the SA$_2$ window is shifted $128$ pixels to the right as shown in Fig. 15. This shift will allow us to check if the self-calibration stitching method will give self-consistent reference results in their common area SA$_1\cap$ SA$_2$.
From Fig. 16, it is not difficult to notice that the stitched SUT shapes using different SAs are very close to each other with only $0.12$ nm RMS height difference on the $180$ mm $\times$ $10$ mm surface area. On the other hand, the estimated reference maps in Fig. 17 also indicate the self-consistency of the stitching method as the reference maps in the common area of the two SAs are similar to each other with only $0.10$ nm RMS difference, and this difference is dominated by the second order terms.
By removing the 2nd order terms, the difference inside the common area drops to $0.03$ nm RMS, which shows excellent self-consistency of the proposed method. Since the terms up to the 2nd order are not involved in the estimation as mentioned in Eq. (15), the small curvature difference in the common area between the two estimated reference maps is due to the different height variations in the non-overlapping regions of these two SAs which result in different curvature values in each SA.
6. Discussion
We have proposed two types of 2D self-calibration stitching algorithms (the CS and SC algorithms) for stitching interferometry. In this work, the calculation is fully based on pixel-to-pixel relations, since the stitching steps are exact integer pixel numbers. However, the proposed algorithms can be easily extended and applied to an interpolation representation of the reference surface with a linear combination of basis functions when the stitching steps are not integer pixels and the interpolation is unavoidable. In that case, the relation matrices in the CS-P and CS-S algorithms are no longer that sparse which could yield heavier computation issues.
The merits and limitations of the proposed CS-P, CS-S, and SC algorithms are discussed below.
In pixel-relation-based stitching strategy, the CS-P algorithm directly estimate the reference and stitched height at the same time. Although it can provide a direct least squares solution with LSMR() function for data with a reasonable size shown in the experiment, it needs to pay careful attentions to the tolerances in LSMR() and the weights of the extended equations to get accurate stitching and estimation results.
Different from the CS-P algorithm, the CS-S and SC algorithms separate the reference estimation and height stitching, which reduces the computational dimension of the problem. The CS-S algorithm estimates the reference first and then uses it for shape stitching. In CS-S algorithm flow, the reference estimation has less influence from the stitched SUT result, and it has a better efficiency in utilizing the estimated reference comparing to the SC algorithm. The merit of SC algorithm is its capability of keeping the size of the relation matrix $\mathbf {D^s}$ in Eq. (6) fixed as $M\times M$, which does not depend on the subset number $N$. Therefore, it has a unique advantage in dealing with stitching with large number of subsets. It is worth to note that the computation in assembling the matrix $\mathbf {D^s}$ can be time consuming for a limited computing resources, if $M$ becomes a huge number.
7. Conclusion
In our work, the repeatable high-order additive systematic errors can be estimated from the redundant data acquired in $x$- and $y$-directions stepping with big overlaps of neighboring subsets. The periodic errors due to the uniform steps in data acquisition are carefully studied. Randomized nonuniform step sizes are suggested to “break” the artificial periodicity during the data acquisition. Regularization on curvature and twist are introduced in the algorithms to potentially preserve the measured reference curvature and twist. The proposed algorithms have been cross-checked and verified using both simulation and experimental data. Further study on the memory cost and speeding up the computation will be our future work to resolve this practical issue.
Funding
Office of Science (DE-SC0012704).
Acknowledgments
The authors want to thank Dennis Kuhne from the Research & Development Group of NSLS-II for the mechanical machining.
This research used resources of the National Synchrotron Light Source II, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Brookhaven National Laboratory under Contract No. DE-SC0012704.
References
1. P. de Groot, “Principles of interference microscopy for the measurement of surface topography,” Adv. Opt. Photonics 7(1), 1–65 (2015). [CrossRef]
2. L. Deck and P. de Groot, “High-speed noncontact profiler based on scanning white-light interferometry,” Appl. Opt. 33(31), 7334–7338 (1994). [CrossRef]
3. P. Murphy, G. Forbes, J. Fleig, P. Dumas, and M. Tricard, “Stitching interferometry: a flexible solution for surface metrology,” Opt. Photonics News 14(5), 38–43 (2003). [CrossRef]
4. J. Fleig, P. Dumas, P. E. Murphy, and G. W. Forbes, “An automated subaperture stitching interferometer workstation for spherical and aspherical surfaces,” in Advanced Characterization Techniques for Optics, Semiconductors, and Nanotechnologies, vol. 5188 (International Society for Optics and Photonics, 2003), pp. 296–307.
5. M. Tricard, A. Kulawiec, M. Bauer, G. DeVries, J. Fleig, G. Forbes, D. Miladinovich, and P. Murphy, “Subaperture stitching interferometry of high-departure aspheres by incorporating a variable optical null,” CIRP Ann. 59(1), 547–550 (2010). [CrossRef]
6. H. Mimura, H. Yumoto, S. Matsuyama, K. Yamamura, Y. Sano, K. Ueno, K. Endo, Y. Mori, M. Yabashi, K. Tamasaku, Y. Nishino, T. Ishikawa, and K. Yamauchi, “Relative angle determinable stitching interferometry for hard x-ray reflective optics,” Rev. Sci. Instrum. 76(4), 045102 (2005). [CrossRef]
7. S. Chen, S. Li, and Y. Dai, “Iterative algorithm for subaperture stitching interferometry for general surfaces,” J. Opt. Soc. Am. A 22(9), 1929–1936 (2005). [CrossRef]
8. A. Vivo, B. Lantelme, R. Baker, and R. Barrett, “Stitching methods at the European Synchrotron Radiation Facility (ESRF),” Rev. Sci. Instrum. 87(5), 051908 (2016). [CrossRef]
9. L. Huang, J. Xue, B. Gao, and M. Idir, “One-dimensional angular-measurement-based stitching interferometry,” Opt. Express 26(8), 9882–9892 (2018). [CrossRef]
10. L. Huang, M. Idir, C. Zuo, T. Wang, K. Tayabaly, and E. Lippmann, “Two-dimensional stitching interferometry based on tilt measurement,” Opt. Express 26(18), 23278–23286 (2018). [CrossRef]
11. C. J. Evans, R. J. Hocken, and W. T. Estler, “Self-calibration: Reversal, redundancy, error separation, and ’absolute testing’,” CIRP Ann. 45(2), 617–634 (1996). [CrossRef]
12. A. Vivo, R. Barrett, and F. Perrin, “Stitching techniques for measuring x-ray synchrotron mirror topography,” Rev. Sci. Instrum. 90(2), 021710 (2019). [CrossRef]
13. F. Polack and M. Thomasset, “Determination and compensation of the “reference surface” from redundant sets of surface measurements,” Nucl. Instrum. Methods Phys. Res., Sect. A 710, 67–71 (2013). The 4th international workshop on Metrology for X-ray Optics, Mirror Design, and Fabrication. [CrossRef]
14. F. Polack, M. Thomasset, S. Brochet, and D. Dennetiere, “Surface shape determination with a stitching michelson interferometer and accuracy evaluation,” Rev. Sci. Instrum. 90(2), 021708 (2019). [CrossRef]
15. J. Nicolas, M. L. Ng, P. Pedreira, J. Campos, and D. Cocco, “Completeness condition for unambiguous profile reconstruction by sub-aperture stitching,” Opt. Express 26(21), 27212–27220 (2018). [CrossRef]
16. M. Bray, “Stitching interferometry: the practical side of things,” in Optical Manufacturing and Testing VIII, vol. 7426 (International Society for Optics and Photonics, 2009), p. 74260Q.
17. P. Su, J. H. Burge, and R. E. Parks, “Application of maximum likelihood reconstruction of subaperture data for measurement of large flat mirrors,” Appl. Opt. 49(1), 21–31 (2010). [CrossRef]
18. S. Chen, Y. Dai, S. Li, X. Peng, and J. Wang, “Error reductions for stitching test of large optical flats,” Opt. Laser Technol. 44(5), 1543–1550 (2012). [CrossRef]
19. S. Chen, W. Liao, Y. Dai, and S. Li, “Self-calibrated subaperture stitching test of hyper-hemispheres using latitude and longitude coordinates,” Appl. Opt. 51(17), 3817–3825 (2012). [CrossRef]
20. D. Golini, G. Forbes, and P. Murphy, “Method for self-calibrated sub-aperture stitching for surface figure measurement,” (2005). US Patent 6,956,657.
21. M. Vannoni and G. Molesini, “Iterative algorithm for three flat test,” Opt. Express 15(11), 6809–6816 (2007). [CrossRef]
22. M. Vannoni and G. Molesini, “Absolute planarity with three-flat test: an iterative approach with zernike polynomials,” Opt. Express 16(1), 340–354 (2008). [CrossRef]
23. L. Huang, T. Wang, K. Tayabaly, D. Kuhne, W. Xu, W. Xu, M. Vescovi, and M. Idir, “Stitching interferometry for synchrotron mirror metrology at National Synchrotron Light Source II (NSLS-II),” Opt. Lasers Eng. 124, 105795 (2020). [CrossRef]
24. D. C.-L. Fong and M. Saunders, “Lsmr: An iterative algorithm for sparse least-squares problems,” SIAM J. on Sci. Comput. 33(5), 2950–2971 (2011). [CrossRef]
25. D. C.-L. Fong and M. Saunders, “LSMR: Sparse Equations and Least Squares,” http://web.stanford.edu/group/SOL/software/lsmr/.