Pseudo-Zernike polynomials

Last updated June 29, 2021

In mathematics, pseudo-Zernike polynomials are well known and widely used in the analysis of optical systems. They are also widely used in image analysis as shape descriptors.

Definition

They are an orthogonal set of complex-valued polynomials defined as

V_{nm}(x,y)=R_{nm}(x,y)e^{jm\arctan({\frac {y}{x}})},

where $x^{2}+y^{2}\leq 1,n\geq 0,|m|\leq n$ and orthogonality on the unit disk is given as

\int _{0}^{2\pi }\int _{0}^{1}r[V_{nl}(r\cos \theta ,r\sin \theta )]^{*}\times V_{mk}(r\cos \theta ,r\sin \theta )\,dr\,d\theta ={\frac {\pi }{n+1}}\delta _{mn}\delta _{kl},

where the star means complex conjugation, and $r^{2}=x^{2}+y^{2}$ , $x=r\cos \theta$ , $y=r\sin \theta$ are the standard transformations between polar and Cartesian coordinates.

The radial polynomials $R_{nm}$ are defined as^[1]

$R_{nm}(r)=\sum _{s=0}^{n-|m|}D_{n,|m|,s}\ r^{n-s}$

with integer coefficients

D_{n,|m|,s}=(-1)^{s}{\frac {(2n+1-s)!}{s!(n-|m|-s)!(n+|m|-s+1)!}}.

Examples

Examples are:

$R_{0,0}=1$

$R_{1,0}=-2+3r$

$R_{1,1}=r$

$R_{2,0}=3+10r^{2}-12r$

$R_{2,1}=5r^{2}-4r$

$R_{2,2}=r^{2}$

$R_{3,0}=-4+35r^{3}-60r^{2}+30r$

$R_{3,1}=21r^{3}-30r^{2}+10r$

$R_{3,2}=7r^{3}-6r^{2}$

$R_{3,3}=r^{3}$

$R_{4,0}=5+126r^{4}-280r^{3}+210r^{2}-60r$

$R_{4,1}=84r^{4}-168r^{3}+105r^{2}-20r$

$R_{4,2}=36r^{4}-56r^{3}+21r^{2}$

$R_{4,3}=9r^{4}-8r^{3}$

$R_{4,4}=r^{4}$

$R_{5,0}=-6+462r^{5}-1260r^{4}+1260r^{3}-560r^{2}+105r$

$R_{5,1}=330r^{5}-840r^{4}+756r^{3}-280r^{2}+35r$

$R_{5,2}=165r^{5}-360r^{4}+252r^{3}-56r^{2}$

$R_{5,3}=55r^{5}-90r^{4}+36r^{3}$

$R_{5,4}=11r^{5}-10r^{4}$

$R_{5,5}=r^{5}$

Moments

The pseudo-Zernike Moments (PZM) of order $n$ and repetition $l$ are defined as

A_{nl}={\frac {n+1}{\pi }}\int _{0}^{2\pi }\int _{0}^{1}[V_{nl}(r\cos \theta ,r\sin \theta )]^{*}f(r\cos \theta ,r\sin \theta )r\,dr\,d\theta ,

where $n=0,\ldots \infty$ , and $l$ takes on positive and negative integer values subject to $|l|\leq n$ .

The image function can be reconstructed by expansion of the pseudo-Zernike coefficients on the unit disk as

f(x,y)=\sum _{n=0}^{\infty }\sum _{l=-n}^{+n}A_{nl}V_{nl}(x,y).

Pseudo-Zernike moments are derived from conventional Zernike moments and shown to be more robust and less sensitive to image noise than the Zernike moments.^[1]

Related Research Articles

In mathematics, the arithmetic–geometric mean (AGM) of two positive real numbers $x$ and $y$ is defined as follows:

Laplaces equation Second order partial differential equation

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace who first studied its properties. This is often written as

In mathematics, an $n$ -sphere is a topological space that is homeomorphic to a standard $n$ -sphere, which is the set of points in $(n + 1)$ -dimensional Euclidean space that are situated at a constant distance $r$ from a fixed point, called the center. It is the generalization of an ordinary sphere in the ordinary three-dimensional space. The "radius" of a sphere is the constant distance of its points to the center. When the sphere has unit radius, it is usual to call it the unit $n$ -sphere or simply the $n$ -sphere for brevity. In terms of the standard norm, the $n$ -sphere is defined as

In physical science and mathematics, Legendre polynomials are a system of complete and orthogonal polynomials, with a vast number of mathematical properties, and numerous applications. They can be defined in many ways, and the various definitions highlight different aspects as well as suggest generalizations and connections to different mathematical structures and physical and numerical applications.

The Chebyshev polynomials are two sequences of polynomials related to the sine and cosine functions, notated as $and . They can be defined several ways that have the same end result; in this article the polynomials are defined by starting with trigonometric functions:$

In mathematics and physical science, spherical harmonics are special functions defined on the surface of a sphere. They are often employed in solving partial differential equations in many scientific fields.

In mathematics, the Clausen function, introduced by Thomas Clausen (1832), is a transcendental, special function of a single variable. It can variously be expressed in the form of a definite integral, a trigonometric series, and various other special functions. It is intimately connected with the polylogarithm, inverse tangent integral, polygamma function, Riemann zeta function, Dirichlet eta function, and Dirichlet beta function.

In mathematics, physics and engineering, the sinc function, denoted by $sinc(x)$ , has two forms, normalized and unnormalized.

In mathematics, Buffon's needle problem is a question first posed in the 18th century by Georges-Louis Leclerc, Comte de Buffon:

In mathematics, the associated Legendre polynomials are the canonical solutions of the general Legendre equation

In mathematics, the Zernike polynomials are a sequence of polynomials that are orthogonal on the unit disk. Named after optical physicist Frits Zernike, winner of the 1953 Nobel Prize in Physics and the inventor of phase-contrast microscopy, they play important roles in various optics branches such as beam optics and imaging.

In mathematics, a change of variables is a basic technique used to simplify problems in which the original variables are replaced with functions of other variables. The intent is that when expressed in new variables, the problem may become simpler, or equivalent to a better understood problem.

Multiple integral Generalization of definite integrals to functions of multiple variables

In mathematics, a multiple integral is a definite integral of a function of several real variables, for instance, $f (x, y)$ or $f (x, y, z)$ . Integrals of a function of two variables over a region in $are called double integrals, and integrals of a function of three variables over a region in are called triple integrals . For multiple integrals of a single-variable function, see the Cauchy formula for repeated integration.$

Spherical multipole moments are the coefficients in a series expansion of a potential that varies inversely with the distance R to a source, i.e., as 1/R. Examples of such potentials are the electric potential, the magnetic potential and the gravitational potential.

A pendulum is a body suspended from a fixed support so that it swings freely back and forth under the influence of gravity. When a pendulum is displaced sideways from its resting, equilibrium position, it is subject to a restoring force due to gravity that will accelerate it back toward the equilibrium position. When released, the restoring force acting on the pendulum's mass causes it to oscillate about the equilibrium position, swinging it back and forth. The mathematics of pendulums are in general quite complicated. Simplifying assumptions can be made, which in the case of a simple pendulum allow the equations of motion to be solved analytically for small-angle oscillations.

Clenshaw–Curtis quadrature and Fejér quadrature are methods for numerical integration, or "quadrature", that are based on an expansion of the integrand in terms of Chebyshev polynomials. Equivalently, they employ a change of variables $and use a discrete cosine transform (DCT) approximation for the cosine series. Besides having fast-converging accuracy comparable to Gaussian quadrature rules, Clenshaw-Curtis quadrature naturally leads to nested quadrature rules, which is important for both adaptive quadrature and multidimensional quadrature (cubature).$

In physics and mathematics, the solid harmonics are solutions of the Laplace equation in spherical polar coordinates, assumed to be (smooth) functions $. There are two kinds: the regular solid harmonics, which vanish at the origin and the irregular solid harmonics, which are singular at the origin. Both sets of functions play an important role in potential theory, and are obtained by rescaling spherical harmonics appropriately:$

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when the diffraction pattern is viewed at a long distance from the diffracting object, and also when it is viewed at the focal plane of an imaging lens.

References

1 2 Teh, C.-H.; Chin, R. (1988). "On image analysis by the methods of moments". IEEE Transactions on Pattern Analysis and Machine Intelligence. 10 (4): 496–513. doi:10.1109/34.3913.

Belkasim, S.; Ahmadi, M.; Shridhar, M. (1996). "Efficient algorithm for the fast computation of zernike moments". Journal of the Franklin Institute. 333 (4): 577–581. doi:10.1016/0016-0032(96)00017-8.
Haddadnia, J.; Ahmadi, M.; Faez, K. (2003). "An efficient feature extraction method with pseudo-zernike moment in rbf neural network-based human face recognition system". EURASIP Journal on Applied Signal Processing. 2003 (9): 890–901. Bibcode:2003EJASP2003..146H. doi: 10.1155/S1110865703305128 .
T.-W. Lin; Y.-F. Chou (2003). A comparative study of zernike moments. Proceedings of the IEEE/WIC International Conference on Web Intelligence. pp. 516–519. doi:10.1109/WI.2003.1241255. ISBN 0-7695-1932-6.
Chong, C.-W.; Raveendran, P.; Mukundan, R. (2003). "The scale invariants of pseudo-Zernike moments" (PDF). Pattern Anal. Applic. 6 (3): 176–184. doi:10.1007/s10044-002-0183-5.
Chong, Chee-Way; Mukundan, R.; Raveendran, P. (2003). "An Efficient Algorithm for Fast Computation of Pseudo-Zernike Moments" (PDF). Int. J. Pattern Recogn. Artif. Int. 17 (6): 1011–1023. doi:10.1142/S0218001403002769. hdl: 10092/448 .
Shutler, Jamie (1992). "Complex Zernike Moments".

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[Teh88-1] 1 2 Teh, C.-H.; Chin, R. (1988). "On image analysis by the methods of moments". IEEE Transactions on Pattern Analysis and Machine Intelligence. 10 (4): 496–513. doi:10.1109/34.3913.

[1]