Wielomiany Zernikego

Wielomiany Zernikego są zbiorem wielomianów ortogonalnych wewnątrz koła jednostkowego wprowadzonych przez Fritsa Zernike.

Definicja

Wielomiany Zernikego zdefiniowane są w postaci zespolonej:

V n ± m ( ρ , θ ) = R n m ( ρ ) exp ( ± j m θ ) , {\displaystyle V_{n}^{\pm m}(\rho ,\theta )=R_{n}^{m}(\rho )\,\exp(\pm jm\theta ),}

gdzie:

n , m {\displaystyle n,m} są liczbami naturalnymi takimi, że 0 m n , {\displaystyle 0\leqslant m\leqslant n,} oraz n m {\displaystyle n-m} jest parzyste,
ρ , θ {\displaystyle \rho ,\theta } współrzędnymi biegunowymi punktu (odpowiednio długością promienia wodzącego i wartością kąta skierowanego).
R n m ( ρ ) {\displaystyle R_{n}^{m}(\rho )} jest wielomianem radialnym postaci:
R n m ( ρ ) = s = 0 n m 2 ( 1 ) s ( n s ) ! s ! ( n + m 2 s ) ! ( n m 2 s ) ! ρ n 2 s . {\displaystyle R_{n}^{m}(\rho )=\sum _{s=0}^{\frac {n-m}{2}}(-1)^{s}{\frac {(n-s)!}{s!({\frac {n+m}{2}}-s)!({\frac {n-m}{2}}-s)!}}\rho ^{n-2s}.}

Czasami spotyka się również definicję wielomianów Zernikego w postaci rzeczywistej. Wyróżnia się parzyste i nieparzyste wielomiany Zernikego

Z n m ( ρ , θ ) = R n m ( ρ ) cos ( m θ ) {\displaystyle Z_{n}^{m}(\rho ,\theta )=R_{n}^{m}(\rho )\,\cos(m\,\theta )} – wielomian parzysty,
Z n m ( ρ , θ ) = R n m ( ρ ) sin ( m θ ) {\displaystyle Z_{n}^{-m}(\rho ,\theta )=R_{n}^{m}(\rho )\,\sin(m\,\theta )} – wielomian nieparzysty.

Przykłady

Kolejne wielomiany Zernike mają rozwinięcie

V 0 0 {\displaystyle V_{0}^{0}} 1 {\displaystyle 1}
V 1 1 {\displaystyle V_{1}^{1}} ρ exp ( j θ ) {\displaystyle \rho \exp(j\theta )}
V 2 0 {\displaystyle V_{2}^{0}} 2 ρ 2 1 {\displaystyle 2\rho ^{2}-1}
V 2 2 {\displaystyle V_{2}^{2}} ρ 2   exp ( j 2 θ ) {\displaystyle \rho ^{2}\ \exp(j2\theta )}
V 3 1 {\displaystyle V_{3}^{1}} ( 3 ρ 3 2 ρ )   exp ( j θ ) {\displaystyle (3\rho ^{3}-2\rho )\ \exp(j\theta )}
V 3 3 {\displaystyle V_{3}^{3}} ρ 3   exp ( j 3 θ ) {\displaystyle \rho ^{3}\ \exp(j3\theta )}
V 4 0 {\displaystyle V_{4}^{0}} 6 ρ 4 6 ρ 2 + 1 {\displaystyle 6\rho ^{4}-6\rho ^{2}+1}
V 4 2 {\displaystyle V_{4}^{2}} ( 4 ρ 4 3 ρ 2 )   exp ( j 2 θ ) {\displaystyle (4\rho ^{4}-3\rho ^{2})\ \exp(j2\theta )}
V 4 4 {\displaystyle V_{4}^{4}} ρ 4   exp ( j 4 θ ) {\displaystyle \rho ^{4}\ \exp(j4\theta )}

Mapy jasności niektórych wielomianów Zernikego:

V 2 0 {\displaystyle V_{2}^{0}} V 2 2 {\displaystyle V_{2}^{2}} V 3 1 {\displaystyle V_{3}^{1}} V 3 3 {\displaystyle V_{3}^{3}} V 4 0 {\displaystyle V_{4}^{0}}
Część rzeczywista
Część urojona

Własności

Wielomiany radialne są ortogonalne:

0 1 ρ R n m ( ρ )   R n m ( ρ )   d ρ = 1 2 ( n + 1 ) δ n n {\displaystyle \int \limits _{0}^{1}\rho R_{n}^{m}(\rho )\ R_{n'}^{m}(\rho )\ d\rho ={\frac {1}{2(n+1)}}\delta _{nn'}}

gdzie δ n n {\displaystyle \delta _{nn'}} oznacza deltę Kroneckera. Podobnie, ortogonalność zachodzi dla wielomianów Zernikego:

ρ [ V n m ( ρ , θ ) ]   V p q ( ρ , θ )   d ρ d θ = π n + 1 δ n m δ p q {\displaystyle \iint \rho [V_{n}^{m}(\rho ,\theta )]^{*}\ V_{p}^{q}(\rho ,\theta )\ d\rho d\theta ={\frac {\pi }{n+1}}\delta _{nm}\delta _{pq}}

Wielomiany te posiadają również własność rotacyjną

V n m ( ρ , θ + α ) = V n m ( ρ , θ ) exp ( j m α ) {\displaystyle V_{n}^{m}(\rho ,\theta +\alpha )=V_{n}^{m}(\rho ,\theta )\exp(jm\alpha )}

co oznacza, że ich moduł jest niezależny od obrotu:

| V n m ( ρ , θ + α ) | = | V n m ( ρ , θ ) | . {\displaystyle |V_{n}^{m}(\rho ,\theta +\alpha )|=|V_{n}^{m}(\rho ,\theta )|.}

Sprzężenie wielomianu Zernikego ma wartość:

( V n m ( ρ , θ ) ) = V n m ( ρ , θ ) {\displaystyle (V_{n}^{m}(\rho ,\theta ))^{*}=V_{n}^{-m}(\rho ,\theta )}

Zastosowanie

W optyce, wielomiany Zernikego stosuje się do opisu aberracji soczewek.

Wielomiany Zernikego znalazły też zastosowanie w cyfrowym przetwarzaniu obrazów, do dekompozycji obrazów na tzw. momenty Zernikego.

Zobacz też

Linki zewnętrzne

  • Eric W.E.W. Weisstein Eric W.E.W., Zernike Polynomial, [w:] MathWorld, Wolfram Research [dostęp 2020-12-13]  (ang.).
  • James C. Wyant: „Zenike Polyniomials”. optics.arizona.edu. [zarchiwizowane z tego adresu (2006-09-06)]. (en.)