Aproximación a los
sistemas ópticos reales
Aproximación a los sistemas ópticos reales
- Las lentes reales no se comportan de modo paraxial.
- Presentan aberraciones.
- Las aberraciones hacen que el sistema ya no se comporte
completamente como un Sistema Lineal e Invariante.
- Se puede abordar el problema reinterpretando algunas de las
ecuaciones vistas.
- En el caso ideal:
[ λf′ ( )]
π
tL (x, y) ∝ exp −j x 2
+ y 2
Aproximación a los sistemas ópticos reales
La teoría paraxial viene de considerar que los ángulos de incidencia y
refracción en un dioptrio son pequeños. Esto permite aproximar las
funciones trigonométricas por el primer término de su desarrollo en
serie de Taylor.
po pi
R n1R(so + R)sin θ = n2 R(si − R)
θ P
S sin θ po pi
so si
n1 n2 Aproximación
sinθ = θ −θ 3
θ 5
θ
+ − +
7
sinθ ≈ θ 3! 5! ⋅⋅⋅7!
Fórmula del n1 + n2 = n2 − n1
fabricante de lentes so si R
Aproximación a los sistemas ópticos reales
En ocasiones (ángulos implicados mayores de 10º) la aproximación no
es válida y deben añadirse más términos a la ecuación
θ3
sin(θ) = θ −
3!
2 2
2s0 ( s0 R ) 2si ( R si )
n1 n2 n2 − n1 n1 1 1 n2 1 1
+ = + h2 + + −
s0 si R
Los términos de corrección se conocen con el nombre de aberraciones
de tercer orden o aberraciones de Seidel.
Por supuesto, pueden cogerse más elementos de la serie e ir afinando
el cálculo.
Hoy en día es posible realizar cálculos numéricos con bastante
precisión, por lo que es posible extender la serie hasta órdenes
arbitrarios.
Aproximación a los sistemas ópticos reales
Se llama aberración a la desviación de la luz del comportamiento
paraxial:
A
Sistema perfecto B
Sistema real
O O’ O O’
La diferencia de camino óptico punto a punto entre el frente de ondas
ideal (paraxial) y el ideal se conoce con el nombre de aberración de
onda:
L = nAB
Aproximación a los sistemas ópticos reales
Recordemos que las diferencias de camino óptico pueden
caracterizarse mediante una función de transmitancia de fase:
[ λ ]
2π
w(x, y) = exp j W(x, y)
Podríamos caracterizar una lente con aberraciones a partir de una
transmitancia en fase genérica.
No obstante, recordemos que los sistemas de imagen con aberraciones
presentan distorsiones en función de la posición del objeto, es decir,
no son sistemas lineales e invariantes.
Opciones:
a) Abandonamos el formalismo de Fourier 👎
b) Tratamos de adaptar el formalismo actual a la nueva situación 👍 👍
Aproximación a los sistemas ópticos reales
Cuando planteamos el modelo de modelo de caja negra:
Aproximación a los sistemas ópticos reales
Las aberraciones se introducen como correcciones a los sistemas
ideales.
Siguiendo esta línea, podemos asumir que los sistemas ópticos son
Sistemas Lineales e Invariantes con correcciones de la función de
entrada.
TODAS las actuaciones sobre el sistema que modifiquen la linealidad
e invariancia del sistema serán consideradas aparte.
Dado que tenemos una caja negra es indiferente en qué punto del
cálculo introducimos las variaciones.
En general, las correcciones debidas a aberraciones ópticas se agrupan
todas en un término que actúa a la entrada del sistema óptico.
Es decir, interpretamos las aberraciones como una función de
transmitancia que modifica la distribución de entrada.
Aproximación a los sistemas ópticos reales
- Recordemos la lente paraxial infinita
[ λf′ ( )]
π
tL (x, y) = exp −j x 2
+ y 2
- La lente paraxial con apertura
[ λf′ ( )]
π
tL (x, y) = A(x) ⋅ exp −j x 2
+ y 2
- A(x) la hemos usado como una función de apertura de la lente
(circular o cuadrada).
- En general, podemos definir una función Pupila Generalizada en la
que se tengan en cuenta no sólo la apertura sino otros efectos como
correcciones al camino óptico (aberraciones)
Aproximación a los sistemas ópticos reales
Por tanto, definimos la función Pupila Generalizada como:
[ λ ]
2π
P (x, y) = A (x, y) exp j W(x, y)
Dicha función no incluye el término de convergencia, de modo que,
en general la transmitancia de una lente se formulará:
[ λ ] [ λf′ ( )]
2π π
tL (x, y) ∝ A(x, y) exp j W(x, y) exp −j x 2
+ y 2
Donde encontramos especificados los términos de amplitud
(apertura/ transparencia), correcciones de camino óptico
(aberraciones) y convergencia.
Aproximación a los sistemas ópticos reales
Manteniendo el formalismo de Fourier:
- La PSF (respuesta impulsional) vendrá caracterizada por la
respuesta del sistema a la pupila generalizada.
- La imagen final será una convolución entre el objeto y la
transformada de Fourier de la función pupila generalizada, que
incluirá:
Términos de apertura
Términos de aberraciones
En las siguientes transparencias se muestran los casos de la PSF de un
sistema óptico (un ojo sano) en función de la apertura, considerando
un sistema ideal (solo efecto de apertura) o real (incluyendo) las
aberraciones.
Aproximación a los sistemas ópticos reales
PSF ojo humano sin aberraciones para diferentes diámetros pupilares:
1 mm 2 mm 3 mm 4 mm
5 mm 6 mm 7 mm
Aproximación a los sistemas ópticos reales
PSF ojo humano con aberraciones para diferentes diámetros
pupilares:
1 mm 2 mm 3 mm 4 mm
5 mm 6 mm 7 mm
Aproximación a los sistemas ópticos reales
Dado que seguimos manteniendo el formalismo de Sistemas Lineales
e Invariantes podemos construir las imágenes mediante convolución
del objeto con la PSF:
Aproximación a los sistemas ópticos reales
En la aberraciones de Seidel se establecía una clasificación en función
de la apariencia de la PSF o del origen de las aberraciones.
Resulta conveniente descomponer la función de Aberración en
términos elementales de modo que sea posible estudiar y cuantificar
simetrías y distorsiones sobre la imagen.
Dado que, en general, las aperturas de las lentes son circulares, se
optó por usar unos polinomios definidos sobre una base circular, los
polinomios de Zernike.
Los polinomios de Zernike describen la superficie de la onda, una vez
descontado el término de focalización. Aunque son los más usados y
existen propuestas normalizadas, tienen problemas para definir
variaciones de alta frecuencia en el frente de ondas.
https://en.wikipedia.org/wiki/Zernike_polynomials
Polinomios de Zernike
En general, los polinomios de Zernike se definen en coordenadas
polares, aunque pueden escribirse en cartesianas si fuera conveniente.
W (ρ, θ) = Cnm Znm (f, θ)
∑
Aberración Coeficiente Polinomio de
de onda de Zernike Zernike
Ejemplo: astigmatismo
Z2 (ρ, θ)
2 2
= ρ cos (2θ)
m: frecuencia angular
n: orden radial
Polinomios de Zernike
Polinomios de Zernike:
Aberraciones de
bajo orden
Aberraciones de
alto orden
Polinomios de Zernike
Estudiemos los diferentes polinomios para entender su significado.
- Los tres primeros polinomios no se consideran aberraciones
propiamente dichas:
El primer polinomio establece una altura de referencia. Es una
constante que establece el plano del “suelo”
El 2 y 3, tampoco son aberraciones propiamente dichas ya que
indican la inclinación del “suelo”. Informan de efectos
prismáticos (dirección general del haz de luz)
En general para el estudio de los sistemas ópticos estos tres se
obvian y no entran el descripción.
Polinomios de Zernike
Polinomios de Zernike (coord. Cartesianas)
Polinomios de Zernike
La función de aberración de ondas se construye mediante la suma
ponderada de los diferentes polinomios.
𝑧! = 1.1
Frente de onda medido
𝑊 𝑥, 𝑦 𝑧" = 0.1 𝑧# = 0.3
+
𝑧$ = 0.25 𝑧& = −0.27 𝑧% = 0.01
= + +
𝑧' = 0.01 𝑧) = 0.01 𝑧( = 0.24 𝑧!* = 0.04
+ + +
+ + + +
𝑊 𝑥, 𝑦 = & 𝑧! ( 𝑍! 𝑥, 𝑦
𝑧!! = −0.01 𝑧!" = 0.01 𝑧!# = 013 𝑧!$ = −0.03 𝑧!& = 0.01
𝑧! (𝜇𝑚)
Aproximación a los sistemas ópticos reales
Veamos los polinomios siguientes:
m
n -5 -4 -3 -2 -1 0 1 2 3 4 5
trefoil coma coma trefoil
quadrafoil 2nd astigmatism spherical 2nd astigmatism quadrafoil
pentafoil 2nd trefoil 2nd coma 2nd coma 2nd trefoil pentafoil
5
Polinomios de Zernike
Consideremos los polinomios de orden radial 2, en la primera fila de
la figura anterior.
- (2,0) Desenfoque: Recordemos que en la definición de
transmitancia se separaba el efecto de convergencia del de las
aberraciones. Este polinomio indica que el haz no focaliza en el
punto teórico. Un haz convergente está definido por un factor de
fase esférico. Este polinomio añade o elimina un factor esférico a la
esfera original.
Polinomios de Zernike
Consideremos los polinomios de orden radial 2, en la primera fila de
la figura anterior.
- (2,2) y (2,-2) Astigmatismos: En un haz de luz con astigmatismo ,
en cualquier plano de propagación, se distingue una dirección en la
que el haz tiene más curvatura y una dirección perpendicular con
menos curvatura. Si por ejemplo el eje más potente estuviera
alineado con el eje x y menos potente con el eje y, matemáticamente
es como si las componentes x tuviesen una focal y las y otra. El
efecto es el haz converge antes en una dirección que en otra.
En la figura posterior se ilustra la propagación de un haz
astigmático. Vemos que tiene forma ligeramente elíptica de modo
que se cierra (focaliza) antes en una dirección que en otra. Las
focales direccionales se conocen como “focales de Sturm”
Polinomios de Zernike
Consideremos los polinomios de orden radial 2, en la primera fila de
la figura anterior.
- (2,2) y (2,-2) Astigmatismos:
Polinomios de Zernike
Uno de los campos en los que el estudio de aberraciones es más
importante es en la instrumentación óptica, astronomía y en
oftalmología.
Las tres aberraciones anteriormente descritas (de bajo orden) son las
más comunes y fáciles de solucionar. En óptica oftálmica, son las que
corrigen las lentes normales.
Las aberraciones de alto orden son más delicadas, ya que pueden
introducir asimetrías en la PSF produciendo distorsiones importantes.
En una transparencia anterior se mostraba la PSF de un ojo, donde se
aprecia el efecto de las aberraciones de alto orden.
Aunque, en general, las aberraciones tienden a minimizarse para
mejorar la calidad de un sistema óptico, en ocasiones se diseñan
elementos con aberraciones determinadas con el fin de producir un
efecto determinado.
Polinomios de Zernike
En la animación se muestra la propagación de un haz de luz a partir
de una córnea con una deformación corneal. Se aprecia la existencia
de astigmatismo y otras aberraciones de alto orden (coma)
Polinomios de Zernike
Para analizar las aberraciones se utilizan dos estrategias
fundamentales:
- Análisis directo del frente de onda (amplitud y fase)
Estas técnicas usan la información del frente de ondas para obtener
toda la información posible para luego descomponer el frente en una
suma de polinomios de Zernike.
Se utilizan métodos de aproximación iterativos: se introducen como
parámetros los patrones de luz de entrada y de salida y se va
modificando la transmitancia de acuerdo a ciertos parámetros hasta
que el sistema converge. El algoritmo más utilizado para ello es el de
Gerchberg-Saxton.
Este método se usa, principalmente en el diseño de elementos ópticos
difractivos y holográficos.
Polinomios de Zernike
- Análisis directo del frente de onda (amplitud y fase)
Otra técnica consiste en utilizar loops cerrados con óptica adaptativa.
Los sistemas de óptica adaptativa son sistemas (espejos o sistemas
ópticos) que actúan localmente sobre el camino óptico, introduciendo
más o menos retraso en la onda.
El objetivo es “construir” un negativo del frente de onda para
compensar las aberraciones. A partir del negativo se puede obtener la
información del frente de ondas.
Este método es ampliamente utilizado tanto en óptica oftálmica como
en astronomía.
Polinomios de Zernike
- Análisis directo del frente de onda (amplitud y fase)
Polinomios de Zernike
- Análisis directo del frente de onda (amplitud y fase)
Polinomios de Zernike
- Análisis directo del frente de onda (amplitud y fase)
Polinomios de Zernike
- Análisis directo del frente de onda: a partir del análisis de la
amplitud y de la fase.
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
Estos dispositivos miden directamente el gradiente de la función de
aberración de onda de forma local.
Localmente, una porción del frente de ondas puede considerarse
como un frente plano con un rayo asociado perpendicular a dicho
frente local.
Con una matriz de microlentes se focaliza una porción muy pequeña
del frente de ondas. A esa escala, el frente será cuasiplano y
focalizarán todos los frentes en el mismo plano. La posición lateral
del punto focal está relacionada con la inclinación local del haz, y este
con el gradiente.
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
Matriz de microlentes
Dy(x1, y1)
Dy(x1, y2)
Frente de ondas
¶W ( x, y ) Dx( x, y )
aberrado =-
¶x f
Dy(x1, y3)
¶W ( x, y ) Dy ( x, y )
=-
¶y f
Dy(x1, y4)
Longitud focal f’
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
Plano imagen
𝛿𝑦!
tan 𝛼! = −
𝑓"#
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
El frente de ondas es una función bidimensional
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
El frente de ondas es una función bidimensional
https://www.researchgate.net/publication/262951605_Real-
time_correlation_reference_update_for_astronomical_adaptive_optics/figures?lo=1
Polinomios de Zernike
- Uso de aberrómetros (Hartmann-Shack)
Este método es ampliamente utilizado tanto en astronomía como en
óptica oftálmica e instrumental.
Aunque proporciona un análisis rápido de las aberraciones del frente
de onda, desprecia mucha información, ya que sólo toma en
consideración la posición del spot local, y no su calidad. Además,
puede haber problemas por cruce de spots.
En aplicaciones reales se suelen usar combinaciones de diferentes
técnicas, tratando de maximizar las ventajas de cada una y minimizar
las desventajas.
Métricas de calidad de imagen
Métricas de calidad de imagen
Hemos visto que en los sistemas formadores de imágenes aparecen
varios factores limitantes que degradan la calidad final de la imagen.
Se hace necesario definir diferentes métricas que evalúen la calidad de
la imagen a partir de parámetros objetivos.
Antes de definir las métricas más utilizadas, recordaremos algunos
conceptos de interés:
Respuesta impulsional, h(x): Es la función respuesta a un impulso
puntual (delta) proporcionado por un Sistema Lineal e Invariante a
Traslaciones.
h (x2 − u, y2 − v) = S {δ (x1 − u, y1 − v)}
+∞
∬−∞
g2 = g1 ⊗ h = g1 (u, v)h (x − u, y − v) dudv
Métricas de calidad de imagen
Point Spread Function (PSF) o Función de dispersión de punto. Es el
nombre que recibe la respuesta impulsional en óptica. En general, el
término PSF hace referencia a la respuesta impulsional de un sistema
formador de imágenes.
El estudio de su forma y extensión proporciona información del
sistema óptico.
La PSF de referencia es la función de Airy, ya que en general, los
sistemas ópticos usan aperturas circulares
https://en.wikipedia.org/wiki/Point_spread_function
Métricas de calidad de imagen
Point Spread Function (PSF) con Aberración Esférica
Arriba: AE negativa; Centro: no AE; Abajo AE positiva
https://en.wikipedia.org/wiki/Point_spread_function
Métricas de calidad de imagen
Point Spread Function (PSF)
https://hubblesite.org/
Métricas de calidad de imagen
https://www.nasa.gov/content/hubbles-mirror-flaw
Aplicaciones: calidad de imagen
Función de transferencia, H. Es la transformada de Fourier de la
respuesta impulsional.
Proporciona información acerca del contenido frecuencial que el
sistema transfiere, es decir, la ventana de transmisión y modulación
de las diferentes frecuencias.
ℱ [h (x, y)] = H (ξ, η)
Dado que en los sistemas ópticos clásicos no se registra la fase, en
ocasiones se define la Función de Transferencia en Amplitud, como la
Transformada de Fourier de la amplitud de la respuesta impulsional.
Métricas de calidad de imagen
Función de transferencia óptica (OTF)
Dado que, en general trabajamos con intensidades, se define la
función de transferencia óptica, ℋ como la TF de la intensidad de la
PSF.
ℱ(Io) = ℱ(Ii) × ℋ
En general, cuando hablamos de transferencias de intensidad,
entendemos que la transferencia total se produce con valor 1. Por ello,
el valor de la OTF suele darse normalizado:
ℱ [ | h(x, y) | ]
2
ℋ(ξ, η) =
∬−∞ | h (x, y) |2 dxdy
+∞
Métricas de calidad de imagen
Función de transferencia de modulación (MTF): es el módulo de la
función de transferencia óptica
. Dado que, en general trabajamos con intensidades y no con
variables complejas, definimos la MTF a partir de la OTF normalizada
como
MTF = | ℋ(ξ, η) |
De manera equivalente, se define la función de transferencia en fase
(PTF) como la fase de la OTF:
OTF (u , v ) = OTF × e = MTF (u , v ) × e
jf jPTF (u , v )
Métricas de calidad de imagen
La Función de Transferencia de Modulación no sólo informa acerca de las
frecuencias que pasan o no pasan a través del sistema, sino de su atenuación en
intensidad, tal y cómo se muestra en la figura:
bajas medias altas
objecto:
contraste
100%
image
1
contraste
0
frecuencia espacial
Métricas de calidad de imagen
La fase de la OTF indica la fase relativa de una frecuencia dada. Si consideramos una
componente aislada de Fourier (una red sinusoidal), la fase es su posición relativa
respecto al origen.
objecto
imagen
180
phase shift
-180
frecuencia espacial
Métricas de calidad de imagen
La Función de Transferencia de Modulación no sólo informa acerca
de las frecuencias que pasan o no pasan a través del sistema, sino de
su atenuación en intensidad.
En la imagen a continuación se muestra la MTF de un sistema óptico,
tanto en intensidad (centro) como un perfil radial (izda.). A la derecha
se muestra la imagen de un test de resolución radial
La fila superior muestra el caso para un sistema perfectamente
enfocado, mientras que en el caso inferior muestra un caso fuera de
foco.
Nótese que la frecuencia de corte de la MTF es la misma en los dos
casos, pero en en el caso inferior se produce una atenuación
importante de las frecuencias a partir de 200 ciclos/grado
Métricas de calidad de imagen
https://en.wikipedia.org/wiki/Optical_transfer_function
Métricas de calidad de imagen
En la siguiente transparencia se muestra la diferencia entre la OTF y
la MTF para la misma imagen.
Fijémonos que en la OTF hay algunos tramos de frecuencia negativos
que, obviamente, en la MTF no se aprecia al tomar valores absolutos.
Una inversión del signo en la OTF significa que la fase se ha invertido
(la fase de la onda ha cambiado una cantidad> 𝜋). Si representamos
visualmente una onda de contraste, esto significa que las bandas se
han desplazado medio período y hay una inversión de contraste.
Podemos apreciar la inversión de contraste en las zonas marcadas en
el tests de resolución.
Esto sifnifica que, pese a que la sencillez de manejo de la MTF es a
costa de perder información.
Métricas de calidad de imagen
https://en.wikipedia.org/wiki/Optical_transfer_function
Métricas de calidad de imagen
La MTF de referencia es la que se obtiene para una abertura circular
básica.
Recordemos que en un sistema formador de imágenes, la respuesta
impulsional venía dada por:
Métricas de calidad de imagen
La MTF de referencia es la que se obtiene para una abertura circular
básica.
é æ u ö u æ u ö ú
2ù
2ê Rapertura
OTF (u ) = arccosçç ÷÷ - 1 - çç ÷÷ uc =
pê è 2uc ø 2uc è 2uc ø ú s′
ë û
u / uc
*Introduction to Fourier Optics 2 ed. (Joseph W. Goodman) pp. 143
Métricas de calidad de imagen
La curva que marca la OTF de una apertura es el caso óptimo. Cualquier
aberración en el sistema provocará una bajada de la curva.
Desde el punto de vista teórico, las aberraciones no recortan la frecuencia límite,
pero desde el punto de vista práctico, puede considerase que cuando la OTF
toma valores muy bajos la frecuencia correspondiente no se transmite por el
sistema
Métricas de calidad de imagen
Métricas de calidad de imagen
Razón de Strehl: Proporciona un valor númérico para la calidad de
un sistema óptico afectado por aberraciones.
Aunque se puede definir de muchos modos, la más sencilla es el
cociente el máximo de la PSF de un sistema con aberraciones y el
máximo de la PSF sólo limitada por difracción. Equivalentemente, se
define como el cociente entre el volumen de la OTF real de un sistema
y el volumen de la OTF del sistema en ausencia de aberraciones (sólo
limitado por difracción):
+∞
PSFA(0,0) ∬−∞ ℋA(ξ, η)dξdη
SR = ≡ +∞
PSFD(0,0) ∬−∞ ℋD(ξ, η)dξdη
Su utilidad reside en que condensa en un único número la calidad de
un sistema óptico, aunque la información que se puede deducir de
este parámetro es muy limitada
Métricas de calidad de imagen
https://www.olympus-ims.com/en/microscope/terms/strehl_ratio/ Aplanático=libre de aberraciones
