Polinomios de Zernike

Permiten desarrollar en serie funciones definidas sobre un disco unidad

En matemáticas, los polinomios de Zernike son una secuencia de polinomios que son ortogonales en el disco unidad. Fueron nombrados en honor del físico óptico Frits Zernike, ganador del Premio Nobel de física en 1953 e inventor del microscopio de contraste de fases. Estos polinomios juegan un papel importante en la modelización del comportamiento de haces de luz en un sistema óptico.[1][2]

Los primeros 21 polinomios de Zernike, ordenados verticalmente por grado radial y horizontalmente por grado azimutal

Definiciones

editar

Los polinomios de Zernike se distinguen en función de su paridad. Los términos pares se definen como:

 

y los impares como:

 

donde m y n son números enteros no negativos con n ≥ m, φ es el ángulo acimutal, ρ es la distancia radial   y Rmn son los polinomios radiales definidos a continuación. Los polinomios de Zernike tienen la propiedad de estar limitados a un rango de -1 a +1, es decir,  . Los polinomios radiales Rmn se definen como

 

para n-m par, y son idénticamente 0 para n-m impar.

Otras representaciones

editar

Reescribiendo las relaciones de factoriales en la parte radial como productos de coeficientes binomiales, se demuestra que los coeficientes son números enteros:

 .

La notación como términos de funciones hipergeométricas gaussianas es útil para revelar recurrencias, para demostrar casos especiales de los polinomios de Jacobi, o para reducir ecuaciones diferenciales.

 

para n-m par.

El factor   en el polinomio radial   se puede expandir en una base de Bernstein de   para   par o   multiplicado por una función de   para   impar en el rango  . Por lo tanto, el polinomio radial puede expresarse mediante un número finito de polinomios de Bernstein con coeficientes racionales:

 

Índices secuenciales de Noll

editar

Las aplicaciones a menudo implican el uso del álgebra lineal, donde las integrales sobre productos de polinomios de Zernike y algún otro factor se pueden organizar como los elementos de una matriz. Una relación para enumerar las filas y las columnas de estas matrices mediante un solo índice fue introducida por Noll.[3]​ La transformación convencional de los dos índices n y m en un solo índice j mediante la asociación   comienza de la siguiente manera (sucesión A176988 en 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

La regla es que para Z par (con la parte azimutal par m,  ) se obtienen los índices j pares, y para Z impar se obtienen los índices j impares. Dentro de un n dado, los valores más bajos de |m| producen los menores valores de j.

Índices estándar OSA / ANSI

editar

Los polinomios de Zernike de un solo índice utilizan los coeficientes de la Sociedad Óptica Estadounidense [4]​ y del ANSI:

 
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

Índices de Fringe / Zemax

editar

Zemax usa el esquema de indexación de Fringe. Los 20 primeros números de Fringe se enumeran a continuación.[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

Propiedades

editar

Ortogonalidad

editar

La ortogonalidad en la parte radial se expresa como[6]

 

La ortogonalidad en la parte angular está representada por las integrales elementales

 
 
 

donde   (a veces llamado factor de Neumann porque aparece con frecuencia junto con las funciones de Bessel) se define como 2 si   y como 1 si  . El producto de las partes angulares y radiales establece la ortogonalidad de las funciones de Zernike con respecto a ambos índices si se integra en el disco unidad,

 

donde   es el jacobiano del sistema de coordenadas circulares, y donde   y   son pares.

Un valor especial es

 

Transformada de Zernike

editar

Cualquier campo de fase de valor real suficientemente uniforme sobre el disco de la unidad   puede representarse en términos de sus coeficientes de Zernike (impar y par), del mismo modo que las funciones periódicas encuentran una representación ortogonal con la serie de Fourier. Siendo

 

los coeficientes se pueden calcular usando productos internos. En el espacio de las funciones de   en el disco de la unidad, existe un producto interno definido por

 

Los coeficientes de Zernike se pueden expresar de la siguiente manera:

 

Alternativamente, se pueden usar los valores conocidos de la función de fase G en la retí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 retícula del disco unidad. Por lo tanto, los coeficientes también se pueden encontrar resolviendo un sistema lineal, por ejemplo, mediante la inversión de una matriz. Los algoritmos rápidos para calcular la transformación de Zernike directa e inversa utilizan las propiedades de simetría de las funciones trigonométricas, la separabilidad de las partes radiales y acimutales de los polinomios de Zernike y sus simetrías rotacionales.

Simetrías

editar

La paridad con respecto a la reflexión en el eje x es

 

La paridad con respecto al punto de reflexión en el centro de coordenadas es

 

donde   también podría escribirse   porque   es par para los valores relevantes que no tienden a cero. Los polinomios radiales también son pares o impares, según el orden n o m:

 

La periodicidad de las funciones trigonométricas implica invarianza si se rota por múltiplos de   radianes alrededor del centro:

 

Relaciones de recurrencia

editar

Los polinomios de Zernike satisfacen la siguiente relación de recurrencia que no depende ni del grado ni del orden acimutal de los polinomios radiales:[7]

 

De la definición de   se puede ver que   y  . La siguiente relación de recurrencia de tres términos[8]​ permite calcular todos los demás  :

 

La relación anterior es especialmente útil, ya que la derivada de   se puede calcular a partir de dos polinomios de Zernike radiales de grado adyacente:[8]

 

Ejemplos

editar

Polinomios radiales

editar

Los primeros pocos polinomios radiales son:

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Polinomios de Zernike

editar

Se muestran a continuación algunos de los primeros modos de Zernike, con índices OSA/ANSI y e índices únicos de Noll. Están normalizados de tal manera que

 
  Índice
OSA/ANSI
( )
Índice
Noll
( )
Grado
Radial
( )
Grado
Azimutal
( )
  Nombre clásico
   0  1 0  0   Pistón (véase distribución semicircular de Wigner)
   1  3 1 −1   Inclinación (inclinación Y, inclinación vertical)
   2  2 1 +1   Inclinación horizontal
   3  5 2 −2   Astigmatismo oblicuo
   4  4 2  0   Desenfoque (posición longitudinal)
   5  6 2 +2   Astigmatismo vertical
   6  9 3 −3   Lobulado vertical
   7  7 3 −1   Coma vertical
   8  8 3 +1   Coma horizontal
   9 10 3 +3   Lobulado oblicuo
  10 15 4 −4   Cuatrilobulado oblicuo
  11 13 4 −2   Astigmatismo secundario oblicuo
  12 11 4  0   Esférica primaria
  13 12 4 +2   Astigmatismo vertical secundario
  14 14 4 +4   Cuatrilobulado vertical

Aplicaciones

editar

Los polinomios de Zernike son una base definida sobre un área de soporte circular, típicamente los planos de las pupilas en imágenes ópticas clásicas en longitudes de onda visibles e infrarrojas, a través de sistemas de lentes y espejos de diámetro finito. Su principal ventaja procede de las propiedades analíticas simples heredadas de la sencillez de las funciones radiales y de la factorización en funciones radiales y acimutales; esto lleva, por ejemplo, a expresiones de forma cerrada de la transformada de Fourier bidimensional en términos de funciones de Bessel.[9][10]​ Su desventaja, en particular si están involucrados n términos, es la distribución desigual de las líneas nodales sobre el disco unidad, lo que introduce efectos de resonancia cerca del perímetro  , que a menudo conducen a la necesidad de definir otras funciones ortogonales sobre el disco circular.[11][12][13]

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 optometría y oftalmología, los polinomios de Zernike se usan para describir aberraciones de la córnea o del cristalino desde una forma esférica ideal, que da como resultado ametropías.

Se usan comúnmente en óptica adaptativa, donde se pueden emplear para calibrar la distorsión atmosférica. Las aplicaciones habituales para esta propiedad se encuentran en la astronomía visual o infrarroja y en tratamiento de imágenes satelitales.

Otra aplicación de los polinomios de Zernike se encuentra en la teoría extendida de Nijboer-Zernike sobre difracción y aberraciones ópticas.

Los polinomios de Zernike también se usan ampliamente como funciones de base de momentos de imagen. Como los polinomios de Zernike son ortogonales entre sí, los momentos de Zernike pueden representar las propiedades de una imagen sin redundancia ni superposición de información entre los distintos modos. Aunque los momentos de Zernike dependen significativamente del escalado y de la tranlación del objeto en una región de interés, sus magnitudes son independientes del ángulo de rotación del objeto.[14]​ Por lo tanto, pueden utilizarse para extraer propiedades de imágenes que describen la forma características de un objeto. Por ejemplo, los momentos de Zernike se utilizan como descriptores de forma para clasificar e identificar cánceres de mama benigno y maligno en imágenes digitalizadas[15]​ o la superficie de discos vibratorios.[16]​ Los momentos de Zernike también se han usado para cuantificar la forma de las líneas celulares de cáncer de osteosarcoma en el nivel de una sola célula.[17]

Dimensiones más altas

editar

El concepto se traduce a dimensiones mayores D si los multinomios   en coordenadas cartesianas se convierten en coordenadas hiperesféricas,  , multiplicadas por un producto de polinomios de Jacobi de las variables angulares. Por ejemplo, en la dimensión  , las variables angulares son armónicos esféricos. Combinaciones lineales de las potencias   definen una base ortogonal   que satisface

 .

(Téngase en cuenta que un factor   se absorbe aquí en la definición de R, mientras que en   la normalización se elige de forma ligeramente diferente. Esto es en gran medida una cuestión arbitraria, dependiendo de si se desea mantener un conjunto entero de coeficientes o se prefieren fórmulas más estrictas si está involucrada la ortogonalización). La representación explícita es

 

incluso para  , o de lo contrario, idéntico a cero.

Véase también

editar

Referencias

editar
  1. Zernike, F. (1934). «Beugungstheorie des Schneidenverfahrens und Seiner Verbesserten Form, der Phasenkontrastmethode». Physica 1 (8): 689-704. Bibcode:1934Phy.....1..689Z. doi:10.1016/S0031-8914(34)80259-5. 
  2. Born, Max, and Wolf, Emil (1999). Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light (7th edición). Cambridge, UK: Cambridge University Press. p. 986. ISBN 9780521642224. 
  3. Noll, R. J. (1976). «Zernike polynomials and atmospheric turbulence». J. Opt. Soc. Am. 66 (3): 207. Bibcode:1976JOSA...66..207N. doi:10.1364/JOSA.66.000207. 
  4. Thibos, L. N.; Applegate, R. A.; Schwiegerling, J. T.; Webb, R. (2002). «Standards for reporting the optical aberrations of eyes». Journal of Refractive Surgery 18 (5): S652-60. 
  5. Proc SPIE 4771, p.276-286 (2002) doi 10.1117/12.482169
  6. Lakshminarayanan, V.; Fleck, Andre (2011). «Zernike polynomials: a guide». J. Mod. Opt. 58 (7). pp. 545-561. Bibcode:2011JMOp...58..545L. doi:10.1080/09500340.2011.554896. 
  7. 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
  8. a b Kintner, E. C. (1976). «On the mathematical properties of the Zernike Polynomials». Opt. Acta 23 (8): 679-680. Bibcode:1976AcOpt..23..679K. doi:10.1080/713819334. 
  9. Tatulli, E. (2013). «Transformation of Zernike coefficients: a Fourier-based method for scaled, translated, and rotated wavefront apertures». J. Opt. Soc. Am. A 30 (4): 726. Bibcode:2013JOSAA..30..726T. doi:10.1364/JOSAA.30.000726. 
  10. Janssen, A. J. E. M. (2011). «New analytic results for the Zernike Circle Polynomials from a basic result in the Nijboer-Zernike diffraction theory». JEOS:RP 6. doi:10.2971/jeos.2011.11028. 
  11. Barakat, Richard (1980). «Optimum balanced wave-front aberrations for radially symmetric amplitude distributions: Generalizations of Zernike polynomials». J. Opt. Soc. 70 (6): 739-742. Bibcode:1980JOSA...70..739B. doi:10.1364/JOSA.70.000739. 
  12. Janssen, A. J. E. M. (2011). «A generalization of the Zernike circle polynomials for forward and inverse problems in diffraction theory». arXiv:1110.2369. 
  13. Mathar, R. J. (2018). «Orthogonal basis function over the unit circle with the minimax property». arXiv:1802.09518. 
  14. Tahmasbi, A. (2010). An Effective Breast Mass Diagnosis System using Zernike Moments. 17th Iranian Conf. on Biomedical Engineering (ICBME'2010). Isfahán, Irán: Institute of Electrical and Electronics Engineers. pp. 1-4. doi:10.1109/ICBME.2010.5704941. 
  15. Tahmasbi, A.; Saki, F.; Shokouhi, S.B. (2011). «Classification of Benign and Malignant Masses Based on Zernike Moments». Computers in Biology and Medicine 41: 726-735. doi:10.1016/j.compbiomed.2011.06.009. 
  16. Rdzanek, W. P. (2018). «Sound radiation of a vibrating elastically supported circular plate embedded into a flat screen revisited using the Zernike circle polynomials». J. Sound Vibr. 434: 91-125. Bibcode:2018JSV...434...92R. doi:10.1016/j.jsv.2018.07.035. 
  17. Alizadeh, Elaheh; Lyons, Samanthe M; Castle, Jordan M; Prasad, Ashok (2016). «Measuring systematic changes in invasive cancer cell shape using Zernike moments». Integrative Biology 8 (11): 1183-1193. doi:10.1039/C6IB00100A. 

Bibliografía

editar

 [math.NA]. 

Enlaces externos

editar