Advanced
Fast Sub-aperture Stitching Algorithm Using Partial Derivatives
Fast Sub-aperture Stitching Algorithm Using Partial Derivatives
Journal of the Optical Society of Korea. 2015. Feb, 19(1): 84-87
Copyright © 2015, Optical Society of Korea
  • Received : October 10, 2014
  • Accepted : December 12, 2014
  • Published : February 25, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Yiwei Chen
University of Chinese Academy of Sciences, Beijing 100049, China
cyw198788@163.com
Erlong Miao
State Key Laboratory of Applied Optics, Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun, Jilin 130033, China
Yongxin Sui
State Key Laboratory of Applied Optics, Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun, Jilin 130033, China
Huaijiang Yang
State Key Laboratory of Applied Optics, Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun, Jilin 130033, China
Abstract
For large optical elements which are tested by many sub-apertures, it takes too much time for a sub-aperture stitching algorithm to get the stitching result. To solve this problem, we propose a fast sub-aperture stitching algorithm to quickly compensate for piston, tilt, and defocus errors. Moreover, the new algorithm is easy to understand and program. We use partial derivatives of measurement data to separately solve piston, tilt, and defocus errors. First, we show that the new algorithm has a lower time complexity than the currently used algorithm. Although simulation results indicate that the accuracy of the new algorithm is lower than the current algorithm in all 20 simulations, our experimental results validate the algorithm and show it is sufficiently accurate for general use.
Keywords
I. INTRODUCTION
With the development of astronomical optics, space optics [1 - 3] and inertial confinement fusion, large-aperture optical systems have been widely used. For very large optical elements, however, the aperture may be too large for the interferometer to test in one pass. In 1982, C. J. Kim solved this problem with the sub-aperture stitching method [4] . However, location errors are inevitable with this method. Thus, researchers use overlapping areas between the sub-apertures to compensate for location errors [5 - 9] such as piston, tilt, defocus, clocking, and position errors. In most cases, only piston, tilt, and defocus errors are corrected. For large optical elements which are tested by many sub-apertures, it takes much time for a sub-aperture stitching algorithm to get the stitching result. To quickly compensate for piston, tilt, and defocus errors, we develop a fast algorithm, which is easy to understand and program.
In Section 2, we briefly introduce our algorithm. In Section 3, we compare the time complexity of our algorithm with the current algorithm. In Section 4, we study the accuracy of our algorithm using simulation. In Section 5, we discuss an experimental result demonstrating the effectiveness of the proposed algorithm. Finally, we present our conclusions in Section 6.
II. FAST ALGORITHM FOR SUB-APERTURE STITCHING
The current algorithm uses a least-squares method to remove relative piston, tilt and defocus errors directly. In our algorithm, we calculate the partial derivative of data before using a least-squares method to remove relative piston, tilt and defocus errors. The reason for calculating the partial derivative of data is that it makes the least-squares method much easier. Actually, using the least-squares method after calculating the partial derivative of data is just averaging, so it makes the algorithm faster than current algorithm.
For convenience, we only discuss the testing of spherical optical elements in this article; however, with some modifications, this algorithm can be easily applied to the testing of both flat and aspheric optical elements.
One sub-aperture, usually the center sub-aperture, is selected as the reference sub-aperture. The remaining sub-apertures have nonzero piston, tilt, and defocus errors that are relative to the reference sub-aperture. Thus, we compensate for relative piston, tilt and defocus errors [9] as follows:
PPT Slide
Lager Image
where x and y are the coordinates of pixels; z i ( x , y ) is the testing value measured for the i-th sub-aperture; zi * ( x , y ) is the value after compensation; a i and bi are the coefficients of the relative tilt of the i-th sub-aperture in the x and y directions, respectively, ci is the coefficient of the relative piston value of the i-th sub-aperture, and di is coefficient of relative defocus.
Assuming that there are N sub-apertures and that the reference sub-aperture is number 1, current sub-aperture stitching algorithms use the least-squares method to directly calculate the necessary compensations [5 - 9] :
PPT Slide
Lager Image
In our algorithm, we first calculate the second-order partial derivative of Equation 1 with respect to the independent variable x:
PPT Slide
Lager Image
Then, we use the least-squares method to calculate the needed compensation between two adjacent sub-apertures:
PPT Slide
Lager Image
Solving Equation 4, we obtain the following over-determined linear equations:
PPT Slide
Lager Image
As we can see, using a least-squares method after calculating the partial derivative of data is just averaging. Here, N ij is the total number of pixels in the overlapping area of the i-th sub-aperture and the j-th sub-aperture, and N k is the total number of pixels in the overlapping area of 1-st sub-aperture and k-th sub-aperture.
We can easily use the weighted least-squares method to solve the over-determined linear equations by setting the weight equal to the total number of pixels in the corresponding overlapping area.
Using the calculated coefficients of the relative defocus, we compensate the defocus errors:
PPT Slide
Lager Image
where fi ( x , y ) is the data of the i-th sub-aperture after compensating for the defocus error. According to Equation 2 and Equation 6, we obtain
PPT Slide
Lager Image
To compensate for the relative tilt errors, we calculate the first-order partial derivative of Equation 7 with respect to independent variable x and the independent variable y:
PPT Slide
Lager Image
PPT Slide
Lager Image
We use a similar procedure to calculate ai and bi . Then, we compensate for the tilt errors:
PPT Slide
Lager Image
where gi ( x , y ) is the data of i-th sub-aperture after compensating for the tilt and defocus errors.
Finally, we use the same method to compensate for the piston errors and stitch the resulting measurement data to obtain a full-aperture result.
III. TIME COMPLEXITY OF THE ALGORITHM
In this section, we assume that there are only two sub-apertures, where one sub-aperture is the reference sub-aperture, and the overlapping area is a square (n × n pixels).
In the algorithm currently in use, Equation 2 leads to linear equations:
PPT Slide
Lager Image
where ∆( x , y )is the disparity between the two measurements ( z 1 ( x , y )- z 2 ( x , y )).
Solving Equation 11, we obtain the coefficients a , b , c , and d , and we can use Equation 1 to compensate for the errors. In the current algorithm, we use about 12 n 2 multiplication operations and about 18 n 2 addition operations (shown in Table 1 ).
Time complexity
PPT Slide
Lager Image
Time complexity
In our new algorithm, we use differences to calculate partial derivatives. For instance, we find the first-order partial derivative using
PPT Slide
Lager Image
where ∆ x is the distance between two adjacent pixels. In our method, we do not divide ( fi ( x + ∆ x , y ) - fi ( x , y )) + fi ( x , y ) - fi ( x -∆ x , y )) by 2∆ x for each point. For instance, when we obtain
PPT Slide
Lager Image
we calculate it using
PPT Slide
Lager Image
Thus, we use about 5 n 2 multiplication operations and about 12 n 2 addition operations in our algorithm (shown in Table 1 ).
In Table 1 , we can see that our algorithm has a much lower time complexity than the current algorithm.
IV. SIMULATION
In the simulation, a measurement result (501 × 801 pixels, shown in Fig. 1 ) was first divided into two sub-apertures. After dividing the measurement result into sub-apertures, we add random noise ([-0.001, 0.001] wavelength) and piston, tilt, and defocus errors to each sub-aperture. We use both the current algorithm and the new algorithm to calculate the compensators. We repeated the simulation calculations 20 times. In all 20 simulations, the accuracy of the new algorithm is lower than that of the current algorithm. We present one typical result in Table 2 .
PPT Slide
Lager Image
Measurement result used for the simulation. All units are in wavelengths.
Simulation results
PPT Slide
Lager Image
Simulation results
V. EXPERIMENT
In this section, we test a convex spherical surface (caliber of 100 mm and radius of 76.56 mm) using a sub-aperture stitching interferometer (a Fizeau interferometer whose working wavelength is 632.8 nm) with a reference surface whose F-number is 1.5. A real system including a real sample is shown in Fig. 2 . The layout of the sub-apertures is shown in Fig. 3 . The measurement results using our algorithm are shown in Fig. 4(a) . Comparing with the measurement results obtained using a full-aperture interferometer (shown in Fig. 4(b) ), the root mean squared error was 0.003014 λ and the residual error image is shown in Fig. 5 .
PPT Slide
Lager Image
Real system including real sample.
PPT Slide
Lager Image
Layout of the sub-apertures.
PPT Slide
Lager Image
Experimental results: (a) stitched using our algorithm; (b) full-aperture measurement.
PPT Slide
Lager Image
Residual error image.
VI. CONCLUSION
We propose a fast algorithm for sub-aperture stitching to quickly compensate for piston, tilt, and defocus errors. We prove that our algorithm runs more than two times faster than the current algorithm. Although the simulation results demonstrate that the accuracy of the new algorithm is lower than that of the current algorithm in all 20 simulations, the new algorithm is sufficiently accurate for generally use. Using an experimental measurement, we verified the validity of the new algorithm and demonstrated successful sub-aperture stitching.
References
Li X. L. , Xu M. , Ren X. D. , Pei Y. T. 2012 “An optical design of off-axis four-mirror-anastigmatic telescope for remote sensing,” J. Opt. Soc. Korea 16 243 - 246
Li X. L. , Xu M. , Pei Y. T. 2012 “Optical design of an off-axis five-mirror-anastigmatic telescope for near infrared remote sensing,” J. Opt. Soc. Korea 16 343 - 348
Jin H. , Lim J. , Kim Y. , Kim S. 2013 “Optical design of a reflecting telescope for cubesat,” J. Opt. Soc. Korea 17 533 - 537
Kim C. J. 1982 “Polynomial fit of interferograms,” Appl. Opt. 21 4521 - 4525
Otsubo M. , Okada K. , Tsujiuchi J. 1992 “Measurement of large plane surface shape with interferometric aperture synthesis,” Proc. SPIE 1720 444 - 447
Otsubo M. , Okada K. , Tsujiuchi J. 1994 “Measurement of large plane surface shapes by connecting small-aperture,” Opt. Eng. 33 608 - 613
Sjödahl M. , Oreb B. F. 2002 “Stitching interferometric measurement data for inspection of large optical components,” Opt. Eng. 41 403 - 408
Golini D. , Forbes G. W. , Murphy P. E. 2005 “Method for self-calibrated subaperture stitching for surface figure measurement,” U.S. Patent 6956657B
Fleig J. , Dumas P. , Murphy P. E. , Forbes G. W. 2003 “An automated subaperture stitching interferometer workstation for spherical and aspherical surfaces,” Optical Science and Technology, SPIE’s 48th Annual Meeting. International Society for Optics and Photonics 296 - 307