stringtranslate.com

Polinomios de Zernike

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

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 , ganador del Premio Nobel de Física en 1953 e inventor de la microscopía de contraste de fases , y desempeñan papeles importantes en varias ramas de la óptica, como la óptica de haces y la formación de imágenes. [1] [2]

Definiciones

Existen polinomios de Zernike pares e impares . Los polinomios de Zernike pares se definen como

(función par sobre el ángulo azimutal ), y los polinomios de Zernike impares se definen como

(función impar sobre el ángulo acimutal ) donde m y n son números enteros no negativos con n ≥ m ≥ 0 ( m = 0 para polinomios esféricos de Zernike), es el ángulo acimutal , ρ es la distancia radial y 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 se definen como

para un número par de nm , mientras que es 0 para un número impar de nm . Un valor especial es

Otras representaciones

Reescribiendo las razones de los factoriales en la parte radial como productos de binomios se muestra que los coeficientes son números enteros:

.

Una notación como funciones hipergeométricas gaussianas terminales es útil para revelar recurrencias, demostrar que son casos especiales de polinomios de Jacobi , escribir las ecuaciones diferenciales, etc.:

para nm par.

La relación inversa se expande para fijo en

con coeficientes racionales [3]

para incluso .

El factor del polinomio radial puede desarrollarse en una base de Bernstein de para pares o multiplicado por una función de para impares 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

Las aplicaciones a menudo involucran álgebra lineal, donde una integral sobre un producto de polinomios de Zernike y algún otro factor construye los elementos de una matriz. Para enumerar las filas y columnas de estas matrices por un índice único, Noll introdujo una asignación convencional de los dos índices n y l a un índice único j . [4] La tabla de esta asociación comienza de la siguiente manera (secuencia A176988 en la OEIS ).

La regla es la siguiente:

Índices estándar OSA/ANSI

Polinomios de Zernike de índice único OSA [5] y ANSI utilizando:

Índices de Fringe/Universidad de Arizona

El esquema de indexación Fringe se utiliza en software de diseño óptico comercial y pruebas ópticas, por ejemplo, en fotolitografía . [6] [7]

donde es la función sign o signum . Los primeros 20 números marginales se enumeran a continuación.

Índices Wyant

James C. Wyant utiliza el esquema de indexación "Fringe", excepto que comienza en 0 en lugar de 1 (restar 1). [8] Este método se utiliza comúnmente, incluido el software de análisis de interferogramas en los interferómetros Zygo y el software de código abierto DFTFringe.

Fórmula de Rodrigues

Satisfacen la fórmula de Rodrigues

y puede relacionarse con los polinomios de Jacobi como

.

Propiedades

Ortogonalidad

La ortogonalidad en la parte radial se lee [9]

o

La ortogonalidad en la parte angular está representada por el elemental

donde (a veces llamado factor de Neumann porque aparece frecuentemente junto con las funciones de Bessel) se define como 2 si y 1 si . El producto de las partes angular y radial establece la ortogonalidad de las funciones de Zernike con respecto a ambos índices si se integran sobre el disco unitario,

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

Transformación de Zernike

Cualquier campo de fase de valor real suficientemente suave sobre el disco unitario se puede representar en términos de sus coeficientes de Zernike (par e impar), de la misma manera que las funciones periódicas encuentran una representación ortogonal con la serie de Fourier . Tenemos

donde los coeficientes se pueden calcular utilizando productos internos . En el espacio de funciones del disco unitario, hay un producto interno definido por

Los coeficientes de Zernike pueden expresarse de la siguiente manera:

Como alternativa, se pueden utilizar 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 coeficientes desconocidos con (valores conocidos) del polinomio de Zernike en toda la cuadrícula unitaria. Por lo tanto, los coeficientes también se pueden encontrar resolviendo un sistema lineal, por ejemplo, mediante la inversión de matrices. Los algoritmos rápidos para calcular la transformada de Zernike directa e inversa utilizan las propiedades de simetría de las funciones trigonométricas , la separabilidad de las partes radial y azimutal de los polinomios de Zernike y sus simetrías rotacionales.

Simetrías

Las reflexiones de las funciones trigonométricas dan como resultado que la paridad con respecto a la reflexión a lo largo del eje x es

para l ≥ 0,
para l < 0.

Los desplazamientos π de las funciones trigonométricas dan como resultado que la paridad con respecto a la reflexión del punto en el centro de las coordenadas sea

donde también se podría escribir porque los números pares son los únicos casos para obtener polinomios de Zernike que no se desvanecen. (Si n es par, entonces l también es par. Si n es impar, entonces l también es impar). Esta propiedad a veces se usa para categorizar polinomios de Zernike en pares e impares en términos de su dependencia angular. (También es posible agregar otra categoría con l = 0 ya que tiene una propiedad especial de no dependencia angular).

Los polinomios radiales también son pares o impares, dependiendo del orden n o m :

Estas igualdades se ven fácilmente ya que con un m impar (par) solo contiene potencias impares (pares) de ρ (ver ejemplos a continuación).

La periodicidad de las funciones trigonométricas resulta invariante si se rotan en múltiplos de radianes alrededor del centro:

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: [10]

De la definición de it se desprende que y . La siguiente relación de recurrencia de tres términos [11] permite entonces 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 radiales de Zernike de grado adyacente: [11]

La ecuación diferencial de la Función Hipergeométrica Gaussiana es equivalente a

Ejemplos

Polinomios radiales

Los primeros polinomios radiales son:

Polinomios de Zernike

A continuación se muestran los primeros modos de Zernike, en varios índices. Están normalizados de modo que: , que es equivalente a .

Aplicaciones

Resultado de los primeros 21 polinomios de Zernike (como el anterior) introducidos como aberraciones en un haz de superficie plana. El haz se capta mediante una lente, lo que produce una transformada de Fourier cuya intensidad se representa en esta imagen.

Las funciones son una base definida sobre el área de soporte circular, típicamente los planos pupilares 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. 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 transformada de Fourier bidimensional en términos de funciones de Bessel. [12] [13] Su desventaja, en particular si se trata de n altos , es la distribución desigual de líneas nodales sobre el disco unitario, lo que introduce efectos de timbre cerca del perímetro , lo que a menudo conduce a intentos de definir otras funciones ortogonales sobre el disco circular. [14] [15] [16]

En la fabricación óptica de precisión, los polinomios de Zernike se utilizan para caracterizar errores de orden superior observados en 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 las derivadas del polinomio de Zernike promediadas sobre las subaperturas de muestreo. [17] En optometría y oftalmología , los polinomios de Zernike se utilizan para describir aberraciones del frente de onda de la córnea o el 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 se pueden utilizar para caracterizar la distorsión atmosférica . Las aplicaciones obvias para esto son la astronomía IR o visual y las imágenes por satélite .

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

Los polinomios de Zernike se utilizan ampliamente como funciones base de los momentos de imagen . Dado que los polinomios de Zernike son ortogonales entre sí, los momentos de Zernike pueden representar propiedades de una imagen sin redundancia o 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. [18] 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 [19] o la superficie de discos vibratorios. [20] Los momentos de Zernike también se han utilizado para cuantificar la forma de las líneas celulares de cáncer de osteosarcoma a nivel de célula única. [21] Además, los momentos de Zernike se han utilizado para la detección temprana de la enfermedad de Alzheimer mediante la extracción de información discriminativa de las imágenes de RM de la enfermedad de Alzheimer, deterioro cognitivo leve y grupos sanos. [22]

Dimensiones superiores

El concepto se traduce a dimensiones superiores 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. En dimensiones, las variables angulares son armónicos esféricos , por ejemplo. Las combinaciones lineales de las potencias definen una base ortogonal que satisface

.

(Tenga en cuenta que aquí se absorbe un factor 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 de gusto, dependiendo de si uno desea mantener un conjunto entero de coeficientes o prefiere fórmulas más estrictas si está involucrada la ortogonalización). La representación explícita es [3]

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

Véase también

Referencias

  1. ^ Zernike, F. (1934). "Beugungstheorie des Schneidenverfahrens und Seiner Verbesserten Form, der Phasenkontrastmethode". Física . 1 (8): 689–704. Código bibliográfico : 1934Phy......1..689Z. doi :10.1016/S0031-8914(34)80259-5.
  2. ^ Born, Max y Wolf, Emil (1999). Principios de óptica: teoría electromagnética de la propagación, interferencia y difracción de la luz (7.ª ed.). Cambridge, Reino Unido: Cambridge University Press. pág. 986. ISBN 9780521642224.(ver también en Google Books)
  3. ^ ab Mathar, RJ (2009). "Base de Zernike para transformaciones cartesianas". Revista astronómica serbia . 179 (179): 107–120. arXiv : 0809.2368 . Código Bibliográfico :2009SerAJ.179..107M. doi :10.2298/SAJ0979107M. S2CID  115159231.
  4. ^ Noll, RJ (1976). "Polinomios de Zernike y turbulencia atmosférica" ​​(PDF) . J. Opt. Soc. Am . 66 (3): 207. Bibcode :1976JOSA...66..207N. doi :10.1364/JOSA.66.000207.
  5. ^ Thibos, LN; Applegate, RA; Schwiegerling, JT; Webb, R. (2002). "Estándares para informar las aberraciones ópticas de los ojos" (PDF) . Journal of Refractive Surgery . 18 (5): S652-60. doi :10.3928/1081-597X-20020901-30. PMID  12361175.
  6. ^ Loomis, J., "Un programa informático para el análisis de datos interferométricos", Interferogramas ópticos, reducción e interpretación, ASTM STP 666, AH Guenther y DH Liebenberg, Eds., Sociedad Americana de Pruebas y Materiales, 1978, págs. 71–86.
  7. ^ Genberg, VL; Michels, GJ; Doyle, KB (2002). "Ortogonalidad de polinomios de Zernike". Diseño optomecánico e ingeniería 2002. Proc SPIE. Vol. 4771. págs. 276–286. doi :10.1117/12.482169.
  8. ^ Eric P. Goodwin; James C. Wyant (2006). Guía de campo para pruebas ópticas interferométricas . pág. 25. ISBN 0-8194-6510-0.
  9. ^ Lakshminarayanan, V.; Fleck, Andre (2011). "Polinomios de Zernike: una guía". J. Mod. Opt . 58 (7): 545–561. Código Bibliográfico :2011JMOp...58..545L. doi :10.1080/09500340.2011.554896. S2CID  120905947.
  10. ^ Honarvar Shakibaei, Barmak (2013). "Fórmula recursiva para calcular polinomios radiales de Zernike". Opt. Lett . 38 (14): 2487–2489. Bibcode :2013OptL...38.2487H. doi :10.1364/OL.38.002487. PMID  23939089.
  11. ^ ab Kintner, EC (1976). "Sobre las propiedades matemáticas de los polinomios de Zernike". Opt. Acta . 23 (8): 679–680. Código Bibliográfico :1976AcOpt..23..679K. doi :10.1080/713819334.
  12. ^ Tatulli, E. (2013). "Transformación de coeficientes de Zernike: un método basado en Fourier para aperturas de frente de onda escaladas, trasladadas y rotadas". J. Opt. Soc. Am. A . 30 (4): 726–32. arXiv : 1302.7106 . Código Bibliográfico :2013JOSAA..30..726T. doi :10.1364/JOSAA.30.000726. PMID  23595334. S2CID  23491106.
  13. ^ Janssen, AJEM (2011). "Nuevos resultados analíticos para los polinomios del círculo de Zernike a partir de un resultado básico en la teoría de difracción de Nijboer-Zernike". Journal of the European Optical Society: Rapid Publications . 6 : 11028. Bibcode :2011JEOS....6E1028J. doi : 10.2971/jeos.2011.11028 .
  14. ^ Barakat, Richard (1980). "Aberraciones de frente de onda óptimas y equilibradas para distribuciones de amplitud radialmente simétricas: generalizaciones de polinomios de Zernike". J. Opt. Soc. Am . 70 (6): 739–742. Bibcode :1980JOSA...70..739B. doi :10.1364/JOSA.70.000739.
  15. ^ Janssen, AJEM (2011). "Una generalización de los polinomios del círculo de Zernike para problemas directos e inversos en la teoría de la difracción". arXiv : 1110.2369 [math-ph].
  16. ^ Mathar, RJ (2018). "Función base ortogonal sobre el círculo unitario con la propiedad minimax". arXiv : 1802.09518 [math.NA].
  17. ^ Akondi, Vyas; Dubra, Alfredo (22 de junio de 2020). "Gradiente promedio de polinomios de Zernike sobre polígonos". Optics Express . 28 (13): 18876–18886. Bibcode :2020OExpr..2818876A. doi : 10.1364/OE.393223 . ISSN  1094-4087. PMC 7340383 . PMID  32672177. 
  18. ^ Tahmasbi, A. (2010). Un sistema eficaz de diagnóstico de masas mamarias mediante momentos de Zernike . 17.ª Conferencia iraní sobre ingeniería biomédica (ICBME'2010). Isfahán , Irán : IEEE . págs. 1–4. doi :10.1109/ICBME.2010.5704941.
  19. ^ Tahmasbi, A.; Saki, F.; Shokouhi, SB (2011). "Clasificación de masas benignas y malignas según los momentos de Zernike". Computers in Biology and Medicine . 41 (8): 726–735. doi :10.1016/j.compbiomed.2011.06.009. PMID  21722886.
  20. ^ Rdzanek, WP (2018). "Radiación sonora de una placa circular elásticamente soportada vibrante incrustada en una pantalla plana revisada utilizando los polinomios del círculo de Zernike". J. Sound Vib . 434 : 91–125. Bibcode :2018JSV...434...92R. doi :10.1016/j.jsv.2018.07.035. S2CID  125512636.
  21. ^ Alizadeh, Elaheh; Lyons, Samanthe M; Castle, Jordan M; Prasad, Ashok (2016). "Medición de cambios sistemáticos en la forma de células cancerosas invasivas utilizando momentos de Zernike". Biología Integrativa . 8 (11): 1183–1193. doi :10.1039/C6IB00100A. PMID  27735002.
  22. ^ Gorji, HT y J. Haddadnia. "Un nuevo método para el diagnóstico temprano de la enfermedad de Alzheimer basado en el pseudomomento de Zernike obtenido mediante resonancia magnética estructural". Neuroscience 305 (2015): 361–371.

Enlaces externos