Polinomis de Zernike

Els primers 21 polinomis de Zernike, ordenats verticalment per grau radial i horitzontalment per grau azimutal

En matemàtiques, els polinomis de Zernike són una seqüència de polinomis que són ortogonals en el disc unitat. Van ser nomenats en honor del físic òptic Frits Zernike, guanyador del Premi Nobel de física de 1953 i inventor del microscopi de contrast de fases. Aquests polinomis tenen un paper important en la modelització del comportament de feixos de llum en un sistema òptic.[1][2]

Definicions

Els polinomis de Zernike es distingeixen en funció de la seva paritat. Els termes parells es defineixen com:

Z n m ( ρ , φ ) = R n m ( ρ ) cos ( m φ ) {\displaystyle Z_{n}^{m}(\rho ,\varphi )=R_{n}^{m}(\rho )\,\cos(m\,\varphi )\!}

i els imparells com:

Z n m ( ρ , φ ) = R n m ( ρ ) sin ( m φ ) , {\displaystyle Z_{n}^{-m}(\rho ,\varphi )=R_{n}^{m}(\rho )\,\sin(m\,\varphi ),\!}

on m i n són nombres enters no negatius amb n ≥ m, φ és l'angle azimutal, ρ és la distància radial 0 ρ 1 {\displaystyle 0\leq \rho \leq 1} i Rmn són els polinomis radials definits a continuació. Els polinomis de Zernike tenen la propietat d'estar limitats a un rang de -1 a +1, és a dir, | Z n m ( ρ , φ ) | 1 {\displaystyle |Z_{n}^{m}(\rho ,\varphi )|\leq 1} . Els polinomis radials Rmn es defineixen com:

R n m ( ρ ) = k = 0 n m 2 ( 1 ) k ( n k ) ! k ! ( n + m 2 k ) ! ( n m 2 k ) ! ρ n 2 k {\displaystyle R_{n}^{m}(\rho )=\sum _{k=0}^{\tfrac {n-m}{2}}{\frac {(-1)^{k}\,(n-k)!}{k!\left({\tfrac {n+m}{2}}-k\right)!\left({\tfrac {n-m}{2}}-k\right)!}}\;\rho ^{n-2\,k}}

per a n-m parell, i són idènticament 0 per a n-m imparella.

Altres representacions

Reescrivint les relacions dels factorials en la part radial com a productes de coeficients binomials, es demostra que els coeficients són nombres enters:

. R n m ( ρ ) = k = 0 n m 2 ( 1 ) k ( n k k ) ( n 2 k n m 2 k ) ρ n 2 k {\displaystyle R_{n}^{m}(\rho )=\sum _{k=0}^{\tfrac {n-m}{2}}(-1)^{k}{\binom {n-k}{k}}{\binom {n-2k}{{\tfrac {n-m}{2}}-k}}\rho ^{n-2k}}

La notació com a termes de funcions hipergeomètriques gaussianes és útil per revelar recurrències, per demostrar casos especials dels polinomis de Jacobi, o per reduir equacions diferencials.

R n m ( ρ ) = ( n n + m 2 ) ρ n   2 F 1 ( n + m 2 , n m 2 ; n ; ρ 2 ) = ( 1 ) n m 2 ( n + m 2 m ) ρ m   2 F 1 ( 1 + n + m 2 , n m 2 ; 1 + m ; ρ 2 ) {\displaystyle {\begin{aligned}R_{n}^{m}(\rho )&={\binom {n}{\tfrac {n+m}{2}}}\rho ^{n}\ {}_{2}F_{1}\left(-{\tfrac {n+m}{2}},-{\tfrac {n-m}{2}};-n;\rho ^{-2}\right)\\&=(-1)^{\tfrac {n-m}{2}}{\binom {\tfrac {n+m}{2}}{m}}\rho ^{m}\ {}_{2}F_{1}\left(1+{\tfrac {n+m}{2}},-{\tfrac {n-m}{2}};1+m;\rho ^{2}\right)\end{aligned}}}

per a n-m parell.

El factor ρ n 2 k {\displaystyle \rho ^{n-2k}} en el polinomi radial R n m ( ρ ) {\displaystyle R_{n}^{m}(\rho )} es pot expandir en una base de Bernstein de b s , n / 2 ( ρ 2 ) {\displaystyle b_{s,n/2}(\rho ^{2})} per a n {\displaystyle n} parell o ρ {\displaystyle \rho } multiplicat per una funció de b s , ( n 1 ) / 2 ( ρ 2 ) {\displaystyle b_{s,(n-1)/2}(\rho ^{2})} per a n {\displaystyle n} imparell en el rang n / 2 k s n / 2 {\displaystyle \lfloor n/2\rfloor -k\leq s\leq \lfloor n/2\rfloor } . Per tant, el polinomi radial pot expressar-se mitjançant un nombre finit de polinomis de Bernstein amb quocients racionals:

R n m ( ρ ) = 1 ( n / 2 m / 2 ) ρ n mod 2 s = m / 2 n / 2 ( 1 ) n / 2 s ( s m / 2 ) ( ( n + m ) / 2 s + m / 2 ) b s , n / 2 ( ρ 2 ) . {\displaystyle R_{n}^{m}(\rho )={\frac {1}{\binom {\lfloor n/2\rfloor }{\lfloor m/2\rfloor }}}\rho ^{n\mod 2}\sum _{s=\lfloor m/2\rfloor }^{\lfloor n/2\rfloor }(-1)^{\lfloor n/2\rfloor -s}{\binom {s}{\lfloor m/2\rfloor }}{\binom {(n+m)/2}{s+\lceil m/2\rceil }}b_{s,\lfloor n/2\rfloor }(\rho ^{2}).}

Índexs seqüencials de Noll

Les aplicacions sovint impliquen l'ús de l'àlgebra lineal, on les integrals sobre productes de polinomis de Zernike i algun altre factor es poden organitzar com els elements d'una matriu. Una relació per enumerar les files i les columnes d'aquestes matrius mitjançant un sol índex va ser introduïda per Noll.[3] La transformació convencional dels dos índexs n i m en un únic índex j mitjançant l'associació Z n m Z j {\displaystyle Z_{n}^{m}\rightarrow Z_{j}} comença de la següent manera:

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,-2 5,5
j 11 12 13 14 15 16 17 18 19 20

La regla és que per a Z parell (amb la part azimutal parell m, cos ( m φ ) {\displaystyle \cos(m\varphi )} ) s'obtenen els índexs j parells, i per a Z imparella s'obtenen els índexs j imparells. Dins d'un n donat, els valors més baixos de |m| produeixen els menors valors de j.

Índexs estàndard OSA / ANSI

Els polinomis de Zernike d'un sol índex utilitzen els coeficients de la Societat Òptica Nord-americana i del ANSI:[4]

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

Índexs de Fringe / Zemax

Zemax usa l'esquema d'indexació de Fringe. Els 20 primers nombres de Fringe s'enumeren a continuació.[5]

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

Propietats

Ortogonalitat

L'ortogonalitat en la part radial s'expressa com

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

L'ortogonalitat en la part angular està representada per les integrals elementals

0 2 π cos ( m φ ) cos ( m φ ) d φ = ϵ m π δ | m | , | m | , {\displaystyle \int _{0}^{2\pi }\cos(m\varphi )\cos(m'\varphi )\,d\varphi =\epsilon _{m}\pi \delta _{|m|,|m'|},}
0 2 π sin ( m φ ) sin ( m φ ) d φ = ( 1 ) m + m π δ | m | , | m | ; m 0 , {\displaystyle \int _{0}^{2\pi }\sin(m\varphi )\sin(m'\varphi )\,d\varphi =(-1)^{m+m'}\pi \delta _{|m|,|m'|};\quad m\neq 0,}
0 2 π cos ( m φ ) sin ( m φ ) d φ = 0 , {\displaystyle \int _{0}^{2\pi }\cos(m\varphi )\sin(m'\varphi )\,d\varphi =0,}

on ϵ m {\displaystyle \epsilon _{m}} (de vegades anomenat factor de Neumann perquè apareix amb freqüència juntament amb les funcions de Bessel) es defineix com 2 si m = 0 {\displaystyle m=0} i com 1 si m 0 {\displaystyle m\neq 0} . El producte de les parts angulars i radials estableix la ortogonalitat de les funcions de Zernike pel que fa a tots dos índexs si s'integra en el disc unitat,

Z n m ( ρ , φ ) Z n m ( ρ , φ ) d 2 r = ϵ m π 2 n + 2 δ n , n δ m , m , {\displaystyle \int Z_{n}^{m}(\rho ,\varphi )Z_{n'}^{m'}(\rho ,\varphi )\,d^{2}r={\frac {\epsilon _{m}\pi }{2n+2}}\delta _{n,n'}\delta _{m,m'},}

on d 2 r = ρ d ρ d φ {\displaystyle d^{2}r=\rho \,d\rho \,d\varphi } és el jacobià del sistema de coordenades circulars, i on n m {\displaystyle n-m} i n m {\displaystyle n'-m'} són parells.

Un valor especial és

R n m ( 1 ) = 1 , {\displaystyle R_{n}^{m}(1)=1,\,}

Transformada de Zernike

Qualsevol camp de fase de valor real prou uniforme sobre el disc de la unitat G ( ρ , φ ) {\displaystyle G(\rho ,\varphi )} pot representar-se en termes dels seus coeficients de Zernike (imparell i parell), de la mateixa manera que les funcions periòdiques troben una representació ortogonal amb la sèrie de Fourier. Sent:

G ( ρ , φ ) = m , n [ a m , n Z n m ( ρ , φ ) + b m , n Z n m ( ρ , φ ) ] , {\displaystyle G(\rho ,\varphi )=\sum _{m,n}\left[a_{m,n}Z_{n}^{m}(\rho ,\varphi )+b_{m,n}Z_{n}^{-m}(\rho ,\varphi )\right],}

els coeficients es poden calcular usant productes interns. A l'espai de les funcions del L 2 {\displaystyle L^{2}} disc de la unitat, existeix un producte intern definit per

F , G := F ( ρ , φ ) G ( ρ , φ ) ρ d ρ d φ . {\displaystyle \langle F,G\rangle :=\int F(\rho ,\varphi )G(\rho ,\varphi )\rho d\rho d\varphi .}

Els coeficients de Zernike es poden expressar de la següent manera:

a m , n = 2 n + 2 ϵ m π G ( ρ , φ ) , Z n m ( ρ , φ ) , b m , n = 2 n + 2 ϵ m π G ( ρ , φ ) , Z n m ( ρ , φ ) . {\displaystyle {\begin{aligned}a_{m,n}&={\frac {2n+2}{\epsilon _{m}\pi }}\left\langle G(\rho ,\varphi ),Z_{n}^{m}(\rho ,\varphi )\right\rangle ,\\b_{m,n}&={\frac {2n+2}{\epsilon _{m}\pi }}\left\langle G(\rho ,\varphi ),Z_{n}^{-m}(\rho ,\varphi )\right\rangle .\end{aligned}}}

Alternativament, es poden usar els valors coneguts de la funció de fase G en el reticle circular per formar un sistema d'equacions. La funció de fase es recupera mitjançant el producte ponderat del coeficient desconegut amb (valors coneguts) del polinomi de Zernike en el reticle del disc unitat. Per tant, els coeficients també es poden trobar resolent un sistema lineal, per exemple, mitjançant la inversió d'una matriu. Els algorismes ràpids per calcular la transformació de Zernike directa i inversa utilitzen les propietats de simetria de les funcions trigonomètriques, la separabilitat de les parts radials i azimutals dels polinomis de Zernike i les seves simetries rotacionals.

Simetries

La paritat pel que fa a la reflexió en l'eix x és

Z n m ( ρ , φ ) = ( 1 ) m Z n m ( ρ , φ ) . {\displaystyle Z_{n}^{m}(\rho ,\varphi )=(-1)^{m}Z_{n}^{m}(\rho ,-\varphi ).}

La paritat pel que fa al punt de reflexió al centre de coordenades és

Z n m ( ρ , φ ) = ( 1 ) m Z n m ( ρ , φ + π ) , {\displaystyle Z_{n}^{m}(\rho ,\varphi )=(-1)^{m}Z_{n}^{m}(\rho ,\varphi +\pi ),}

on ( 1 ) m {\displaystyle (-1)^{m}} també podria escriure's ( 1 ) n {\displaystyle (-1)^{n}} perquè n m {\displaystyle n-m} és parell per als valors rellevants que no tendeixen a zero. Els polinomis radials també són parells o imparells, segons l'ordre n o m:

R n m ( ρ ) = ( 1 ) n R n m ( ρ ) = ( 1 ) m R n m ( ρ ) . {\displaystyle R_{n}^{m}(\rho )=(-1)^{n}R_{n}^{m}(-\rho )=(-1)^{m}R_{n}^{m}(-\rho ).}

La periodicitat de les funcions trigonomètriques implica constància si és trencada per múltiples de 2 π / m {\displaystyle 2\pi /m} radiants al voltant del centre:

Z n m ( ρ , φ + 2 π k m ) = Z n m ( ρ , φ ) , k = 0 , ± 1 , ± 2 , . {\displaystyle Z_{n}^{m}\left(\rho ,\varphi +{\tfrac {2\pi k}{m}}\right)=Z_{n}^{m}(\rho ,\varphi ),\qquad k=0,\pm 1,\pm 2,\cdots .}

Relacions de recurrència

Els polinomis de Zernike satisfan la següent relació de recurrència que no depèn ni del grau ni de l'ordre azimutal dels polinomis radials:[6]

R n m ( ρ ) + R n 2 m ( ρ ) = ρ [ R n 1 | m 1 | ( ρ ) + R n 1 m + 1 ( ρ ) ]  . {\displaystyle {\begin{aligned}R_{n}^{m}(\rho )+R_{n-2}^{m}(\rho )=\rho \left[R_{n-1}^{\left|m-1\right|}(\rho )+R_{n-1}^{m+1}(\rho )\right]{\text{ .}}\end{aligned}}}

De la definició de R n m {\displaystyle R_{n}^{m}} es pot veure que R m m ( ρ ) = ρ m {\displaystyle R_{m}^{m}(\rho )=\rho ^{m}} i R m + 2 m ( ρ ) = ( ( m + 2 ) ρ 2 ( m + 1 ) ) ρ m {\displaystyle R_{m+2}^{m}(\rho )=((m+2)\rho ^{2}-(m+1))\rho ^{m}} . La següent relació de recurrència de tres termes permet calcular tots els altres R n m ( ρ ) {\displaystyle R_{n}^{m}(\rho )} :[7]

R n m ( ρ ) = 2 ( n 1 ) ( 2 n ( n 2 ) ρ 2 m 2 n ( n 2 ) ) R n 2 m ( ρ ) n ( n + m 2 ) ( n m 2 ) R n 4 m ( ρ ) ( n + m ) ( n m ) ( n 2 )  . {\displaystyle R_{n}^{m}(\rho )={\frac {2(n-1)(2n(n-2)\rho ^{2}-m^{2}-n(n-2))R_{n-2}^{m}(\rho )-n(n+m-2)(n-m-2)R_{n-4}^{m}(\rho )}{(n+m)(n-m)(n-2)}}{\text{ .}}}

La relació anterior és especialment útil, ja que la derivada de R n m {\displaystyle R_{n}^{m}} es pot calcular a partir de dos polinomis de Zernike radials de grau adjacent:[7]

d d ρ R n m ( ρ ) = ( 2 n m ( ρ 2 1 ) + ( n m ) ( m + n ( 2 ρ 2 1 ) ) ) R n m ( ρ ) ( n + m ) ( n m ) R n 2 m ( ρ ) 2 n ρ ( ρ 2 1 )  . {\displaystyle {\frac {\operatorname {d} }{\operatorname {d} \!\rho }}R_{n}^{m}(\rho )={\frac {(2nm(\rho ^{2}-1)+(n-m)(m+n(2\rho ^{2}-1)))R_{n}^{m}(\rho )-(n+m)(n-m)R_{n-2}^{m}(\rho )}{2n\rho (\rho ^{2}-1)}}{\text{ .}}}

Exemples

Polinomis radials

Els primers pocs polinomis radials són:

R 0 0 ( ρ ) = 1 {\displaystyle R_{0}^{0}(\rho )=1\,}
R 1 1 ( ρ ) = ρ {\displaystyle R_{1}^{1}(\rho )=\rho \,}
R 2 0 ( ρ ) = 2 ρ 2 1 {\displaystyle R_{2}^{0}(\rho )=2\rho ^{2}-1\,}
R 2 2 ( ρ ) = ρ 2 {\displaystyle R_{2}^{2}(\rho )=\rho ^{2}\,}
R 3 1 ( ρ ) = 3 ρ 3 2 ρ {\displaystyle R_{3}^{1}(\rho )=3\rho ^{3}-2\rho \,}
R 3 3 ( ρ ) = ρ 3 {\displaystyle R_{3}^{3}(\rho )=\rho ^{3}\,}
R 4 0 ( ρ ) = 6 ρ 4 6 ρ 2 + 1 {\displaystyle R_{4}^{0}(\rho )=6\rho ^{4}-6\rho ^{2}+1\,}
R 4 2 ( ρ ) = 4 ρ 4 3 ρ 2 {\displaystyle R_{4}^{2}(\rho )=4\rho ^{4}-3\rho ^{2}\,}
R 4 4 ( ρ ) = ρ 4 {\displaystyle R_{4}^{4}(\rho )=\rho ^{4}\,}
R 5 1 ( ρ ) = 10 ρ 5 12 ρ 3 + 3 ρ {\displaystyle R_{5}^{1}(\rho )=10\rho ^{5}-12\rho ^{3}+3\rho \,}
R 5 3 ( ρ ) = 5 ρ 5 4 ρ 3 {\displaystyle R_{5}^{3}(\rho )=5\rho ^{5}-4\rho ^{3}\,}
R 5 5 ( ρ ) = ρ 5 {\displaystyle R_{5}^{5}(\rho )=\rho ^{5}\,}
R 6 0 ( ρ ) = 20 ρ 6 30 ρ 4 + 12 ρ 2 1 {\displaystyle R_{6}^{0}(\rho )=20\rho ^{6}-30\rho ^{4}+12\rho ^{2}-1\,}
R 6 2 ( ρ ) = 15 ρ 6 20 ρ 4 + 6 ρ 2 {\displaystyle R_{6}^{2}(\rho )=15\rho ^{6}-20\rho ^{4}+6\rho ^{2}\,}
R 6 4 ( ρ ) = 6 ρ 6 5 ρ 4 {\displaystyle R_{6}^{4}(\rho )=6\rho ^{6}-5\rho ^{4}\,}
R 6 6 ( ρ ) = ρ 6 . {\displaystyle R_{6}^{6}(\rho )=\rho ^{6}.\,}

Polinomis de Zernike

Es mostren a continuació algunes de les primeres maneres de Zernike, amb índexs GOSA/ANSI i índexs únics de Noll. Estan normalitzats de tal manera que

0 2 π 0 1 Z j 2 ρ d ρ d θ = π . {\displaystyle \int _{0}^{2\pi }\int _{0}^{1}Z_{j}^{2}\,\rho \,d\rho \,d\theta =\pi .}
  Índex

OSA/ANSI

( j {\displaystyle j} )

Índex

Noll

( j {\displaystyle j} )

Grau radial

( j {\displaystyle j} )

Grau azimutal

( j {\displaystyle j} )

Z j {\displaystyle Z_{j}} Nom clàssic
Z 0 0 {\displaystyle Z_{0}^{0}}  0  1 0  0 1 {\displaystyle 1} Pistó (vegeu distribució semicircular de Wigner)
Z 1 1 {\displaystyle Z_{1}^{-1}}  1  3 1 −1 2 ρ sin θ {\displaystyle 2\rho \sin \theta } Inclinació (inclinació I, inclinació vertical)
Z 1 1 {\displaystyle Z_{1}^{1}}  2  2 1 +1 2 ρ cos θ {\displaystyle 2\rho \cos \theta } Inclinació horitzontal
Z 2 2 {\displaystyle Z_{2}^{-2}}  3  5 2 −2 6 ρ 2 sin 2 θ {\displaystyle {\sqrt {6}}\rho ^{2}\sin 2\theta } Astigmatisme oblic
Z 2 0 {\displaystyle Z_{2}^{0}}  4  4 2  0 3 ( 2 ρ 2 1 ) {\displaystyle {\sqrt {3}}(2\rho ^{2}-1)} Desenfoc (posició longitudinal)
Z 2 2 {\displaystyle Z_{2}^{2}}  5  6 2 +2 6 ρ 2 cos 2 θ {\displaystyle {\sqrt {6}}\rho ^{2}\cos 2\theta } Astigmatisme vertical
Z 3 3 {\displaystyle Z_{3}^{-3}}  6  9 3 −3 8 ρ 3 sin 3 θ {\displaystyle {\sqrt {8}}\rho ^{3}\sin 3\theta } Lobulat vertical
Z 3 1 {\displaystyle Z_{3}^{-1}}  7  7 3 −1 8 ( 3 ρ 3 2 ρ ) sin θ {\displaystyle {\sqrt {8}}(3\rho ^{3}-2\rho )\sin \theta } Coma vertical
Z 3 1 {\displaystyle Z_{3}^{1}}  8  8 3 +1 8 ( 3 ρ 3 2 ρ ) cos θ {\displaystyle {\sqrt {8}}(3\rho ^{3}-2\rho )\cos \theta } Coma horitzontal
Z 3 3 {\displaystyle Z_{3}^{3}}  9 10 3 +3 8 ρ 3 cos 3 θ {\displaystyle {\sqrt {8}}\rho ^{3}\cos 3\theta } Lobulat oblic
Z 4 4 {\displaystyle Z_{4}^{-4}} 10 15 4 −4 10 ρ 4 sin 4 θ {\displaystyle {\sqrt {10}}\rho ^{4}\sin 4\theta } Cuatrilobulat oblic
Z 4 2 {\displaystyle Z_{4}^{-2}} 11 13 4 −2 10 ( 4 ρ 4 3 ρ 2 ) sin 2 θ {\displaystyle {\sqrt {10}}(4\rho ^{4}-3\rho ^{2})\sin 2\theta } Astigmatisme secundari oblic
Z 4 0 {\displaystyle Z_{4}^{0}} 12 11 4  0 5 ( 6 ρ 4 6 ρ 2 + 1 ) {\displaystyle {\sqrt {5}}(6\rho ^{4}-6\rho ^{2}+1)} Esfèrica primària
Z 4 2 {\displaystyle Z_{4}^{2}} 13 12 4 +2 10 ( 4 ρ 4 3 ρ 2 ) cos 2 θ {\displaystyle {\sqrt {10}}(4\rho ^{4}-3\rho ^{2})\cos 2\theta } Astigmatisme vertical secundari
Z 4 4 {\displaystyle Z_{4}^{4}} 14 14 4 +4 10 ρ 4 cos 4 θ {\displaystyle {\sqrt {10}}\rho ^{4}\cos 4\theta } Cuatrilobulat vertical

Aplicacions

Els polinomis de Zernike són una base definida sobre una àrea de suport circular, típicament els plànols de les pupil·les en imatges òptiques clàssiques en longituds d'ona visibles i infraroges, a través de sistemes de lents i miralls de diàmetre finit. El seu principal avantatge procedeix de les propietats analítiques simples heretades de la senzillesa de les funcions radials i de la factorització en funcions radials i azimutals; això porta, per exemple, a expressions de forma tancada de la transformada de Fourier bidimensional en termes de funcions de Bessel.[8][9] El seu desavantatge, en particular si estan involucrats n termes, és la distribució desigual de les línies nodals sobre el disc unitat, la qual cosa introdueix efectes de ressonància prop del perímetre ρ 1 {\displaystyle \rho \approx 1} , que sovint condueixen a la necessitat de definir altres funcions ortogonals sobre el disc circular.[10]

En la fabricació òptica de precisió, els polinomis de Zernike s'utilitzen per caracteritzar els errors d'ordre superior observats en les anàlisis interferomètriques.

En optometria i oftalmologia, els polinomis de Zernike s'usen per descriure aberracions de la còrnia o del cristal·lí des d'una forma esfèrica ideal, que dona com a resultat ametropies.

S'usen comunament en òptica adaptativa, on es poden emplear per calibrar la distorsió atmosfèrica. Les aplicacions habituals per a aquesta propietat es troben en l'astronomia visual o infraroja i en el tractament d'imatges provinents de satèl·lits.

Una altra aplicació dels polinomis de Zernike es troba en la teoria estesa de Nijboer-Zernike sobre difracció i aberracions òptiques.

Els polinomis de Zernike també s'usen àmpliament com a funcions de base de moments d'imatge. Com els polinomis de Zernike són ortogonals entre sí, els moments de Zernike poden representar les propietats d'una imatge sense redundància ni superposició d'informació entre les diferents maneres. Encara que els moments de Zernike depenen significativament de l'escalat i de la translació de l'objecte en una regió d'interès, les seves magnituds són independents de l'angle de rotació de l'objecte. Per tant, poden utilitzar-se per extreure propietats d'imatges que descriuen la forma característiques d'un objecte. Per exemple, els moments de Zernike s'utilitzen com a descriptores de forma per classificar i identificar càncers de mama benignes i malignes en imatges digitalitzades o en la superfície de discos vibratoris.[11][12] Els moments de Zernike també s'han usat per quantificar la forma de les línies cel·lulars de càncer d'osteosarcoma en el nivell d'una sola cèl·lula.[13]

Dimensions més altes

El concepte es tradueix a dimensions majors D si els multinomis x 1 i x 2 j x D k {\displaystyle x_{1}^{i}x_{2}^{j}\cdots x_{D}^{k}} en coordenades cartesianes es converteixen en coordenades hiperesfèriques, ρ s , s D {\displaystyle \rho ^{s},s\leq D} , multiplicades per un producte de polinomis de Jacobi de les variables angulars. Per exemple, en la dimensió D = 3 {\displaystyle D=3} , les variables angulars són harmònics esfèrics. Combinacions lineals de les potències ρ s {\displaystyle \rho ^{s}} defineixen una base ortogonal R n ( l ) ( ρ ) {\displaystyle R_{n}^{(l)}(\rho )} que satisfà

. 0 1 ρ D 1 R n ( l ) ( ρ ) R n ( l ) ( ρ ) d ρ = δ n , n {\displaystyle \int _{0}^{1}\rho ^{D-1}R_{n}^{(l)}(\rho )R_{n'}^{(l)}(\rho )d\rho =\delta _{n,n'}}

(Tingui's en compte que un factor 2 n + D {\displaystyle {\sqrt {2n+D}}} s'absorbeix aquí en la definició de R, mentre que en D = 2 {\displaystyle D=2} la normalització es tria de forma lleugerament diferent. Això és en gran manera una qüestió arbitrària, depenent de si es desitja mantenir un conjunt sencer de coeficients o es prefereixen fórmules més estrictes si està involucrada la ortogonalització). La representació explícita és

R n ( l ) ( ρ ) = 2 n + D s = 0 n l 2 ( 1 ) s ( n l 2 s ) ( n s 1 + D 2 n l 2 ) ρ n 2 s = ( 1 ) n l 2 2 n + D s = 0 n l 2 ( 1 ) s ( n l 2 s ) ( s 1 + n + l + D 2 n l 2 ) ρ 2 s + l = ( 1 ) n l 2 2 n + D ( n + l + D 2 1 n l 2 ) ρ l   2 F 1 ( n l 2 , n + l + D 2 ; l + D 2 ; ρ 2 ) {\displaystyle {\begin{aligned}R_{n}^{(l)}(\rho )&={\sqrt {2n+D}}\sum _{s=0}^{\tfrac {n-l}{2}}(-1)^{s}{{\tfrac {n-l}{2}} \choose s}{n-s-1+{\tfrac {D}{2}} \choose {\tfrac {n-l}{2}}}\rho ^{n-2s}\\&=(-1)^{\tfrac {n-l}{2}}{\sqrt {2n+D}}\sum _{s=0}^{\tfrac {n-l}{2}}(-1)^{s}{{\tfrac {n-l}{2}} \choose s}{s-1+{\tfrac {n+l+D}{2}} \choose {\tfrac {n-l}{2}}}\rho ^{2s+l}\\&=(-1)^{\tfrac {n-l}{2}}{\sqrt {2n+D}}{{\tfrac {n+l+D}{2}}-1 \choose {\tfrac {n-l}{2}}}\rho ^{l}\ {}_{2}F_{1}\left(-{\tfrac {n-l}{2}},{\tfrac {n+l+D}{2}};l+{\tfrac {D}{2}};\rho ^{2}\right)\end{aligned}}}

fins i tot per n l 0 {\displaystyle n-l\geq 0} , o en cas contrari, idèntic a zero.

Vegeu també

  • Polinomis de Jacobi
  • Teoria de Nijboer-Zernike
  • Polinomis pseudo-Zernike

Referències

  1. Zernike, F. «Beugungstheorie des Schneidenverfahrens und Seiner Verbesserten Form, der Phasenkontrastmethode». Physica, 1, 1934, pàg. 689–704. Bibcode: 1934Phy.....1..689Z. DOI: 10.1016/S0031-8914(34)80259-5.
  2. Born, Max; Wolf, Emil. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light (en anglès). 7a. Cambridge, UK: Cambridge University Press, 1999, p. 986. ISBN 9780521642224. 
  3. Noll, R. J. «Zernike polynomials and atmospheric turbulence». J. Opt. Soc. Am., 66, 3, 1976, pàg. 207. Bibcode: 1976JOSA...66..207N. DOI: 10.1364/JOSA.66.000207.[Enllaç no actiu]
  4. Thibos, L. N.; Applegate, R. A.; Schwiegerling, J. T.; Webb, R. «Standards for reporting the optical aberrations of eyes». Journal of Refractive Surgery, 18, 5, 2002, pàg. S652-60.
  5. Proc SPIE 4771, p.276-286 (2002) doi:10.1117/12.482169
  6. Honarvar Shakibaei Asli, Barmak; Raveendran, Paramesran (July 2013). "Recursive formula to compute Zernike radial polynomials" Opt. Lett. (OSA) 38 (14): 2487–2489. doi:10.1364/OL.38.002487
  7. 7,0 7,1 Kintner, E. C. «On the mathematical properties of the Zernike Polynomials». Opt. Acta, 23, 8, 1976, pàg. 679–680. Bibcode: 1976AcOpt..23..679K. DOI: 10.1080/713819334.
  8. Tatulli, E. «Transformation of Zernike coefficients: a Fourier-based method for scaled, translated, and rotated wavefront apertures». J. Opt. Soc. Am. A, 30, 4, 2013, pàg. 726. Bibcode: 2013JOSAA..30..726T. DOI: 10.1364/JOSAA.30.000726.
  9. Janssen, A. J. E. M. «New analytic results for the Zernike Circle Polynomials from a basic result in the Nijboer-Zernike diffraction theory». JEOS:RP, 6, 2011. DOI: 10.2971/jeos.2011.11028.
  10. Barakat, Richard «Optimum balanced wave-front aberrations for radially symmetric amplitude distributions: Generalizations of Zernike polynomials». J. Opt. Soc., 70, 6, 1980, pàg. 739-742. Bibcode: 1980JOSA...70..739B. DOI: 10.1364/JOSA.70.000739.
  11. Tahmasbi, A.; Saki, F.; Shokouhi, S.B. «Classification of Benign and Malignant Masses Based on Zernike Moments». Computers in Biology and Medicine, 41, 2011, pàg. 726–735. DOI: 10.1016/j.compbiomed.2011.06.009.
  12. Rdzanek, W. P. «Sound radiation of a vibrating elastically supported circular plate embedded into a flat screen revisited using the Zernike circle polynomials». J. Sound Vibr., 434, 2018, pàg. 91-125. Bibcode: 2018JSV...434...92R. DOI: 10.1016/j.jsv.2018.07.035.
  13. Alizadeh, Elaheh; Lyons, Samanthe M; Castle, Jordan M; Prasad, Ashok «Measuring systematic changes in invasive cancer cell shape using Zernike moments». Integrative Biology, 8, 11, 2016, pàg. 1183–1193. DOI: 10.1039/C6IB00100A.

Bibliografia

  • Weisstein, Eric W., «Polinomis de Zernike» a MathWorld (en anglès).
  • Andersen, Torben B. «Efficient and robust recurrence relations for the Zernike circle polynomials and their derivatives in Cartesian coordinates». Opt. Exp., 26, 2018, pàg. 18878-18896. DOI: 10.1364/OE.26.018878.
  • Bhatia, A. B.; Wolf, E. «The Zernike circle polynomials occurring in diffraction theory». Proc. Phys. Soc. B, 65, 11, 1952, pàg. 909–910. Bibcode: 1952PPSB...65..909B. DOI: 10.1088/0370-1301/65/11/112.
  • Callahan, P. G.; De Graef, M. «Precipitate shape fitting and reconstruction by means of 3D Zernike functions». Model. Simul. Mat. Sci. Engin., 20, 2012, pàg. 015003. Bibcode: 2012MSMSE..20a5003C. DOI: 10.1088/0965-0393/20/1/015003.
  • Campbell, C. E. «Matrix method to find a new set of Zernike coefficients form an original set when the aperture radius is changed». J. Opt. Soc. Am. A, 20, 2, 2003, pàg. 209. Bibcode: 2003JOSAA..20..209C. DOI: 10.1364/JOSAA.20.000209.
  • Cerjan, C. «The Zernike-Bessel representation and its application to Hankel transforms». J. Opt. Soc. Am. A, 24, 6, 2007, pàg. 1609. Bibcode: 2007JOSAA..24.1609C. DOI: 10.1364/JOSAA.24.001609.
  • Comastri, S. A.; Perez, L. I.; Perez, G. D.; Martin, G.; Bastida Cerjan, K. «Zernike expansion coefficients: rescaling and decentering for different pupils and evaluation of corneal aberrations». J. Opt. Soc. Am. A, 9, 3, 2007, pàg. 209–221. Bibcode: 2007JOptA...9..209C. DOI: 10.1088/1464-4258/9/3/001.
  • Conforti, G. «Zernike aberration coefficients from Seidel and higher-order power-series coefficients». Opt. Lett., 8, 7, 1983, pàg. 407–408. Bibcode: 1983OptL....8..407C. DOI: 10.1364/OL.8.000407.
  • Dai, G-m.; Mahajan, V. N. «Zernike annular polynomials and atmospheric turbulence». J. Opt. Soc. Am. A, 24, 2007, pàg. 139. Bibcode: 2007JOSAA..24..139D. DOI: 10.1364/JOSAA.24.000139.
  • Dai, G-m. «Scaling Zernike expansion coefficients to smaller pupil sizes: a simpler formula». J. Opt. Soc. Am. A, 23, 3, 2006, pàg. 539. Bibcode: 2006JOSAA..23..539D. DOI: 10.1364/JOSAA.23.000539.
  • Díaz, J. A.; Fernández-Dorado, J.; Pizarro, C.; Arasa, J. «Zernike Coefficients for Concentric, Circular, Scaled Pupils: An Equivalent Expression». Journal of Modern Optics, 56, 1, 2009, pàg. 149–155. Bibcode: 2009JMOp...56..149D. DOI: 10.1080/09500340802531224.
  • «Zernike Coefficients for Concentric, Circular, Scaled Pupils». del Proyecto de Demostraciones Wolfram.
  • Farokhi, Sajad; Shamsuddin, Siti Mariyam; Flusser, Jan; Sheikh, U.U; Khansari, Mohammad «Rotation and noise invariant near-infrared face recognition by means of Zernike moments and spectral regression discriminant analysis». Journal of Electronic Imaging, 22, 1, 2013, pàg. 013030. Bibcode: 2013JEI....22a3030F. DOI: 10.1117/1.JEI.22.1.013030.
  • Gu, J.; Shu, H. Z.; Toumoulin, C.; Luo, L. M. «A novel algorithm for fast computation of Zernike moments». Pattern Recogn., 35, 2002, pàg. 2905–2911. DOI: 10.1016/S0031-3203(01)00194-7.
  • Herrmann, J. «Cross coupling and aliasing in modal wave-front estimation». J. Opt. Soc. Am., 71, 8, 1981, pàg. 989. Bibcode: 1981JOSA...71..989H. DOI: 10.1364/JOSA.71.000989.
  • Hu, P. H.; Stone, J.; Stanley, T. «Application of Zernike polynomials to atmospheric propagation problems». J. Opt. Soc. Am. A, 6, 10, 1989, pàg. 1595. Bibcode: 1989JOSAA...6.1595H. DOI: 10.1364/JOSAA.6.001595.
  • Kintner, E. C. «On the mathematical properties of the Zernike Polynomials». Opt. Acta, 23, 8, 1976, pàg. 679–680. Bibcode: 1976AcOpt..23..679K. DOI: 10.1080/713819334.
  • Lawrence, G. N.; Chow, W. W. «Wave-front tomography by Zernike Polynomial decomposition». Opt. Lett., 9, 7, 1984, pàg. 267. Bibcode: 1984OptL....9..267L. DOI: 10.1364/OL.9.000267.
  • Liu, Haiguang; Morris, Richard J.; Hexemer, A.; Grandison, Scott; Zwart, Peter H. «Computation of small-angle scattering profiles with three-dimensional Zernike polynomials». Acta Crystallogr. A, 68, 2, 2012, pàg. 278–285. DOI: 10.1107/S010876731104788X.
  • Lundström, L.; Unsbo, P. «Transformation of Zernike coefficients: scaled, translated and rotated wavefronts with circular and elliptical pupils». J. Opt. Soc. Am. A, 24, 3, 2007, pàg. 569. Bibcode: 2007JOSAA..24..569L. DOI: 10.1364/JOSAA.24.000569.
  • Mahajan, V. N. «Zernike annular polynomials for imaging systems with annular pupils». J. Opt. Soc. Am., 71, 1981, pàg. 75. Bibcode: 1981JOSA...71...75M. DOI: 10.1364/JOSA.71.000075.
  • . 

math.NA.

  • Mathar, R. J. «Zernike Basis to Cartesian Transformations». Serbian Astronomical Journal, 179, 179, 2009, pàg. 107–120. arXiv: 0809.2368. Bibcode: 2009SerAJ.179..107M. DOI: 10.2298/SAJ0979107M.
  • Prata Jr, A.; Rusch, W. V. T. «Algorithm for computation of Zernike polynomials expansion coefficients». Appl. Opt., 28, 4, 1989, pàg. 749. Bibcode: 1989ApOpt..28..749P. DOI: 10.1364/AO.28.000749.
  • Schwiegerling, J. «Scaling Zernike expansion coefficients to different pupil sizes». J. Opt. Soc. Am. A, 19, 10, 2002, pàg. 1937. Bibcode: 2002JOSAA..19.1937S. DOI: 10.1364/JOSAA.19.001937.
  • Sheppard, C. J. R.; Campbell, S.; Hirschhorn, M. D. «Zernike expansion of separable functions in Cartesian coordinates». Appl. Opt., 43, 20, 2004, pàg. 3963. Bibcode: 2004ApOpt..43.3963S. DOI: 10.1364/AO.43.003963.
  • Shu, H.; Luo, L.; Han, G.; Coatrieux, J.-L. «General method to derive the relationship between two sets of Zernike coefficients corresponding to different aperture sizes». J. Opt. Soc. Am. A, 23, 8, 2006, pàg. 1960. Bibcode: 2006JOSAA..23.1960S. DOI: 10.1364/JOSAA.23.001960.
  • Swantner, W.; Chow, W. W. «Gram-Schmidt orthogonalization of Zernike polynomials for general aperture shapes». Appl. Opt., 33, 10, 1994, pàg. 1832. Bibcode: 1994ApOpt..33.1832S. DOI: 10.1364/AO.33.001832.
  • Tango, W. J. «The circle polynomials of Zernike and their application in optics». Appl. Phys. A, 13, 4, 1977, pàg. 327–332. Bibcode: 1977ApPhy..13..327T. DOI: 10.1007/BF00882606.
  • Tyson, R. K. «Conversion of Zernike aberration coefficients to Seidel and higher-order power series aberration coefficients». Opt. Lett., 7, 6, 1982, pàg. 262. Bibcode: 1982OptL....7..262T. DOI: 10.1364/OL.7.000262.
  • Wang, J. Y.; Silva, D. E. «Wave-front interpretation with Zernike Polynomials». Appl. Opt., 19, 9, 1980, pàg. 1510. Bibcode: 1980ApOpt..19.1510W. DOI: 10.1364/AO.19.001510.
  • Barakat, R. «Optimum balanced wave-front aberrations for radially symmetric amplitude distributions: Generalizations of Zernike polynomials». J. Opt. Soc. Am., 70, 6, 1980, pàg. 739. Bibcode: 1980JOSA...70..739B. DOI: 10.1364/JOSA.70.000739.
  • ten Brummelaar, T. A. «Modeling atmospheric wave aberrations and astronomical instrumentation using the polynomials of Zernike». Opt. Commun., 132, 3–4, 1996, pàg. 329–342. Bibcode: 1996OptCo.132..329T. DOI: 10.1016/0030-4018(96)00407-5.
  • Novotni, M.; Klein, R. «3D Zernike Descriptors for Content Based Shape Retrieval». Proceedings of the 8th ACM Symposium on Solid Modeling and Applications. DOI: 10.1145/781606.781639.
  • Novotni, M.; Klein, R. «Shape retrieval using 3D Zernike descriptors». Computer Aided Design, 36, 11, 2004, pàg. 1047–1062. DOI: 10.1016/j.cad.2004.01.005.
  • Farokhi, Sajad; Shamsuddin, Siti Mariyam; Sheikh, U.U; Flusser, Jan «Near Infrared Face Recognition: A Comparison of Moment-Based Approaches». Lecture Notes in Electrical Engineering. Springer, 291, 1, 2014, pàg. 129–135. DOI: 10.1007/978-981-4585-42-2_15.
  • Farokhi, Sajad; Shamsuddin, Siti Mariyam; Flusser, Jan; Sheikh, U.U; Khansari, Mohammad «Near infrared face recognition by combining Zernike moments and undecimated discrete wavelet transform». Digital Signal Processing, 31, 1, 2014, pàg. 13–27. DOI: 10.1016/j.dsp.2014.04.008.