Study of Zernike polynomials of an elliptical aperture ... - OSA Publishing

2 downloads 0 Views 82KB Size Report
Aug 20, 2013 - 2The Aerospace Corporation, 2310 El Segundo Boulevard, El Segundo, California 90245, USA. *Corresponding author: [email protected].
Study of Zernike polynomials of an elliptical aperture obscured with an elliptical obscuration: comment José Antonio Díaz1,* and Virendra N. Mahajan2 1

Departamento de Óptica, Universidad de Granada, Granada 18071, Spain

2

The Aerospace Corporation, 2310 El Segundo Boulevard, El Segundo, California 90245, USA *Corresponding author: [email protected] Received 11 March 2013; accepted 16 May 2013; posted 11 July 2013 (Doc. ID 186616); published 15 August 2013

Recently, Hasan and Shaker published a set of orthonormal polynomials for an annular elliptical pupil obtained by the Gram–Schmidt orthogonalization of the Zernike circle polynomials [Appl. Opt. 51, 8490 (2012)]. However, the expressions for many of the polynomials are incorrect, apparently due to wrong usage of the Gram–Schmidt orthogonalization process. We provide the correct equations for the orthogonalization process and the expressions for the orthonormal polynomials obtained by applying them. © 2013 Optical Society of America OCIS codes: (110.0110) Imaging systems; (010.7350) Wave-front sensing; (220.1010) Aberrations (global); (120.3180) Interferometry; (220.0220) Optical design and fabrication. http://dx.doi.org/10.1364/AO.52.005962

Recently, Hasan and Shaker [1] published a set of orthonormal polynomials for an annular elliptical pupil obtained by the Gram–Schmidt orthogonalization of the Zernike circle polynomials. However, they did not mention the relevance of these polynomials by way of any practical applications. Without relevance, such work may not have utility, and thus reduce to a mathematical exercise. Upon further examination, we find that the expressions for many of the annular elliptical polynomials given by them are incorrect. This is apparently due to the wrong usage of the Gram–Schmidt orthogonalization process. The annular elliptical polynomials Ej x; y; b; ϵ for an elliptical pupil of aspect ratio b and obscuration ratio ϵ can be obtained by Gram–Schmidt orthonormalizing the Zernike circle polynomials Zj x; y over the annular elliptical pupil recursively according to [2]

1559-128X/13/245962-03$15.00/0 © 2013 Optical Society of America 5962

APPLIED OPTICS / Vol. 52, No. 24 / 20 August 2013

 Ej1  N j1 Zj1 −

j X

 hZj1 Ek iEk ;

(1)

k1

where E1  1; Z Z bp 1 1−x2 1 hZj1 Ek i  dx p Zj1 Ek dy 2 πb1 − ϵ  −1 −b 1−x2  Zϵ Z bp 2 2 ϵ −x − dx Z E dy ; p j1 k −b ϵ2 −x2

−ϵ

(2)

(3)

and N j1 is a normalization constant so that the E-polynomials are orthonormal over the pupil, i.e., Z Z bp 1 1−x2 1 dx hEj Ej0 i  p Ej Ej0 dy πb1 − ϵ2  −1 −b 1−x2  Zϵ Z bp ϵ2 −x2 − dx p Ej Ej0 dy  δjj0 : −ϵ

−b ϵ2 −x2

(4)

Table 1.

Orthonormal Annular Elliptical Polynomials Ex ;y ; b; ϵ for an Elliptical Pupil with an Aspect Ratio b and Obscuration Ratio ϵ

Ej

Ex; y; b; ϵ

E1

1

E2

N2x

E3

N3y

E4

N 4 2ρ2 − 12 1  b2 1  ϵ2 

E5

  3−1  b4 1  ϵ4 1 − 4ρ2  ϵ2  b2 1  ϵ2  N 5 x2 − y2  14 −1  b2 1  ϵ2  − 431  ϵ4  − 2b2 1  4ϵ2  ϵ4   3b4 1  ϵ4 

E6

2N 6 xy

E7

  N7 1  3b2 1  ϵ2  ϵ4  y 6ρ2 − 2 2 1ϵ

E8

  N8 3  b2 1  ϵ2  ϵ4  x 6ρ2 − 2 1  ϵ2

E9

a

E10 b

8 9 −3b2 −1  b2 1  2ϵ2  3ϵ4  3ϵ6  2ϵ8  ϵ10  = 4N 9 < y 3b2 x2 −1 − 2ϵ2 − 6ϵ4 − 2ϵ6 − ϵ8  3b2 1  2ϵ2 − 2ϵ4  2ϵ6  ϵ8  ; A : 2 y −5 − 10ϵ2 − 6ϵ4 − 10ϵ6 − 5ϵ8  3b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8  8 4N 10 < x B :

9 −3b2 −1  b2 1  2ϵ2  3ϵ4  3ϵ6  2ϵ8  ϵ10  = 2 4 6 8 2 2 4 6 8  2ϵ − 2ϵ  2ϵ  ϵ   b 1  2ϵ  6ϵ  2ϵ  ϵ  ; b2 x2 −31  2ϵ2  6ϵ4  2ϵ6  ϵ8   b2 5  10ϵ2  6ϵ4  10ϵ6  5ϵ8  3y2 −31

E11

   1 31  ϵ2  3x2 2  ϵ2  2ϵ4   3b4 ϵ2 x2 3ϵ2 y2  3b4 y2 2  ϵ2  2ϵ4  2b2 x2 − 2y2  N 11 6ρ4  3  2b2  3b4 1  4ϵ2  ϵ4  − − 8 4 1  ϵ2  ϵ4  b2 1  ϵ2  ϵ4 

E12

  3−1  b2 1  ϵ2 1  ϵ4  4N 12 xy x2 − y2  2 4 81  ϵ  ϵ 

E13

−10.2857ρ2  18.1919ρ4  sin 2θ

E14

3.2184 − 20.4282ρ2  31.6338ρ4  13.7867ρ2 − 35.1829ρ4  cos 2θ  8.2278ρ4 cos 4θ

E15

7.5667ρ2 − 18.4404ρ4  sin 2θ  7.3881ρ4 sin 4θ

A  5  10ϵ2  6ϵ4  10ϵ6  5ϵ8 − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8   9b4 1  2ϵ2 − 2ϵ4  2ϵ6  ϵ8 : B  91  2ϵ2 − 2ϵ4  2ϵ6  ϵ8  − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8   b4 5  10ϵ2  6ϵ4  10ϵ6  5ϵ8 :

a b

Here, δjj0 is a Kronecker delta. Letting j0  1, we find that the mean value hEj i of a polynomial is zero (except when j  1). Similarly, letting j  j0 , its mean square value hE2j i is unity. The new polynomials can also be obtained nonrecursively using a matrix approach, as discussed in [3]. The closed-form analytical expressions of the first 12 annular elliptical polynomials obtained by using the Gram–Schmidt algorithm Eqs. (1)–(3) are given in Tables 1 and 2. Of course, these polynomials reduce to (1) unobscured elliptical polynomials [4–6] as ϵ → 0; (2) annular polynomials as b → 1; and (3) circle polynomials as b → 1 and ϵ → 0. Hasan and Shaker use a different numbering of the polynomials, but without any explanation. For example, our E15 is their A13. However, their polynomials A5 and

A9–A11 given in Table 2 are incorrect. Also their normalization coefficients C4, C5, and C7–C9 in Table 3 are incorrect. Hence, all of their orthonormal polynomials Cj Aj are incorrect, except for j  1, 2, 3, 6, and 12. Even the unobscured elliptical polynomials given by them in their Table 1 have errors. For example, their pE2 should be 2x and E4 should be multiplied by 2. We have also included the numerical form of our polynomials E13 , E14 , and E15 for b  0.7 and k  0.5, because the closed-form analytical expressions for them are too long. Comparing with the corresponding numerical results of Hasan and Shaker, given by them as polynomials A12, A15, and A13 in their Table 6, we find that they are incorrect. Further examination shows that their polynomials A6, A10, 20 August 2013 / Vol. 52, No. 24 / APPLIED OPTICS

5963

Table 2.

Normalization Constants N j for the Annular Elliptical Polynomials Given in Table 1

Nj

Expression

N1

1

N2

21  ϵ2 −1∕2

N3

2b−1 1  ϵ2 −1∕2

N4

p 1231  ϵ4  − 2b2 1  4ϵ2  ϵ4   3b4 1  ϵ4 −1∕2

N5

91  ϵ4   9b4 1  ϵ4  − 6b2 1  4ϵ2  ϵ4 1∕2 2b4 −1  ϵ2 2 1  ϵ2  ϵ4 −1∕2

N6

p 2 6b 1  ϵ2  ϵ4 −1∕2

N7

641  ϵ2 1∕2 fb2 5  10ϵ2  6ϵ4  10ϵ6  5ϵ8 − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8   9b4 1  2ϵ2 − 2ϵ4  2ϵ6  ϵ8 g−1∕2

N8

641  ϵ2 1∕2 91  2ϵ2 − 2ϵ4  2ϵ6  ϵ8  − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8   b4 5  10ϵ2  6ϵ4  10ϵ6  5ϵ8 −1∕2 

N9

5

N 12



6ϵ4



10ϵ6

b6 1  2ϵ2  ϵ4 1  ϵ2 1  ϵ4 1  4ϵ2  ϵ4   5ϵ8  9b4 1  2ϵ2 − 2ϵ4  2ϵ6  ϵ8  − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8 

3201  ϵ2  ϵ4 1∕2

  451  b8 1  ϵ2  3ϵ4 − 4ϵ6  3ϵ8  ϵ10  ϵ12  − 60b2  b8 1  3ϵ2  3ϵ4  4ϵ6  3ϵ8  3ϵ10  ϵ12  −1∕2  b4 94  78ϵ2  282ϵ4  712ϵ6  282ϵ8  78ϵ10  94ϵ12 

 2  b 5  10ϵ2  5ϵ4 3  3ϵ4  2ϵ6  ϵ8   6b2 1  2ϵ2  3ϵ4  8ϵ6  3ϵ8  2ϵ10  ϵ12   5b4 1  2ϵ2  3ϵ4  3ϵ8  2ϵ10  ϵ12  −1∕2 2 4 1601  ϵ  ϵ 

and A14 are also incorrect, as may be seen by comparing with our polynomials E5 , E10 , and E12 , respectively, for the above values of b and k. References 1. S. Y. Hasan and A. S. Shaker, “Study of Zernike polynomials of an elliptical aperture obscured with and elliptical obscuration,” Appl. Opt. 51, 8490–8497 (2012). 2. G. A. Korn and T. M. Korn, Mathematical Handbook for Scientists and Engineers (McGraw-Hill, 1968), p. 454.

5964

−1∕2

b4 1  2ϵ2  ϵ4 1  ϵ2 1  ϵ4 1  4ϵ2  ϵ4 −1∕2 ×91  2ϵ2 − 2ϵ4  2ϵ6  ϵ8  − 6b2 1  2ϵ2  6ϵ4  2ϵ6  ϵ8   b4 5  10ϵ2  6ϵ4  10ϵ6  5ϵ8 1∕2

N 10

N 11

10ϵ2

APPLIED OPTICS / Vol. 52, No. 24 / 20 August 2013

3. G.-m. Dai and V. N. Mahajan, “Nonrecursive orthonormal polynomials with matrix formulation,” Opt. Lett. 32, 74–76 (2007). 4. V. N. Mahajan and G.-m. Dai, “Orthonormal polynomials in wavefront analysis: analytical solution,” J. Opt. Soc. Am. A 24, 2994–3016 (2007). 5. V. N. Mahajan, “Orthonormal polynomials in wavefront analysis: analytical solution: errata” J. Opt. Soc. Am. A 29, 1673–1674 (2012). 6. V. N. Mahajan, “Orthonormal polynomials in wavefront analysis,” in Handbook of Optics, V. N. Mahajan and E. V. Stryland, eds., 3rd ed., Vol. II (McGraw-Hill, 2010), pp. 11.3–11.41.