Polinomis de Zernike

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]

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

Definicions modifica

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

 

i els imparells com:

 

on m i n són nombres enters no negatius amb n ≥ m, φ és l'angle azimutal, ρ és la distància radial   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,  . Els polinomis radials Rmn es defineixen com:

 

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

Altres representacions modifica

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

. 

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.

 

per a n-m parell.

El factor   en el polinomi radial   es pot expandir en una base de Bernstein de   per a   parell o   multiplicat per una funció de   per a   imparell en el rang  . Per tant, el polinomi radial pot expressar-se mitjançant un nombre finit de polinomis de Bernstein amb quocients racionals:

 

Índexs seqüencials de Noll modifica

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ó   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,  ) 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 modifica

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 modifica

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 modifica

Ortogonalitat modifica

L'ortogonalitat en la part radial s'expressa com

 

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

 
 
 

on   (de vegades anomenat factor de Neumann perquè apareix amb freqüència juntament amb les funcions de Bessel) es defineix com 2 si   i com 1 si  . 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,

 

on   és el jacobià del sistema de coordenades circulars, i on   i   són parells.

Un valor especial és

 

Transformada de Zernike modifica

Qualsevol camp de fase de valor real prou uniforme sobre el disc de la unitat   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:

 

els coeficients es poden calcular usant productes interns. A l'espai de les funcions del   disc de la unitat, existeix un producte intern definit per

 

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

 

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 modifica

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

 

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

 

on   també podria escriure's   perquè   é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:

 

La periodicitat de les funcions trigonomètriques implica constància si és trencada per múltiples de   radiants al voltant del centre:

 

Relacions de recurrència modifica

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]

 

De la definició de   es pot veure que   i  . La següent relació de recurrència de tres termes permet calcular tots els altres  :[7]

 

La relació anterior és especialment útil, ja que la derivada de   es pot calcular a partir de dos polinomis de Zernike radials de grau adjacent:[7]

 

Exemples modifica

Polinomis radials modifica

Els primers pocs polinomis radials són:

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Polinomis de Zernike modifica

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

 
  Índex

OSA/ANSI

( )

Índex

Noll

( )

Grau radial

( )

Grau azimutal

( )

  Nom clàssic
   0  1 0  0   Pistó (vegeu distribució semicircular de Wigner)
   1  3 1 −1   Inclinació (inclinació I, inclinació vertical)
   2  2 1 +1   Inclinació horitzontal
   3  5 2 −2   Astigmatisme oblic
   4  4 2  0   Desenfoc (posició longitudinal)
   5  6 2 +2   Astigmatisme vertical
   6  9 3 −3   Lobulat vertical
   7  7 3 −1   Coma vertical
   8  8 3 +1   Coma horitzontal
   9 10 3 +3   Lobulat oblic
  10 15 4 −4   Cuatrilobulat oblic
  11 13 4 −2   Astigmatisme secundari oblic
  12 11 4  0   Esfèrica primària
  13 12 4 +2   Astigmatisme vertical secundari
  14 14 4 +4   Cuatrilobulat vertical

Aplicacions modifica

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  , 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 modifica

El concepte es tradueix a dimensions majors D si els multinomis   en coordenades cartesianes es converteixen en coordenades hiperesfèriques,  , multiplicades per un producte de polinomis de Jacobi de les variables angulars. Per exemple, en la dimensió  , les variables angulars són harmònics esfèrics. Combinacions lineals de les potències   defineixen una base ortogonal   que satisfà

. 

(Tingui's en compte que un factor   s'absorbeix aquí en la definició de R, mentre que en   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

 

fins i tot per  , o en cas contrari, idèntic a zero.

Vegeu també modifica

Referències modifica

  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 modifica

math.NA.