World Library  
Flag as Inappropriate
Email this Article

Zernike polynomials

Article Id: WHEBN0002025291
Reproduction Date:

Title: Zernike polynomials  
Author: World Heritage Encyclopedia
Language: English
Subject: Optical aberration, Frits Zernike, Fourier optics, Point spread function, Series expansion, List of special functions and eponyms, Piston (optics), Tilt (optics), Deformable mirror, Pseudo-Zernike polynomials
Collection:
Publisher: World Heritage Encyclopedia
Publication
Date:
 

Zernike polynomials


In mathematics, the Zernike polynomials are a sequence of polynomials that are orthogonal on the unit disk. Named after Nobel Prize winner and optical physicist, and inventor of the phase contrast microscopy, Frits Zernike, they play an important role in beam optics.[1][2]


Definitions

There are even and odd Zernike polynomials. The even ones are defined as

Z^{m}_n(\rho,\varphi) = R^m_n(\rho)\,\cos(m\,\varphi) \!

and the odd ones as

Z^{-m}_n(\rho,\varphi) = R^m_n(\rho)\,\sin(m\,\varphi), \!

where m and n are nonnegative integers with nm, φ is the azimuthal angle, and ρ is the radial distance 0\le\rho\le 1. Zernike polynomials have the property of being limited to a range of -1 to +1, i.e. |Z^{m}_n(\rho,\varphi)| \le 1. The radial polynomials Rmn are defined as

R^m_n(\rho) = \! \sum_{k=0}^{(n-m)/2} \!\!\! \frac{(-1)^k\,(n-k)!}{k!\,((n+m)/2-k)!\,((n-m)/2-k)!} \;\rho^{n-2\,k}

for nm even, and are identically 0 for nm odd.

Other Representations

Rewriting the ratios of factorials in the radial part as products of binomials shows that the coefficients are integer numbers:

R_n^m(\rho)=\sum_{k=0}^{(n-m)/2}(-1)^k \binom{n-k}{k} \binom{n-2k}{(n-m)/2-k} \rho^{n-2k}.

A notation as terminating Gaussian Hypergeometric Functions is useful to reveal recurrences, to demonstrate that they are special cases of Jacobi polynomials, to write down the differential equations, etc.:

R_n^m(\rho)= \binom{n}{(n+m)/2}

\rho^n {}_2F_{1}\left(-\frac{n+m}{2},-\frac{n-m}{2};-n;\rho^{-2}\right)

= (-1)^{(n+m)/2}\binom{(n+m)/2}{(n-m)/2}\rho^m

{}_2F_{1}\left(1+n,1-\frac{n-m}{2};1+\frac{n+m}{2};\rho^2\right) for nm even.

Noll's sequential indices

Applications often involve linear algebra, where integrals over products of Zernike polynomials and some other factor build the matrix elements. To enumerate the rows and columns of these matrices by a single index, a conventional mapping of the two indices n and m to a single index j has been introduced by Noll.OEIS)

n,m 0,0 1,1 1,-1 2,0 2,-2 2,2 3,-1 3,1 3,-3 3,3
j 1 2 3 4 5 6 7 8 9 10
n,m 4,0 4,2 4,-2 4,4 4,-4 5,1 5,-1 5,3 5,-3 5,5
j 11 12 13 14 15 16 17 18 19 20

The rule is that the even Z (with even azimuthal part m, \cos(m\varphi)) obtain even indices j, the odd Z odd indices j. Within a given n, lower values of m obtain lower j.

Properties

Orthogonality

The orthogonality in the radial part reads

\int_0^1 \rho \sqrt{2n+2}R_n^m(\rho)\,\sqrt{2n'+2}R_{n'}^{m}(\rho)d\rho = \delta_{n,n'}.

Orthogonality in the angular part is represented by

\int_0^{2\pi} \cos(m\varphi)\cos(m'\varphi)d\varphi=\epsilon_m\pi\delta_{|m|,|m'|},
\int_0^{2\pi} \sin(m\varphi)\sin(m'\varphi)d\varphi=(-1)^{m+m'}\pi\delta_{|m|,|m'|};\quad m\neq 0,
\int_0^{2\pi} \cos(m\varphi)\sin(m'\varphi)d\varphi=0,

where \epsilon_m (sometimes called the Neumann factor because it frequently appears in conjunction with Bessel functions) is defined as 2 if m=0 and 1 if m\neq 0. The product of the angular and radial parts establishes the orthogonality of the Zernike functions with respect to both indices if integrated over the unit disk,

\int Z_n^m(\rho,\varphi)Z_{n'}^{m'}(\rho,\varphi)d^2r

=\frac{\epsilon_m\pi}{2n+2}\delta_{n,n'}\delta_{m,m'}, where d^2r=\rho\,d\rho\,d\varphi is the Jacobian of the circular coordinate system, and where n-m and n'-m' are both even.

A special value is

R_n^m(1)=1.

Zernike Transform

Any sufficiently smooth real-valued phase field over the unit disk G(\rho,\varphi) can be represented in terms of its Zernike coefficients (odd and even), just as much as periodic functions find an orthogonal representation with the Fourier series. We have

G(\rho,\varphi) = \sum_{m,n}\left[ a_{m,n} Z^{m}_n(\rho,\varphi) + b_{m,n} Z^{-m}_n(\rho,\varphi) \right],

where the coefficients can be calculated using inner products. On the space of L^2 functions on the unit disk we have an inner product defined by

\langle F, G \rangle := \int F(\rho,\varphi)G(\rho,\varphi)\rho d\rho d\varphi.

The Zernike coefficients can then be expressed as follows:

a_{m,n} = \frac{2n+2}{\epsilon_m\pi}\langle G(\rho,\varphi),Z^{m}_n(\rho,\varphi) \rangle,\quad b_{m,n} = \frac{2n+2}{\epsilon_m\pi}\langle G(\rho,\varphi),Z^{-m}_n(\rho,\varphi) \rangle.

Alternatively, we can use the known values of phase function G on the circular grid to form a system of equations. The phase function is retrieved by the unknown-coefficient weighted product with (known values) of Zernike polynomial across the unit grid. Hence coefficients can also be found by solving linear system, for instance by matrix inversion. Fast algorithms to calculate the forward and inverse Zernike transform use symmetry properties of trigonometric functions, separability of radial and azimuthal parts of Zernike polynomials, and their rotational symmetries.

Symmetries

The parity with respect to reflection along the x axis is

Z_n^{m}(\rho,\varphi)=(-1)^m Z_n^{m}(\rho,-\varphi).

The parity with respect to point reflection at the center of coordinates is

Z_n^m(\rho,\varphi) = (-1)^m Z_n^m(\rho,\varphi+\pi),

where (-1)^m could as well be written (-1)^n because n-m is even for the relevant, non-vanishing values. The radial polynomials are also either even or odd, depending on order n or m:

R_n^m(\rho)=(-1)^n R_n^m(-\rho)=(-1)^m R_n^m(-\rho).

The periodicity of the trigonometric functions implies invariance if rotated by multiples of 2\pi/m radian around the center:

Z_n^m(\rho,\varphi+2\pi k/m)=Z_n^m(\rho,\varphi),\quad k= 0, \pm 1,\pm 2,\ldots.

Examples

Radial polynomials

The first few radial polynomials are:

R^0_0(\rho) = 1 \,
R^1_1(\rho) = \rho \,
R^0_2(\rho) = 2\rho^2 - 1 \,
R^2_2(\rho) = \rho^2 \,
R^1_3(\rho) = 3\rho^3 - 2\rho \,
R^3_3(\rho) = \rho^3 \,
R^0_4(\rho) = 6\rho^4 - 6\rho^2 + 1 \,
R^2_4(\rho) = 4\rho^4 - 3\rho^2 \,
R^4_4(\rho) = \rho^4 \,
R^1_5(\rho) = 10\rho^5 - 12\rho^3 + 3\rho \,
R^3_5(\rho) = 5\rho^5 - 4\rho^3 \,
R^5_5(\rho) = \rho^5 \,
R^0_6(\rho) = 20\rho^6 - 30\rho^4 + 12\rho^2 - 1 \,
R^2_6(\rho) = 15\rho^6 - 20\rho^4 + 6\rho^2 \,
R^4_6(\rho) = 6\rho^6 - 5\rho^4 \,
R^6_6(\rho) = \rho^6. \,

Zernike polynomials

The first few Zernike modes, ordered by Noll index j are[3]

Noll index (j) Radial degree (n) Azimuthal degree (m) Z_j Classical name
1 0 0 1 Piston
2 1 1 2 \rho \cos \theta Tip (lateral position) (X-Tilt)
3 1 -1 2 \rho \sin \theta Tilt (lateral position) (Y-Tilt)
4 2 0 \sqrt{3} (2 \rho^2 - 1) Defocus (longitudinal position)
5 2 -2 \sqrt{6} \rho^2 \sin 2 \theta Astigmatism
6 2 2 \sqrt{6} \rho^2 \cos 2 \theta Astigmatism
7 3 -1 \sqrt{8} (3 \rho^3 - 2\rho) \sin \theta Coma
8 3 1 \sqrt{8} (3 \rho^3 - 2\rho) \cos \theta Coma
9 3 -3 \sqrt{8} \rho^3 \sin 3 \theta Trefoil
10 3 3 \sqrt{8} \rho^3 \cos 3 \theta Trefoil
11 4 0 \sqrt{5} (6 \rho^4 - 6 \rho^2 +1) Third order spherical
12 4 2 \sqrt{10} (4 \rho^4 - 3\rho^2) \cos 2 \theta
13 4 -2 \sqrt{10} (4 \rho^4 - 3\rho^2) \sin 2 \theta
14 4 4 \sqrt{10} \rho^4 \cos 4 \theta
15 4 -4 \sqrt{10} \rho^4 \sin 4 \theta

Applications

The functions are a basis defined over the circular support area, typically the pupil planes in classical optical imaging at visible and infrared wavelengths through systems of lenses and mirrors of finite diameter. Their advantages are the simple analytical properties inherited from the simplicity of the radial functions and the factorization in radial and azimuthal functions; this leads, for example, to closed form expressions of the two-dimensional Fourier transform in terms of Bessel Functions. Their disadvantage, in particular if high n are involved, is the unequal distribution of nodal lines over the unit disk, which introduces ringing effects near the perimeter \rho\approx 1, which often leads attempts to define other orthogonal functions over the circular disk.

In precision optical manufacturing, Zernike polynomials are used to characterize higher-order errors observed in interferometric analyses, in order to achieve desired system performance.

In optometry and ophthalmology Zernike polynomials are used to describe aberrations of the cornea or lens from an ideal spherical shape, which result in refraction errors.

They are commonly used in adaptive optics where they can be used to effectively cancel out atmospheric distortion. Obvious applications for this are IR or visual astronomy and Satellite imagery. For example, one of the Zernike terms (for m = 0, n = 2) is called 'de-focus'. By coupling the output from this term to a control system, an automatic focus can be implemented.

Another application of the Zernike polynomials is found in the Extended Nijboer-Zernike (ENZ) theory of diffraction and aberrations.

Zernike polynomials are widely used as basis functions of image moments. Since Zernike polynomials are orthogonal to each other, Zernike moments can represent properties of an image with no redundancy or overlap of information between the moments. Although Zernike moments are significantly dependent on the scaling and the translation of the object in an ROI, their magnitudes are independent of the rotation angle of the object.[4] Thus, they can be utilized to extract features from images that describe the shape characteristics of an object. For instance, Zernike moments are utilized as shape descriptors to classify benign and malignant breast masses.[5][6]

Higher Dimensions

The concept translates to higher dimensions D if multinomials x_1^ix_2^j\cdots x_D^k in Cartesian coordinates are converted to hyperspherical coordinates, \rho^s, s\le D, multiplied by a product of Jacobi Polynomials of the angular variables. In D=3 dimensions, the angular variables are Spherical harmonics, for example. Linear combinations of the powers \rho^s define an orthogonal basis R_n^{(l)}(\rho) satisfying

\int_0^1 \rho^{D-1}R_n^{(l)}(\rho)R_{n'}^{(l)}(\rho)d\rho = \delta_{n,n'}.

(Note that a factor \sqrt{2n+D} is absorbed in the definition of R here, whereas in D=2 the normalization is chosen slightly differently. This is largely a matter of taste, depending on whether one wishes to maintain an integer set of coefficients or prefers tighter formulas if the orthogonalization is involved.) The explicit representation is

R_n^{(l)}(\rho) = \sqrt{2n+D}\sum_{s=0}^{(n-l)/2}

(-1)^s{(n-l)/2 \choose s}{n-s-1+D/2 \choose (n-l)/2}\rho^{n-2s}

=(-1)^{(n-l)/2}\sqrt{2n+D}\sum_{s=0}^{(n-l)/2}

(-1)^s{(n-l)/2 \choose s} {s-1+(n+l+D)/2 \choose (n-l)/2} \rho^{2s+l}

=(-1)^{(n-l)/2}\sqrt{2n+D}{ (D+n+l)/2-1 \choose (n-l)/2}\rho^l

{}_2F_1( -(n-l)/2,(n+l+D)/2;l+D/2;\rho^2)

for even n-l\ge 0, else identical to zero.

See also

References

  • MathWorld.
  • from The Wolfram Demonstrations Project.
  • Mathar, R. J. (2007). "Third Order Newton's Method for Zernike Polynomial Zeros". 2007arXiv0705.1329M.

External links

  • The Extended Nijboer-Zernike website.
  • MATLAB code for fast calculation of Zernike moments
  • Python/NumPy library for calculating Zernike polynomials
  • Telescope Optics
  • Example: using WolframAlpha to plot Zernike Polynomials
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
 
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
 
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.
 


Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.