Polinomios de Zernike

En matemáticas, los polinomios de Zernike son una secuencia de polinomios que son ortogonales en el disco unitario. Llevan el nombre del físico óptico Frits Zernike, premio Nobel de Física de 1953 e inventor de la microscopía de contraste de fases, y desempeñan un papel importante en diversas ramas de la óptica, como la óptica de haces y la formación de imágenes.
Definiciones
Hay polinomios de Zernike pares e impares. Los polinomios pares de Zernike se definen como
- Znm()*** *** ,φ φ )=Rnm()*** *** )# ()mφ φ ){displaystyle Z_{n} {m}(rhovarphi)=R_{n} {m}rho),cos(m,varphi)!}
(incluso la función sobre el ángulo azimutal φ φ {displaystyle varphi }), y los extraños polinomios Zernike se definen como
- Zn− − m()*** *** ,φ φ )=Rnm()*** *** )pecado ()mφ φ ),{displaystyle Z_{n} {-m}(rhovarphi)=R_{n} {m}rho),sin(m,varphi),!}
(función anida sobre el ángulo azimutal φ φ {displaystyle varphi }Donde m y n son enteros no negativos con ≥ m ≥ 0 ()m = 0 para polinomios esféricos Zernike), φ φ {displaystyle varphi } es el ángulo azimutal, *** es la distancia radial 0≤ ≤ *** *** ≤ ≤ 1{displaystyle 0leq rho leq 1}, y Rnm{displaystyle R_{n} {m} son los polinomios radiales definidos a continuación. Los polinomios Zernike tienen la propiedad de ser limitados a una gama de −1 a +1, es decir. SilencioZnm()*** *** ,φ φ )Silencio≤ ≤ 1{displaystyle TENZ_{n} {rhovarphi) 1}. Los polinomios radiales Rnm{displaystyle R_{n} {m} se definen como
- Rnm()*** *** )=. . k=0n− − m2()− − 1)k()n− − k)!k!()n+m2− − k)!()n− − m2− − k)!*** *** n− − 2k{displaystyle R_{n} {m}(rho)=sum ¿Por qué? {n-m}{2}{frac {(-1)^{k},(n-k)}{k!left({tfrac] {n+m}{2}-kright) {n-m}{2}-kright)};rho ^{n-2k}
para un número par de n − m, mientras que es 0 para un número impar de n − m. Un valor especial es
- Rnm()1)=1.{displaystyle R_{n} {m}(1)=1.}
Otras representaciones
Reescribir las razones de factoriales en la parte radial como productos de binomios muestra que los coeficientes son números enteros:
- Rnm()*** *** )=. . k=0n− − m2()− − 1)k()n− − kk)()n− − 2kn− − m2− − k)*** *** n− − 2k{displaystyle R_{n} {m}(rho)=sum ¿Por qué? {n-m}{2}(-1)}{k}{binom {n-k}{binom} {n-2k}{tfrac} {n-m}{2}-k}rho ^{n-2k}.
Una notación como terminación de funciones hipergeométricas gaussianas es útil para revelar recurrencias, demostrar que son casos especiales de polinomios de Jacobi, escribir las ecuaciones diferenciales, etc.:
- Rnm()*** *** )=()nn+m2)*** *** n 2F1()− − n+m2,− − n− − m2;− − n;*** *** − − 2)=()− − 1)n− − m2()n+m2m)*** *** m 2F1()1+n+m2,− − n− − m2;1+m;*** *** 2){fnMicrosoft Sans Serif} {fn} {fn} {fnMicrosoft Sans Serif} {n+m}{2}rho ^{n} {fnK} {cHFF} {cHFF}cHFF}cHFF}ccHFF}cHFF}ccHFF} {n+m}{2},-{tfrac {n-m}{2};-n;rho ################################################################################################################################################################################################################################################################ {n-m}{2}{binom {tfrac} {n+m}{m}}rho ^{m} {}_{2}F_{1}left(1+{tfrac {n+m}{2},-{tfrac {n-m}{2};1+m;rho ^{2}end{aligned}}}
para n − m par.
La relación inversa se expande *** *** j{displaystyle rho ^{j} fija m≤ ≤ j{displaystyle mleq j} en Rnm()*** *** ){displaystyle ¿Qué?
- *** *** j=. . n↑ ↑ m()mod2)jhj,n,mRnm()*** *** ){displaystyle rho ^{j}=sum _{nequiv m{pmod {2}} {j}h_{j,n,m}R_{n} {m} {rho)}
con coeficientes racionales hj,n,m{displaystyle h_{j,n,m}
- hj,n,m=n+11+j+n2()()j− − m)/2()n− − m)/2)()()j+n)/2()n− − m)/2){displaystyle h_{j,n,m}={ {n+1}{1+{frac {fn} {fn} {fn}} {fn}}} {fn}}} {fn}}} {fn}}}} {fn}}}}} {fn}}}}}}} {fn}}} {fn}}}}}}} {fnf}}}}} {f}}}}}}}}}}}}}}}}}}}}} {f} {f}}}}}} {f}}}}} {f}}} {f}} { {binom {} {binom}{} {binom}{} {nse}{} {nse} {nse}}}}}}}}} {benom {binom} {nse]} {cH00}}}}} {benom} {benom}}} {benom}}}}} {
para incluso j− − m=0,2,4,... ... {displaystyle J-m=0,2,4,ldots.
El factor *** *** n− − 2k{displaystyle rho ^{n-2k} en el polinomio radial Rnm()*** *** ){displaystyle ¿Qué? puede ampliarse en una base de Bernstein bs,n/2()*** *** 2){displaystyle b_{s,n/2}(rho ^{2}} para incluso n{displaystyle n} o *** *** {displaystyle rho } tiempos una función bs,()n− − 1)/2()*** *** 2){displaystyle b_{s,(n-1)/2}(rho ^{2} por extraño n{displaystyle n} en el rango ⌊ ⌊ n/2⌋ ⌋ − − k≤ ≤ s≤ ≤ ⌊ ⌊ n/2⌋ ⌋ {displaystyle lfloor n/2rfloor -kleq sleq lfloor n/2rfloor }. El polinomio radial puede ser expresado por un número finito de los polinomios de Bernstein con coeficientes racionales:
- Rnm()*** *** )=1()⌊ ⌊ n/2⌋ ⌋ ⌊ ⌊ m/2⌋ ⌋ )*** *** nmod2. . s=⌊ ⌊ m/2⌋ ⌋ ⌊ ⌊ n/2⌋ ⌋ ()− − 1)⌊ ⌊ n/2⌋ ⌋ − − s()s⌊ ⌊ m/2⌋ ⌋ )()()n+m)/2s+⌈ ⌈ m/2⌉ ⌉ )bs,⌊ ⌊ n/2⌋ ⌋ ()*** *** 2).{displaystyle ¿Qué? - ¿Qué? _{s=lfloor m/2rfloor }{lfloor n/2rfloor }(-1)^{lfloor n/2rfloor - Sí. {fnMicrosoft} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} m/2rfloor # {binom {(n+m)/2}{s+lceil m/2rceil }b_{s,lfloor n/2rfloor }(rho ^{2}). }
Índices secuenciales de Noll
Las aplicaciones a menudo implican álgebra lineal, donde un integral sobre un producto de polinomios Zernike y algún otro factor construye un elemento matriz.
Para enumerar las filas y columnas de estas matrices por un solo índice, un mapeo convencional de los dos índices n y l a un único índice j ha sido presentado por Noll. El cuadro de esta asociación Znl→ → Zj{displaystyle Z_{n}{l}rightarrow Z_{j} comienza como sigue (secuencia) A176988 en el OEIS).
0land nequiv {0,1}{pmod {4}};\0,&lj=n()n+1)2+SilenciolSilencio+{}0,l■0∧ ∧ n↑ ↑ {}0,1}()mod4);0,lc)0∧ ∧ n↑ ↑ {}2,3}()mod4);1,l≥ ≥ 0∧ ∧ n↑ ↑ {}2,3}()mod4);1,l≤ ≤ 0∧ ∧ n↑ ↑ {}0,1}()mod4).{fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fn} {fnMicrosoft Sans Serif} {fn0} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft}fnMicrosoft} {fnMicrosoft}f}fnMicrosoft}f}fnMicrosoft}fnMicrosoft}fnMicrosoft ##f}fnMicrosoft #########fnMicrosoft ## ## #################################### ## ######## ## ############## #########################0land nequiv {0,1}{pmod {4}};\0,&l
n,l | 0,0 | 1.1 | 1, 1 - 1 | 2,0 | 2,2 a 2 | 2,2 | 3,1 a 1 | 3,1 | 3, 3 a 3 | 3,3 |
---|---|---|---|---|---|---|---|---|---|---|
j | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
n,l | 4.0 | 4,2 | 4, 2 | 4.4 | 4,4 a 4 | 5,1 | 5,1 a | 5,3 | 5, 3, 3 | 5,5 |
j | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 |
La regla es la siguiente.
- Los polinomios incluso Zernike Z (con partes azimutales # ()mφ φ ){displaystyle cos(mvarphi)}, donde m=l{displaystyle m=l} como l{displaystyle l} es un número positivo) obtener incluso índices j.
- Lo extraño Z obtiene (con partes azimutales impares pecado ()mφ φ ){displaystyle sin(mvarphi)}, donde m=SilenciolSilencio{displaystyle m=leftvert lrightvert} como l{displaystyle l} es un número negativo) índices impares j.
- Dentro de un n, un bajo SilenciolSilencio{displaystyle leftvert lrightvert} resultados en un menor j.
Índices estándar OSA/ANSI
AOS y polinomios de Zernike de índice único ANSI usando:
- j=n()n+2)+l2{displaystyle j={frac {n+2)+l}{2}} {f}}
n,l | 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,l | 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 |
Índices Fringe/Universidad de Arizona
El esquema de indexación Fringe se utiliza en software comercial de diseño óptico y pruebas ópticas en, por ejemplo, fotolitografía.
j=()1+n+SilenciolSilencio2)2− − 2SilenciolSilencio+1− − Sgn l2{displaystyle j=left(1+{frac {fnfn} {2}}derecha)} {2}}}}
Donde Sgn l{displaystyle operatorname {sgn} l} es la función Signum o Signum. Los primeros 20 números de franja se enumeran a continuación.
n,l | 0,0 | 1.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,l | 3, 3 | 4,2 | 4, 2 | 5,1 | 5,1 a | 6,0 | 4.4 | 4, 4 | 5,3 | 5,-3 |
j | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 |
Índices de Wyant
James C. Wyant utiliza el estilo "Fringe" esquema de indexación excepto que comienza en 0 en lugar de 1 (resta 1). Este método se utiliza habitualmente, incluido el software de análisis de interferogramas en los interferómetros Zygo y el software de código abierto DFTFringe.
Propiedades
Ortogonalidad
La ortogonalidad en la parte radial dice
- ∫ ∫ 012n+2Rnm()*** *** )2n.+2Rn.m()*** *** )*** *** d*** *** =δ δ n,n.{displaystyle int _{0}{1}{sqrt {2n+2}R_{n} {m}(rho),{sqrt {2n'+2}R_{n'} {m}(rho),rho drho =delta _{n'n}}}}} {rho)}
o
∫ ∫ 10Rnm()*** *** )Rn.m()*** *** )*** *** d*** *** =δ δ n,n.2n+2.{displaystyle {underset {0}{ncip}{nMithop {int} ¿Qué? {deleta {n,{n}} {2n+2}}
La ortogonalidad en la parte angular está representada por la elemental
- ∫ ∫ 02π π # ()mφ φ )# ()m.φ φ )dφ φ =ε ε mπ π δ δ m,m.,{displaystyle int _{0}^{2pi }cos(mvarphi)cos(m'varphi),dvarphi =epsilon ¿Por qué?
- ∫ ∫ 02π π pecado ()mφ φ )pecado ()m.φ φ )dφ φ =π π δ δ m,m.;mل ل 0,{displaystyle int _{0}^{2pi}sin(mvarphi)sin(m'varphi),dvarphi =pi delta _{m,m'};quad mneq 0,}
- ∫ ∫ 02π π # ()mφ φ )pecado ()m.φ φ )dφ φ =0,{displaystyle int _{0}{2pi }cos(mvarphi)sin(m'varphi),dvarphi =0,}
Donde ε ε m{displaystyle epsilon _{m} (a veces llamado factor Neumann porque aparece con frecuencia en conjunto con las funciones de Bessel) se define como 2 si m=0{displaystyle m=0} y 1 si mل ل 0{displaystyle mneq 0}. El producto de las piezas angulares y radiales establece la ortogonalidad de las funciones de Zernike con respecto a ambos índices si se integra sobre el disco de unidad,
- ∫ ∫ Znl()*** *** ,φ φ )Zn.l.()*** *** ,φ φ )d2r=ε ε lπ π 2n+2δ δ n,n.δ δ l,l.,{displaystyle int Z_{n} {rhovarphi)Z_{n'}(rhovarphi),d^{2}r={frac {epsilon _{l}pi }{2n+2}delta _{n,n'}delta _{l,l'}}} {d}} {d}}} {d}}} {d}}}}}}}}}}d}d}d}d}d}d}d}d}cccccccccccccccccccccncnccncccncccnccnccccncncccncccncccncncccn
Donde d2r=*** *** d*** *** dφ φ {displaystyle d^{2}r=rho ,drho ,dvarphi } es el Jacobiano del sistema de coordenadas circulares, y donde n− − l{displaystyle n-l} y n.− − l.{displaystyle n'-l'} ambos son incluso.
Transformación de Zernike
Cualquier campo de fase de valor real suficientemente suave sobre el disco de unidad G()*** *** ,φ φ ){displaystyle G(rhovarphi)} puede ser representado en términos de sus coeficientes Zernike (odd e incluso), así como las funciones periódicas encuentran una representación ortogonal con la serie Fourier. Tenemos
- G()*** *** ,φ φ )=. . m,n[am,nZnm()*** *** ,φ φ )+bm,nZn− − m()*** *** ,φ φ )],{displaystyle G(rhovarphi)=sum _{m,n}left[a_{m,n}Z_{n}{m}(rhovarphi)+b_{m,n}Z_{n}{-m}(rhovarphi)right],}
donde los coeficientes se pueden calcular utilizando productos internos. En el espacio de L2{displaystyle L^{2} funciones en el disco de unidad, hay un producto interior definido por
- . . F,G. . :=∫ ∫ F()*** *** ,φ φ )G()*** *** ,φ φ )*** *** d*** *** dφ φ .{displaystyle langle F,Grangle:=int F(rhovarphi)G(rhovarphi)rho drho dvarphi.}
Los coeficientes de Zernike se pueden expresar de la siguiente manera:
- am,n=2n+2ε ε mπ π .G()*** *** ,φ φ ),Znm()*** *** ,φ φ ).,bm,n=2n+2ε ε mπ π .G()*** *** ,φ φ ),Zn− − m()*** *** ,φ φ )..{displaystyle {begin{aligned}a_{m,n} {2n+2}{epsilon _{m}pi }leftlangle G(rhovarphi),Z_{n}^{m}(rhovarphi)rightrangle\b_{m,n} limit={frac] {2n+2}{epsilon _{m}i}leftlangle G(rhovarphi),Z_{n}{-m}(rhovarphi)rightrangle.
Alternativamente, se pueden usar los valores conocidos de la función de fase G en la cuadrícula circular para formar un sistema de ecuaciones. La función de fase se recupera mediante el producto ponderado de coeficiente desconocido con (valores conocidos) del polinomio de Zernike en la cuadrícula de unidades. Por lo tanto, los coeficientes también se pueden encontrar resolviendo un sistema lineal, por ejemplo mediante inversión de matrices. Los algoritmos rápidos para calcular la transformada de Zernike directa e inversa utilizan propiedades de simetría de funciones trigonométricas, separabilidad de partes radiales y azimutales de polinomios de Zernike y sus simetrías rotacionales.
Simetrías
Las reflexiones de funciones trigonométricas dan como resultado que la paridad con respecto a la reflexión a lo largo del eje x es
- Znl()*** *** ,φ φ )=Znl()*** *** ,− − φ φ ){displaystyle Z_{n} {rhovarphi)=Z_{n} {l}(rho-varphi)} para l ≥ 0,
- Znl()*** *** ,φ φ )=− − Znl()*** *** ,− − φ φ ){displaystyle Z_{n} {rhovarphi)=-Z_{n}(rho-varphi)} para l 0.
Los desplazamientos π de funciones trigonométricas dan como resultado que la paridad con respecto a la reflexión puntual en el centro de las coordenadas es
- Znl()*** *** ,φ φ )=()− − 1)lZnl()*** *** ,φ φ +π π ),{displaystyle ¿Por qué?
Donde ()− − 1)l{displaystyle (-1)} bien podría ser escrito ()− − 1)n{displaystyle (-1)} porque n− − l{displaystyle n-l} como incluso los números son sólo casos para obtener polinomios Zernike no-vanishing. (Si) n es incluso entonces l también lo es. Si n es extraño, entonces l también es extraño.) Esta propiedad se utiliza a veces para clasificar polinomios Zernike en polinomios uniformes y extraños en términos de su dependencia angular. (también es posible añadir otra categoría con l = 0 ya que tiene una propiedad especial de ninguna dependencia angular.)
- Polinomios anulares incluso zernike: polinomios Zernike con incluso l así Znl()*** *** ,φ φ )=Znl()*** *** ,φ φ +π π ).{displaystyle Z_{n} {rhovarphi)=Z_{n}(rhovarphi +pi).}
- Polinomios Zernike anularmente extraños: polinomios Zernike con extraños l así Znl()*** *** ,φ φ )=− − Znl()*** *** ,φ φ +π π ).{displaystyle Z_{n} {rhovarphi)=-Z_{n}(rhovarphi +pi).}
Los polinomios radiales también son pares o impares, dependiendo del orden n o m:
- Rnm()*** *** )=()− − 1)nRnm()− − *** *** )=()− − 1)mRnm()− − *** *** ).{displaystyle ¿Qué?
Estas igualdades se ven fácilmente desde Rnm()*** *** ){displaystyle ¿Qué? con un extraño (incluso) m contiene sólo poderes extraños (incluso) *** (véanse ejemplos de Rnm()*** *** ){displaystyle ¿Qué? infra).
La periodicidad de las funciones trigonométricas resulta en invariancia si rota por múltiplos de 2π π /l{displaystyle 2pi /l} radian alrededor del centro:
- Znl()*** *** ,φ φ +2π π kl)=Znl()*** *** ,φ φ ),k=0,± ± 1,± ± 2,⋯ ⋯ .{displaystyle Z_{n}{l}left(rhovarphi +{tfrac {2pi k}{l}right)=Z_{n}{l}(rhovarphi),qquad k=0,pm 1,cdots.}
Relaciones de recurrencia
Los polinomios de Zernike satisfacen la siguiente relación de recurrencia que no depende ni del grado ni del orden azimutal de los polinomios radiales:
- Rnm()*** *** )+Rn− − 2m()*** *** )=*** *** [Rn− − 1Silenciom− − 1Silencio()*** *** )+Rn− − 1m+1()*** *** )].{fnMicrosoft Sans Serif} {fn0}fn1} {fn2}rho)=rho left[R_{n-1}}{left perpetuam-1right alive}(rho)+R_{n-1}{}{}{}{}{}{}{}{}{}{}{}{}{}{}}}{}{}{}}}}}}{}{}{}}}{}}}}{}}}}{}}}}{}}{}{}}}}}}{}}}}}}{}}}}}}}}}}}}}}}{}}}}}}}}}}}{}}}}}}{}}}}}}}}}}}}{}}}}}}}}}}}}}}}}}}}}}}}}}}}}}}{}}}}}}{}}}}}}}}}}}}}}}
De la definición Rnm{displaystyle R_{n} {m} se puede ver que Rmm()*** *** )=*** *** m{displaystyle ¿Qué? y Rm+2m()*** *** )=()()m+2)*** *** 2− − ()m+1))*** *** m{displaystyle (m+2)rho ^{2}-(m+1))rho. La siguiente relación de recurrencia a tres plazos permite calcular todo lo demás Rnm()*** *** ){displaystyle ¿Qué?:
- Rnm()*** *** )=2()n− − 1)()2n()n− − 2)*** *** 2− − m2− − n()n− − 2))Rn− − 2m()*** *** )− − n()n+m− − 2)()n− − m− − 2)Rn− − 4m()*** *** )()n+m)()n− − m)()n− − 2).{displaystyle R_{n}{m}(rho)={frac {2(n-1)(2n(n-2)rho ¿Qué?
La relación anterior es especialmente útil ya que el derivado de Rnm{displaystyle R_{n} {m} se puede calcular a partir de dos radiales Polinomios Zernike de grado adyacente:
- dd*** *** Rnm()*** *** )=()2nm()*** *** 2− − 1)+()n− − m)()m+n()2*** *** 2− − 1)))Rnm()*** *** )− − ()n+m)()n− − m)Rn− − 2m()*** *** )2n*** *** ()*** *** 2− − 1).{displaystyle {frac {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif}}}n}n}n} {n}n} {n}}n} {n}}n} {n}n} {n}n} {n}n} {n}}n} {n}}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n} {nn}n} {n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n}n
Ejemplos
Polinomios radiales
Los primeros polinomios radiales son:
- R00()*** *** )=1{displaystyle ¿Qué?
- R11()*** *** )=*** *** {displaystyle R_{1} {rho]=rho ,}
- R20()*** *** )=2*** *** 2− − 1{displaystyle ¿Qué?
- R22()*** *** )=*** *** 2{displaystyle R_{2} {rho]=rho ^{2},}
- R31()*** *** )=3*** *** 3− − 2*** *** {displaystyle ¿Qué?
- R33()*** *** )=*** *** 3{displaystyle R_{3} {rho]=rho ^{3},}
- R40()*** *** )=6*** *** 4− − 6*** *** 2+1{displaystyle R_{4} {0}=6rho ^{4}-6rho ^{2}+1,}
- R42()*** *** )=4*** *** 4− − 3*** *** 2{displaystyle ¿Qué?
- R44()*** *** )=*** *** 4{displaystyle R_{4} {rho]=rho ^{4},}
- R51()*** *** )=10*** *** 5− − 12*** *** 3+3*** *** {displaystyle ¿Qué?
- R53()*** *** )=5*** *** 5− − 4*** *** 3{displaystyle R_{5}{3}(rho)=5rho ^{5}-4rho ^{3},}
- R55()*** *** )=*** *** 5{displaystyle R_{5} {rho)=rho ^{5},}
- R60()*** *** )=20*** *** 6− − 30*** *** 4+12*** *** 2− − 1{displaystyle R_{6} {0}=20rho ^{6}-30rho ^{4}+12rho ^{2}-1}
- R62()*** *** )=15*** *** 6− − 20*** *** 4+6*** *** 2{displaystyle ¿Qué?
- R64()*** *** )=6*** *** 6− − 5*** *** 4{displaystyle ¿Qué?
- R66()*** *** )=*** *** 6.{displaystyle R_{6} {rho)=rho ^{6}
Polinomios de Zernike
Los primeros pocos modos Zernike, en varios índices, se muestran a continuación. Se normalizan de tal manera que: ∫ ∫ 02π π ∫ ∫ 01Z2⋅ ⋅ *** *** d*** *** dφ φ =π π {displaystyle int _{0}{2pi}int ¿Por qué?, que equivale a Var ()Z)círculo de unidad=1{displaystyle operatorname {fnMicrosoft Sans Serif}=1}.
Znl{displaystyle Z_{n} {}} {fn}} | OSA/ANSI índice ()j{displaystyle j}) | Noll índice ()j{displaystyle j}) | Wyant índice ()j{displaystyle j}) | Fringe/UA índice ()j{displaystyle j}) | Radial grado ()n{displaystyle n}) | Azimuthal grado ()l{displaystyle l}) | Zj{displaystyle Z_{j} | Nombre clásico |
---|---|---|---|---|---|---|---|---|
Z00{displaystyle Z_{0} {0} | 00 | 01 | 00 | 01 | 0 | 00 | 1{displaystyle 1} | Piston (ver, distribución de semicírculos Wigner) |
Z1− − 1{displaystyle Z_{1} {-1} | 01 | 03 | 02 | 03 | 1 | −1 | 2*** *** pecado φ φ {displaystyle 2rho sin phi } | Tilt (Y-Tilt, inclinación vertical) |
Z11{displaystyle Z_{1} {1}} | 02 | 02 | 01 | 02 | 1 | + 1 | 2*** *** # φ φ {displaystyle 2rho cos phi } | Tilt (X-Tilt, inclinación horizontal) |
Z2− − 2{displaystyle Z_{2} {-2} | 03 | 05 | 05 | 06 | 2 | −2 | 6*** *** 2pecado 2φ φ {displaystyle {sqrt {6}rho ^{2}sin 2phi } | Astigmatismo oblicuo |
Z20{displaystyle Z_{2} {0} | 04 | 04 | 03 | 04 | 2 | 00 | 3()2*** *** 2− − 1){displaystyle {sqrt {3}(2rho ^{2}-1)} | Defocus (Posición longitudinal) |
Z22{displaystyle Z_{2} {2}} | 05 | 06 | 04 | 05 | 2 | +2 | 6*** *** 2# 2φ φ {displaystyle {sqrt {6}rho ^{2}cos 2phi} | Astigmatismo vertical |
Z3− − 3{displaystyle Z_{3} {-3} | 06 | 09 | 10 | 11 | 3 | −3 | 8*** *** 3pecado 3φ φ {displaystyle {sqrt {8}rho ^{3}sin 3phi } | Trefil vertical |
Z3− − 1{displaystyle Z_{3} {-1} | 07 | 07 | 07 | 08 | 3 | −1 | 8()3*** *** 3− − 2*** *** )pecado φ φ {displaystyle {sqrt {8}}(3rho ^{3}-2rho)sin phi } | coma vertical |
Z31{displaystyle Z_{3} {1} | 08 | 08 | 06 | 07 | 3 | + 1 | 8()3*** *** 3− − 2*** *** )# φ φ {displaystyle {sqrt {8}}(3rho ^{3}-2rho)cos phi } | coma horizontal |
Z33{displaystyle Z_{3} {3} | 09 | 10 | 09 | 10 | 3 | +3 | 8*** *** 3# 3φ φ {displaystyle {sqrt {8}rho }cos 3phi } | Trefoil oblicuo |
Z4− − 4{displaystyle Z_{4} {-4} | 10 | 15 | 17 | 18 | 4 | −4 | 10*** *** 4pecado 4φ φ {displaystyle {sqrt}rho ^{4}sin 4phi } | Quadrafoil oblicuo |
Z4− − 2{displaystyle Z_{4} {-2} | 11 | 13 | 12 | 13 | 4 | −2 | 10()4*** *** 4− − 3*** *** 2)pecado 2φ φ {displaystyle {sqrt {10}(4rho ^{4}-3rho ^{2}sin 2phi } | Astigmatismo secundario oblicuo |
Z40{displaystyle Z_{4} {0} | 12 | 11 | 08 | 09 | 4 | 00 | 5()6*** *** 4− − 6*** *** 2+1){displaystyle {sqrt {5}(6rho ^{4}-6rho ^{2}+1)} | Esférica primaria |
Z42{displaystyle Z_{4} {2} | 13 | 12 | 11 | 12 | 4 | +2 | 10()4*** *** 4− − 3*** *** 2)# 2φ φ {displaystyle {sqrt {10}(4rho ^{4}-3rho ^{2}cos 2phi } | Astigmatismo secundario vertical |
Z44{displaystyle Z_{4} {4}} | 14 | 14 | 16 | 17 | 4 | +4 | 10*** *** 4# 4φ φ {displaystyle {sqrt}rho ^{4}cos 4phi } | Cuadrícula vertical |
Aplicaciones

Las funciones son una base definida sobre el área de soporte circular, típicamente los planos del alumno en imágenes ópticas clásicas a longitudes de onda visibles e infrarrojas a través de sistemas de lentes y espejos de diámetro finito. Sus ventajas son las propiedades analíticas simples heredadas de la simplicidad de las funciones radiales y la factorización en funciones radiales y azimutales; esto conduce, por ejemplo, a expresiones de forma cerrada de la transformación bidimensional de Fourier en términos de funciones de Bessel. Su desventaja, en particular si es alta n están involucrados, es la distribución desigual de líneas nodal sobre el disco unitario, que introduce efectos de anillo cerca del perímetro *** *** . . 1{displaystyle rho approx 1}, que a menudo conduce intentos de definir otras funciones ortogonales sobre el disco circular.
En la fabricación óptica de precisión, los polinomios de Zernike se utilizan para caracterizar los errores de orden superior observados en los análisis interferométricos. En sensores de pendiente de frente de onda como el Shack-Hartmann, los coeficientes de Zernike del frente de onda se pueden obtener ajustando las pendientes medidas con derivadas polinómicas de Zernike promediadas sobre las subaperturas de muestreo. En optometría y oftalmología, los polinomios de Zernike se utilizan para describir aberraciones del frente de onda de la córnea o del cristalino a partir de una forma esférica ideal, que dan lugar a errores de refracción. También se utilizan habitualmente en óptica adaptativa, donde pueden utilizarse para caracterizar la distorsión atmosférica. Las aplicaciones obvias para esto son la astronomía visual o IR y las imágenes de satélite.
Otra aplicación de los polinomios de Zernike se encuentra en la teoría extendida de difracción y aberraciones de Nijboer-Zernike.
Los polinomios de Zernike se utilizan ampliamente como funciones base de momentos de imágenes. Dado que los polinomios de Zernike son ortogonales entre sí, los momentos de Zernike pueden representar propiedades de una imagen sin redundancia ni superposición de información entre los momentos. Aunque los momentos de Zernike dependen significativamente de la escala y la traslación del objeto en una región de interés (ROI), sus magnitudes son independientes del ángulo de rotación del objeto. Por lo tanto, se pueden utilizar para extraer características de imágenes que describen las características de forma de un objeto. Por ejemplo, los momentos de Zernike se utilizan como descriptores de forma para clasificar masas mamarias benignas y malignas o la superficie de discos vibrantes. Los Momentos Zernike también se han utilizado para cuantificar la forma de líneas celulares de cáncer de osteosarcoma a nivel de células individuales. Además, los Momentos Zernike se han utilizado para la detección temprana de la enfermedad de Alzheimer extrayendo información discriminativa de las imágenes de resonancia magnética de la enfermedad de Alzheimer, el deterioro cognitivo leve y los grupos sanos.
Dimensiones superiores
El concepto se traduce en dimensiones superiores D si multinomios x1ix2j⋯ ⋯ xDk{displaystyle x_{1} {i}x_{2} {j}cdots # en las coordenadas cartesianas se convierten en coordenadas hiperesféricas, *** *** s,s≤ ≤ D{displaystyle rho ^{s},sleq D}, multiplicado por un producto de polinomios Jacobi de las variables angulares. In D=3{displaystyle D=3} dimensiones, las variables angulares son armónicas esféricas, por ejemplo. Combinaciones lineales de los poderes *** *** s{displaystyle rho ^{s} definir una base ortogonal Rn()l)()*** *** ){displaystyle R_{n} {rho)} satisfacción
- ∫ ∫ 01*** *** D− − 1Rn()l)()*** *** )Rn.()l)()*** *** )d*** *** =δ δ n,n.{displaystyle int ¿Por qué?.
(Nota que un factor 2n+D{displaystyle {sqrt {2n+D}} se absorbe en la definición de R aquí, mientras que en D=2{displaystyle D=2} la normalización es elegida ligeramente diferente. Esto es en gran parte una cuestión de gusto, dependiendo de si uno desea mantener un conjunto entero de coeficientes o prefiere fórmulas más estrictas si la ortogonalización está implicada.) La representación explícita es
- Rn()l)()*** *** )=2n+D. . s=0n− − l2()− − 1)s()n− − l2s)()n− − s− − 1+D2n− − l2)*** *** n− − 2s=()− − 1)n− − l22n+D. . s=0n− − l2()− − 1)s()n− − l2s)()s− − 1+n+l+D2n− − l2)*** *** 2s+l=()− − 1)n− − l22n+D()n+l+D2− − 1n− − l2)*** *** l 2F1()− − n− − l2,n+l+D2;l+D2;*** *** 2){displaystyle {begin{aligned}R_{n} {l)} {rho)} {sqrt {2n+D}}sum ¿Por qué? {cHFF} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft} {fnMicrosoft} {fnMicrosoft Sans} {fnMicrosoft Sans}} {fnMicrosoft}}}}} {f}}}}}}}}}}} {f}}}}}} {f}}} {f}}}}}}} {cf}}}}}}}}}} {f}}}}}}}}} {f}}}}}} {f} {f}} {f}}}}}} {f}} {cf}}}}f}}} {cf}}}}}}}}}}}}}}}}}}}}}}}}}}}} {ccfn {n-l}{2} {n-s-1+{tfrac {D}{2} choose {tfrac} {n-l}{2}rho ################################################################################################################################################################################################################################################################ {n-l}{2}{sqrt {2n+D}sum ¿Por qué? {cHFF} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft Sans Serif} {fnMicrosoft} {fnMicrosoft} {fnMicrosoft Sans} {fnMicrosoft Sans}} {fnMicrosoft}}}}} {f}}}}}}}}}}} {f}}}}}} {f}}} {f}}}}}}} {cf}}}}}}}}}} {f}}}}}}}}} {f}}}}}} {f} {f}} {f}}}}}} {f}} {cf}}}}f}}} {cf}}}}}}}}}}}}}}}}}}}}}}}}}}}} {ccfn {n-l}{2} {S-1+{tfrac {n+l+D}{2} choose {tfrac} {n-l}{2}rho ################################################################################################################################################################################################################################################################ {n-l}{2}{sqrt {2n+D}{tfrac {n+l+D}{2}-1 choose {tfrac} {n-l}{2}}rho ^{l} {fnK} {cHFF} {cHFF}cHFF}cHFF}ccHFF}cHFF}ccHFF} {n-l}{2},{tfrac {n+l+D}{2};l+{tfrac {D}}}rho ^{2}right)end{aligned}}}
para incluso n− − l≥ ≥ 0{displaystyle n-lgeq 0}, más idéntico a cero.
Contenido relacionado
Conjunto vacío
Precisión y exactitud
Historia de la lógica