Instrumentation, Techniques, and Measurement

Optical surface evaluation by correlating bi-Ronchigram images

[+] Author Affiliations
Alberto Cordero-Dávila

Benemérita Universidad Autónoma de Puebla, Facultad de Ciencias Físico Matemáticas, Avenida San Claudio y Río Verde s/n, Col. San Manuel, Puebla Pue. C. P. 72570, México

Jorge González-García

Universidad Tecnológica de la Mixteca, Instituto de Física y Matemáticas, km 2.5 Carretera a Acatlima, Huajuapan de León, Oaxaca C. P. 69000, México

Opt. Eng. 54(3), 034108 (Mar 16, 2015). doi:10.1117/1.OE.54.3.034108
History: Received November 22, 2014; Accepted February 26, 2015
Text Size: A A A

Open Access Open Access

Abstract.  By correlating bi-Ronchigram images, surface errors without any supposed symmetry are evaluated. No approximation and interference orders are needed. First, only one experimental bi-Ronchigram (Ronchigram with a square grid) image was recorded. Second, given a surface parameter set (curvature radius, conic constant, and/or symmetric and asymmetric deformation coefficients), a bi-Ronchigram image is simulated and correlated with the bi-Ronchigram experimental image. Third, genetic algorithms are used to find the parameters for which the correlation coefficient reaches its maximum value. Finally, the parameters of the experimental surface are estimated. Evaluations of curvature radius, conic constant, and error functions for reflecting surfaces will be shown.

Figures in this Article

In optical shops, experimental and simulated Ronchigrams1,2 are visually compared in order to qualitatively evaluate reflecting surfaces. On the other hand, to quantitatively evaluate the aberrations of any asymmetric surface (or wavefront), two crossed Ronchigrams have to be recorded.3,4 From each one, a component of the transversal aberration function is calculated.57 After an analytical integration is applied3,4 and the optical path differences (OPD) function is calculated, it is possible to record only one pattern (bi-Ronchigram) if a square grid is used instead of the classical Ronchi ruling.8 In this latter case, a numerical8 or analytical integration9 can be applied to estimate the OPD function. To evaluate phase function from a pattern of fringes,10 the correlation coefficient between the experimental and simulated interferogram is maximized by using genetic algorithms. Recently, a similar procedure has been used to evaluate Ronchigrams.11

What follows in this paper is the description of a procedure to estimate parameters and error function of any kind of mirror. By using the experimental setup shown in Fig. 1, we recorded the experimental bi-Ronchigram image, which was correlated with several simulated images of bi-Ronchigrams. Before any bi-Ronchigram was recorded, the ruling was aligned in order to guarantee the parallelism between (1) ruling lines and CCD array, (2) the normal of the CCD and ruling plane normal, and (3) the normal of the ruling plane and optical axis of the mirror. For each theoretical surface used for simulation, the error function was a different one until a correlation coefficient was maximized. As a result, the parameters of the experimental surface were estimated.

Graphic Jump Location
Fig. 1
F1 :

(a) Experimental setup of Ronchi test and (b) CCD camera, Ronchi tester (λ=660nm) and square grid.

In Sec. 2, the problem is stated. In Sec. 3, the algorithm to simulate bi-Ronchigram images is described. Algorithms used to normalize experimental bi-Ronchigram images and to correlate simulated and experimental bi-Ronchigrams are described in Sec. 4. Experimental results are reported in Sec. 5. Finally, our conclusions are reported in Sec. 6.

Any surface can be described by means of a vector of parameters, a (curvature radius, conic constant, symmetric and asymmetric deformation constants, etc.). The experimental surface is defined by an unknown experimental vector of parameters, aex, and the surface used to simulate a bi-Ronchigram image is defined by its parameters vector, asi. By correlating the experimental and simulated bi-Ronchigram images, a correlation coefficient, C(aex,asi), is obtained. We used genetic algorithms to find the asi* value for which the correlation coefficient, C(aex,asi*), is maximized. We assume that the best estimator of aex is asi* and then the experimental surface parameters are evaluated. It is important to point out that no symmetry assumption, approximation, and interference orders are needed to apply this procedure.

Figure 2 shows the geometry of the Ronchi test with a square grid. An incident ray comes from a point source light, (α,β,γ), to a point, (x,y,z), on the surface. The reflected ray crosses12 the square grid plane, z=z0, at the point coordinates (x0,y0) given by Display Formula

x0=x+(z0z){(xα)(1zx2+zy2)2zx[zy(yβ)+(γz)](γz)(1zx2zy2)+2[zx(xα)+zy(yβ)]},(1)
Display Formula
y0=y+(z0z){(yβ)(1zx2+zy2)2zy[zx(xα)+(γz)](γz)(1zx2zy2)+2[zx(xα)+zy(yβ)]},(2)
where zx and zy represent partial derivatives (with respect to x and y) of the surface sagitta, z(x,y).

Graphic Jump Location
Fig. 2
F2 :

Geometry of Ronchi test with square grid.

The latter is given as the sum of the ideal sagitta, zid=(x,y), plus surface error, w, i.e., Display Formula

z(x,y)=zid(x,y)+w(x,y).(3)

In this application, w is the surface error measured along the z axis.

After the incident ray with unitary irradiance crosses the grid, its irradiance, I, can be calculated approximately by Display Formula

I=[1+cos(2πx0/d)2][1+cos(2πy0/d)2],(4)
where d is the period of the square grid along the X or Y direction. Even though Eq. (4) is an approximation, it is, however, very useful to match the experimental and simulated bi-Ronchigram images.13 The bi-Ronchigram image is obtained by assigning many values to points (x,y,z) on the mirror and then calculating I(x,y).

In order to evaluate the correlation coefficient, the point coordinates (x,y) where the simulated irradiance is calculated have to correspond to the point coordinates where the irradiance of the bi-Ronchigram image is recorded.

It is important to point out that Eqs. (1)–(3) can be applied to Hartmann’s test with no changes, i.e., from a mathematical point of view, there is no difference between calculations for the Ronchigram and the hartmanngrams since the locations of the filtering and observation planes depend only on two parameters. In fact, a common mathematical model has been established for Ronchi’s and Hartmann’s tests. Even from the physical point of view, the Ronchigram and Hartmanngram correspond to virtual and real patterns, respectively.12

Normalization of Bi-Ronchigram Images

The file of the experimental bi-Ronchigram image is given by the triad set: (xi, yi, Ii), i=1,2,NP, where NP is the number of pixels of the CCD and Ii is the measured irradiance at the pixel center coordinates (xi, yi).

First of all, by applying the least-squares method to data points of the bi-Ronchigram border, the center and radius of the bi-Ronchigram image are evaluated.14 And then a new triad set (xci, yci, IEci), i=1,2,,NPC, is evaluated, where (xci, yci) are the NPC point coordinates with respect to the center of the experimental bi-Ronchigram image and within a unitary circle.

Second, the Zernike polynomials are used to normalize the experimental bi-Ronchigram image,15 see Fig. 3, i.e., to calculate the set of triads (xECNi, yECNi, IECNi), see Fig. 4, where IECNi represents the normalized irradiances at the i’th point. In Fig. 5, the irradiance plots are shown along the x axis of the bi-Ronchigram before and after the normalization algorithm had been applied.

Graphic Jump Location
Fig. 3
F3 :

Experimental bi-Ronchigram image for mirror of 14 cm diameter and 60.5 cm curvature radius.

Graphic Jump Location
Fig. 4
F4 :

Normalized experimental bi-Ronchigram image for mirror of 14 cm diameter and 60.5 cm curvature radius.

Graphic Jump Location
Fig. 5
F5 :

Bi-Ronchigram irradiance plots, along the x axis, (a) before and (b) after the normalization algorithm was applied.

Correlating Simulated and Experimental Bi-Ronchigram Images

As it has been pointed out, by using Eqs. (1)–(4), at the point coordinates (xECNi, yECNi), the irradiances, ISi, of the simulated bi-Ronchigram are calculated. Then the triad sets (xECNi, yECNi, ISi) are obtained.

For simulations of bi-Ronchigram images, the following two functions will be used. The sagitta function Display Formula

f=cρ21+1(k+1)c2ρ2,(5)
where ρ2=x2+y2. And the second equation is the known Kingslake polynomial Display Formula
g=txy+tyx+df(x2+y2)+as(x2+3y2)+coy(x2+y2)+sp(x2+y2)2.(6)

For the examples of Secs. 5.1 and 5.2, zid=f and w=g, i.e., the ideal surface is a conic mirror and the error function is described by the Kingslake function.

However, for the examples of Sec. 5.3, zid=0 and w=f. This means that the ideal surface is a plane one and the error function is a conic mirror, i.e., we want to know the curvature and conic constant of the mirror under test.

Digital image correlation then becomes a task of comparing the sets of triads (xECNi, yECNi, ISi) and (xECNi, yECNi, ICNi). The typical correlation function that measures how well a set matches is Display Formula

C=i=1Npc[(IECNiI¯ECN)(ISNiI¯SN)][i=1Npc(IECNiI¯ECN)2i=1Npc(ISNiI¯SN)2]1/2,(7)
where I¯SN and I¯ECN represent the mean values of ISNi and IECNi sets, respectively.

Figure 3 shows the experimental bi-Ronchigram obtained for a spherical mirror (under test) of 14 cm diameter and 60.50±0.05cm curvature radius, r. The latter was measured by locating the Ronchi ruling position for which the field on the mirror is either totally bright or dark.4 The (on axis) point source and the square grid of 42holes/in. (i.e., a period of 0.61 mm) were located at a distance 58.8 cm away from the mirror vertex.

Two different ideal surfaces were assumed: a spherical (Sec. 5.1) and a parabolic (Sec. 5.2) surface.

Spherical Mirror as Ideal Surface

We substituted k=0.0 and c=1/60.5cm into Eq. (5), and, by using genetic algorithms, calculated the coefficients of Eq. (6) for which the correlation coefficient of Eq. (7) reached its maximum value. We found that the correlation coefficient is maximized to a value of 0.9264 if coefficients are given by tx=3.73±0.05λ, ty=0.56±0.05λ, df=0.31±0.05λ, sp=0.00025±0.00005λ, c0=0.0000±0.0025λ, and as=0.03720±0.00025λ.

The tilt coefficients are different from zero since the dot of that zero order is decentered along the X and Y directions, as can be seen in Fig. 3. The spherical aberration value indicates a symmetrical error value of 3.6×105cm at the edge of the mirror. Finally, the astigmatic coefficient indicates that the mirror is not an axisymmetric surface. This can be seen from the experimental bi-Ronchigram shown in Fig. 3. The period of the dots along the X direction (13 dots) is different from the period along the Y direction (12 dots). The reproduced bi-Ronchigram is shown in Fig. 6.

Graphic Jump Location
Fig. 6
F6 :

Reproduced bi-Ronchigram from an ideal spherical mirror.

It is important to draw attention to the square dots of the simulated bi-Ronchigram (see Fig. 6), compared to the circular dots of the experimental bi-Ronchigram (see Fig. 3).

Parabolic Mirror as Ideal Surface

We substituted k=1.0 and c=1/60.5cm into Eq. (5). In this case, a maximal correlation coefficient of 0.9261 was reached for the coefficients tx=3.71±0.05λ, ty=0.56±0.05λ, df=0.32±0.05λ, sp=0.01148±0.00005λ, c0=0.0000±0.0025λ, and as=0.03700±0.00025λ. As can be seen, all of the coefficients are reproduced with the exception of the spherical aberration coefficient. This result was to be expected since the ideal surface is now a parabolic mirror. In addition, from the spherical aberration coefficient, a surface error of 1.40×103cm at the border of the mirror can be calculated. This result can be compared with the sagittae difference between the spherical and parabolic mirrors at their border (1.36×103cm).

The uncertainties written after each estimated coefficients correspond to the interval of seeking the maximum value of the correlation coefficient.

Evaluation of Curvature Radius and Conic Constant of Mirrors

Two mirrors of 20 and 14 cm diameters were tested in order to estimate, using our method, their curvature, c, and conic constant, k.

For the 20 cm mirror diameter, the point source and the squared grid (42holes/in.) were located at a distance 157.7 cm away from the mirror’s vertex. The experimental bi-Ronchigram obtained is shown in Fig. 7. After our algorithm was applied, we obtained r=160.171cm and k=0.981 for the paraxial curvature radius and the conic constant, respectively. The maximum value obtained for the correlation coefficient was 0.8113. The technician reported a paraxial curvature radius of 160.10±0.05cm and a conic constant k=1.0 (a parabolic mirror). As can be seen, the difference between nominal values and those estimated with our method is 2% of the conic constant and 0.05% of the paraxial curvature radius.

Graphic Jump Location
Fig. 7
F7 :

Experimental bi-Ronchigram of a parabolic mirror.

The second mirror with a 14 cm diameter with an internal hole of 2.7 cm diameter was tested with a squared grid of 33.02holes/in. located 57.2 cm away from the mirror’s vertex. The experimental bi-Ronchigram obtained is shown in Fig. 8(a). In Fig. 8(b), a plot of the irradiance along the x axis is shown. We applied our algorithm and we estimated the values of r=60.04cm and k=0.703 for the paraxial curvature radius and conic constant, respectively. The obtained maximum correlation coefficient was 0.7570. The technician in 16 reported a paraxial curvature radius of 60.10±0.05cm and a conic constant k=0.72 with our elliptical mirror. As can be seen, the difference between nominal values and those estimated with our method is 2.4% of the conic constant and 0.1% of the paraxial curvature radius.

Graphic Jump Location
Fig. 8
F8 :

(a) An experimental bi-Ronchigram image for an elliptical mirror and (b) experimental bi-Ronchigram irradiance along x axis.

A method to evaluate surface parameters and errors of surfaces, based on the correlation of bi-Ronchigram images, was presented. For this method, the experimental setup of the Ronchi test was used, and only a square grid instead of a Ronchi ruling was used. Thus, it is not necessary to purchase any extra hardware to implement this procedure. Added to this, no rotator is required to record two Ronchigrams. The required software is the one used for (1) interferograms normalization and (2) image correlations. It is important to point out that no symmetry assumption is required to apply this method.

De Vany  A. S., “Interpreting wave-front and glass-error slopes in an interferogram,” Appl. Opt.. 19, , 173  (1980). 0003-6935 CrossRef
De Vany  A. S., “Patterns correlation of interferograms and Ronchigrams,” Appl. Opt.. 20, , A40  (1981). 0003-6935 CrossRef
Cornejo  A., “Ronchi test,” in Optical Shop Testing. , Malacara  D., Ed., pp. 317 –360,  John Wiley and Sons ,  New York  (2007).
Cornejo-Rodríguez  A., and Cordero-Dávila  A., “Wavefront slope measurements in optical shop testing,” in Handbook of Optical Engineering. , Malacara  D., and Thompson  B. J., Eds., pp. 53 –71,  Marcel Dekker ,  New York  (2001).
Malacara-Doblado  D., and Ghozeil  I., “Hartman, Hartmann-Shack, and other screen tests,” in Optical Shop Testing. , Malacara  D., Ed., pp. 361 –397,  John Wiley and Sons ,  New York  (2007).
Takeda  M., , Ina  H., and Kobayashi  S., “Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry,” J. Opt. Soc. Am.. 72, , 156 –160 (1982). 0030-3941 CrossRef
Servin  M., , Malacara  D., and Cuevas  F. J., “Direct-phase detection of modulated Ronchi rulings using a phase-locked loop,” Opt. Eng.. 33, (4 ), 1193 –1199 (1994). 0091-3286 CrossRef
Cordero-Dávila  A.  et al., “Ronchi test with a square grid,” Appl. Opt.. 37, , 672 –675 (1998). 0003-6935 CrossRef
Cordero-Dávila  A.  et al., “Only one fitting for bironchigrams,” Appl. Opt.. 40, , 5600 –5609 (2001). 0003-6935 CrossRef
Vázquez-Montiel  S., , Sánchez-Escobar  J. J., and Fuentes  O., “Obtaining the phase of an interferogram by use of an evolution strategy: part I,” Appl. Opt.. 41, , 3448 –3452 (2002). 0003-6935 CrossRef
Aguirre-Aguirre  D.  et al., “Algorithm for ronchigram recovery with random aberrations coefficients,” Opt. Eng.. 52, (5 ), 053606  (2013). 0091-3286 CrossRef
Cordero-Dávila  A., , Cornejo-Rodríguez  A., and Núñez  O. C., “Ronchi and Hartmann tests with the same mathematical theory,” Appl. Opt.. 31, , 2370 –2376 (1992). 0003-6935 CrossRef
Cordero-Dávila  A., , Kantum-Montiel  J. R., and González-García  J., “Ronchigram simulations for free-form concave reflective surfaces,” Optik. 124, , 4892 –4895 (2013). 0030-4026 CrossRef
Cordero-Dávila  A., , Cardona-Núñez  O., and Cornejo-Rodríguez  A., “Least squares estimators for the center and radius of circular patterns,” Appl. Opt.. 32, , 5683 –5685 (1993). 0003-6935 CrossRef
Juarez-Salazar  R.  et al., “Generalized phase-shifting interferometry by parameter estimation with the least squares method,” Opt. Lasers Eng.. 51, , 626 –632 (2013). 0143-8166 CrossRef
Cordero-Dávila  A.  et al., “Local and global surface errors evaluation using Ronchi test, without both approximation and integration,” Appl. Opt.. 50, , 4817 –4823 (2011). 0003-6935 CrossRef

Alberto Cordero-Dávila studied his doctorate degrees in INAOE at Puebla. His has been a member of the National System of Researchers since 1989. He has coauthored several research and dissemination articles. He has been working for the past 43 years as a research lecturer at the Benemérita Universidad Autónoma de Puebla, Puebla, México. He has also been doing research on design, construction and testing of optical systems and is now a specialist in polishing techniques.

Jorge González-García received his master’s degree in optical instrumentation at BUAP and his PhD degree in optical polishing surfaces. He has coauthored studies on optical design, adaptive optics, mechanical design, optical test, and optical polishing surfaces. Right now, he is a member of the National System of Researchers (SNI). He has been working for the past 11 years as a research lecturer at the Technological University of the Mixteca UTM Oaxaca, México.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Alberto Cordero-Dávila and Jorge González-García
"Optical surface evaluation by correlating bi-Ronchigram images", Opt. Eng. 54(3), 034108 (Mar 16, 2015). ; http://dx.doi.org/10.1117/1.OE.54.3.034108


Figures

Graphic Jump Location
Fig. 1
F1 :

(a) Experimental setup of Ronchi test and (b) CCD camera, Ronchi tester (λ=660nm) and square grid.

Graphic Jump Location
Fig. 2
F2 :

Geometry of Ronchi test with square grid.

Graphic Jump Location
Fig. 5
F5 :

Bi-Ronchigram irradiance plots, along the x axis, (a) before and (b) after the normalization algorithm was applied.

Graphic Jump Location
Fig. 4
F4 :

Normalized experimental bi-Ronchigram image for mirror of 14 cm diameter and 60.5 cm curvature radius.

Graphic Jump Location
Fig. 3
F3 :

Experimental bi-Ronchigram image for mirror of 14 cm diameter and 60.5 cm curvature radius.

Graphic Jump Location
Fig. 6
F6 :

Reproduced bi-Ronchigram from an ideal spherical mirror.

Graphic Jump Location
Fig. 7
F7 :

Experimental bi-Ronchigram of a parabolic mirror.

Graphic Jump Location
Fig. 8
F8 :

(a) An experimental bi-Ronchigram image for an elliptical mirror and (b) experimental bi-Ronchigram irradiance along x axis.

Tables

References

De Vany  A. S., “Interpreting wave-front and glass-error slopes in an interferogram,” Appl. Opt.. 19, , 173  (1980). 0003-6935 CrossRef
De Vany  A. S., “Patterns correlation of interferograms and Ronchigrams,” Appl. Opt.. 20, , A40  (1981). 0003-6935 CrossRef
Cornejo  A., “Ronchi test,” in Optical Shop Testing. , Malacara  D., Ed., pp. 317 –360,  John Wiley and Sons ,  New York  (2007).
Cornejo-Rodríguez  A., and Cordero-Dávila  A., “Wavefront slope measurements in optical shop testing,” in Handbook of Optical Engineering. , Malacara  D., and Thompson  B. J., Eds., pp. 53 –71,  Marcel Dekker ,  New York  (2001).
Malacara-Doblado  D., and Ghozeil  I., “Hartman, Hartmann-Shack, and other screen tests,” in Optical Shop Testing. , Malacara  D., Ed., pp. 361 –397,  John Wiley and Sons ,  New York  (2007).
Takeda  M., , Ina  H., and Kobayashi  S., “Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry,” J. Opt. Soc. Am.. 72, , 156 –160 (1982). 0030-3941 CrossRef
Servin  M., , Malacara  D., and Cuevas  F. J., “Direct-phase detection of modulated Ronchi rulings using a phase-locked loop,” Opt. Eng.. 33, (4 ), 1193 –1199 (1994). 0091-3286 CrossRef
Cordero-Dávila  A.  et al., “Ronchi test with a square grid,” Appl. Opt.. 37, , 672 –675 (1998). 0003-6935 CrossRef
Cordero-Dávila  A.  et al., “Only one fitting for bironchigrams,” Appl. Opt.. 40, , 5600 –5609 (2001). 0003-6935 CrossRef
Vázquez-Montiel  S., , Sánchez-Escobar  J. J., and Fuentes  O., “Obtaining the phase of an interferogram by use of an evolution strategy: part I,” Appl. Opt.. 41, , 3448 –3452 (2002). 0003-6935 CrossRef
Aguirre-Aguirre  D.  et al., “Algorithm for ronchigram recovery with random aberrations coefficients,” Opt. Eng.. 52, (5 ), 053606  (2013). 0091-3286 CrossRef
Cordero-Dávila  A., , Cornejo-Rodríguez  A., and Núñez  O. C., “Ronchi and Hartmann tests with the same mathematical theory,” Appl. Opt.. 31, , 2370 –2376 (1992). 0003-6935 CrossRef
Cordero-Dávila  A., , Kantum-Montiel  J. R., and González-García  J., “Ronchigram simulations for free-form concave reflective surfaces,” Optik. 124, , 4892 –4895 (2013). 0030-4026 CrossRef
Cordero-Dávila  A., , Cardona-Núñez  O., and Cornejo-Rodríguez  A., “Least squares estimators for the center and radius of circular patterns,” Appl. Opt.. 32, , 5683 –5685 (1993). 0003-6935 CrossRef
Juarez-Salazar  R.  et al., “Generalized phase-shifting interferometry by parameter estimation with the least squares method,” Opt. Lasers Eng.. 51, , 626 –632 (2013). 0143-8166 CrossRef
Cordero-Dávila  A.  et al., “Local and global surface errors evaluation using Ronchi test, without both approximation and integration,” Appl. Opt.. 50, , 4817 –4823 (2011). 0003-6935 CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

PubMed Articles
The coherent radiation fraction of low-emittance synchrotrons. J Synchrotron Radiat 2017;24(Pt 3):566-575.
Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.