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 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]

Definiciones

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

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

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

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

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

Ortogonalidad

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

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

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

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

Polinomios radiales

Los primeros pocos polinomios radiales son:

Polinomios de Zernike

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

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

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

Referencias

  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

 [math.NA]. 

Enlaces externos

Read other articles:

画像提供依頼:大崎公園の画像提供をお願いします。(2021年10月) おおさきかみじまちょう 大崎上島町 海と島の歴史資料館 大崎上島町旗 大崎上島町章 国 日本地方 中国地方、山陽地方中国・四国地方都道府県 広島県郡 豊田郡市町村コード 34431-1法人番号 8000020344311 面積 43.11km2総人口 6,703人 [編集](推計人口、2023年10月1日)人口密度 155人/km2隣接自治体 (海上で隣

Fabriekspand Cuypers-Stoltzenberg (nu Stedelijk Museum Roermond) Het atelier Cuypers-Stoltzenberg was een Nederlandse onderneming in Roermond (1852-1947), gespecialiseerd in het maken van beeldhouwkunst en meubelstukken voor rooms-katholieke kerken. Geschiedenis Rond 1838 richtte textielfabrikant François Stoltzenberg (1805-1875) in Roermond het atelier F. Stoltzenberg op, waar goudborduurwerk werd gemaakt voor vaandels en kerkelijke gewaden. Na 1850 werd het atelier Stoltzenberg uitgebreid ...

عبد الباسط عبد العليم حجاب معلومات شخصية تاريخ الميلاد سنة 1906  تاريخ الوفاة 1988 الحياة العملية المهنة كاتب  تعديل مصدري - تعديل   يفتقر محتوى هذه المقالة إلى الاستشهاد بمصادر. فضلاً، ساهم في تطوير هذه المقالة من خلال إضافة مصادر موثوق بها. أي معلومات غير موثقة يمكن ال...

2000 television film by Roger Young This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed.Find sources: The Thin Blue Lie – news · newspapers · books · scholar · JSTOR (October 2023) (Learn how and when to remove this template message) The Thin Blue LieWritten byDaniel HelfgottDirected byRoger YoungStarringRob MorrowRandy Quaid...

Artikel ini sebatang kara, artinya tidak ada artikel lain yang memiliki pranala balik ke halaman ini.Bantulah menambah pranala ke artikel ini dari artikel yang berhubungan atau coba peralatan pencari pranala.Tag ini diberikan pada Oktober 2022. Artikel ini bukan mengenai Mockbuster. Film kembar adalah film dengan alur yang sama atau sangat mirip yang diproduksi atau dirilis pada saat yang sama oleh dua studio film berbeda.[1] Fenomena tersebut dapat timbul dari dua perusahaan produksi...

Sebuah cluster Xserve kecil dengan Xserve RAID. Apple Computer mengenalkan Xserve pertama kalinya di bulan Mei 2002. Produk ini adalah sebuah komputer kelas server berukuran 1U. Xserve dijual dengan konfigurasi satu atau dua prosesor PowerPC G4 berkecepatan 1.33 GHz. Server ini bisa menyimpan sampai empat hard-drive, sebuah optical drive dan juga dilengkapi dengan video card buatan ATI. Di bulan Maret 2003 Apple merilis Xserve Cluster node yang hanya bisa menyimpan satu hard-drive dan ti...

Ця стаття має кілька недоліків. Будь ласка, допоможіть удосконалити її або обговоріть ці проблеми на сторінці обговорення. Ця стаття містить текст, що не відповідає енциклопедичному стилю. Будь ласка, допоможіть удосконалити цю статтю, погодивши стиль викладу зі стилісти

この記事は検証可能な参考文献や出典が全く示されていないか、不十分です。出典を追加して記事の信頼性向上にご協力ください。(このテンプレートの使い方)出典検索?: ブレーメンのマルクト広場の市庁舎とローラント像 – ニュース · 書籍 · スカラー · CiNii · J-STAGE · NDL · dlib.jp · ジャパンサーチ · TWL(2010年10月) ブレ...

Calendar year Millennium: 2nd millennium Centuries: 17th century 18th century 19th century Decades: 1730s 1740s 1750s 1760s 1770s Years: 1747 1748 1749 1750 1751 1752 1753 1750 by topic Arts and science Archaeology Architecture Art Literature Poetry Music Science Countries Austria Canada Denmark France Great Britain Ireland Norway Russia Scotland Spain Sweden Lists of leaders State leaders Colonial governors Religious leaders Birth and death categories Births Deaths Establ...

International trips made by Narendra Modi while Prime Minister This article is part of a series aboutNarendra Modi Prime Minister of IndiaIncumbent Electoral history Public image Awards and honours Bibliography Chief Minister of Gujarat 2002 2007 2012 Gujarat Council of Ministers First Second Third Fourth Premiership 2014 campaign Achhe Din Aane Waale Hain 2019 campaign Main Bhi Chowkidar Swearing-in 2014 2019 Union Council of Ministers First Second Lok Sabha Sixteenth Seventeenth Timeline Ma...

For another convent with the name, see Ascension Convent (Kyiv). Katholikon of Ascension Convent (1580s), from an early 19th-century drawing. Ascension Convent, known as the Starodevichy Convent or Old Maidens' Convent until 1817 (Russian: Вознесенский монастырь, Voznesensky monastyr'), was an Orthodox nunnery in the Moscow Kremlin which contained the burials of grand princesses, tsarinas, and other noble ladies from the Muscovite royal court. It was destroyed in 1929...

Período Dinástico Arcaico Localización geográfica aproximada DatosCronología Edad del BronceLocalización Mesopotamia[editar datos en Wikidata] Período Dinástico Arcaico, período protodinástico de Mesopotamia o período presargónico fue la etapa de la historia de Mesopotamia comprendida entre 2900 a. c. y 2334 a. c.: esto es, en el período sumerio. Está marcada por el desplazamiento de las poblaciones rurales hacia las ciudades; lo cual se traducirá ...

S. JanakiInformasi latar belakangNama lainJanakiammaLahir23 April 1938 (umur 85)Repalle, Guntur, Andhra Pradesh, IndiaGenreNyanyian playback, Klasik IndiaPekerjaanPenyanyiInstrumenVokalTahun aktif1957–sekarang Sishtla Sreeramamurthy Janaki (ಕನ್ನಡ: ಎಸ್. ಜಾನಕಿ) (kelahiran 23 April 1938), yang umumnya dikenal sebagai S. Janaki, adalah seorang penyanyi playback India. Hingga 2007[update] Janaki telah menyanyikan sekitar 40,000 lagu, kebanyakan dalam bahasa...

هذه المقالة يتيمة إذ تصل إليها مقالات أخرى قليلة جدًا. فضلًا، ساعد بإضافة وصلة إليها في مقالات متعلقة بها. (أبريل 2019) ستانلي جورج براون معلومات شخصية الميلاد 8 ديسمبر 1907  تاريخ الوفاة 20 يناير 1986 (78 سنة)   مواطنة نيجيريا  الحياة العملية المهنة جراح  اللغات الإنجليزية...

UNESCO World Heritage Site in Catalonia, Spain Eglésia de Santa Maria de l'AssumpcióUNESCO World Heritage SiteLocationCóll, Vall de Boí, Alta Ribagorça, Province of Lleida, Catalonia, SpainPart ofCatalan Romanesque Churches of the Vall de BoíCriteriaCultural: (ii), (iv)Reference988-005Inscription2000 (24th Session)Area0.46 ha (1.1 acres)Coordinates42°28′20″N 0°46′22″E / 42.47222°N 0.77278°E / 42.47222; 0.77278Location of Santa Maria, Cóll i...

Type of roller coaster Stand-up roller coasterShockwave, a former stand-up coaster that featured a zero-g roll inversion.StatusIn ProductionFirst manufactured1984No. of installations21ManufacturersTOGOArrow DynamicsBolliger & MabillardIntaminFirst retrofit1982 A stand-up roller coaster is a roller coaster where passengers aboard a train stand throughout the course of the ride. The first manufacturer to employ the format was TOGO, a Japanese company that converted two traditional roller co...

Day to commemorate war dead Not to be confused with Remembrance Day. This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed.Find sources: Remembrance Sunday – news · newspapers · books · scholar · JSTOR (November 2021) (Learn how and when to remove this template message) Remembrance SundayThe poppy is worn around the time of Rem...

Kadabanahalli—  Làng  — Quốc gia Ấn ĐộBangKarnatakaHuyệnHassanThành phố/thị trấnSakleshpurMúi giờIST (UTC+5:30) Kadabanahalli là một làng thuộc tehsil Sakleshpur, huyện Hassan, bang Karnataka, Ấn Độ.[1][2] Tham khảo ^ “Danh sách làng thuộc bang Karnataka”. Bộ Nhà ở Ấn Độ. Truy cập ngày 27 tháng 10 năm 2012. ^ “Làng Kadabanahalli”. Onefivenine.com. Truy cập ngày 27 tháng...

Albanian game showThis article's factual accuracy is disputed. Relevant discussion may be found on the talk page. Please help to ensure that disputed statements are reliably sourced. (September 2023) (Learn how and when to remove this template message) Kush do të bëhet milioner?GenreQuiz showPresented byVeton Ibrahimi (2008-2009) Fehmi Ferati (2009-2010) Agron Llakaj (2011-2012) Enkel Demi (2014)ComposersKeith Strachan Matthew StrachanCountry of originKosovo AlbaniaOriginal languageAlbanian...

 TV series or program AnimaliaCreated byGraeme Base (book)Based onAnimaliaby illustrator Graeme BaseDeveloped by Graeme Base Tom Ruegger Ewan Burnett Doug MacLeod Robyn Base Bruce D. Johnson Written by Sherri Stoner John P. McCann Deanna Oliver Tom Ruegger Nicholas Hollander Mark Seidenberg John K. Ludin John Loy Directed byDavid ScottVoices of Brooke Anderson Katie Leigh Robert Mark Klein Kate Higgins Dean O'Gorman Chris Hobbs Peta Johnson Theme music composer Graeme Base Yuri Worontsch...