Curso Fourier Uba
Curso Fourier Uba
Curso Fourier Uba
Claudio Iemmi
OPTICA DE FOURIER
INDICE
i
Óptica de Fourier Claudio Iemmi
ii
Óptica de Fourier Claudio Iemmi
VIII. HOLOGRAFÍA
Introducción 207
Proceso de síntesis y reconstrucción de un holograma 209
Disposición de Gabor 213
Disposición de Leith – Upatnieks 215
Hologramas bi y tri dimensionales 217
Hologramas planos 218
Holograma de Fresnel 219
Holograma de Fraunhofer 224
Holograma de Fourier 224
Hologramas rainbow 228
Hologramas de fase 231
Eficiencia de hologramas planos 235
Factores que afectan la resolución de una imagen holográfica 236
Hologramas de volumen 241
Franjas en 3D 242
Difracción por una red 3D 244
Sensibilidad a la orientación y a la longitud de onda de las redes 3D 252
Hologramas de volumen de transmisión y reflexión 257
iii
Óptica de Fourier Claudio Iemmi
X. DIGITALIZACIÓN DE IMÁGENES
Digitalización 285
Teorema del muestreo 286
Transformada de Fourier discreta de señales digitales 289
Hologramas generados por computadora (CGH) 293
Holografía digital 297
iv
Óptica de Fourier C. Iemmi
1
Óptica de Fourier C. Iemmi
Fue Kirchhoff en 1882 quién mediante fundamentos matemáticos logra demostrar que
estas hipótesis son resultados que surgen naturalmente de la aplicación de la teoría
ondulatoria de la luz, es decir a partir de las ecuaciones de Maxwell (1864). Sin embargo
Kirchhoff también hace dos suposiciones que se demostrarán incompatibles.
Sommerfeld en 1894 modifica lo planteado por Kirchhoff y logra así evitar esta
incompatibilidad. Surge así la teoría de difracción llamada de Rayleigh – Sommerfeld.
Si bien tanto la teoría de Kirchhoff como la de Rayleigh – Sommerfeld se mostraron muy
precisas en muchos casos, cabe destacar que no funcionan en otros ya que realizan
algunas aproximaciones y simplificaciones de importancia. La más importante de ellas es
que la luz es tratada como un fenómeno escalar. Es decir se supone que el campo
eléctrico (en este caso interesa el campo eléctrico y no el magnético ya que es el vector
óptico) tiene una sola componente ó, si tuviese dos, la otra puede ser tratada de igual
forma. Esta aproximación desprecia el hecho de que las distintas componentes del campo
eléctrico están acopladas por las ecuaciones de Maxwell y no pueden tratarse de forma
independiente. Afortunadamente esta teoría da resultados muy precisos si se cumple que:
Está claro que estas condiciones no se cumplen en problemas de gran importancia como
por ejemplo difracción por nano-estructuras, difracción cónica,etc. En estos casos debe
utilizarse una teoría vectorial rigurosa para obtener resultados precisos. Estos problemas
quedan excluidos de este curso. Esta teoría también falla en sistemas formadores de
imágenes si se trabaja con ángulos muy grandes. Analicemos por ejemplo la Figura 1.
Sobre el lado izquierdo de la lente no hay interferencia ya que los campos E1 y E2 son
2
Óptica de Fourier C. Iemmi
90º
E1
Figura 1: Sobre el lado izquierdo de la lente los campos son perpendiculares y por lo tanto no
interfieren, en cambio sí lo hacen sobre el lado derecho
1 B
X E Ley de Faraday
c t
4 1 D
X H J Ley de Ampere
c c t
J E B H D E
3
Óptica de Fourier C. Iemmi
1
. D . E ; X H X B
1 B 4 E 4
X E ; X B J ; . E ; . B 0
c t c c t
Para llegar a la ecuación de ondas podemos, por ejemplo, tomar rotor en ambos
miembros de la primera ecuación, tenemos así que
1 XB
X X E . E 2 E
c t
4 E
J
c t
. E E
2 1 c
c t
2 E 4 J
2 E 4
c 2 2t c 2 t
4
Óptica de Fourier C. Iemmi
fuentes S = 0 y recordando que nosotros vamos a trabajar con campos escalares una
solución posible de la ecuación es:
E r , t E0 e
i k r t
Si E0 Cte esta expresión corresponderá a una onda plana, si E0 Cte será una onda
r
esférica y si E0 Cte describirá a una onda cilíndrica.
r
Dado que la ecuación de ondas es lineal y de coeficientes reales, la solución debe ser
real, con lo cual la expresión que tenga sentido físico será E r , t e E0 e
i k r t
.
2
E
1 2E
0 con
E r ,t E r ei t
v 2
t
2
k E 0
2 2
(1)
términos del valor del campo E y de su derivada normal, sobre una superficie cerrada
arbitraria S que rodea a dicho punto. (Esta solución había sido propuesta con
anterioridad por Helmholtz para acústica).
Para encontrar esta solución utiliza el Teorema de Green que establece que si U1 y U 2
5
Óptica de Fourier C. Iemmi
U1 U 2 U 2 U1 dV U1 U 2 U 2 U1 n dS
2 2
V S
P1 S
r10
P0
V
2 U k 2 U 0 U 2 U U k 2 U 0
1 1 2 1 2 1
2 U k 2 U 0 U 2 U U k 2 U 0
2 2 1 2 1 2
U1 U 2 U 2 U1 n dS 0 (2)
S
ikr10
Podemos elegir U1 E y U G P , P e con
2 1 0 r10
1
r10 x0 x1 y0 y1 z0 z1
2 2 2 2
6
Óptica de Fourier C. Iemmi
debemos tener presente que esta función no representa una onda real sino que es una
herramienta matemática para resolver el problema.
Ahora bien, si observamos el dibujo vemos que si queremos llevar P1 hasta P0 aparece
una divergencia, entonces lo que se debe hacer es rodear a P0 con una esfera de radio
P1 S’
r10
P0
V V
S
V’
Entonces
S
lim
S'
0
S
. Vamos a calcular la segunda integral teniendo en
cuenta que el gradiente se toma sobre la variable 1 de modo que 1 , ,
1x y1 z1
eikr10 1 eikr10
1G 1 ik rˆ10
r10 r10 r10
1 eikr10 eikr10
S r10 r10
E ik ˆ
r10
r
E n dS
10
U
Debemos recordar que n U con lo cual la ecuación anterior queda
n
7
Óptica de Fourier C. Iemmi
eikr10 1 E
E ik dS pero sobre S r10
S
r10 r10 n
eik 1 E
S
E ik
n
dS dado que 0 y recordando que E y G y
0 0
(3)
Veamos ahora como aplicar el teorema integral al caso de difracción por una pantalla
plana, infinita y opaca, con una abertura.
8
Óptica de Fourier C. Iemmi
SR
SP
P0
r10
P1 R
fuente
real
Ahora S = SP + + SR
S SP SR
de tomar n hacia adentro (al revés que en el Goodman) es para que coincida con los z
positivos sobre la abertura. Por otra parte se puede ver también la razón de tomar ondas
emergentes de P1 en vez de P0 ya que ahora sobre esas onditas coincidirán con las
fuentes virtuales secundarias del Principio de Huygens – Fresnel.
Vamos a calcular primero la integral (3) en la superficie SR.
2 1 1
De acuerdo a las hipótesis iniciales R k ; además R se
R R
hace tender a . Así la integral queda:
9
Óptica de Fourier C. Iemmi
1 exp ikR E
SR
4
SR
R
n
ikE n ˆ
r10 dS R
Ahora bien, uno estaría tentado a decir que dado que G y E decaen cuando
Otro argumento sería que, dado que la radiación se propaga a la velocidad finita c, basta
tomar SR tan grande que las ondas aún no hayan llegado; pero existe el inconveniente
que trabajamos con ondas monocromáticas y por lo tanto estas ondas, por definición,
existieron siempre.
Obviamente en la práctica esto siempre es válido ya que en una experiencia siempre hay
alguien que enciende la fuente, es decir estrictamente hablando, aunque usemos un láser
las ondas son cuasi-monocromáticas.
E
lim R ikE n rˆ10 0
R
n
Esta es la condición de radiación de Sommerfeld y es válida para campos que se atenúen
Ejercicio 1:
Probar que si la fuente real emite ondas esféricas, se verifica la condición de radiación
Dado que la perturbación que ilumina la apertura es una onda esférica, o una
combinación lineal de ondas esféricas, esta condición siempre se cumple. Podemos pasar
ahora a resolver las integrales sobre la pantalla.
10
Óptica de Fourier C. Iemmi
E
Sobre la superficie Sp el campo E y son idénticamente cero.
n
Esta aproximación tampoco es estrictamente cierta ya que siempre hay una zona de
penumbra cercana a la abertura. Por otra parte este requisito lleva a una inconsistencia,
tal como lo veremos más adelante. Estas condiciones a menudo son conocidas como
condiciones de contorno de Kirchhoff.
Así las integrales sobre y Sp se reducen a la integral sobre
1 exp ikr10 E 1
E P0 r10 n r10 10 dS
ik E n ˆ
r
4
1
Como antes r10 k
r10
exp ikr10 E
E P0
1
4
r10 n ik E cos n rˆ10 dS (4)
1 e 21
ikr
eikr21 1
E P A 1 E A ik r 21 nuevamente k
1 r21 r r21 r21
21
11
Óptica de Fourier C. Iemmi
P0
r10
P2 r21 P1
fuente
real
Figura 5: Difracción por una pantalla plana iluminada con una onda esférica
(5)
exp ikr10
E P0 E ' P1 dS con
r10
E ' P1
A exp ikr21 cos n, rˆ21 cos n, r10
i r21 2
12
Óptica de Fourier C. Iemmi
Vemos que podemos interpretar el campo en P0 como el provisto por infinitas fuentes
puntuales virtuales, ubicadas en los distintos P1 de la abertura. Es decir llegamos a la
formulación que intuitivamente habían hecho Huygens y Fresnel. Es más, vemos que la
eikr21
onda primaria que llega a P1 tiene amplitud y fase A pero la onda emergente está
r21
desfasada con un adelanto de 90º y además está atenuada por el factor de oblicuidad
[cos, cos], que vale 1 como máximo, y por la cantidad 1 . Es decir que surgen de forma
natural todos los requerimientos que Fresnel había impuesto en forma arbitraria para
ajustar sus mediciones.
Obviamente estas fuentes virtuales no tienen un real significado físico, cabe entonces
preguntarse como surge la alteración del campo al pasar por la abertura. Young (1802)
sugirió que el campo observado correspondía a la interacción entre el campo de la onda
incidente, directamente transmitido, con la onda difractada en los bordes de la abertura.
Posteriores investigaciones confirmaron este punto de vista.
CLASE 2
Si bien los resultados hallados por Kirchhoff coinciden muy bien con las experiencias, se
presenta una inconsistencia en una de las hipótesis que utiliza. La misma está originada
en el hecho de pedir condiciones simultáneas sobre el campo y sus derivadas.
Existe un teorema que dice que si una función potencial y su derivada normal se anulan
simultáneamente a lo largo de un segmento de curva, entonces dicha función debe
anularse en todo el plano (ver por ejemplo Classical Electrodynamics de Jackson).
Evidentemente esto no es lo que sucede, por lo tanto alguna hipótesis está mal planteada.
Tales inconsistencias fueron solucionadas por Sommerfeld quien eliminó la necesidad de
E
imponer condiciones sobre E y . Para ello eligió funciones de Green que fuesen
n
G
combinaciones lineales de la anteriormente utilizada y tales que, o bien G ó se
n
13
Óptica de Fourier C. Iemmi
E
anulasen sobre SP+ . De esta forma sólo se debe pedir condiciones sobre E ó en
n
forma separada.
Cabe destacar que el hecho que una teoría sea autoconsistente y la otra no, no significa
que una sea más precisa que la otra.
Sommerfeld plantea como soluciones de G
r’10 r10
P’0 P0
1G
exp ik r10 1
ik rˆ10
exp ik r10' ik 1 r '
10
r10 r10 r10' r10'
14
Óptica de Fourier C. Iemmi
E G
1
E P0 G n E n dS
4
SP
1 exp ikr10 1
E P0 E P1 ik n r10 dS como k
1
2
r10 r10 r10
exp ikr10
E P0
1
i
E P1
r10
cos n, r10 dS (6)
Análogamente si usamos G+
1 exp ikr10
E P0 n 1E dS
2
r10 (7)
Comparemos los resultados que se obtienen para una fuente puntual iluminando
E P A
1
eikr21
r21
Para Kirchhoff habíamos obtenido la expresión (5)
15
Óptica de Fourier C. Iemmi
A exp ik r10 r21 cos n, r 21 cos n, r10
E P0
dS
i r21r10 2
Para G_ la expresión (6) toma la forma
exp ikr21 1
n 1 E A ik n r 21
r21 r21
Vemos que las expresiones varían sólo en el factor de oblicuidad y que para el caso de
fuente y punto de observación lejanos, estos coinciden.
Se puede ver, además, que si cambiamos el punto fuente por el punto campo el resultado
no varía. Esto se conoce como principio de reciprocidad
Con anterioridad supusimos que las fuentes eran monocromáticas, veamos cómo se
modifica esta formulación cuando ello no sucede. (Lo hacemos para R – S pero para
Kirchhoff es semejante).
Para ondas monocromáticas teníamos que
16
Óptica de Fourier C. Iemmi
Dado que las ecuaciones de Maxwell son lineales, podemos sintetizar una onda
policromática como la suma de ondas monocromáticas, así
2 2
por ejemplo utilizamos la expresión (6) y, teniendo en cuenta que, k
c
exp 2 i r10
i
EPol P0 , t E0 P1 ,
c
r10
c
exp 2 i t cos n, r10 d dS
Dividiendo y multiplicando por 2 tenemos que
EPol P0 , t
1 r10 cos n, r10
0 1
2 c 2 i E P , exp 2
i t c d dS
r10
EPol P1 , t 2 i E0 P1 , exp 2 i t d por lo tanto tendremos que:
t
EPol P0 , t
1
r10 cos n, r10
2 c
E Pol 1
P , t dS (8)
t c r10
17
Óptica de Fourier C. Iemmi
Ejercicio 2:
PRINCIPIO DE BABINET
Supongamos que queremos calcular el campo producido por una pantalla con una
abertura de determinada forma en un punto P0 . Según lo visto, tenemos que:
exp ikr10
E P0 E P1
i r10
cos n, r10 dS E P1 f r10 dS
E P1 E P1 1 C P1 E P1 C P1
complementarias, por lo tanto ’ y ’’ son aberturas complementarias. Vale decir que la
distribución de campo que produce una pantalla será igual a la suma de los campos
producidos por ’ y ’’ . Esto es E E' E'' .
Por ejemplo:
18
Óptica de Fourier C. Iemmi
Comentarios:
Es evidente que si una de las pantallas es opaca, por ejemplo ’, entonces ’’
campo producido por la pantalla ’ será igual al producido por la pantalla ’’
I P0 E P0
2
ambos casos dado que
S’
L
π
19
Óptica de Fourier C. Iemmi
Aproximaciones iniciales
Hasta ahora el tratamiento, si bien escalar, fue totalmente general. Vamos a ver algunas
aproximaciones adicionales que permiten calcular figuras de difracción con métodos
matemáticos menos complejos.
Habíamos visto que para ondas monocromáticas, y según la teoría de Rayleigh –
Sommerfeld, la expresión para el campo difractado en un punto P0 era
2
exp i r10
1
E P0 E P1
i
r10
cos n, r10 dS
(6)
variable en amplitud. Por ejemplo si r10 cambia en una longitud de onda, la fase varía en
2 mientras que la amplitud casi no cambia. Así deberemos tener cuidado con qué clase
20
Óptica de Fourier C. Iemmi
Y1
Y0
r10 P0
P1
n 1 0
1 X0
X1
π
Figura 9: Geometría empleada para analizar la difracción por una pantalla plana
Esta es la aproximación paraxial ya que los ángulos que intervienen son pequeños. Si
observamos el dibujo vemos que
cos
z
; 0 1
z 2
2
2
1
cos
1
1
2 2 z
1
z
21
Óptica de Fourier C. Iemmi
Así si deseamos trabajar con un error porcentual del 1% 1.41 0.1 8º , esto es
podemos trabajar con un cono de 16º. Si en cambio admitimos un error del 5% podemos
trabajar con un cono de 36º.
1
2
y0 y1
2 2
r10 x0 x1 z2 z 1 u con
2 2
x x y y
u 0 1 0 1 1
z z
1 1 1 u u 2 1 1
1 , basta tomar para la amplitud
r10 z 1 u z 2 8
r10 z
Veamos ahora que sucede con las fases. Como dijimos, la fase es una función de
variación rápida por lo que tendremos que tomar órdenes superiores en los desarrollos en
serie. Del tipo de aproximación que tomemos surgirán los tratamientos de Fresnel y
Fraunhofer.
Aproximación de Fresnel
Primero desarrollemos en serie r10
u u2
r10 z 1 Vamos a quedarnos, en principio, con el desarrollo hasta el
2 8
término lineal, despreciando el cuadrático y superiores. Veamos en qué condiciones esto
es válido.
En realidad estamos aproximando una onda esférica eik r10 por superficies cuadráticas
r10
exp ik z ik 2
exp x0 x1 y0 y1
2
z 2z
22
Óptica de Fourier C. Iemmi
2
E x1 , y1 exp ik f x 1, y 1 dx1 dy1 ; k
1
f x1 , y1 x0 x1 y0 y1
2 2
;
2z
lenta. Veamos si para la función f x1, y1 existen dos puntos críticos xc , yc tales que
sus derivadas parciales se anulen, esto es:
f x1 , y1 1
x0 x1 0 xc x1 x0
x1 xc , yc z xc , yc
f x1 , y1 1
y0 y1 0 yc y1 y0
y1 xc , yc z xc , yc
Vemos que para f x1, y1 existen tales puntos críticos y que en ellos su derivada es
23
Óptica de Fourier C. Iemmi
xc , yc
Se podría argumentar que también en los extremos de f x1, y1 la derivada es nula, pero
la función es tan fluctuante y apretada que en los extremos la derivada prácticamente no
está definida. Entonces vemos que el término lineal en u, esto es cuadrático en x, y ,
contribuye sólo en las inmediaciones de x1 x0 y y1 y0 es decir en aquellos puntos que
se hallan fuera de la sombra geométrica. Veamos esto en la figura siguiente
exp ik z k 2
E P0 E x1 , y1 exp i x0 x1 y0 y1 dx1 dy1
2
i z 2z
(9)
Aproximación de Fraunhofer
La ecuación anterior la podemos escribir como
24
Óptica de Fourier C. Iemmi
1
exp ik z k
E x0 , y0 exp i x02 y02 .
i z 2z
k k
. E x1 , y1 exp i x0 x1 y0 y1 exp i x12 y12 dx1 dy1
z 2z
2
i k x12 y12
E x1, y1
El término 2 es la transformada de Fourier de
2z
a menos del factor e
.
difracción de Fraunhofer. Esto es, la figura de difracción de Fraunhofer, a menos del factor
1 es la transformada de Fourier del campo en la abertura . Ahora bien, el factor 1, a
menos de la atenuación 1/z es básicamente un factor global de fase. Dado que lo que se
registra es la intensidad (I=E.E*) esa fase se elimina y la figura que se obtiene coincidirá
con el módulo de la transformada de Fourier.
Veremos que esta no es la única manera de obtener E . Por ahora recordemos que la
clase pasada dijimos que si en lugar de tener una abertura, la misma tiene un factor de
transmisión t x1, y1 y por otra parte el campo que llega a la misma es E x1, y1 ,
entonces el campo inmediatamente emergente será E ' x1, y1 E x1, y1 t x1, y1 y ese
Veremos más adelante que una lente convergente de distancia focal f es justamente el
i k x12 y12
filtro que introduce el desfasaje e
2f
.
25
Óptica de Fourier C. Iemmi
FUNCIONES PERIÓDICAS
Para comenzar veamos que sucede con una función periódica arbitraria. Supongamos
que sumamos dos funciones armónicas (senos o cosenos) de distinta amplitud y período.
Vemos que tenemos una función periódica pero no senoidal. Es decir que si tenemos una
función arbitraria periódica, la podríamos sintetizar mediante la suma de senos y cosenos
con amplitudes, períodos y fases relativas adecuadas. Justamente esto es en lo que se
basa la teoría desarrollada por el físico francés Jean Baptiste Joseph Fourier (1768-1830).
Dicho de otra forma, Fourier elige una base de senos y cosenos para expandir una
función periódica. También pueden elegirse otras bases, por ejemplo la constituida por
funciones rectángulo, que conduce a las transformadas de Hadamard.
Volviendo al caso que nos ocupa, el teorema de Fourier establece que si f x es una
función de período espacial 0 con k0 2 (si fuese una función temporal en período
0
sería T y k0 2 ) entonces
T0
a0
f x
an cos nk0 x bn sin nk0 x
2 n1
(10)
26
Óptica de Fourier C. Iemmi
donde cos k0 x ; s en k0 x constituyen las componentes fundamentales. Luego, dado
Ahora debemos encontrar las expresiones para los an y bn . Para ello recurramos a la
0
sin nk0 x cos mk0 x dx 0
0
0
0 1 si n m
cos nk0 x cos mk0 x dx nm nm
0 2 0 si n m
0
0
sin nk0 x sin mk0 x dx nm
0 2
nm es la delta de Kronecker
Entonces multipliquemos (10) por cos mk0 x e integremos
0
a0 0
f x cos mk0 x dx cos mk0 x dx
0 2 0
0 0
an cos mk0 x cos nk0 x dx bn cos mk0 x s en nk0 x dx
n 1 0 0
0 0
am 0 2
f x cos mk0 x dx am f x cos mk0 x dx (11)
0 2 0 0
m=0,1,2…
0
f x dx equivale al
La expresión (11) también es válida para a0 y el término a0 1
2 0 0
27
Óptica de Fourier C. Iemmi
Multiplicando (10) por sen mk0 x e integrando, obtenemos
0
2
bm f x sin mk0 x dx (12)
0 0
m=1,2…
Cabe destacar que la integración puede realizarse sobre cualquier intervalo espacial igual
a 0 . En general entre x ' y x ' 0
Vimos que una función periódica f x se puede sintetizar como una sumatoria infinita
a0
N
a0
f x N x
2
an cos nk0 x bn sin nk0 x 2
n1 n1
0
0 N x
2
N2 x
0
Puede demostrarse que:
0
0 f x dx
2
N2 x 0
a02 1 N 2 2
an bn
4 2 n1
28
Óptica de Fourier C. Iemmi
Con anterioridad dijimos que f x era una función periódica arbitraria, sin embargo debe
valores finitos.
Supongamos que en x1 hay una discontinuidad. Entonces la serie de Fourier en x1 toma
f x1 f x1
el valor f x1
2
x1
0
iii) 0 f x dx ; f x debe ser de módulo integrable ya que es proporcional a a 0
a0 f x
2
Sea f x periódica y tal que f x . Si además f ' x es también
n1 x
continua por tramos entonces
f ' x nk0 an sin nk0 x bn cos nk0 x
n1
Existen algunas consideraciones de simetría que deben realizarse antes de iniciar el
desarrollo de la función a fin de ahorrar esfuerzos en el cálculo de los coeficientes an y
29
Óptica de Fourier C. Iemmi
0 f x cos nk0 x dx
2
4
f x f x
a
f x 0 an cos n k0 x an
2 n1 0
bn 0
IMPAR an 0
f x f x f x bn s en n k0 x
0
f x s en nk x dx
2
4
n 1 bn
0 0
0
1/2 ONDA 0
0 f x cos 2n 1k0 x dx
2
4
a2n1
0
f x f x
0 a
2 n 1 cos 2n 1k0 x
0
2 f x
100
0/2
80
n 1 b2 n 1 s en 2n 1k0 x
0
f x s en 2n 1k x dx
2
4
60 b2 n1
Este
0 0
0
40 Oeste
1/4 20
ONDA PAR Norte
0
0
f x f x
2
0 f x cos 2n 1k0 x dx
4
8
a2n1
f x a2n1 cos 2n 1k0 x
0
1er 2do 3er 4to 0
trim. trim. f x f x
trim. trim. n 1
b2 n1 0
1/4 ONDA IMPAR 0 a2n1 0
f x f x
2
f x b2n1 s en 2n 1k0 x 0
f x s en 2n 1k x dx
4
f x f x
8
n 1 b2 n1
0 0
0
30
Óptica de Fourier C. Iemmi
Forma compleja
Otra forma de representar las series de Fourier es la forma compleja. Sabemos que:
f x
Cne
ik0 nx
Cn
an ibn
n (13) con 2
Cn
an ibn
2
Usando las propiedades de ortogonalidad también se puede llegar (usar las definiciones
de an y bn ) a que
0
Cn 1 0 f x e ink0 x
dx (14)
Si f x Cn Cn
a
Función rectángulo
1
x0
x0 -a/2 x0 +a/2
31
Óptica de Fourier C. Iemmi
x x0 1
0
a 2
x x0 1 x x0 1
rect
a 2 a 2
x x0 1
1
a 2
Función escalón
a>0 a<0
1 1
x0 x0
0
x x0 0
a
x x0 1
step x x0 0 a 1 o 1
a 2
1
x x0 0
a
Función signo
a>0 a<0
1
1
x0 x0
-1 -1
32
Óptica de Fourier C. Iemmi
1
x x0 0
a
x x0
sign 0 x x0 0 a 1 o 1
a
x x0 0
1 a
x x0 x x0
sign 2 step 1
a a
Función rampa
0 x x0
x x0 1
ramp x x0
a x x0
a
x0 x0+a
1
tan
a
Función triángulo
1
x x0
x x0 0 1
tri
a
a x x0 x x0
1 1 x0-a x0+a
a a x0
área a
33
Óptica de Fourier C. Iemmi
Función sinc
sinc 1 sinc2
1
0.5
x x
sen 0
x x0 a
sinc x0 nos dice dónde está centrada y |a| es el área
a x x0
a
Función gaussiana
a
2
x x0
x x0
ae
-1/2 gaus 2b 2
a.e = a. 0.607
b
2
x x0
ae
2b 2 dx ab 2
x0 x0+b
1
a, x0 y b pertenecen a los reales. El valor de la integral es 1 si a , en cuyo caso
b 2
la función gaussiana es la función densidad de una variable aleatoria con distribución
normal de media x0 y varianza b.
34
Óptica de Fourier C. Iemmi
comb x x n
n
xn
x n
0 x n
-3 -2 -1 0 1 2 3
f x x n dx f n
x x0
comb
a
a x x0 na
n
-2 - 0 2 3
35
Óptica de Fourier C. Iemmi
2
x n ' cos nk x dx cos nn ' k
2 2 2
an
2
cos nk x
1 2
f x
n1
1
x
2. Veamos otro ejemplo f x 1 0 x
0
Conviene redefinir la función para buscar una simetría
1 x 1/2
g x 0 x
2
Ahora tenemos una función impar g(x) = -g(-x) an = 0 0
-1/2
1 x 2
4 cos nk x 2 2
bn sen nk x dx 0 sen nk x
4 x
dx
0 2 2nk 0
1 x cos nk x 1 sen nk x 1
n1
4
1 2 1n 1 1 1 1
n
2nk nk n2k 2 0 n n n
1 1
sen nkx f x sen nkx
1
g x
n1 n 2 n1 n
d d
A x
3. Veamos ahora el siguiente ejemplo f x 2 2
0
36
Óptica de Fourier C. Iemmi
A
-d/2 d/2
La función es par bn = 0;
d
2
2 2 Ad
a0
d
Adx
2
d sin nk d
d
2
2 A d 2dA 2
an A cos nkx dx sin nkx 2
2 2A
sin nk
d nk d n 2 nk d
2
2 2
sin n d
Ad 2 Ad cos 2n x
f x
n1 n d
Vemos que los coeficientes an consisten en la clásica función sinc que, si recordamos de
óptica, está relacionada con la difracción de Fraunhofer por una ranura. Habíamos visto
con anterioridad que la difracción en campo lejano correspondía a la transformada de
Fourier de la abertura. Luego volveremos sobre este tema.
Sea una función continua por tramos, con discontinuidades en los puntos .
Además su derivada está definida en todas partes excepto en dichas
las discontinuidades
f(x)
a1
a2
x1 x2
g(x)
37
Óptica de Fourier C. Iemmi
CLASE 4
FUNCIONES NO PERIÓDICAS
Hasta ahora habíamos analizado como representar funciones periódicas, veamos cómo
se puede hacer un desarrollo para aquellas que no lo son. Para ello comencemos con la
siguiente función que vimos con anterioridad
0
A d d
A x
f x 2 2
0
-d/2 d/2
0 nk0 d
2 sin
1
f x exp ink0 x dx
Ad 2
Cn
0 0 nk0 d
0 2
2
Cn
Ad/0
k0
nk0
38
Óptica de Fourier C. Iemmi
Pensemos en un análogo óptico en el que tenemos una red de difracción constituida por
bandas oscuras y transparentes como la de la figura
Los órdenes difractados por esa red estarán separados de acuerdo a la frecuencia
espacial de la misma, cuanto más alta más separados. La cantidad de luz que va a cada
orden estará modulada por la campana de difracción, que corresponderá a la de una sola
ranura. A medida que vamos disminuyendo esa frecuencia espacial, esto es aumentando
el período, dichos órdenes se irán juntando. En el límite nos quedaríamos con la luz
difractada por una sola ranura que coincide con la envolvente que antes modulaba los
órdenes. Ahora los órdenes están tan juntos que constituyen un continuo dentro de esa
campana.
Volvamos al caso que nos ocupa. Para ver que sucede en este límite escribamos
39
Óptica de Fourier C. Iemmi
0
x Cn exp i nk0 x 2k f x '
2
f
n n
0
exp i nk0 x ' dx ' exp i nk0 x
2
f x
1
f x ' exp i k x ' dx ' exp i k x dk
2 (15)
F k
1
f x F k exp i k x dk ; F k f x exp i k x dx (16)
2
i) f x debe ser tal que f x dx es decir, debe ser absolutamente integrable.
adelante para poder usar funciones tales como cos kx ; eik x es tomar
x
f x f x rect donde f x es la función que no es integrable.
a
t
t t0
Por ejemplo, un voltaje V t V0 cos wt rect
t
encendido t0 apagado
40
Óptica de Fourier C. Iemmi
Las funciones con sentido físico son aquellas que están acotadas espacial y
temporalmente, es decir, cuando escribimos coseno, en realidad estamos usando una
función recortada que dura en forma limitada, pero nosotros medimos un t ' t de
manera de despreciar los efectos de borde. Por simplicidad no se trabaja
matemáticamente con funciones limitadas, de modo que muchas veces i) se viola.
iii) f x debe tener un número finito de máximos y mínimos en cualquier rectángulo finito
(o cubo si es tridimensional)
Vimos que F k f x
f x e
ik x
dx . Esto también puede expresarse como
f x f x cos kx dx i
f x sin kx dx C
f x i
S
f x
f x F k
REAL Y PAR REAL Y PAR
IMAGINARIA Y PAR IMAGINARIA Y PAR
REAL E IMPAR IMAGINARIA E IMPAR
IMAGINARIA E IMPAR REAL E IMPAR
41
Óptica de Fourier C. Iemmi
a
f x b g x a F k bG k
k
f x F k f ax
1
F
a a
f x F k
f x a eik a F k
f x, y dx dy F kx , k y dkx dk y
2 2
f x F k
f x F k
Si f x Re F k F k
t x, y f x g y t x, y f x g y
1
F k
1
f x exp i k x dx exp i k x ' dk
2
1 1
f x
2
exp
ik x x ' dk dx f
x
2
exp ik x x ' dk x x '
42
Óptica de Fourier C. Iemmi
Ejemplos de transformadas
1. f x x x0 es la delta de Dirac
f
x f x e
ik x
dx x x0 eik x dx eik x0
f(x) |F(k)|
x0 x k
e ei k0x
i k0 x
e i ei
2. f x cos k0 x ik0 x
e eik0 x
2 2 2
i i
f
x e 2 eik0 x e
2
eik0 x
43
Óptica de Fourier C. Iemmi
eik0 x eik0 x eik x dx ei k k0 x dx 2
k k0
Análogamente
eik0 x eik0 x eik x dx ei k k0 x dx 2
k k0
i
ei
f
x e 2 2 k k0 2
2 k k0
Si 0 es la transformada del coseno
k k0 k k0
cos
k x
0
2
2 2
f(x) F(k)
-k0 k0
Aca se debe tener cuidado en cómo se grafican las deltas ya que son imaginarias.
44
Óptica de Fourier C. Iemmi
f(x)
R
Si f
x 2
k k0 k k0
2 2
f(x) R
Convolución y correlación
Vamos a definir ahora dos funciones que resulta importantísimo conocer cuando uno
trabaja con transformadas de Fourier. Por ahora sólo veremos su interpretación
matemática, para luego aplicarlas en problemas de difracción y sistemas lineales.
Correlación
45
Óptica de Fourier C. Iemmi
q fg x f x g x f g x d
*
(17)
Si f x = g x se denomina autocorrelación.
Examinemos el proceso que realiza esta operación en forma gráfica ya que no resulta
f()
f rect
4 -2 0 2
46
Óptica de Fourier C. Iemmi
f() 1 f(-x)
-4 -3 -2 -1 0 1 2 3 4
x
4
qff(x)
3
Convolución
C fg x f x g x f g x d (18)
3 f()
f ramp step
3/ 2 1 2 ramp
1 step
1 2 3
47
Óptica de Fourier g() C. Iemmi
1
1
g rect
4
-1 0 1 2 3
g(-)
1
1
g rect
4
-3 -2 -1 0 1
2
f()
g(-+x) 1
-4 -3 -2 -1 0 1 2 3 4
cfg(x)
48
Óptica de Fourier C. Iemmi
Debemos tener en cuenta que las condiciones para que q f g x y c f g x existan son
f x g x g x f x conmutatividad
f x g x t x f x g x t x asociatividad
f x g x t x f x g x f x t x propiedad distributiva
f x g x f x g x
* * *
conjugación
f x a x b f x a b
f x x f x
f x x a f x a
f x f x f x f x
f x g x g x f x en general no es conmutativa
f x , g x f x g x g x f x
f x , g x ; g x g x f x g x f x g x
f x g* x f x g x
f x g x F k G k
f x g * x F k G* k F k
G k
f x g x F k G k Teorema de la convolución
f x
g x F k G* k
si G k f x
g x F k G k y si además
49
Óptica de Fourier C. Iemmi
G k G k f x
g x F k G k
f x f x F k F k
2
si es par
f x f x F k F k
2
si
Volvamos ahora a los temas de óptica que nos ocupan, pero tratemos de comparar las
expresiones anteriormente obtenidas con las funciones matemáticas que hemos
estudiado.
La operación convolución en dos dimensiones puede expresarse como
C fg x, y f x, y g x, y f , g x , y d d
exp ik z
E x0 , y0 E x1, y1, z 0 k 2
exp i x0 x1 y0 y1 dx1 dy1
2
i z 2z
h
E x0 , y0 E h
Volveremos sobre esta expresión cuando veamos sistemas lineales.
50
Óptica de Fourier C. Iemmi
exp ik z k x
exp i x02 y02 E x1 , y1 exp i 2 0 x1 0 y1 dx1 dy1
y
E x0 , y0
i z 2z z z
2 2
f x, y f x, y exp i k x x k y y dx dy con kx ; ky
x y
x0 1
fx
z x
Si definimos las frecuencias espaciales como entonces
f y 1
0
y
z y
A fin de aclarar un poco más estos conceptos veamos algunos diagramas de difracción
básicos y de suma importancia.
exp ik z k 2
E P0 E x1 , y1 exp i x0 x1 y0 y1 dx1 dy1
2
i z 2z
51
Óptica de Fourier C. Iemmi
Abertura rectangular
x1 y
E x1, y1 rect rect 1
a b
a b
exp ik z 2 k 2 k
2
2
E P0
i z a
exp i x0 x1 dx1
2z 2z
exp i y0 y1 dy1
2 b
2
z ; v k y y dy dv z
u
k
z x1 x0 dx1 du
k z 1 0 1 k
a k a b k b
u1 x0 ; v1 y0
2 z 2 2 z 2
a k a b k b
u2 x0 ; v2 y0
2 z 2 2 z 2
Entonces
exp ik z z u 2 2
v2
2
E x0 , y0 exp i u du exp i v dv
i z k u 2 v 2
1 1
u'
Es conveniente llevar estas integrales a la forma 0 , esto es
0 u 0 v u1 u2 v1 v2
E x0 , y0
eik z 2
2
eik z
2i u 0 v 0
2i 0 0
1
1
0 0
52
Óptica de Fourier C. Iemmi
w w
C w cos w '2 dw ' ; S w sin w '2 dw '
0 2 0 2
un un
2
exp i u du cos u 2 i sin u 2 du C
2
un i S un n 1, 2
0 0 2 2
Entonces
2 2 2 2
I x0 , y0 C u2 C u1 S u2 S u1 C v2 C v1 S v2 S v1
1
4
53
Óptica de Fourier C. Iemmi
S(w)
S(w1)
Z(w1) - Z(w0)
0.5
S(w0)
w0 es el largo de la
Z(w1) curva
Z(w0)
-0.5
C(w)
0.5 C(w1) C(w0)
-0.5
w w
Z w C w i S w cos w '2 dw ' i sin w '2 dw '
0 2 0 2
2 2
dC dS
dl dC dS
2 2 2
dw
2
cos2 w2 sin 2 w2 dw2 dw2
dw dw 2 2
dl dw
54
Óptica de Fourier C. Iemmi
graficarse
Z w1 Z w0 C w1 C w0 i S w1 S w0 . Pero esta
x , y0 no es más que Z w1 Z w0
2
describe la intensidad I 0 .
k a k b
u1 x0 ; v1 y0
z 2 z 2 b
k a k b
u2 x0 ; v2 y0
z 2 z 2
k k
Si hacemos u2 u1 a ; v 2 v 1 b encontramos expresiones que
z z
son independientes del punto de observación, sólo dependen de la geometría de la
abertura y son las que determinan el largo de la curva sobre la espiral (no su
posición). Por ejemplo
x0 = 0 u2 = -u1 x0 = -a/2 u1 = 0
Z(u2)-Z (u1)
Z(u2)-Z(u1)
55
Óptica de Fourier C. Iemmi
k
valiendo la aproximación de Fesnel. En este caso es muy grande y podemos
z
tomar
a a
x0 2 x0 2
u1 ; u2
x a x a
0
2 0
2
b b
y0 y0
v1 2 ; v 2
2
y b y b
0
2 0
2
1 a 1 a
x0 x0
2 ; C u S u 2
C u1 S u1 2 2 2 2
1 x a 1 x a
2 0
2 2 0
2
1 b 1 b
y0 y0
C v1 S v1 2 2 ; C v S v 2
2 2
2
1 y b 1 y b
2 0
2 2 0
2
exp i k z 1 1 1 1 1 1 1 1
E x0 , y0 x a ; y b i i
0 2 02 2i 2 2 2 2 2 2 2 2
exp i k z x y
2i exp i k z rect 0 rect 0
2i a b
56
Óptica de Fourier C. Iemmi
E x0 , y0 0
x0 a ; y0 b
2 2
exp ik z k k
E x0 , y0 exp i x02 y02 E x1 , y1 exp i x12 y12 .
i z 2z 2z
*
k
.exp i x0 x1 y0 y1 dx1 dy1
z
vemos que podemos interpretarla como la transformada de Fourier del término señalado
con la llave, o sea del campo en la abertura por un factor de fase cuadrático. Para más
detalles puede verse el trabajo [Appl. Opt. 38, 7085 (1999)].
En esta aproximación veremos los diagramas generados por una red cuya amplitud varía
sinusoidalmente, una cuya fase varía sinusoidalmente y una abertura circular.
E0(x,y)
E x, y
t x, x e
i x , y
t x, x
E0 x, y
E(x,y)
t(x,y)
57
Óptica de Fourier C. Iemmi
1 m y
t x, y cos 2 f 0 x rect rect
x
2 2 a b
espacial de la red. Las funciones rectángulo dan cuenta del tamaño finito de la red y
0<m 1 es la modulación.
t(x,y)
m 1/2
0 = 1/f0
Supongamos que sobre la red incide normalmente una onda plana, monocromática de
A e i k n
A k
n
58
Óptica de Fourier C. Iemmi
1 m m x y
t x, y exp i 2 f 0 x exp i 2 f 0 x rect rect
2 4 4 a b
En el producto de los rectángulos no hace falta aplicar la convolución dado que son
variables separables. Calculemos entonces las transformadas
1 1
exp i 2 f x x exp i 2 f y y dx dy f x , f y
1
2 2 2
ik x dx 2 k pero existe una propiedad de las deltas que
conviene recordar que
e
x 2 k 2 2 f f
1
establece que cte x
cte
m m
f x f0 , f y
m
4
Análogamente
m m
4 exp i 2 f 0 x 4 f x f 0 , f y
59
Óptica de Fourier C. Iemmi
x y x
y
rect rect rect exp i 2 f x x dx rect exp i 2 f y y dy
a b a b
a b
exp i 2 f x x dx exp i 2 f y y dy
2 2
a b
2 2
a b
a b
2 2
2 sin 2 f y 2
b b
sin 2 f x a 2
a
cos 2 f x
2
cos 2 f y
i i
x x y y
2 f 2 f a
2 f 2 f
2
x a x y y
2
2 b b
2
sin f x a sin f y b
a b sinc f x a sinc f y b
fx fy
Así pues
1
t x, y f x , f y f x f 0 , f y f x f 0 , f y
m m
2 4 4
a b sinc f x a sinc f y b
sinc f y b sinc f x a sinc f x f 0 a sinc f x f 0 a
ab m m
2 2 2
sin px
Teniendo en cuenta que sinc px tiene el máximo principal en x0
px
, grafiquemos E x0 , y0
1
y el primer cero cuando px x
p
60
Óptica de Fourier C. Iemmi
f0 f0
fy
Si a y b tienden a infinito entonces las funciones sinc tienden a ser una delta. Por otra
2
parte dado que f 0 la intensidad será
a
a 2b 2 m2 m2
I E.E *
4 2 z2
sinc 2
b f y sinc 2
a f x sinc 2
a f x f 0 sinc 2
a f x f 0
4 4
ya que los términos cruzados se pueden despreciar por estar las tres funciones sinc
centradas en puntos muy alejados. Vemos entonces que una red cuya transmisión en
amplitud varía sinusoidalmente da lugar al orden 0 y sólo al 1. Obviamente siempre
que d0>> que es una de nuestras hipótesis de trabajo.
Examinemos este resultado que es muy importante. Habíamos visto que las
frecuencias espaciales que aparecían en la expresión de la transformada de Fourier
eran
x0
fx x0 f x z
z
f y0 y0 f y z
y z
61
Óptica de Fourier C. Iemmi
Por otra parte acabamos de ver que un objeto con frecuencia espacial f 0 en el eje x
x0 x0 x0
sin f 0 x0 f 0 z
z z z
Vale decir que cada frecuencia espacial presente en el objeto contribuirá con un orden
en el plano transformado, cuya posición estará unívocamente relacionada con su
frecuencia espacial y la orientación de la misma. A la inversa, si nosotros montamos la
experiencia de Young de forma tal que la separación entre fuentes sea x0, a una
distancia z obtendremos una figura de interferencia cosenoidal con frecuencia
x0
f0 .
z
Este ejemplo que acabamos de ver es muy importante porque, de acuerdo a lo que
estudiamos de Fourier, una función arbitraria (por ejemplo una imagen en una
diapositiva) la podemos descomponer en senos y cosenos de distintas frecuencias y
orientaciones, por lo tanto podemos imaginarla como conformada por una
superposición de redes sinusoidales de amplitudes, frecuencias y orientaciones
distintas; ahora ya sabemos cómo se comporta cada una de ellas en el plano
transformado.
CLASE 6
Es muy común decir que una red sinusoidal produce sólo tres órdenes difractados, el 0, el
1 y el -1. Debe remarcarse que esto es cierto si la red es de amplitud. Para estudiar qué
sucede si lo que varía sinusoidalmente es la fase, veamos este segundo ejemplo.
62
Óptica de Fourier C. Iemmi
1 x y
t x, y exp i k m cos 2 f 0 x rect rect
2 a b
Una forma de hacer esto es, por ejemplo, modelando un sustrato transparente con
una forma de superficie . Otra, es elaborar una placa de
un material transparente pero cuyo índice de refracción varíe de forma sinusoidal
(cuando estudiemos materiales de registro veremos esto en más detalle). El término
exp i a cos J n a exp i n exp i n
n 2
Donde Jn a es la función J de Bessel de primera especie y de orden n.
a . b .sinc a f x .sinc b f y
n
t x, y a . b .sinc b f y J n exp i sinc a f x n f 0
n 2 2
63
Óptica de Fourier C. Iemmi
ancho 2/a. Vemos que la separación entre órdenes así como su ancho es una
cuestión geométrica ya que es igual en la red de amplitud, sin embargo acá aparecen
órdenes superiores y la amplitud de cada uno de ellos está determinada por la función
n
i
2
J n e 2 . Suponiendo nuevamente que f 0
2 a
2
2
a .b a
I x, y sinc b f y . J n .sinc2 x0 n zf 0
z n 2 z
Con este tipo de redes es posible distribuir más energía en los órdenes laterales que
en el central, según el elegido.
3. Abertura circular
Y1
Y0, fy
r1 X0, fx
X1
π
fx , f y difieren de x0 , y0 en un factor 1
z
64
Óptica de Fourier C. Iemmi
f x2 f y2
r1 x12 y12
fy
y1 arctan
arctan fx
x1
r0 x02 y02 z
r 1 r1 a
t x1 , y1 t r1 circ 1
a 0 r1 a
x1 r1 cos ; f x cos
y1 r1 sin ; f y sin ; dx1 dy1 r1 dr1 d
2
r1
t r1 circ a exp i 2 r cos cos sin sin r dr d
1 1 1
0 0
2
r1
circ r1 exp i 2 r1 cos d dr1
0 a 0
65
Óptica de Fourier C. Iemmi
2 2
in 1
Jn x exp i x cos exp i n d para n 0 J 0 x exp i x cos d
2 0
2 0
x 2 r1
Luego si identificamos tenemos que d d y los límites de
integración que para eran de 0 a 2 pasan a ser para de a 2 .
Por otra parte, teniendo en cuenta que el signo de la exponencial no tiene
importancia por ser J0 una función par y que integrar entre 0 y 2 ó entre y
2 es igual por ser una función de período 2 . Entonces podemos escribir
r
t r1 2 circ 1 r1 J 0 2 r1 dr1
0 a
a
t r1 2 r1 J 0 2 r1 dr1
0
2 2 a
t u u J 0 u du
2
2
0
pero existe una propiedad de las funciones J de Bessel que establece que
2 a
d 2
u n 1 J n1 u u n 1 J n u t u u J1 u
2
2
du
0
Entonces
exp ik z kr02 za 2 a r0
E x0 , y0 exp i J1
i z 2 z 0
r z
66
Óptica de Fourier C. Iemmi
r0
donde se tuvo en cuenta que . La intensidad vendrá dada por
z
2
k ar0
J 1 z a2k
2
I r0 I 0 2 ; I0
k ar 0 2z
z
Se elige agruparlo de esta forma para que la expresión entre corchetes tienda a 1
cuando r0 tiende a 0. Veamos algunas características de esta figura conocida
I(r0)/I0
0.0175 0.0042
67
Óptica de Fourier C. Iemmi
z = fob
L
estrellas Oc
Ob
kar0 z
1.22 r01min 1.22
z 2a
68
Óptica de Fourier C. Iemmi
Nx partículas
2 L + S Z 5
Con anterioridad habíamos visto que la integral de difracción de Fresnel venía dada por
exp ik z k
E x0 , y0 exp i x02 y02 .
i z 2z
k k
. E x1 , y1 exp i x0 x1 y0 y1 exp i x12 y12 dx1 dy1
z 2z
2
y que el término 2 era la transformada de Fourier de E x1, y1 a menos del factor
i k x12 y12
e 2z
. También vimos que una de las posibilidades para eliminarlo era pedir
z
k 2 2
x y
2 1 1 Max
esto es, estar en las condiciones de difracción de Fraunhofer. La otra
69
Óptica de Fourier C. Iemmi
i k x12 y12
consistía en introducir un factor de fase e
2z
. Vamos a estudiar ahora que
dicho factor de fase es justamente el introducido por una lente convergente.
Cuando un haz de luz atraviesa un elemento constituido por un medio transparente (en el
caso de las lentes generalmente es vidrio), su amplitud permanece casi inalterada pero su
distribución de fase puede variar fuertemente, dependiendo de la forma de dicho
elemento. Consideremos el siguiente dibujo
Y
(y)
Luz
n=1 n n=1
0 Z
y k n y 1 0 y k 0 k n 1 y
k n 1 y es el desfasaje que sufre el rayo con respecto a otro que hubiese viajado
forma del elemento transparente. Veamos entonces la expresión de y para una lente
esférica
70
Óptica de Fourier C. Iemmi
A B C
R1
y
zB Z
(R12-y2)1/2
zC
Vamos a introducir una convención de signos: la luz incide desde la izquierda y las
superficies convexas que encuentra ( ( ) tienen radio positivo, mientras que las
tiene radio positivo y la segunda negativo. Vamos a separar entonces la lente en dos
partes. Sobre la primer superficie tendremos:
z R R 2 y 2 y2
B 1 1
; S 1 y zc zb 01 R1 1 1 2
zC 01 R1
y2
S 2 y 02 R2 1 1 2
R2
Esto, para simplificar, lo hemos planteado en una dimensión pero suponemos que tiene
simetría de revolución (en el caso de tratarse de una lente cilíndrica, convergente,
obtendríamos la transformada de Fourier en una sola dimensión). Entonces
71
Óptica de Fourier C. Iemmi
x2 y 2 x2 y 2
x, y x, y x, y 01 02 R1 1 1
S1 S2
R2 1 1
R12 R22
x2 y 2 x2 y 2
x,y pequeños. En este caso 1 1 . Así resulta
R2 2R2
x2 y 2 1 1
x, y 0
2 R1 R2
En consecuencia la fase introducida por la lente será:
x2 y 2 1 1 x2 y 2 1 1
x, y k 0 k n 1 0 k n 0 k
n 1
2 R1 R2 2 R1 R2
Recordemos la fórmula del constructor de lentes que relaciona las propiedades de la lente
1 1 1
(n, R1 y R2) ; n 1 . Con lo cual la expresión para la fase queda
f R1 R2
k x2 y 2
x, y k n 0
2f
Cabe destacar que a pesar de la elección de la forma () de la lente adoptada para esta
deducción, la expresión de la fase es válida para cualquier combinación de radios,
siempre que se tenga en cuenta cuales son positivos y cuales negativos. Según su
geometría tenemos que:
Tienen f >0
72
Óptica de Fourier C. Iemmi
Tienen f <0
k
tL x, y P x, y exp i k n 0 exp i x 2 y 2
(19)
2f
A i k z i 2z x y2
k 2
A ik r
esférica e era reemplazada por su aproximación cuadrática e e
r
z
Consideremos que sobre una lente incide normalmente una onda plana de amplitud
k
E ' x, y 1 t x, y exp i x 2 y 2
2f
Vemos que este término puede considerarse como una aproximación cuadrática a una
onda esférica que, si f >0 converge al plano focal detrás de la lente (justamente f >0
corresponde a lentes convergentes) y si f <0 la onda diverge de una fuente virtual situada
a una distancia f delante de la lente (en este caso la lente es divergente)
73
Óptica de Fourier C. Iemmi
f f
Lente convergente Lente divergente
CLASE 7
x x'
74
Óptica de Fourier C. Iemmi
o L d
xo xL xd
z
do d
Primero debemos efectuar una propagación libre desde el objeto situado en el plano o
hasta el plano de la lente L . Entonces, según la convención adoptada deberemos hacer
los siguientes reemplazos en la integral
z do
exp i k do k xL2 k xo2 k xL xo
x xo Ex
exp i E x exp i exp i dxo
i do
L o
x ' x 2do 2d o do
L
k xL2
E x E x t L xL E x P xL exp i
L
L
L
2f
Ahora haremos nuevamente una propagación libre desde el plano de la lente L hasta un
plano genérico d situado a una distancia d de la lente.
75
Óptica de Fourier C. Iemmi
exp i k d k xd2
E xd exp i .
i d 2d
exp i k d o k xL2 k xo2 k xL xo
o o
exp i E x exp i exp i dx
. i d o 2 d o 2 d o d o
k x2 k x2 kx x
.P xL exp i L exp i L exp i L d dxL
2f 2d d
2
b2
exp i a x b x dx a exp i 4 exp i 4a
k1 1 1 xo xd
Con lo cual si tomamos a ; b k ; P xL 1
2 d do f do d
tenemos que:
76
Óptica de Fourier C. Iemmi
2 xf
Recordando que k y que la frecuencia espacial f X podemos escribir
f f
fX f
E xo E x exp i 2 f
o Xf
xo dxo
Con lo cual llegamos al resultado final, que en dos dimensiones toma la forma
exp i k f d o k x 2f y 2f d
E xf , yf exp i 1
o
i f 2f f (20
f X f , fY f o o
. E x , y )
Vemos entonces que el campo en el plano focal de la lente es, a menos de un factor
global de fase, proporcional a la transformada de Fourier del campo que emerge de la
1
transparencia. Como lo que se mide, en general, es la intensidad I E E
f2
2
el factor global de fase se elimina. De todos modos vemos que si ubicamos el objeto en
77
Óptica de Fourier C. Iemmi
Otro caso particular a considerar es cuando el objeto está contra la lente. Basta tomar
para ello do 0 .
Para llegar a la ecuación (20) hemos supuesto P xL , yL 1 despreciando el tamaño
Yo
Xo YL
do
Yf
o XL
f
Xf
L
xf , yf
f
78
Óptica de Fourier C. Iemmi
Dado que el punto x f , yf se halla en el plano focal, la luz que arriba a dicho punto
provendrá de un haz paralelo cuya dirección está especificada por los ángulos
xf yf
; . Sin embargo de todos los rayos que viajen paralelos a dicha dirección
f f
sólo algunos entrarán en la lente, debido a su tamaño finito. Los que entren serán
justamente los que están contenidos en la proyección hacia atrás, en la dirección ;,
de la pupila. Esto nos delimita en el plano objeto o cuáles son los puntos que
exp i k f d o k x 2f y 2f d
E xf , yf exp i 1
o
i f 2f f
do x f do y f
E xo , yo P xo
(21)
, yo
f X f , fY f
f f
Analicemos ahora el caso del objeto ubicado a una distancia d o del plano focal de la lente
79
Óptica de Fourier C. Iemmi
L o f
xL xo Xf
z
do
objeto incide una onda esférica que converge sobre el plano focal f .
En este caso también debemos estudiar la propagación por tramos. Para ello analicemos,
al igual que antes, el caso unidimensional para simplificar los cálculos y luego
extendemos el resultado final a dos dimensiones. Inmediatamente a la salida de la lente el
campo será
k xL2
E x A P xL exp i
L
2f
emergente de la lente con el plano o. Esta área se representa matemáticamente por la
proyección de la pupila P xL a través del cono de rayos que cae sobre el objeto. Esto
x f x x
es, sobre el plano objeto tenemos una función pupila efectiva P o ya que L o .
d0 f do
x
Por ejemplo si la lente es circular y de radio R entonces P xL circ L . Ahora
R
bien, la proyección de este círculo sobre el plano o corresponderá a otro círculo de radio
80
Óptica de Fourier C. Iemmi
x x f xo f
P ' xo circ o
Rd o circ o
con lo cual P . Si el área iluminada
f Rd o Rd o do
f
x f
es mayor que la extensión del objeto entonces se toma P o 1.
d0
Luego de este análisis, hagamos la propagación libre desde la lente hasta el objeto.
Entonces ahora, de acuerdo a la nomenclatura utilizada, tenemos que:
z f do
exp i k f d o k xo2
x xL Ex
exp i
i f d o 2 f d o
o
x' x
o
x f k xL2 k xL2 k xL xo
AP o exp i exp i exp i dxL
do 2 f
2 f d
o f do
exp i a x 2
b x
dx
a
exp i
4
exp
b2
i
4a
kdo k xo
a ; b
2 f f do
y tomando tenemos que
f do
81
Óptica de Fourier C. Iemmi
k x2 x f
E xo exp i k f do
f
exp i o AP o exp i
ido 2d o do 4
z do
exp i k d o k x 2f
x xo
x ' x
E xf
i d o
exp i
2d
o
exp i k f d o exp i 4
f
id o
A.
f
x f k xo2 k xo2 k x f xo
t xo P o exp i exp i exp i dxo
do 2 d o 2 d o d o
k x 2f xo f
E x f exp i k f exp i
4
f
i 2 do2
A exp i
2d
o
fX f t xo P
do
k x 2f y 2f
E x f , y f exp i k f
Af
exp i
i d o2 2d o
(22)
xo f yo f
.
f X f ,Y f t xo , yo P ,
o do
d
obtenemos el mismo resultado ya que ambos corresponden al caso del objeto contra la
lente.
82
Óptica de Fourier C. Iemmi
Sin embargo cuando trabajamos con el objeto delante de la lente obtuvimos frecuencias
xf yf
espaciales fX ; fY para cualquier posición del objeto. Al contrario, para el
f f f f
xf yf
objeto situado detrás de la lente, encontramos que fX ; fY por lo tanto
f do f do
el tamaño de la transformada depende de la posición del objeto.
Esta propiedad puede resultar muy útil en las aplicaciones de filtrado espacial, tal como
veremos más adelante.
s o L T
Xs xo xL xT
S
z
do dT
ds
83
Óptica de Fourier C. Iemmi
k x2 k xo2 k xo xs k xL xo k xL2
. t xo exp i o exp i exp i exp i exp i dxo
2d o 2 d s d o d s do do 2 f
k xo2 d s k xo xs
t xo exp i exp i
2d o d s d o d s do
1 1 1 2 xo xT
exp i k xL exp i k xL dxL dxo
2d o 2dT 2 f d o dT
b2
i a x 2
b x i i
4a ; a k 1 1 1 ; b k xo xT
e dx
a
e 4 e
2d o 2dT
2f do dT
84
Óptica de Fourier C. Iemmi
k1
1 1
exp i
4
exp i
2 1
k xT2
1 1
2 d T
2 d o dT f o
d d T f
ds 1 x2
. t xo exp i k
o
2d o d s d o 2d 2 1 1 1
o
d o dT f
x dx
.exp i k
x x
s
T
o
d s do 1 1 1
o
dT d o
o
d d f
T
k 1 1 x2
E xT exp i T
1 1 1
2 dT dT
od d f
T
ds 1 x2 .
t xo exp i k
o
2d o d o d s 2d 2 1 1 1
o
d d f
o T
85
Óptica de Fourier C. Iemmi
xs xT x dx
exp i k o
do d s 1 1 1
o
dT d o
od d f
T
1 1 1 1 1 1
esto es . Reemplacemos esta expresión en la de E xT
dT d s f dT d s f
k d s f f d s d s d o 2
E xT exp i xT
2d s f f d s f d o
f
x T 1
d x dx
. t xo exp i k
x
s
o o
s
do d s
f 1 o
d
d s
Vemos que, aunque con una expresión complicada para las frecuencias espaciales, la
Analicemos por último el caso en el que el objeto está detrás de la lente y la fuente se
encuentra a una distancia finita ds de ella.
86
Óptica de Fourier C. Iemmi
s L o T
Xs xL xo XT
z
do
ds dT
A continuación realicemos una propagación libre justo hasta después del objeto
k xs2 k xo2
E x exp i
exp i
2 dT d o
o
2d s
k 1 1 k xL2 k xs xL k xo xL
. t xo exp i xL2 exp i exp i exp i dxL
s
2 d f 2 d T d
o d s d T d
o
87
Óptica de Fourier C. Iemmi
Acá nuevamente debemos compatibilizar desde donde medimos las distancias con la
convención de signos de la fórmula de las lentes, de modo tal que cambiaremos d s por
k x2 k 1 1 kx x
E xT exp i T t xo exp i xo2 exp i o T
2d o 2 dT d o d o do
k 1 1 1 2 x xo
exp i xL exp i k s xL dxL dxo
d d d
2 d s f dT d o s T o
2
b2
exp i a x b x dx a exp i 4 exp i 4a
1 1 1 xs xo
a k ; b k
2d s 2 dT do 2 f
y tomaremos
d s d T d o
88
Óptica de Fourier C. Iemmi
k do x x xo
a ; b k s s
2 dT dT do d
T f dT d o
k xT2 k 1 d 2 2
E xT exp i exp i dT d o T2 xs .
2d o 2 dT d o d o f d o f
x x xd
. t xo exp i k T s s T xo dxo
do do f do
Hasta ahora todo el análisis que hemos hecho ha sido para iluminación monocromática,
esto es, coherente. Sin embargo sabemos que ninguna fuente es estrictamente
monocromática (ni siquiera un láser) y por otra parte todas las fuentes presentan un cierto
grado de coherencia (aún el sol).
Si bien el procesado óptico coherente es elegante desde el punto de vista matemático,
experimentalmente presenta algunas desventajas. Por ejemplo el ruido coherente
(speckle) afecta la calidad de la imagen final, es necesario trabajar con elementos
perfectamente limpios y en un ambiente libre de polvo ya que cualquier partícula causará
difracción, limita el tipo de fuentes luminosas que pueden emplearse, admite sólo ciertos
objetos de entrada (quedan excluidos los objetos autoluminosos, y aquellos que
introduzcan fases espurias), la luz empleada es monocromática por lo tanto se pierde la
información relacionada con el color del objeto, etc. Sin embargo no todo tipo de
procesado puede llevarse a cabo con iluminación incoherente.
89
Óptica de Fourier C. Iemmi
Qué significa que una fuente sea coherente, incoherente o parcialmente coherente? Que
significa coherencia temporal y coherencia espacial? Algunos de estos conceptos los
estudiaremos en esta capítulo ya que son necesarios para analizar el fenómeno de
formación de imágenes y todos los procesos relacionados desde un punto de vista más
completo.
Nosotros veremos solamente los fundamentos de la Teoría de la Coherencia Parcial ya
que es sumamente extensa, de hecho toda una rama de la óptica se dedica a realizar
estos estudios. Un tratamiento más extenso puede encontrarse en Statistical Optics de
J.Goodman o en Principles of Optics de M.Born y E.Wolf y aún más profundo en Optical
Coherence and Quantum Optics e Introduction to the Theory of Coherence and
Polarization of Light, ambos de E.Wolf.
Las primeras investigaciones sobre este tema fueron realizadas por Verdet en 1869, quien
encontró que bajo ciertas condiciones fuentes comúnmente consideradas incoherentes
(como por ejemplo el sol) producían interferencia al realizar una experiencia tipo Young.
Sin embargo fue Michelson en 1890 el primero en establecer la relación existente entre el
contraste de las franjas de interferencia y la distribución de intensidad de una fuente
extensa.
Algunos de los desarrollos más importantes sobre la teoría de coherencia parcial fueron
realizados por Van Cittert en 1934 y Zernike en 1938. Ellos determinaron cuál es el grado
de coherencia entre dos puntos cualesquiera de una pantalla iluminada por una fuente
extensa. Otro avance importante lo realizó Wolf en 1957 con la introducción de la función
coherencia mútua.
Antes de comenzar con el desarrollo de la teoría conviene aclarar algunos conceptos. A
menudo resulta útil dividir los efectos de coherencia en dos clasificaciones: coherencia
temporal y espacial. La primera se relaciona con el ancho de banda finito de la fuente y la
segunda con el tamaño de la misma. Esta división resulta un tanto artificial ya que ambas
características suelen estar, como veremos, íntimamente relacionadas.
Por ejemplo a veces se habla de una fuente “extensa monocromática” para describir a
una fuente extensa que emite con una cierta longitud de onda media y ancho de
90
Óptica de Fourier C. Iemmi
emite una onda cuya diferencia de fase con la emitida por otro punto varía al azar, esto
es, tal diferencia no se mantiene constante. Por definición una onda monocromática existe
desde t por lo que no tiene sentido hablar de diferencias de fase variables en el
tiempo. Por ejemplo tomemos un láser como una aproximación a una onda
monocromática (teniendo claro que no lo es) y para simular una fuente extensa
hagámoslo incidir sobre un vidrio esmerilado. Si el difusor permanece quieto dicha fuente
no es extensa ya que la diferencia de fase entre dos puntos cualesquiera permanece
constante. Será en todo caso un frente de ondas con retardos distintos por zonas, pero
fijos. Si en cambio tal difusor se hace rotar, tenemos entonces una buena aproximación a
r2
donde 1 y 2 son las fases iniciales. S2
*
El campo total será E E1 E 2 y la intensidad vendrá dada por I EE donde el
01 02
I E 2 E 2 E 01 E 02 exp i k r1 r2 1 2 exp i k r1 r2 1 2
E 2 E 2 2 E 01 E 02 cos k r1 r2 1 2
01 02
(23)
91
Óptica de Fourier C. Iemmi
M1
M2
S
BS
92
Óptica de Fourier C. Iemmi
r1
S
P
r2
diferencia de camino entre SS1 y SS2 por lo tanto podremos incorporarlo a la diferencia
de caminos total r .
X X’
D. S1
P
S
D
S’
S2 D.
L L’
Para estudiar el efecto que tiene sobre la figura de interferencia el tamaño de la fuente
comencemos analizando que sucede si tenemos dos fuentes puntuales primarias S y S’,
de igual intensidad.
Supongamos también que ambas aberturas S1 y S2 , tienen la misma intensidad
93
Óptica de Fourier C. Iemmi
I1 E01
2
I 2 E02
2
I 0 . Si reemplazamos en la ecuación (23) obtenemos que cada
fuente primaria producirá sobre el plano X’ un diagrama de franjas descripto por
k
I 4 I 0 cos 2 r
2
r SS2 S2 P SS1 S1P S2 P S1P D
Para S’
r ' S ' S2 S2 P S ' S1 S1P D D
un caso (para S) cae en 0 y en el otro (para S’) en . Así pues, tendremos dos
sistemas de franjas uno para S y otro para S’ que se suman en intensidad, porque S y S’
son incoherentes, y que se hallan desplazados provocando una caída del contraste.
k k
I 4 I 0 cos 2 D cos 2 D
2 2
kD
dI A cos 2 d
2
Así la intensidad total será
94
Óptica de Fourier C. Iemmi
max
2 kD
I A cos d max
max 2
sin k D max
Amax 1 cos k D
k D max
Hasta ahora nosotros trabajamos con campos monocromáticos y adoptamos, por ser útil
matemáticamente, la notación compleja.
E R t E t
Dado que ahora comenzaremos a trabajar con campos no monocromáticos deseamos
encontrar una notación compleja para representarlos, es decir, buscamos una
generalización de lo hecho para campos monocromáticos.
95
Óptica de Fourier C. Iemmi
0
E R t exp 2 i t d exp 2 i t d
0
exp 2 i t d exp 2 i t d
0 0
Pero es la antitransformada de E R
t , o sea E R t exp 2 i t dt
Entonces
*
0 exp 2 i t d 0 *
exp 2 i t d exp 2 i t d
0
Con lo cual
E t 2 exp 2 i t d 2 exp 2 i t d
R
0 0
Si ahora llamamos E t 2 exp 2 i t d podemos escribir un campo
I
0
policromático complejo como
E t E R
t i E t 2 exp 2 i t d
I
96
Óptica de Fourier C. Iemmi
Consideremos una fuente extensa con un ancho de banda finito y dos puntos
arbitrarios en el espacio S1 y S2 . En dichos puntos la perturbación óptica a tiempo t es
d S2
r2
E2(t)
P
S1
r1
E1(t)
originado por la superposición del que a tiempo t t1 fue emitido por S1 y el que a
tiempo t t2 fue emitido por S2. Los factores K1 y K2 están relacionados con las
97
Óptica de Fourier C. Iemmi
EPol P, t
1
r1 cos n, r1
E Pol 1
P , t dS
2 c t c r1
Si suponemos que los pinholes son muy pequeños y que el campo sobre su superficie se
mantiene constante, entonces podemos sacarlo de la integral de modo que
EPol P, t
1
cos n, r1
E P , t r1
2 c t Pol 1
dS
r1 c
Tengamos en cuenta esta expresión para más adelante y recordemos que el factor entre
r r
EPol P1 , t 1 P1 , exp 2 i t 1 d se reduce a
c c
r r
E P1 , t 1 P1 , exp 2 i t 1 y por lo tanto
c c
r r1 2 r
E P1 , t 1 2 i P1 , exp 2 i t E P1 , t 1
t c c i c
con lo cual
E P, t
cos n, r1
r
dS E P1 , t 1
i c r1 c
Vemos que justamente el factor entre corchetes es el propagador K1 que en este caso es
un número imaginario puro.
Volviendo al tema que nos ocupaba, la intensidad en el punto P producida por los dos
pinholes será I P EP t EP t
*
98
Óptica de Fourier C. Iemmi
Supondremos ahora, como en general sucede, que los campos son estacionarios, es
decir que no alteran su naturaleza estadística con el tiempo. Así pues, el promedio
temporal es independiente de cualquier momento que se escoja.
Entonces en el primer término de la ecuación anterior haremos un corrimiento del origen
I P K1 E1 t E1* t K 2 E2 t E2* t
2 2
Ahora bien, vamos a definir la función coherencia mutua entre los campos luminosos en
S1 y S2 como
1
12 E1 t E2* t Lim E1 t E2* t dt (24)
2
99
Óptica de Fourier C. Iemmi
I Sn la irradiancia en Sn.
Recordando que para fuentes de ancho de banda no muy grande K1 y K2 son números
imaginarios puros, podemos escribir
I P I1 P I 2 P 2 K1 K2 12
Comentario
Vale la pena aclarar que si la fuente tiene un ancho de banda extenso, entonces de
acuerdo a lo visto
d d
EP t K1 E1 t t1 K 2 E2 t t2 K1 ; K 2
dt dt
Entonces
d d
I P EP t EP* t K1 E1 t t1 E1* t t1
2
dt dt
d d d d
E2 t t2 E2* t t2 2 K1K 2* E1 t t1 E2* t t2
2
K2
dt dt dt dt
Se llega a que
2
I P I1 P I 2 P 2 K1 K 2 2 12
100
Óptica de Fourier C. Iemmi
I P I1 P I 2 P 2 K1 I S1 K 2 I S 2 12
I S 1 I S 2
I1 P I 2 P
12
12
11 0 22 0
(25)
12 12 exp i ARG 12 Luego veremos que es conveniente expresar
2
ARG 12 12 2 12 con 2 r2 r1
Resulta así
CLASE 9
1
2
2
1
1
Lim E1 t E1* t dt Lim E2 t E2* t dt
2 2
101
Óptica de Fourier C. Iemmi
con lo cual
12 12
2
1 12 1
11 0 22 0 11 0 22 0
corresponde a una fuente incoherente. El máximo peso de dicho término ocurre cuando
I P I1 P I 2 P 2 I1 P I 2 P 12 cos 12
2
Recordando que 2 r2 r1 , vemos que este factor tiene en cuenta las
diferencias de camino entre las fuentes S1 y S2 con el punto P, es decir cumple el papel
cumple el papel de 1 2 .
I max I min
Por otra parte vimos que el contraste de las franjas se definía como C
I max I min
102
Óptica de Fourier C. Iemmi
I max e I min se dan cuando cos 1 respectivamente. Con lo cual tenemos que
I1 P I 2 P
C2 12
I1 P I 2 P
término de bias dado por I1 P I 2 P sobre el que se halla un patrón de franjas con
un período determinado por modulado por una variación mucho mas lenta dada por
12 y por 12 .
S1 S1 r1
r1’ r1
S r1’
r2
r2’ r2 S
S2 r2’ S2
Franja central en r1=r2 r1’=r2’ 1=2 Franja central desplazada r1r2 12
103
Óptica de Fourier C. Iemmi
S1
r1’ r1
S
r2
r2’
S2
P
Dado que r1’= r2’ las perturbaciones ópticas en S1 y S2 serán idénticas. Entonces el
estudio del diagrama de franjas estará relacionado con los efectos de coherencia temporal
y la función coherencia mutua será igual a la autocoherencia, esto es
Esta función, así expresada, lo que hace es comparar el campo en un instante t con el
para ver si se superponen en P los trenes de onda secundarios generados a partir del
mismo tren de ondas primario, debemos comparar como eran los campos en S1 y S2 con
esa diferencia de tiempos.
Una situación análoga ocurre en un interferómetro por división de amplitud como el
104
Óptica de Fourier C. Iemmi
E1 t E1* t
11 2
exp i t exp i t
E1
1
1
Lim
exp i t exp i t dt Lim
0
exp i dt
0
fuente.
Analicemos que sucede ahora con la coherencia espacial. Para ello retornemos al
ejemplo de la experiencia de Young donde la fuente luminosa empleada es extensa y de
ancho de banda muy angosto.
S1
S r1
P
r2
S2
evaluación de 12 0 y 12 0 , a partir del contraste de las franjas, nos permitirá tener
idea de la coherencia espacial de la fuente. Estas funciones nos permiten evaluar cuán
parecidos son los campos en S1 y S2 en el mismo instante de tiempo.
105
Óptica de Fourier C. Iemmi
1
12 E1 t E t Lim
*
2 1 exp 2 i t d
2
2
0
*
2 2 ' exp 2 i ' t d ' dt
0
1
12 Lim 0 41 ' exp 2 i exp 2 i ' t dt d d '
*
2
2
1 2* '
12 4 Lim exp 2 i ' d d '
0
2
1 2*
12 4 Lim exp 2 i d
0
2
G12
12 4 G12 exp 2 i d
0
106
Óptica de Fourier C. Iemmi
Dado que 11 0 y 22 0 son reales, la fase de 12 será igual a la fase de 12
. Entonces
12 12 exp i 12 2 4 G12 exp 2 i d
0
Luego
12 exp i 12 4 G12 exp 2 i d
0
Ahora bien, la longitud de coherencia, esto es el largo del tren de ondas, viene dado por
r1 r2 1
coherencia entonces 1 . En consecuencia
C
12 exp i 12 4 G12 d 12 0 12 0 exp i 12 0
0
107
Óptica de Fourier C. Iemmi
12 0 J12
12 12 0 factor complejo de coherencia
11 0 22 0 J11 J 22 2
1 1
2
Este teorema es muy importante y permite establecer una relación entre el tamaño de la
fuente utilizada para iluminar un interferómetro y el contraste de las franjas de
interferencia obtenidas.
Supongamos que tenemos una fuente extensa y cuasi-monocromática que ilumina una
plantalla sobre la que se ubican los pinholes que actúan como fuentes secundarias S1 y
S2 en una experiencia de Young. Lo que queremos establecer es qué área sobre dicha
pantalla sería iluminada coherentemente de modo que si se ubicasen los pinholes dentro
de ella producirían sobre un plano de observación franjas de interferencia con buen
contraste.
Para analizar esta situación veamos el siguiente dibujo y consideremos que la fuente
está conformada por emisores elementales de área dm cuyas dimensiones son tales que
d m
Y
dm Rm2 S2
X
Rm1
S1
108
Óptica de Fourier C. Iemmi
Vamos a denotar Em1 t al campo en S1 producido por el elemento dm, así resulta que
Dado que estamos interesados en estudiar cuál es el efecto que produce la extensión
espacial de la fuente sobre el contraste de las franjas, analizaremos la función intensidad
mutua J12
tanto el campo Em t irradiado por cada uno de ellos corresponderá al de una onda
esférica. Recordemos que una onda de este tipo la podemos describir matemáticamente
A 2 2
como exp i k r t , o teniendo en cuenta que k ; 2 y
r c
que podemos agrupar la amplitud y la fase inicial en un cierto A' , tomará la forma
A' r
exp i 2 t . Con lo cual podemos escribir
r c
R n
exp i 2 t m1
R n c
Em1 t Am t m1
c R m1
R m1 n
tiempo que le demanda en llegar a la radiación desde ese elemento a S1 es por lo
c
109
Óptica de Fourier C. Iemmi
R n
tanto la amplitud y fase inicial están contempladas en Am t m1 , luego a S1 llega
c
esa amplitud dividida por R m1 y con la fase adicional exp i k R m1 t .
Reemplacemos estas expresiones en el promedio anterior
exp i 2 R m1 R m 2
n
R n R n
J12
c
Am t m1 Am* t m 2
m c c R m1 R m 2
Dado que los campos son estacionarios podemos efectuar un corrimiento en el tiempo sin
afectar la estadística, por lo que podemos escribir
exp i 2 R m1 R m 2
n
n
J12 Am t Am* t R m 2 R m1
c
m c R m1 R m 2
n
Como estamos bajo la hipótesis de campos cuasi-monocromáticos ( Rm 2 Rm1 t ),c
c
cada elemento de la fuente emite trenes de onda lo suficientemente largos como para que
luz proveniente del mismo tren llegue a las fuentes secundarias S1 y S2. Quiere decir
entonces que la luz proveniente de cada elemento y que llega a S1 y S2 tiene la misma
fase inicial. En consecuencia podemos escribir
exp i 2 R m1 R m 2
n
J12 Am t Am* t c
m R m1 R m 2
Ahora debemos pasar al continuo y para ello vamos a definir la densidad de irradiancia de
Am t Am* t
la fuente como I m ,m . Con lo cual tenemos que
d m
110
Óptica de Fourier C. Iemmi
exp i 2 R m1 R m 2
n
J12 I m ,m c d
m
m R m1 R m 2
Ahora podemos hacer tender el área del elemento a cero y reemplazar la sumatoria por
una integral, de modo que
exp i 2
n
c
R1 R 2
d d
R1 R1 ,
J12 I , con
R1 R 2 R 2 R 2 ,
Veamos ahora que aproximaciones geométricas podemos hacer en esta integral. Acá, tal
como hicimos en la integral de difracción, se deberán tomar distintas aproximaciones para
la amplitud y para la fase.
2 2
paraxial R 2 en consecuencia podemos escribir
x1,2 y1,2
2
x y1,2
2 2
1 1
R 2
1,2
R1,2 En la amplitud aproximaremos
2R R1 R 2 R
En la fase tomaremos
x y1 x2 2 y2 2
2 2
R1 R 2 R 1 R
2R 2R
x 2
1 y12 x22 y22 2 x1 x2 2 y1 y2
2R
111
Óptica de Fourier C. Iemmi
Así tenemos que podemos escribir el factor complejo de coherencia (cuyo módulo
gobernaba el contraste de las franjas de interferencia para el caso cuasi-monocromático)
exp i
n 2
Rc
x1 y1
2
x2
2
y 2
2
12
R 2 I1 I 2
n
I , exp 2i x1 x2 y1 y2 d d
Rc
I , I ,
I1 J11 d d d d
R12 R2
I , I ,
I 2 J 22 d d d d
R22 R2
Obtenemos una expresión final para el factor complejo de coherencia dada por
k
x12 y12 x22 y22
f X 1,2, fY 1,2
exp i I ,
2R
12
I , d d (27)
f X1,2
x1 x2 ; fY1,2
y1 y2
donde
R R
Luego de esta deducción encontramos un resultado muy importante: para fuentes cuasi-
monocromáticas el contraste de las franjas de interferencia es igual a la transformada de
Fourier normalizada de la distribución de intensidad de la fuente.
112
Óptica de Fourier C. Iemmi
Un caso particular en donde vimos esta propiedad con anterioridad fue cuando
analizamos la experiencia de Young con una fuente extensa lineal, al comienzo de este
capítulo. Allí habíamos encontrado que el contraste era igual a la función sinc
Veamos un ejemplo. Tomemos una fuente circular, uniforme de radio h.
I0 2 2 h
I ,
0 2 2 h
Entonces
hk
2 J1
R
12 x1 x2 y1 y2
2 2
con
hk
R
Vemos que es el radio de un círculo tomando como origen uno de los pinholes. Vamos
a graficar la función que representa a 12
1 |12|
0.88
1 kh /R
Si suponemos, como se hace habitualmente, que una buena coherencia corresponde a un
contraste de 0.88 o mejor, entonces
k h R R R
1 0.16
R kh 2 h h
h
Ahora bien sin 0.16
R sin
113
Óptica de Fourier C. Iemmi
Así el área circular máxima que puede ser iluminada coherentemente por una fuente
circular tiene un radio max 0.16 O sea, existe un área para la cual la
sin
iluminación es coherente por más grande que sea la fuente. Es notable como en este
caso el contraste luego de pasar por un mínimo vuelve a crecer a pesar de aumentar el
tamaño de la fuente. Esta experiencia fue propuesta por Thompson y Wolf (JOSA 47, 895,
1957) y luego repetida numerosas veces. A continuación se muestran los resultados
obtenidos para una fuente extensa con forma de ranura cuyo ancho varía de 0 a 1500 m.
Esta fuente ilumina un par de ranuras de 40 m de ancho separadas 250 m entre sí y
ubicadas a 28 cm de dicha fuente. En el grafico se ilustra el contraste teórico y los puntos
medidos experimentalmente. Las fotografías corresponden a los casos de 400, 700 y 900
m respectivamente. Se puede distinguir la oscilación en el módulo del contraste pero
además es notorio el cambio de fase pasando la franja central de máximo, en el caso
correspondiente a 400 m, a mínimo en el caso correspondiente a 900 m, lo que
concuerda con la fase de la función sinc que describe la transformada de Fourier de una
ranura.
114
Óptica de Fourier C. Iemmi
CLASE 10
FÓRMULA DE HOPKINS
R n
exp i 2 t m1,2
R n c
Em1,2 t Am t m1,2
c R m1,2
exp ikR m1,2
Ahora deberemos reemplazar el factor por una función más general
R m1,2
h , , x, y, que tendrá en cuenta la transmisión del medio. Esto es, el campo en x,y
será para una perturbación monocromática
115
Óptica de Fourier C. Iemmi
E x, y E , h , , x, y, d d
exp ikR
h , , x, y, i donde R es la distancia del punto , al x, y y
R
donde hemos despreciado el factor de oblicuidad.
exp ikR m1,2
Así, generalizando, puede intuirse que el factor debe ser reemplazado por
R m1,2
J12 2 I , h , , x1 , y1 , h* , , x2 , y2 , d d
En consecuencia
2
12 I , h , , x , y , h , , x , y , d d
*
I x1 , y1 I x2 , y2
1 1 2 2
U , , x1 , y1 , i I , h , , x1 , y1 ,
U , , x2 , y2 , i I , h , , x2 , y2 ,
Con lo cual
1
12 U , , x , y , U , , x , y , d d
*
I x1 , y1 I x2 , y2
1 1 2 2
J12 U , , x1 , y1 , U , , x2 , y2 , d d
* (28)
116
Óptica de Fourier C. Iemmi
Estas son las fórmulas de Hopkins y podrían interpretarse como el factor complejo de
coherencia y la intensidad mutua entre dos puntos sobre el plano sobre el que arriban
equivalentes ubicadas en los distintos puntos , de que emiten con una frecuencia
.
Y
S2
r2
P2
X
S1
P1
r1
’
117
Óptica de Fourier C. Iemmi
J P1P2 U , , P1 U * , , P2 d d
Ahora bien, dado que U es una perturbación monocromática podemos aplicar para su
propagación de a ’ la ecuación de Huygens – Fresnel, de modo que
U , , P1 U , , x1 , y1
exp i kr1 q dx dy
1 1 1
r1
1
q1 cos nˆ, rˆ1 ´. Análogamente
i
U , , P2 U , , x2 , y2
exp i kr2 q dx2 dy2
2
r2
exp i k r1 r2 *
J P1P2 q1 q2 U , , x1 , y1 U , , x2 , y2 d d dx1 dy1 dx2 dy2
*
r1 r2
exp i k r1 r2 (29)
J P1P2 q1 q2* J S 1S 2 dx1 dy1 dx2 dy2
r1 r2
118
Óptica de Fourier C. Iemmi
exp i k r1 r2
I P I x1 , y1 I x2 , y2 S 1S 2 q1 q2* dx1 dy1 dx2 dy2
r1 r2
Esta integral expresa a la intensidad en un punto P como la suma de las contribuciones
de cada par de puntos S1 y S2 de una superficie . La contribución de cada par de
elementos depende de las intensidades sobre S1 y S2 y es pesada por el factor grado de
coherencia S S
1 2
. Podría pensarse que esta fórmula es el equivalente a Huygens –
exp i k r
el factor q por la apropiada función de transmisión h x, y, P con lo cual
r
SISTEMAS LINEALES
119
Óptica de Fourier C. Iemmi
L a g E1 x, y b g E 2 x, y a L g E1 x, y b L g E 2 x, y
a g S1 x ', y ' b g S 2 x ', y '
g E x, y g x , y x x , y y dx
E o o o o o dyo
g S x ', y ' L g E xo , yo x xo , y yo dxo dyo
g S x ', y ' g x , y L x x , y y dx
E o o o o o dyo
120
Óptica de Fourier C. Iemmi
respuesta del sistema en el punto x ', y ' del espacio de salida para la función
h x ', y ', xo , yo L x xo , y yo
La función h es llamada la respuesta al impulso del sistema. Así la relación entre la señal
de entrada y la de salida será
g S x ', y ' g x , y h x ', y ', x , y dx
E o o o o o dyo
Una subclase de sistemas lineales son los denominados invariantes. En óptica un sistema
sólo depende de las distancias x ' xo ; y ' yo . Para tales sistemas tenemos que
121
Óptica de Fourier C. Iemmi
eje óptico
g S x ', y ' g x , y h x ' x , y ' y dx
E o o o o o dyo (30)
gS
gE
h
esto es GS f x , f y GE f x , f y H f x , f y
(31)
La función H es llamada función transferencia del sistema y da cuenta de los efectos que
este causa en el dominio de las frecuencias.
122
Óptica de Fourier C. Iemmi
Airy centrado en el punto imagen ideal correspondiente. Así el plano imagen estaba
cubierto por una distribución de diagramas de Airy más o menos superpuestos.
Claramente esta línea de razonamiento es la que se sigue cuando escribimos la función
de entrada como una suma de funciones y encontrábamos como actuaba el sistema
sobre cada una de ellas.
Alternativamente Ernst Abbe (1840 – 1905) en 1873 propuso que un objeto actuaba como
una superposición de redes de difracción de distintas frecuencias y orientaciones,
emitiendo luz en las direcciones correspondientes a sus órdenes. Si el sistema formador
de imágenes no tiene un tamaño infinito entonces no podrá captar toda la luz que emerge
del objeto. Actúa así como un filtro pasa-bajos que rechaza frecuencias espaciales
mayores que un valor dado. En este caso la imagen no corresponde exactamente al
objeto sino que más bien se relacionará con un objeto ficticio cuya transformada de
Fourier coincida con la brindada por la lente. Los órdenes difractados presentes en el
plano transformado se pueden pensar como fuentes secundarias con amplitudes, fases y
posiciones necesarias como para sintetizar la imagen.
Esta visión coincide con la que hemos estudiado anteriormente y con lo expresado en la
H f x , f y se están teniendo en cuenta los efectos del sistema sobre las frecuencias
123
Óptica de Fourier C. Iemmi
o L i
xo xL xi
Eo EL- EL+ Ei
z
do di
El objeto está iluminado por una onda plana, monocromática y la distribución de campo
se forma una imagen del mismo con una distribución de campo Ei xi , yi . Dado que el
fenómeno de propagación de ondas es lineal podemos escribir:
Ei xi , yi h x , y , x , y E x , y dx
i i o o o o o o dyo
expresión para h xi , yi , xo , yo .
Para hallar esta función tomemos como objeto al emisor puntual ubicado en xo , yo . De
ese punto emergerán ondas esféricas que van a incidir sobre la lente luego, en la
aproximación paraxial, el campo justo delante de la misma será
1 ik 2
EL xL , yL xL xo yL yo
2
exp
i do 2d o
Al pasar la lente tendremos
ik 2
EL xL , yL EL xL , yL P xL , yL exp xL yL2
2f
124
Óptica de Fourier C. Iemmi
ik 2 2 ik 2 2
h xi , yi , xo , yo
1
2 di d o
exp xi yi exp xo yo
i
2 d o
2 d
i k 1 1 1 2
P x , y exp
L L 2 do di f L L
x y 2
x x y y
exp i k o i xL o i yL dxL dyL
d o di d o di
ik 2
Veamos cómo se puede simplificar esta expresión. El factor exp xi yi2 sobre
2d i
el plano imagen en general podrá eliminarse dado que todos los detectores (ojo, films,
desaparece.
ik 2
Por el contrario el factor exp xo yo2 depende de las coordenadas del punto
2d o
objeto de modo que cuando debamos resolver la integral Ei hE o dxo dyo
deberemos integrar sobre esas variables. Sin embargo, tal como hemos dicho con
anterioridad, la formación de imágenes es un fenómeno local y por lo tanto en la
125
Óptica de Fourier C. Iemmi
Luego si en un entorno de esa región la fase no cambia más que una fracción de radián
ik 2 2
i k xi2 yi2
se puede aproximar exp o o
x y exp 2 . Ahora que dicho
M 2
2d o 2d o M
factor está escrito en términos de las coordenadas de la imagen podemos aplicar el
mismo argumento que hemos empleado anteriormente para despreciar el factor
previamente analizado.
Por último debemos tener en cuenta que si estamos en el plano imagen las distancias do
1 1 1
y di cumplen con la ecuación gaussiana de las lentes . Así la expresión para
d o di f
h xi , yi , xo , yo queda
1
xo xi yo yi
h xi , yi , xo , yo 2 L L do di L do di yL dxL dyL
P x , y exp i k x
di d o
di
Debemos considerar que el aumento M también puede escribirse como M con lo
do
cual obtenemos
1 i 2
h xi , yi , xo , yo P x , y exp
x Mx x y My y
dxL dyL
2 di d o L L
d i
i o L i o L
126
Óptica de Fourier C. Iemmi
pupila P xL , yL tienda a 1.
h xi , yi , xo , yo
M exp i 2 x Mx x y My y
i o L i o L dxL dyL
1 xi y
M xi Mxo , yi Myo xo , i yo
M M M
1 xi yi 1 xi yi
Ei xi , yi M M x , M y o o Eo xo , yo dxo dyo M Eo M , M
1 x y
invertido EGi xi , yi Eo i , i .
M M M
Habiendo hecho este paso intermedio, continuemos con nuestro análisis. Para ello,
127
Óptica de Fourier C. Iemmi
h xi , yi , xo , yo M P d x , d y exp i 2 x x x y y y
i L i L i o L i o L dxL dyL
De esta forma logramos expresar a la función respuesta al impulso de forma tal que sea
1 xo yo
Ei xi , yi h x x , y y M E
i o i o o M , M dxo dyo
Ei xi , yi h xi , yi EGi xi , yi
(32)
Para la óptica geométrica cada punto objeto brinda un punto imagen perfecto. Al tener
una lente de tamaño finito, existe difracción que ocasiona que ese punto imagen deba ser
reemplazado por un patrón más o menos extenso que corresponderá a la figura de
difracción de una fuente puntual a través de la lente. Esto es, tomamos cada punto objeto
por separado y vemos cómo lo afecta la lente, luego reemplazamos el punto imagen ideal
por esa respuesta individual. Claramente esta línea de razonamiento es la planteada por
Rayleigh.
128
Óptica de Fourier C. Iemmi
CLASE 11
exp i k f d o k x 2f y 2f d
E xf , yf exp i 1
o
i f 2f f
2
. Eo xo , yo exp i x f xo y f yo dxo dyo
f
f X f , fY f Eo xo , yo
o L f i
xo xL xf xi
f z
do di
Para calcular el campo sobre el plano imagen debemos hacer una propagación libre
desde el plano focal
129
Óptica de Fourier C. Iemmi
exp i k f d o k x 2f y 2f d
.exp i
k
xi x f yi y f
exp i 1
o
di f i f 2f f
2
Eo xo , yo exp i x f xo y f yo dxo dyo dx f dy f
f
1 1 1
Teniendo en cuenta que en formación de imágenes se cumple que , resulta
d o di f
xf yf
Pero fX f ; fY f con lo cual dx f f df X f ; dy f f dfYf .
f f
Por otra parte si en el factor exp i
k
di f
xi x f yi y f multiplicamos y dividimos el
2 f
exponente por f tenemos que exp i
di f
xi f X f yi fY f . Recordando
1 1 1 f di 1 do f
nuevamente que entonces do .
d o di f di f M di di f
2
Luego podemos escribir la exponencial como exp i
M
xi f X f yi fY f que al
reemplazarla en la integral da
130
Óptica de Fourier C. Iemmi
Vemos que si del plano objeto al plano focal de la lente obtenemos la transformada de
Fourier de la entrada, del plano focal al plano imagen se origina la antitransformada y
obtenemos la imagen
k xi2 yi2 1 x y
Ei xi , yi exp i k d o di exp i Eo i , i
2 di f M M M
Hasta el momento tratamos casos en los que el sistema óptico era iluminado por fuentes
puntuales, monocromáticas. Analicemos ahora qué sucede si utilizamos fuentes que no lo
son.
Cuando estudiamos la propagación de la intensidad mutua entre dos puntos (S1,S2) y
(P1,P2) a través de sistemas ópticos, habíamos obtenido que:
131
Óptica de Fourier C. Iemmi
forma
I P hx , y , x , y h x
o o i i
*
o ', yo ', xi , yi J xo , yo , xo ', yo ' dxo dyo dxo ' dyo '
Cabe aclarar que los límites de las integrales son ya que estamos realizando el
proceso de integración sobre todo un plano y eventualmente los límites del objeto estarán
tenidos en cuenta en su factor de transmisión. Vamos a suponer que el sistema es
isoplanático, luego tenemos que
I P hx x , y y h x x
i o i o
*
i o ', yi yo ' J xo , yo , xo ', yo ' dxo dyo dxo ' dyo '
Ahora bien, sabemos que J xo , yo , xo ', yo ' Eo xo , yo Eo* xo ', yo ' donde Eo es el
campo a la salida del objeto. A este último podemos escribirlo como el campo que llega al
objeto por el factor de transmisión del mismo, esto es:
132
Óptica de Fourier C. Iemmi
Ahora que tenemos esto expresado en términos del factor complejo de coherencia vamos
a analizar la forma que toma la ecuación de la intensidad en un punto imagen para los
casos límites.
Si la fuente de iluminación es espacialmente coherente, esto es, tenemos una fuente
puntual, entonces de acuerdo a lo que habíamos visto cuando estudiamos la teoría de la
t xo , yo t * xo ', yo ' J xo , yo J xo ', yo ' dxo dyo dxo ' dyo '
h x x , y y t x , y
i o i o o o J xo , yo dxo dyo
h x x ', yi yo ' t * xo ', yo ' J xo ', yo ' dxo ' dyo '
*
i o
2
h x x , y y t x , y xo , yo dxo dyo
i o i o o o J
Recordemos que Eo xo , yo t xo , yo J xo , yo
luego
2
I P h xi xo , yi yo E xo , yo dxo dyo h xi , yi E xi , yi
2
Ei xi , yi h xi , yi EGi xi , yi
I P Ei xi , yi Ei* xi , yi Ei xi , yi h xi , yi EGi xi , yi
2 2
133
Óptica de Fourier C. Iemmi
con el otro con lo cual xo , yo , xo ', yo ' xo xo ', yo yo ' . Luego la intensidad
I P hx x , y y t x , y J xo , yo dxo dyo
2 2
i o i o o o
hx x , y y E xo , yo dxo dyo h xi , yi I xi , yi
2 2 2
i o i o
I P h xi , yi I Gi xi , yi
2
I P hx x , y y h x x
i o i o
*
i o ', yi yo '
Supongamos que el sistema óptico con el que trabajamos no está formado por una sola
lente delgada sino por varias, incluso algunas pueden ser gruesas, otras divergentes, etc.
Vamos a pensar dicho sistema como una caja negra cuyas propiedades pueden
especificarse completamente a partir de las propiedades terminales del conjunto, esto es,
conociendo la pupila de entrada (apertura real o efectiva que limita el cono de luz entrante
al sistema) y la pupila de salida (apertura real o efectiva que limita el cono de salida).
134
Óptica de Fourier C. Iemmi
Yo caja
negra
Xo
objeto
pupila de Yi
salida
o
imagen
pupila de
entrada Xi
i
Ei xi , yi h xi , yi EGi xi , yi
En tanto que cuando era iluminado incoherentemente su respuesta era lineal en
intensidad
I xi , yi h xi , yi I Gi xi , yi
2
135
Óptica de Fourier C. Iemmi
h xi , yi P d x i P , di yP exp i 2 xi xP yi yP dxP dyP
Dijimos que cuando el objeto está iluminado en forma coherente el campo a la salida es
E xi , yi h x x , y y E x , y dx
i o i o G o o o dyo
Dado el tamaño finito de las lentes, las frecuencias espaciales en el objeto no estarán
todas presentes en la imagen. Claramente la imagen ideal que se obtendría si no hubiese
difracción sí contendría todas las frecuencias espaciales y su espectro vendría dado por
f , f E x , y exp i 2 f x
G x y G o o x o f y yo dxo dyo
f , f E x , y exp i 2 f x f y dx dy
i x y i i i x i y i i i
Ahora bien, de acuerdo a las propiedades que vimos para la convolución tendremos que
Ei xi , yi h xi , yi EGi xi , yi i f x , f y H f x , f y G f x , f y
136
Óptica de Fourier C. Iemmi
h xi , yi
H fx , f y h x , y exp i 2 f x f y dx dy
i i x i y i i i
H fx , f y P di xP , di yP P di f x , di f y
f x , f y xi , yi
Si el sistema está limitado por difracción, esto es no existen aberraciones, la función
pupila es real y sus valores son 0 ó 1. Luego en estos tipos de sistemas y bajo iluminación
coherente las frecuencias espaciales que pasan no sufren atenuación o distorsión de
fase, sólo existe una frecuencia de corte que impone un límite a la resolución del sistema
óptico. Frecuencias superiores directamente no pasan.
Veamos el ejemplo de la pupila circular:
x2 y 2
P xP , yP circ P P
donde R es el radio de la pupila de salida
R
xP d i x P
de acuerdo al cambio de variables hecho anteriormente tenemos que
yP di yP
d x2 y 2
Luego la función pupila toma la forma P di xP , di yP circ con lo
i P P
R
cual la función transferencia coherente es
d f x2 f y2
H f x , f y P di f x , di f y circ
i
R
137
Óptica de Fourier C. Iemmi
R
Existe una frecuencia de corte fc que establece cual es el valor de la máxima
di
frecuencia espacial que podrá ser encontrada en la imagen
R
fc = R/di = sin /
di
imagen
pupila de
salida
Al igual que con iluminación coherente, debido al tamaño finito de las lentes, habrá una
degradación de las frecuencias espaciales.
I x , y exp i 2 f x
G o o x o f y yo dxo dyo
G f , f
x y
I x , y dx
G o o o dyo
138
Óptica de Fourier C. Iemmi
I x , y exp i 2 f x
i o o x o f y yo dxo dyo
i f , f
x y
I x , y dx
i o o o dyo
i f , f f , f f , f
x y x y G x y
donde f , f es
x y la función transferencia incoherente, comúnmente llamada OTF
(optical transfer function). Su expresión viene dada por:
h xi xo , yi yo exp i 2 f x xo f y yo dxo dyo
2
f , f
x y
h xi xo , yi yo dxo dyo
2
Hay dos razones por las que se toman los espectros normalizados, una menos importante
es porque resultará más simple matemáticamente para futuros cálculos, la otra es que la
intensidad nunca toma valores negativos y en general siempre está presente un fondo
constante. Dado que la calidad visual de una imagen depende fuertemente del contraste,
esto es de las intensidades relativas, conviene analizar el espectro normalizado por ese
fondo.
Vimos que tanto la función transferencia coherente como la incoherente dependen de
h xi , yi , por lo tanto podemos esperar que exista una relación entre ambas.
Analicemos cuál es. Por un lado tenemos que
H fx , f y h x , y
fx , f y i i ; por el otro hemos encontrado la expresión:
h x , y 2
f x , f y i i
f , f
x y
h x , y 2
pero esto lo podemos escribir como
f x 0, f y 0 i i
139
Óptica de Fourier C. Iemmi
h x , y 2 h x , y h* xi , yi h x , y h* xi , yi
f x , f y i i
fx , f y i i fx , f y i i fx , f y
h xi , yi h* xi , yi h xi , yi * h x , y
fx , f y fx , f y fx , f y f x , f y i i
H fx , f y H fx , f y
f , f
x y
H f x' , f y' df x' df y'
2
donde f , f son variables mudas de integración que denominamos de esta forma sólo
'
x
'
y
fx f
f x f x' ; f y f y' y . Tenemos entonces que
2 2
fx f f f
H f x , f y y H * f x x , f y y df x df y
f , f
2 2 2 2
x y
H f x , f y df x df y
2
0 f , f 1
x y
140
Óptica de Fourier C. Iemmi
0,0 1
f , f f , f
x y
*
x y
CLASE 12
1 m cos 2 f o xo
Ii xi , yi h xi , yi I Gi xi , yi .
2
Io y recordemos que
2
x
1 m cos 2 f o o
1
M
entonces
En este caso tenemos que I Gi
M 2
2
I i xi h xi xo I Gi xo dxo h xo I Gi xi xo dxo
2 2
1
xi xo dxo
1 m fo
I i xi 2 h xo dxo h xo cos 2
2 2
M 2 2 M
1 m
f f
2 2
h xo dx o h x o cos 2 o xi cos 2 o xo dxo
1 2 2 M M
2
M m fo fo
2 o
2
h x sin 2 xi sin 2 xo dxo
M M
141
Óptica de Fourier C. Iemmi
1 m fo fo fo
cos 2
M i x exp[ ]
1 2 2 M M
I i xi 2 h xo dxo
2
M m sin 2 f o x f o exp[ f o ]
2 M i M M
1
1 m fo fo fo
I i xi 2
2
M o o 2 2
h x dx M
cos
2
M i x M
1 m 1 m
I max I min 2 2 m
Co
I max I min 1 m 1 m
2 2
El contraste de la imagen ideal también tendrá este valor. Sin embargo veamos que
sucede con la imagen real
fo fo
1 m M 1 m M
I max I min 2 2 fo
Ci m M
I max I min fo fo
1 m M 1 m M
2 2
fo fo
Luego dado que M 1 el contraste disminuye. Por otra parte da
M
cuenta de un corrimiento o distorsión en la fase.
142
Óptica de Fourier C. Iemmi
di f x di f y * di f x d f
P f x ,fy P fx ,fy i y df x df y
f , f
2 2 2 2
x y
P f x , f y df x df y
2
Si el sistema está libre de aberraciones la función pupila es real y el hecho de que valga 0
di f x di f y di f x di f y
P f x , f y P f x , f y df x df y
x y
2 2 2 2
f , f
P f x , f y df x df y
Esta ecuación representa el cociente entre dos áreas. El numerador corresponde al valor
del área de superposición de dos funciones pupila desplazadas; una de ellas está
di f x di f y
centrada en , y la otra se halla centrada en el punto simétrico opuesto
2 2
143
Óptica de Fourier C. Iemmi
di f x di f y
, . El denominador es simplemente el área total de la pupila.
2 2
Hagamos un gráfico para ilustrar esto
fy fy
fx fx
numerador denominador
f , f .
x y
Es interesante poner el resultado que se obtiene para el caso de la pupila circular de radio
R, dado que en general las lentes son de esa forma. Puede demostrarse que
1 d i f d f di f
2
2 2R
cos i 1 f
f 2R 2R 2 R di
0 afuera
2R
Vemos que la frecuencia de corte incoherente f ci 2 f cc es dos veces la
di
frecuencia de corte coherente que habíamos obtenido con anterioridad. Sin embargo debe
destacarse que en iluminación coherente las frecuencias pasaban o no lo hacían, pero
había un corte brusco. En el caso incoherente, si bien el rango se extendió al doble, existe
una atenuación gradual de las frecuencias. Un gráfico 3D y un corte transversal de esta
OTF puede verse en la siguiente figura
144
Óptica de Fourier C. Iemmi
(f)
(f)
1
1
-fci
fci
0 fci
Cuando trabajábamos con sistemas limitados por difracción suponíamos que a la pupila
de entrada llegaba una onda esférica y de la pupila de salida emergía una onda esférica.
Si ahora tomamos en cuenta las aberraciones del sistema óptico, podemos pensar que la
pupila de salida es iluminada por una onda esférica pero que en dicha apertura hay
ubicada una placa transparente con una distribución de espesores tal que introduce
145
Óptica de Fourier C. Iemmi
pupila generalizada.
Veamos ahora cómo influyen las aberraciones en las funciones transferencia.
df x df y
f , f
área de superposición f x , f y
x y
df x df y
área de total 0,0
di f x di f y di f x di f y
exp ik
W f x
2
, f y
2
W x
f
2
, f y
2
df x df y
área de f x , f y
f , f superposición
x y
df x df y
área de total 0,0
146
Óptica de Fourier C. Iemmi
Un ejemplo relativamente sencillo y que ilustra muy bien lo que veníamos discutiendo es
el caso de la aberración originada por desenfoque en un sistema con pupila cuadrada.
Supongamos entonces que tenemos una pupila cuadrada de lado a y que el error de
1 1 1
desenfoque viene dado por .
d o di f
Habíamos visto, cuando calculamos la función respuesta al impulso, que
i k 1 1 1 2
h xi , yi , xo , yo 2
1
di d o P x P , y P exp x P y 2
P
2 d o di f
i 2
exp xi Mxo xP yi Myo yP dxP dyP
di
1 1 1
Luego en el plano imagen tomábamos 0 . Con lo cual ahora podemos
d o di f
escribir la pupila generalizada que describe este sistema desenfocado como
147
Óptica de Fourier C. Iemmi
Así para una apertura cuadrada el máximo error de fase sobre un borde de la apertura
a2 a2
será k siendo el máximo error de camino óptico. Este último número da
8 8
una idea de cuan desenfocado está el sistema.
tenemos que
fy 8 f x 8 f fy
f x , f y tri 2 ffx
fx
tri sinc 1 sinc
y
1
cc 2 f cc 2 f cc 2 f cc 2 f cc 2 f cc
a
donde f cc .
2 di
1
/ = 0
0.25
1 0.5
0.75
fx /2 fcc
0.2 0.4 0.6 0.8 1
Se ve en esta figura que la función toma valores negativos, lo que indica una inversión del
contraste.
Para el caso de las pupilas circulares, que es bastante más complicado desde el punto de
148
Óptica de Fourier C. Iemmi
Acá también se ve que hay una inversión del contraste, fenómeno que se ilustra en las
siguientes fotografías
Este es un clásico test, que se usa para evaluar sistemas ópticos, en el cual las
frecuencias espaciales aumentan a medida que nos acercamos al centro. En este ejemplo
se puede ver cómo cuando el sistema está desenfocado aparece la mencionada inversión
de contraste.
Aquellos interesados en leer un poco más acerca de cómo se evalúa la MTF de un
sistema óptico puede ir a
http://www.optikos.com/resource_files/pdfs/how_to_measure_mtf.pdf
Hemos visto que la función transferencia incoherente poseía una frecuencia de corte igual
al doble de la coherente, lo que nos tentaría presuponer que este tipo de iluminación da
invariablemente imágenes de mejor calidad que la coherente. Esta conclusión no siempre
es válida ya que estamos comparando frecuencias de corte que no son estrictamente
comparables. La fcc determina la componente de máxima frecuencia en la amplitud de la
imagen, mientras que la fci se refiere a la máxima frecuencia en la intensidad de la misma.
149
Óptica de Fourier C. Iemmi
I i h EG h* EG* h EG h* EG*
H f f EG H f f EG H f G f H f G f
* *
Esto es, obtenemos la autocorrelación del espectro del campo geométrico que está
multiplicado por la función transferencia coherente.
2
Por otra parte, para el caso incoherente habíamos visto que I i h I G con lo cual
I i h h* EG EG* H f H * f G f G* f
H f H f G f G f
2
Las transmitancias de ambos objetos en amplitud vienen dadas por:
150
Óptica de Fourier C. Iemmi
Objeto A: t A xo , yo cos 2 f o xo
Objeto B: tB xo , yo cos 2 fo xo
Supongamos ahora que nuestro sistema óptico está limitado por difracción y tiene una
x
pupila cuadrada de lado a P x p rect p .
a
Comencemos comparando, en una tabla, los dos tipos de iluminación para el objeto A
COHERENTE INCOHERENTE
G GG
1/2 1/2
1/4
-fo fx -2fo fx
fo 2fo
H HH
1
1
fx fx
151
Óptica de Fourier C. Iemmi
[Ii] [Ii]
1/2 1/2
1/4 1/4
<1/4 <1/4
-2fo 2fo fx -2fo 2fo fx
Vemos que para el objeto A resulta una imagen con mejor contraste aquella brindada por
la iluminación coherente.
Estudiemos ahora que sucede con el objeto B. Para ello realicemos nuevamente la tabla
comparativa con los dos tipos de iluminación teniendo en cuenta que la transformada de
medida que se alejan del centro y todas ellas separadas a una distancia 2 f o . Acá nos
COHERENTE INCOHERENTE
G GG
fx fx
-2fcc -fcc fcc 2fcc -2fcc -fcc fcc 2fcc
152
Óptica de Fourier C. Iemmi
H HH
1
H
fx fcc
[Ii]
fx -2fo 2fo fx
En este caso vemos que la iluminación incoherente resulta mejor que la coherente ya que
en esta última se pierde toda la información. Este ejemplo muestra claramente que lo que
denominaríamos como la iluminación más adecuada depende de cada caso particular.
Como vimos este es un criterio basado en el análisis de los espectros de frecuencia,
estudiemos ahora otro criterio posible.
153
Óptica de Fourier C. Iemmi
2
En el caso de iluminación incoherente habíamos visto que I i h I G . En este caso la
2
imagen geométrica corresponde a una delta mientras que h es justamente el diagrama
de Airy (suponemos obviamente una pupila circular). La convolución da como resultado el
diagrama centrado en las dos deltas, luego las dos fuentes puntuales están justamente
resueltas cuando en el plano imagen la separación entre dichas funciones delta sea
1.22 di
asumiendo que D es el diámetro de la lente y di la distancia al plano
D
imagen
I i(x i)
0.19 h
X
xi D
Vamos a hacer un cambio de variables tal que llamamos de esta forma el
di
gráfico queda
I i(x)
-0.61 0.61 X
154
Óptica de Fourier C. Iemmi
Ic X 2 e 2
x 0.61 x 0.61
donde es fase relativa entre los dos puntos. En el siguiente gráfico puede verse que
cuando los puntos están en cuadratura ( = 90º) la resolución es equivalente al caso
incoherente, en cambio para = 0º la resolución empeora y para = 180º mejora
notablemente.
I c(x)
=0º
=180º
=90º
Otros efectos
155
Óptica de Fourier C. Iemmi
Bordes abruptos
Supongamos que parte del objeto viene descripta por una función escalón, en el gráfico
contiguo se ilustra el comportamiento distinto entre un tipo de iluminación y la otra. Con
iluminación coherente aparece un efecto conocido como ringing debido a la difracción en
el borde, además se produce un corrimiento de la sombra geométrica
I(x)
ideal
coherente
incoherente
Speckle
Es la apariencia granulosa de una imagen formada con iluminación altamente coherente.
Su presencia se hace más evidente cuando la imagen corresponde a un objeto difusor, en
este caso podemos pensar que cada punto objeto se comporta como una fuente puntual
que emite ondas esféricas que al ser coherentes con las producidas por los vecinos
interfieren. Luego volveremos sobre este tema.
Salvo que se utilice esta estructura para un fin muy especial, por ejemplo interferometría,
en general resulta un efecto no deseado que degrada la calidad de la imagen.
156
Óptica de Fourier C. Iemmi
157
Óptica de Fourier Claudio Iemmi
157
Óptica de Fourier Claudio Iemmi
o atenúan frecuencias de toda una banda de la señal de entrada, así si el ruido tiene
una frecuencia coincidente, por ejemplo, con parte del espectro del cantante, cuando
este es eliminado también se cambiará el timbre de la voz.
En óptica el primero en modificar intencionalmente las frecuencias espaciales de una
imagen fue Abbe en 1893 y luego Porter en 1906. A continuación describiremos estas
sencillas experiencias que muestran claramente el proceso de formación de imágenes
y su relación con el análisis espectral.
Abbe y Porter utilizaron como objeto una red cuadriculada de alambre y la iluminaron
coherentemente en un dispositivo como el esquematizado en la figura
o F i
158
Óptica de Fourier Claudio Iemmi
Claramente los filtros usados en estas experiencias son los más elementales en donde
la luz pasa o no pasa. Veamos ahora que sucede con un filtro que modifique la fase en
vez de la amplitud.
sea que su transmisión en campo viene dada por t x, y exp i x, y .
Esta condición es necesaria para que exista una relación lineal entre la intensidad en
el plano final y la fase.
Supongamos que en un microscopio común observamos este objeto, iluminándolo con
una fuente coherente. Consideremos además, por simplicidad, que la amplitud del
campo que ilumina al objeto es 1 y que la pupila del sistema es muy grande con lo cual
h . Entonces
I ic h EG 1 i x, y 1 2
2 2
159
Óptica de Fourier Claudio Iemmi
Si, como dijimos las variaciones de fase son pequeñas el término 2 es despreciable y
en consecuencia obtenemos una distribución de intensidad uniforme, como era de
esperar.
Si observamos la ecuación anterior vemos que existen dos términos: uno que
corresponde a la luz no difractada, y otro que contiene la información de la distribución
de fase, y por lo tanto describe a la luz difractada, i x, y . Este último está
desfasado en 90º con respecto al fondo, ya que uno es real puro y el otro imaginario
puro, físicamente esto significa que ambas ondas no interfieren entre sí. Zernike se dio
cuenta de este hecho y propuso agregar un desfasaje adicional de 90º o 270º de
manera que sí pudiesen hacerlo y lograr de esta forma un cambio de intensidad. Así
pensó en introducir, en el plano transformado, un filtro de fase que retardara al orden
cero con respecto a la luz difractada. Una manera de hacerlo es por ejemplo
evaporando una sustancia transparente, de espesor adecuado, sobre un sustrato de
vidrio. Un esquema del dispositivo puede verse en la figura siguiente
luz difractada
vidrio
lente
muestra
plano transformado
i
2
I ic e 2 i x, y i 1 x, y 2 1 2
160
Óptica de Fourier Claudio Iemmi
i 3
2
I ic e 2 i x, y i 1 x, y 2 1 2
que corresponde también a una relación lineal pero con el contraste cambiado.
Cabe destacar que en los microscopios de contraste de fase comerciales la fuente
tiene forma de anillo (en vez de ser puntual) y el filtro correspondiente también ajusta a
esa forma. Debemos recordar que el plano transformado es el conjugado de la fuente
y por lo tanto allí se obtiene una imagen de la misma.
Quienes estén interesados en este tema pueden encontrar información en el sitio
http://www.microscopyu.com/articles/phasecontrast/phasemicroscopy.html
O F xI I
f f f f f yI
focal anterior de la lente LT. El campo a la salida del objeto viene dado por Eo xo , yo
y constituye la señal de entrada. Luego, sobre el plano focal posterior F de la lente
161
Óptica de Fourier Claudio Iemmi
xF yF
transmisión compleja de este filtro viene dada por t F xF , yF C2 H , . Así
f f
la distribución de campo detrás del filtro será proporcional a H . Una segunda lente
LA, cuyo plano focal anterior coincide con el plano focal posterior de LT, brinda la
transformada de esta distribución de campo sobre su plano focal posterior I. Esto es,
2
I xi , yi K E , h x , y d d
i i
Vemos la analogía con el resultado obtenido para formación de imágenes, sólo que
antes h representaba la función respuesta al impulso del sistema óptico y ahora
representa la respuesta al impulso del filtro. Este resultado es razonable ya que
cuando estudiamos formación de imágenes el factor limitante para la luz de entrada
era la pupila y justamente h era la transformada de Fourier de dicha función. Ahora el
objeto que limitará la entrada de luz al sistema será el filtro y en consecuencia h es su
transformada de Fourier. En este razonamiento estamos asumiendo que el sistema
óptico no es el responsable del recorte de frecuencias espaciales, aunque en realidad
la resolución del mismo debe tenerse en cuenta a la hora de diseñar un procesador en
función del ancho de banda con el que se va a trabajar.
Tal como habíamos visto con anterioridad, el objeto algunas veces conviene ubicarlo
detrás de la lente LT de forma tal que el aumento de su transformada pueda regularse
y no esté determinada por la distancia focal de esta lente. De esta forma resulta
menos crítico el ajuste entre y la escala con la que se realizó el filtro.
La configuración 4f es clásica porque, como vimos con anterioridad, no introduce fases
espurias y permite obtener la transformada de Fourier exacta.
También existen otras configuraciones más sencillas en las que se emplean menos
elementos ópticos. Veamos por ejemplo la esquematizada en la siguiente figura
162
Óptica de Fourier Claudio Iemmi
sobre el plano conjugado del mismo I. La distancia 2f a la que se encuentra el objeto
de la lente no es esencial sólo se toma así si se desea tener aumento unitario en el
plano final. Lo único importante es que O y I sean planos conjugados.
yO yF
LC L
xO xF
S
O F xI I
f 2f f f yI
S xO xF
O F I
xI
2f 2f 2f yI
En esta disposición se utiliza una primera lente para colectar la luz de la fuente puntual
y para realizar la transformada de Fourier de la transparencia que se halla después de
ella sobre el plano O. El plano de frecuencias F se ubica justo después de la
segunda lente y es, tal como explicamos anteriormente, el plano conjugado de la
fuente. A la vez esta segunda lente antitransforma el espectro sobre el plano de salida
que es el plano conjugado del plano objeto. Este montaje en particular presenta la
ventaja de que hay menos recorte de frecuencias espaciales que en el 4f y que al igual
que en este no se agrega ninguna fase adicional ya que el factor de fase cuadrática
que se origina por propagación desde el plano objeto al transformado es cancelado
por el factor de transmisión de la lente. La desventaja es que es un montaje más largo.
Al igual que en los otros casos hay montajes similares en donde el objeto no esta justo
detrás de la primera lente ni la disposición es 2f-2f.
163
Óptica de Fourier Claudio Iemmi
O P
fo
fLc
Ejercicio 9: Vamos a suponer que se utiliza un láser de He-Ne ( = 632.8 nm) que
emite en un modo transversal TEM00 de perfil gaussiano cuya distribución de
intensidad viene dada por
2 PT 2r 2
I r exp 2 donde P I r 2 r dr es la potencia total y R el radio
R2 R
T
0
donde la intensidad cae a e-2.
i) Calcule la distribución de intensidad en el plano focal del expansor
ii) Suponiendo que la distancia focal del objetivo es fo = 20 mm y que R, antes de
entrar al mismo es 1 mm, calcule cuanto vale R en el plano focal del expansor. Este es
el diámetro que se elige para el pinhole
164
Óptica de Fourier Claudio Iemmi
En esta sección vamos a ver algunos ejemplos que son “históricos” pero ilustran muy
bien el proceso de filtrado.
En la introducción a este tema hemos visto dos tipos de filtros, uno de amplitud (Abbe-
Porter) y otro de fase (Zernike). En general el procesado de una señal requerirá de
filtros compuestos en los que tanto la amplitud como la fase sean alteradas
simultáneamente. En un principio esto se lograba por la combinación de un filtro de
amplitud y un filtro transparente que, mediante adecuados cambios de espesor,
modificaba la fase. Esta técnica a menudo es bastante limitada ya que es muy difícil
de implementar salvo algunos casos sencillos.
Para ubicarnos históricamente en el tema, cabe destacar que recién a comienzos de
los ´50 el procesado óptico comienza a tener auge. Esto se debe fundamentalmente a
que en esos años se hizo evidente la conexión entre óptica y la teoría de la
información mediante publicaciones tales como “Fourier treatment of optical
processes” [JOSA 42, 127 (1952) y “Optics and Communication Theory” [JOSA 43,
229 (1953) de P.Elias et. al y “Spatial filtering in optics” de E.O’Neil [IRE Trans.
Inform. Theory IT-2, 56 (1956). También en la parte experimental A.Maréchal hizo
grandes aportes principalmente en el mejoramiento de imágenes fotográficas con
filtros compuestos. Algunos años antes, en 1949, D. Gabor [Proc. Roy. Soc. A197,
454, (1949) había propuesto un nuevo proceso a partir del cual era posible reconstruir
una imagen sin lentes. El lo denominó reconstrucción del frente de ondas ya que
permitía modificar la amplitud y fase de un haz de luz; posteriormente fue conocido
como holografía. Este proceso no fue tomado en un principio con mucho entusiasmo
debido a la dificultad de su implementación, pero el advenimiento del láser en 1960
junto con una modificación de la técnica original sugerida por E.Leith y J.Upatnieks
[JOSA 52, 1123 (1962) provocaron una verdadera revolución en este campo. La
síntesis de filtros complejos avanzó rápidamente a partir de que estos pudieron
realizarse de forma holográfica.
Por ahora analizaremos unos ejemplos sencillos y después de estudiar holografía
volveremos sobre el tema para tratar filtros más complicados.
Como vimos anteriormente, el filtro tiene en general una transmisión compleja en
campo dada por:
x y x y x y
tF xF , yF C H F , F H F , F exp i F , F
f f f f f f
165
Óptica de Fourier Claudio Iemmi
x y
que debe cumplir las condiciones H F , F 1 ya que el filtro no puede
f f
xF yF
amplificar, y 0 , 2 .
f f
E x, y
entrada E x, y , esto es . Así en el plano de salida del procesador óptico
x
E x, y
1
f x , f y H f x , f y f ,f
buscamos
x donde x y es la
transformada de Fourier de la entrada y H f x , f y es la transmisión compleja del filtro.
O sea f , f H f , f
x y x y es la distribución de campo a la salida del plano
transformado.
2 i f x f y y
Ahora bien, tenemos que E x, y f , f
x e y
x
df x df y
Luego
E x, y
2 i f x f y y
x
f ,xf 2 i f
y e x
x
df x df y 1
f x , f y 2 i f x
i
debe crecer linealmente con la frecuencia mientras que la fase debe ser e 2 para
3
i
frecuencias positivas y e 2 para negativas. Por lo tanto debemos hacer una
transparencia fotográfica con máxima transmisión para las altas frecuencias y que
disminuya linealmente a medida que estas bajan. El filtro de fase se puede realizar
mediante dos evaporaciones de fluoruro de magnesio de espesor conveniente sobre
un sustrato de vidrio plano.
166
Óptica de Fourier Claudio Iemmi
Más adelante veremos que esta operación se puede llevar a cabo en forma más
sencilla con un filtro holográfico.
El filtrado inverso fue unas de las técnicas introducidas por Maréchal para mejorar
imágenes empobrecidas por defectos en el registro de las mismas. Esto es imágenes
movidas, fuera de foco, etc.
H fx , f y
C
. Notemos que en general esta función será físicamente no
D fx , f y
realizable, por ejemplo para las frecuencias espaciales en las que D f x , f y 0 no
tendremos solución. Para comprender un poco mejor este proceso analicemos los
siguientes gráficos
A
f1 f2 f3 fx
tm B
fx
C
fx
D
fx
167
Óptica de Fourier Claudio Iemmi
Supongamos, por simplicidad, que D f x , f y está compuesta por tramos rectos (ver
Fig. A), esto es, cada tramo será de la forma a f x b . Luego H f x , f y deben ser
hipérbolas de la forma H f x , f y C
a fx b
(ver Fig. B). Los tramos negativos deben
Ahora bien, analicemos la Fig. D. Dado que H max f x , f y 1 entonces de la ecuación
tm Dmax
1 podremos despejar las frecuencias que establecen los intervalos,
a fx b
f , f D f , f H f , f f , f D f , f 1
x y x y x y x y x y
Veamos el ejemplo de una fotografía movida. En este caso cada punto ideal de la
transparencia, en vez de serlo, es un segmento de recta. Luego, si suponemos que el
movimiento se produjo en la dirección X la función de entrada que lo describa
individualmente será
168
Óptica de Fourier Claudio Iemmi
1 1
1 x x x
d x 2 2
0 afuera
x
A sin f x
plano transformado tendremos la función D f x d x 2
que es la
x
fx
2
que debemos filtrar. En las figuras siguientes se esquematizan el espectro del defecto
y el filtro H f x , f y (módulo y fase) que permite compensarlo.
1 D(fx)/A
f1 f2 f3 f4 fx
|H(fx)|
1
tm
f1 f2 f3 f4 fx
H(fx)
f1 f2 f3 f4 fx
169
Óptica de Fourier Claudio Iemmi
CLASE 14
PROCESADORES ÓPTICOS INCOHERENTES
Vamos a analizar dos tipos de procesadores, uno con fuente extensa y otro con fuente
“puntual” (obviamente el término puntual en este caso no tiene las mismas
implicancias que en el caso coherente).
Los sistemas procesadores basados e óptica geométrica utilizan de una u otra forma
la proyección de imágenes (en general se denomina shadow casting o image casting).
Veamos uno en particular que resulta interesante ya que permite implementar la
operación de convolución entre dos funciones. Para ello analicemos la siguiente figura
yF
L1
S T1(-x, -y) T2(x, y) L2 (xs, ys)
xF
F
(-xs, -ys)
f d f
170
Óptica de Fourier Claudio Iemmi
Una fuente extensa S está ubicada en el plano focal anterior de la lente L1.
Inmediatamente después de dicha lente se ubica una transparencia cuya transmitancia
en intensidad es T1(-x, -y). A distancia d de la misma e inmediatamente delante de la
lente L2 se ubica una segunda transparencia con transmitancia T2(x, y). La intensidad
final se registra sobre el plano focal posterior de esta segunda lente.
Para entender como funciona el sistema consideremos un punto particular de la fuente
(-xs, -ys). Los rayos provenientes de ese punto emergerán de L1, y en consecuencia de
T1 , en forma paralela ya que la fuente se encuentra en el plano focal de esta primera
lente. La distribución de intensidad que llega a la segunda transparencia es
d d
proporcional a T1 x xs , y ys . Después de pasar por T2 los rayos se
f f
focalizan en el punto (xs, ys) del plano focal de la segunda lente. Sobre este plano final
la distribución de intensidad será
d d
I xF xs , yF ys K T1 x xs , y ys T2 x, y dx dy
f f
171
Óptica de Fourier Claudio Iemmi
obtener una fuente lo más estrecha y brillante posible. En estos casos el diámetro del
diafragma es mucho mayor que el tamaño de los pinholes utilizados en iluminación
coherente, además su función es diferente ya que se utilizan con el fin de limitar la
extensión de la fuente y no para limpiar el haz.
Veamos un caso donde se aprovecha la policromaticidad de la iluminación empleada.
El pseudocoloreado es una técnica que se emplea para resaltar determinada
información del objeto mediante una codificación de la misma en colores falsos. Por
ejemplo supongamos que queremos resaltar la distribución de frecuencias espaciales
en un objeto. En general el espectro de frecuencias de una imagen posee una
distribución bastante compleja. Las bajas frecuencias corresponden a los detalles
burdos y las altas a los detalles más finos. Si en el plano transformado ubicamos un
filtro tal como el esquematizado en la figura, podremos obtener a la salida una imagen
coloreada según el contenido de frecuencias de cada zona.
Obviamente la imagen anterior es sólo ilustrativa ya que los bordes bien definidos de
los objetos más grandes también poseen cierto contenido de altas frecuencias y en
consecuencia los colores no serán tan puros sino que habrá distinto tipo de mezclas.
Vamos a ver como serían los canales reales y su combinación
172
Óptica de Fourier Claudio Iemmi
S yO yF
diafragma LC G(xo, yo) LT LA
banda
Hn
ancha
xO H1 xF
F xI I
L
f t(xo, yo) O f f f f yI
173
Óptica de Fourier Claudio Iemmi
esta forma sobre el objeto de transmisión t(xo, yo) incide un haz de luz paralelo. Dicho
objeto a su vez se encuentra en contacto con una red de difracción cuya transmisión
f , f , C t xo , yo 1 cos 2 po yo exp 2 i xFf xo yFf yo dxo dyo
xF yF
x f , y f , C T x f , y f C T x f , y f p
F F F F F F xF y F
C3T ,
po
f f
1 2 o
orden 1
orden 0
orden -1
174
Óptica de Fourier Claudio Iemmi
Basta quedarse con uno de los órdenes difractados y bloquear los otros.
Ahora bien, para cada longitud de onda estamos en condiciones de realizar un filtrado
x y
coherente, esto es, para n deberemos construir un filtro H F , F que debemos
f f
centrar en yF n n f po
yF
yF3 yF
H3
yF2
H2
yF1
H1
xF
espaciales p propias del objeto t(xo, yo) . Cada filtro actúa sobre una extensión tal que
donde
max
yFn n f po yF ; min
yFn n f po yF
red objeto
175
Óptica de Fourier Claudio Iemmi
4 po p
n nmax nmin n
po2 p
2
4p
Como en general po p entonces n n
po
Cada operación de filtrado en particular requerirá una coherencia temporal distinta y se
deberá estudiar específicamente. Una forma de regular dicho ancho de banda es
eligiendo en forma adecuada la frecuencia espacial de la red empleada.
Si en vez de un espectro continuo tenemos uno discreto como el del dibujo, puede
calcularse cuál debe ser la frecuencia de la red para tener en cada filtro un ancho de
banda igual al ancho espectral de la línea.
Ahora bien, dado que este filtrado espacial es realizado por separado, en distintas
bandas espectrales, las señales filtradas serán mutuamente incoherentes y se
sumarán en intensidad. Esto es, el procesador parcialmente coherente es capaz de
procesar la señal en campo (como si fuera coherente) para determinadas bandas
espectrales, luego, al ser cada una de ellas incoherente con las demás, las respuestas
individuales se suman en intensidad y es posible eliminar el molesto ruido coherente.
Entonces para cada canal tenemos
Ein xi , yi , n t xi , yi h xi , yi , n
Luego
N
I xi , yi t x , y h x , y ,
2
n i i i i n
n 1
176
Óptica de Fourier Claudio Iemmi
x y
donde h xi , yi , n es la función respuesta al impulso de H F , F
f f
Hasta el momento describimos los filtros y las señales de entrada mediante su función
transmisión pero no mencionamos cómo implementarlos experimentalmente. Esto se
puede llevar a cabo por distintos métodos y utilizando diferentes materiales de registro
o medios de representación. A la hora de elegir los mismos deben analizarse sus
características principales tales como: resolución, eficiencia de difracción, sensibilidad
espectral, energía requerida para el grabado de la información, si constituye un medio
que permite modular sólo la fase o la amplitud, o ambas, si es reutilizable o no, si
permite trabajar a tiempo real, su precio, etc.
De la gran variedad de elementos que permiten modular el frente de ondas acá sólo
analizaremos algunos pocos de ellos y comenzaremos haciéndolo por los materiales
fotográficos que, si bien están siendo cada vez menos empleados, constituyen un
punto de partida para comprender el proceso de registro y representación de la
información.
MATERIALES FOTOGRÁFICOS
Fueron tradicionalmente los más utilizados por su bajo precio, fácil procesado, amplio
rango espectral y sensibilidad. Se pueden elegir dentro de un gran rango de
resoluciones por ejemplo:
Película o placa Resolución (l/mm) Sensibilidad (Asas)
Kodak Tri-X (film) aprox. 50 400
Agfa Copex-pan (film) aprox. 300 25
Holotest 8E-75 (placa y film) > 3000 0.015
Para que tengamos una idea comparativa podemos mencionar que una cámara CCD
tiene aproximadamente una resolución de 100 l/mm. Ahora bien, una cosa es registrar
la imagen y otra poder representarla con igual resolución espacial.
El material fotográfico puede utilizarse tanto para representar la señal de entrada,
como para sintetizar el filtro o como medio de registro de la señal de salida. Algunas
177
Óptica de Fourier Claudio Iemmi
de las desventajas que presentan es que no son reutilizables, no pueden usarse como
medios dinámicos y no son lineales con la exposición.
La película fotográfica virgen está constituida por una gran cantidad de granos de
haluro de plata suspendidos en un soporte de gelatina, el cual a la vez se halla
adherido a un sustrato de acetato (film) o a una de vidrio (placa).
sustrato
capa antihalo
(puede no tenerla)
178
Óptica de Fourier Claudio Iemmi
I transmitida en x, y
T x, y
I incidente en x, y
donde < > significa un promedio local, esto es, un promedio sobre un área grande con
respecto al tamaño de un grano del film pero pequeña comparada con un cambio
significativo de la transmitancia T.
En 1890 F.Hurter y V.Driffield publicaron un trabajo (J.Soc.Chem.Indust. 9, 445, 1890)
1
en donde demuestran que log log T es proporcional a la masa de granos de
T
plata por unidad de área de la transparencia revelada. Definieron a esa cantidad como
la densidad fotográfica D.
T D
hombro
(shoulder) saturación
0.01 2
región lineal
0.1 1 talón
(toe)
tan =
1 0
Log W
densidad del
soporte + velo Do
Energías por debajo del talón no impresionan el negativo. Aquellas por arriba del
hombro tampoco producen un cambio de densidad ya que la película está saturada.
En fotografía convencional se trabaja generalmente en la región lineal, la pendiente de
esta región se conoce como el de la película. Aquellas con un grande (2,3) son las
179
Óptica de Fourier Claudio Iemmi
de alto contraste, mientras que las que tienen un chico (<1) son de bajo contraste.
Esta región lineal corresponde al rango dinámico de la película. Un determinado valor
de está influenciado por tres factores: tipo de emulsión, tipo de revelado y tiempo de
revelado. Cuanto más tiempo es revelado el film, más aumenta el .
Qué tipo de este parámetro se desea elegir depende, en parte, de en qué tipo de
procesador se utilizará la película.
Consideremos a continuación las propiedades de una película bajo distintos tipos de
iluminación.
D n log W Do n log I o . Do
donde n es la pendiente de la curva y n denota que estamos trabajando con un
nuevo film arribe una intensidad I1.Tn luego la transmisión de esta película, que será
un positivo, viene dada por
p . .
Tp K p I1.Tn K p .K n p .I1 p .I o p n K .I o p n
180
Óptica de Fourier Claudio Iemmi
uniforme uniforme
Io I1 I2
negativo
positivo
Tn
Tp
negativo
positivo
revelado revelado
I Io
En la figura anterior se esquematiza el proceso. Ahora vemos que mediante dos pasos
es posible obtener una fotografía con un factor de transmisión lineal con respecto a la
tenga una relación lineal con la intensidad I o a la que se expuso ó con el campo Eo
respectivo.
toma en cuenta las variaciones de fase introducidas por la película. Estas se deben
principalmente a dos motivos. Uno es que el sustrato donde se deposita la gelatina no
es plano desde el punto de vista óptico, el otro es que el espesor de la emulsión no es
constante ya que el mismo varía con la densidad de granos de plata. Así en general
181
Óptica de Fourier Claudio Iemmi
plano óptico
líquido film
n p n p
1 1
n p
Luego t x, y K .I o
2 2
K . Eo .E
2 *
o 2 K´. Eo
Si lo que se busca es que el campo transmitido por t(x,y) guarde una relación lineal
general se elige n 12 ; p 4
t x, y K´ Eo x, y K´I o x, y es posible
2
Cuando estamos interesados en que
obtener esta relación en una sola exposición para un pequeño rango donde
t
1
tb
0.5
0
Wmin Wmax W
Wb
182
Óptica de Fourier Claudio Iemmi
t x, y tb W x, y Wb
t x, y tb ´ I o x, y I b
Obviamente los materiales fotográficos, tal como los hemos estudiado hasta el
momento, nos permiten sintetizar elementos (objeto de entrada, filtros, hologramas
etc.) de amplitud. Esto no quiere decir que no pueda codificarse en amplitud
información de fase, tal como veremos después.
La transmisión de dichos elementos será del tipo
T x, y
transmisión .exp x, y con x, y cte
t x, y
Veremos a continuación que a partir de una película fotográfica es posible obtener un
elemento de fase mediante un proceso llamado blanqueado.
CLASE 15
Fotografías blanqueadas
Los medios de registro cuya transmisión es del tipo t x, y exp i x, y , aparte
de ser utilizados cuando se busca un elemento que tenga esta característica, muchas
veces se emplean por poseer una eficiencia de difracción mucho más alta que los de
amplitud. Esto, por ejemplo, es muy deseable en hologramas, redes, etc. Más
adelante, cuando estudiemos holografía, nos adentraremos un poco más en el tema
183
Óptica de Fourier Claudio Iemmi
de las eficiencias, por ahora sólo veamos en qué consisten estos medios y cómo es
posible transformar una película fotográfica en un modulador de fase.
Existen básicamente dos tipos de blanqueo, los denominados de rehalogenación o
directos y los solventes o inversos.
Veamos primero el blanqueo directo. En este tipo de proceso la película es revelada y
fijada, como se hace habitualmente, luego se sumerge en un baño donde la plata
metálica es reconvertida en una sal de plata transparente con diferente índice de
refracción que la gelatina. Así pues las diferencias de intensidad registradas en la
película se convierten en diferencias de camino óptico. Esto se deben a dos razones:
la primera es el cambio debido a la diferencia de índices de refracción entre la gelatina
y la sal de plata; la segunda es por una modulación en la superficie que se puede
acentuar según el tipo de drogas empleadas. Tal modulación de superficie, en general,
disminuye mucho para frecuencias espaciales mayores a las 100 l/mm. En el caso de
hologramas, donde se registra la luz proveniente de objetos difusores, además es
registrado el speckle provocado por la iluminación coherente. El speckle, usualmente,
es un ruido no deseado con frecuencias espaciales menores a 100 l/mm, que en el
caso de los blanqueos directos es amplificado por sumarse las diferencias de camino
producidas por el cambio de índice de refracción y la modulación en superficie. Un
baño blanqueador típico es el ferrocianuro de potasio [K 3 Fe (CN) 6] que se utiliza por
proveer alta eficiencia y ser relativamente estable con el transcurso del tiempo. Existen
otras drogas más eficientes pero las películas así tratadas luego ennegrecen con el
tiempo.
Hablamos recién del efecto no deseado que produce el speckle, ya que da origen a
una gran cantidad de luz dispersada (scattering) cuando se ilumina la película
fotográfica una vez sometida al proceso total. Para evitar este efecto es conveniente
usar el blanqueo inverso o solvente.
En este proceso la película es expuesta y luego revelada pero no fijada. Así en la
gelatina hay zonas con plata metálica (las que recibieron luz) y zonas con haluro de
plata. Luego del revelado se introduce la película en un baño que disuelve la plata
metálica, así ahora los haluros de plata constituyen la sal transparente con distinto
índice de refracción que la gelatina. En este proceso la modulación en superficie no se
suma al cambio de índices por lo que el efecto que produce la luz dispersada
disminuye notoriamente. Un blanqueado de este tipo es por ejemplo el Kodak R-9 que
es una combinación de dicromato de potasio y ácido sulfúrico.
Para lograr una buena eficiencia con cualquiera de los procesos mencionados es
necesario que las películas estén sobreexpuestas, aproximadamente 10 veces más
184
Óptica de Fourier Claudio Iemmi
que para un elemento de amplitud (para la placa 8E75 que mencionamos antes sería
100 J/cm2). Claramente estos procesos no son lineales
En el siguiente dibujo se esquematizan las etapas de los dos tipos de blanqueo
analizados.
exposición exposición
secado secado
sustrato sustrato
Las gelatinas dicromatadas es otro medio sensible a la luz que una vez revelado
permite obtener un elemento de fase. Poseen algunas características muy deseables
para realizar hologramas de fase: alta resolución, baja absorción y dispersión de la luz,
altísima eficiencia de difracción (casi el 90%), etc. Sin embargo este material no se
vende comercialmente y es muy difícil manejar todos los parámetros que dan lugar a
un registro de alta calidad. Su resolución es prácticamente a nivel molecular,
superando las 5000 l/mm. El proceso físico-químico que durante la exposición produce
el registro de la imagen aún no está completamente comprendido; es más muchos
autores opinan cosas totalmente opuestas. Sin embargo un modelo más o menos
aceptado es el siguiente:
En el proceso de sensibilización un sustrato cubierto con gelatina (similar al de las
placas fotográficas) se sumerge en una solución de dicromato de amonio. Luego se
extrae y se deja secar en una atmósfera saturada con amoníaco para que el dicromato
no cristalice sobre la superficie. Cuando la placa es expuesta en el rango violeta-azul
185
Óptica de Fourier Claudio Iemmi
+6
de longitudes de onda (luego la sensibilidad baja muchísimo) iones de Cr
+3
hexavalentes pasan a ser Cr trivalentes formando encadenamientos cruzados de
moléculas de gelatina. Estos encadenamientos endurecen la gelatina en forma
diferencial siendo más dura en las regiones más expuestas y menos dura en las
menos expuestas. Cabe destacar que dependiendo de la humedad ambiente y de la
temperatura esta reacción puede también producirse en ausencia de luz y por lo tanto
arruinar la placa. Una vez expuesto, el film se lava en agua corriente a temperatura
controlada, esto produce un hinchamiento de la capa de gelatina. Luego se extrae y se
sumerge en isopropanol, lo que produce una rápida deshidratación y revela la imagen
latente. Obviamente esta placa será muy sensible a la humedad por lo que debe
sellarse con una resina que la proteja.
Los primeros hologramas así realizados partían de una gelatina sin endurecer, con lo
cual en el revelado las zonas no expuestas eran eliminadas dando lugar a una
modulación de superficie. Posteriormente se utilizaron gelatinas pre-endurecidas de tal
forma que la exposición y el revelado originaban zonas de diferente dureza y, en
consecuencia con distinto índice de refracción pero ninguna zona era eliminada. Esto
dio lugar a hologramas de mejor calidad.
Dijimos que este producto no se vende y se debe fabricar en el laboratorio. La forma
más sofisticada consiste en preparar uno mismo la capa de gelatina sobre el soporte
deseado y con el espesor adecuado. Otra más sencilla pero no tan efectiva consiste
en pasar una placa fotográfica por fijador, sin revelar, de forma tal que sean eliminados
todos los haluros de plata y quede solamente el soporte con la gelatina, luego esta se
sensibiliza con dicromato de amonio.
Todo este proceso es extremadamente complicado y lleno de secretos que las
publicaciones habitualmente no cuentan, de esta forma resulta poco factible que uno
por si solo llegue a un final feliz. Por otra parte se necesita un láser que emita en el
rango de longitudes de onda antes mencionado y con una potencia adecuada ya que
el material es muy poco sensible, unas mil veces menos que una placa holográfica.
Cabe destacar que el producto final sólo consta de zonas de distinta dureza de
gelatina.
Una técnica alternativa y mucho más sencilla es la que permite combinar la
relativamente alta sensibilidad y amplia respuesta espectral de las placas holográficas,
con la alta eficiencia y relación señal ruido de las gelatinas dicromatadas. En esta
técnica, denominada Silver Halide Sensitized Gelatin, se parte de una placa
holográfica, se expone, se revela y se realza un blanqueado inverso que contiene
dicromato de amonio. Este blanqueador reacciona con la plata metálica y convierte,
como en las gelatinas dicromatadas, iones de Cr +6 en Cr +3 y estos a su vez dan lugar
186
Óptica de Fourier Claudio Iemmi
a zonas de endurecimiento diferencial. Luego la placa se pasa por fijador para eliminar
los haluros de plata que no fueron expuestos y posteriormente se lava y se deshidrata
con isopropílico. Obtenemos así nuevamente un objeto de fase donde las diferencias
de camino son producidas por las diferencias de índice de refracción entre la gelatina
endurecida y la menos endurecida. En este caso la resolución está determinada por la
placa holográfica empleada y no a nivel molecular como en las gelatinas
dicromatadas, sin embargo su resolución es alta ya que la misma puede alcanzar las
3000 l/mm.
PHOTORESIST
photoresist
metal plástico
187
Óptica de Fourier Claudio Iemmi
1.0
Cambio de espesor ()
0.8
0.6
0.4
0.2
0
W (mJ/cm2)
50 100 150
Hasta ahora los materiales analizados no son ni reutilizables ni aptos para trabajar en
tiempo real. Veamos a continuación un ejemplo de material regrabable.
FOTOTERMOPLÁSTICOS
188
Óptica de Fourier Claudio Iemmi
In O2 vidrio
separador
V termoplástico
In O2 fotoconductor
vidrio
+++++++++
1. Primera carga termoplástico
fotoconductor
-------------
2. Exposición
+++++++++
termoplástico
---- ---- fotoconductor
-----
En una tercera instancia, y para lograr una distribución no uniforme de cargas sobre la
superficie del termoplástico, el dispositivo se vuelve a cargar como en el primer paso.
Como resultado, una cantidad adicional de cargas se depositan en las zonas
expuestas, en forma proporcional al potencial reducido y, en consecuencia, a la
exposición recibida.
+++ +++
+++++++++
termoplástico
---- ---- fotoconductor
3. Segunda carga -------------
Así sobre la superficie del termoplástico tenemos una distribución de cargas que será
proporcional a la cantidad de luz recibida.
En un cuarto paso el termoplástico se calienta (aproximadamente 60 ºC), aplicándole
una corriente alterna a las superficies conductoras hasta su ablandamiento. Así el
termoplástico sufre una deformación superficial originada por la presencia no uniforme
de cargas y se afina donde recibió luz. Luego se deja enfriar a temperatura ambiente y
el registro se torna estable.
189
Óptica de Fourier Claudio Iemmi
4. Revelado
termoplástico
fotoconductor
CRISTALES FOTORREFRACTIVOS
Los cristales fotorrefractivos se caracterizan por sufrir una variación del índice de
refracción ante la iluminación con una longitud de onda adecuada. Dependiendo de las
condiciones experimentales (potencia, longitud de onda de los haces de registro y del
haz de lectura), y uso que se les va a dar, se deberá elegir convenientemente el cristal
fotorrefractivo a utilizar. El efecto fotorrefractivo ha sido observado en una variedad de
materiales electro-ópticos (BaTiO3, KNbO3, LiNbO3, GaAs, Bi12SiO20, Bi12GeO20, etc.).
Existen cristales, como el LiNbO3, que mantienen el registro holográfico durante
mucho tiempo (años). Ellos son utilizados como memorias holográficas, dado que los
mismos poseen una gran capacidad para almacenar información (en un cristal de 1
cm3 pueden almacenarse 10000 hologramas). Requieren láseres de gran potencia
para ser grabados y luego de un proceso adecuado pueden leerse sin ser dañados.
Otros, como los pertenecientes a la familia de los selenites (Bi12SiO20, Bi12GeO20) son
mucho más sensibles pero no son aptos para almacenar información. Se los utiliza en
cambio como medio holográfico dinámico. Para ello es preciso que tanto el registro
como el borrado de la información se realicen en tiempos cortos, de manera de no
aumentar el tiempo de procesado de las señales.
Entre los cristales fotorrefractivos utilizados para procesado dinámico el Bi 12SiO20
(BSO) es uno de los más empleados. Veamos a continuación algunas de sus
características y propiedades más relevantes.
Estos cristales son sensibles en el rango verde-azul del espectro (rango en el cual
suelen usarse los haces de escritura) y muy poco sensibles en el rango del rojo (color
que se emplea entonces para la lectura, evitando el borrado de la información
190
Óptica de Fourier Claudio Iemmi
Poseer trampas de cargas, esto es, centros parcialmente ocupados con cargas
eléctricas que pueden fotoexcitarse. Esto permite disponer de niveles donores
y aceptores de electrones y huecos dando lugar a la migración de cargas
cuando estas son fotoexcitadas.
I(x,y)
n(x,y)
BSO
191
Óptica de Fourier Claudio Iemmi
e-
banda de
conducción
h
donores
aceptores
banda de
valencia
192
Óptica de Fourier Claudio Iemmi
Dependiendo del corte del cristal (según el corte cambia su birrefringencia) este
material puede utilizarse para distintos procesos. Uno de los cortes es apropiado para
registrar hologramas en tiempo real mediante mezcla de cuatro ondas. Como
mencionamos, si para escritura y lectura se utiliza la misma longitud de onda, el
proceso se llama mezcla degenerada y si las longitudes de onda son distintas la
mezcla es no degenerada. También cambiando la dirección del campo eléctrico
aplicado puede utilizarse para amplificación óptica mediante la mezcla de dos ondas.
En este caso un haz de alta energía pero sin información le transfiere su potencia a
uno de baja energía pero que contiene la información que se desea amplificar.
En otro de los cortes posibles se utiliza para construir un dispositivo llamado PROM
(Pockels read-out optical modulator) cuyo funcionamiento veremos a continuación.
electrodos
transparentes
luz roja
luz azul
espejo dicroico
BSO
aislante
(parylene)
V
El espesor del cristal oscila entre 0.2 a 1 mm y el mismo está cubierto por una capa
aislante de parylene de 3 a 10 m. A la vez sobre la capa aislante se evaporan
193
Óptica de Fourier Claudio Iemmi
VC Vo exp KW
donde K es una constante positiva y W es la energía de exposición.
La lectura se realiza iluminando uniformemente con luz roja, por ejemplo un láser de
He-Ne, linealmente polarizado. Esta luz prácticamente no altera la conducción del
cristal, esto es no genera cargas libres, ya que el mismo es aproximadamente 200
veces menos sensible al rojo que al azul.
Dijimos al comienzo que el BSO presentaba efecto Pockels, esto es que presentaba
una birrefringencia proporcional al campo eléctrico aplicado. Luego en las zonas
originalmente oscuras habrá máxima diferencia de índices de refracción y en aquellas
previamente iluminadas con luz azul será mínima.
El voltaje aplicado y el espesor del cristal se eligen de forma tal que las zonas oscuras
actúen como láminas de media onda, girando así 90º la polarización del haz rojo. Esto
es si el haz del láser de He-Ne entró polarizado verticalmente, luego de reflejarse en el
espejo dicroico y pasar nuevamente por el cristal, saldrá polarizado horizontalmente de
las zonas oscuras y verticalmente de las claras. Esta codificación de niveles de
iluminación en polarizaciones puede reconvertirse en cambios de intensidad si se hace
pasar el haz emergente por un polaroid.
Si por ejemplo el polaroid se pone con su eje vertical entonces la amplitud del campo a
la salida será
VC
E Eo cos
2 Vo
donde Eo es la amplitud de la luz roja de lectura y se elige Vo = V(/2)
Supongamos que originalmente iluminamos con luz azul la parte superior del cristal,
entonces tendremos la siguiente respuesta
194
Óptica de Fourier Claudio Iemmi
BSO
Eo
polarizador
PROM polarizador
expansor - colimador
Láser He-Ne
ESPEJOS DEFORMABLES
195
Óptica de Fourier Claudio Iemmi
de los primeros el dispositivo más avanzado fue desarrollado por Texas Instruments
(DMD, Digita Mirror Device) y está orientado a sistemas de proyección de imágenes de
alta definición. Su tamaño y resolución depende del modelo pero es de
aproximadamente 1 ó 2 cm2 y contiene más de 1 millón de micro-espejos móviles,
cada uno de 10 a 7.6 m de lado. Su factor de llenado es del 90%. La eficiencia de
difracción es de 86% y puede trabajar a 4KHz con señales binarias o 120 Hz a 8 bits.
Los espejos pueden orientarse, mediante fuerzas electrostáticas, en ángulos que van
de +12º (estado ON) a -12º (estado OFF). En el estado ON el haz de iluminación es
redirigido a una lente de proyección, mientras que en el estado OFF la luz es enviada
fuera de dicha lente. Los niveles de gris se logran regulando cuanto tiempo del ciclo de
barrido cada espejo permanece en un estado ON u OFF.
En la figura se puede observar el dispositivo y un esquema de su funcionamiento
-Vbias
Los espejos están montados sobre una estructura que está conectada, por medio de
dos barras de torsión que actúan como bisagras, a dos puntos de soporte. El espejo
se halla a un potencial negativo –Vbias y debajo del mismo se hallan dos electrodos
de orientación que pueden cargarse alternativamente, uno u otro, a un potencial
positivo +Va de forma que se establece una fuerza atractora que permite deflectar al
196
Óptica de Fourier Claudio Iemmi
espejo deformable
computadora
197
Óptica de Fourier Claudio Iemmi
posicionadores posicionadores
imanes
estructura de
electrodos
bobinas
198
Óptica de Fourier Claudio Iemmi
Información suplementaria sobre este tipo de moduladores del frente de ondas y óptica
adaptativa puede obtenerse, por ejemplo, en
http://www.ti.com/lit/an/dlpa008a/dlpa008a.pdf
http://www.ti.com/lsds/ti/dlp/dlp-getting-started.page
http://www.okotech.com/images/pdfs/catwww3.pdf
http://www.thorlabs.com/NewGroupPage9.cfm?ObjectGroup_ID=3258
píxel
pantalla
199
Óptica de Fourier Claudio Iemmi
El primer polarizador se halla orientado de modo que su eje coincida con la orientación
de las moléculas de la primera cara. Las moléculas de cristal son birrefringentes y van
haciendo que la polarización de la luz incidente rote a medida que avanza en el cristal.
Cuando no hay campo aplicado entre electrodos, al emerger la luz por la cara opuesta
de la celda, rotó aproximadamente 90 grados. El eje del segundo polarizador se halla
cruzado al primero y en consecuencia la intensidad transmitida es máxima. Cuando
hay aplicado un voltaje máximo entre electrodos, las moléculas se orientan paralelas al
campo aplicado, su birrefringencia disminuye y la luz atraviesa la celda sin rotar su
polarización de modo que a la salida tenemos intensidad mínima. Así podemos lograr
toda una gama de atenuaciones de la luz en función del voltaje aplicado.
polarizador
vidrio
ITO +
capa de
alineado
analizador
Sin embargo con esta disposición no es posible alcanzar una modulación desacoplada
de amplitud y fase, aunque para los fines comerciales para los que habitualmente se
usan las pantallas no tiene importancia. Para que un dispositivo de este tipo sea útil en
procesado óptico se debe lograr que actúe ya sea como modulador puro de fase o
puro de amplitud. Esto es posible si además de los polarizadores se utilizan entre ellos
y los vidrios, láminas retardadoras que permitan generar estados de polarización
elíptica adecuados a partir de los cuales se logran obtener dichas modulaciones.
200
Óptica de Fourier Claudio Iemmi
1
Intensidad y Fase (u. de 2 pi)
0.9
Intensidad
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1 Fase
0
0
20
40
60
80
100
120
140
160
180
200
220
240
-0.1
Nivel de gris
Modulación de amplitud
Intensidad y Fase (u de 2 pi)
1
0.9 Fase
Intensidad
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
20
40
60
80
100
120
140
160
180
200
220
240
255
Nivel de gris
Modulación de fase
201
Óptica de Fourier Claudio Iemmi
mediante una LUT se le envía el nivel de gris necesario para obtener el desfasaje
deseado. Es decir tenemos un elemento de fase programable píxel a píxel con el cual
pueden obtenerse respuestas que usualmente no se consiguen con otro tipo de
dispositivos y además puede modificarse a velocidad de video.
Veamos dos ejemplos en los que este dispositivo actúa como lente difractiva.
En el primero de ellos se combinaron varias lentes difractivas de diferente distancia
focal, de manera de mejorar algunas de las características del sistema. Primero
veamos como se obtiene una lente difractiva simple en una LCTV. No es más que la
distribución de fase cuadrática que describe una lente común, calculada módulo 2 pi y
enviada a través de una PC a la LCTV. En este ejemplo vamos a ver como aumentar
la profundidad de enfoque (DOF) mediante una lente compuesta (esto tiene un gran
interés en litografía donde es importante tener una gran DOF para obtener una imagen
estable sobre el photoresist). La lente difractiva simple tiene una distancia focal de 100
cm. Diseñamos una lente compuesta conformada por 33 lentes simples, cuyas
distancias focales iban de 92 a 108 cm con incrementos de 0.5 cm. Se vio que la mejor
manera de obtener esta lente compuesta era eligiendo en forma random píxeles de
cada una de las lentes simples [C.Iemmi, J.Campos, J.Escalera, O.López-Coronado,
R.Gimeno, M.Yzuel Depth of focus increase by multiplexing programmable diffractive
lenses. Opt. Exp. 14, 10207 (2006)].
202
Óptica de Fourier Claudio Iemmi
Lente simple
Z = 0 cm Z = 10cm
Lente compuesta
El plano de mejor enfoque para la lente simple corresponde a Z=0. A Z=+- 10cm la
imagen esta muy deteriorada; vemos en cambio que para la lente compuesta la
imagen mantiene su nitidez en todo el rango para la que fue programada.
Otro ejemplo es aquel en el que mediante un sistema de moduladores se puede
realizar un zoom programable, anamórfico. Mediante la combinación de lentes
cilíndricas con distintas orientaciones y distintas distancias focales, es posible
magnificar, demagnificar y distorsionar la imagen de un objeto. A continuación se
muestran algunos resultados [C.Iemmi, J.Campos, Anamorphic zoom system based on
liquid crystal displays J. Eur. Opt. Soc. – Rapid Pub. 4, 9029 (2009)]
.
203
Óptica de Fourier Claudio Iemmi
204
Óptica de Fourier Claudio Iemmi
PD
L
T2 SS T1
DM
sample
Se muestran los resultados en los que se programaron fases lineales que permiten
desplazar de forma arbitraria el haz a través del plano donde esta el objeto bajo
estudio. Además permite generar barridos lineales (programando una lente cilíndrica)
o estudiar que sucede a distintas profundidades (programando lentes de distintas
distancias focales) [M.Capeluto, C. LaMela, C.Iemmi, M.Marconi, Scanning mechanism
based on a programmable liquid crystal display, Opt. Comm. 232, 107, (2004)].
En los últimos años ha sido desarrollada una nueva tecnología de LCDs , dando lugar
a un dispositivo denominado LCoS (liquid cristal on silicon). Las pantallas LCoS están
basadas en el desarrollo de la tecnología CMOS, la cual permite obtener a un precio
relativamente accesible un dispositivo que posee muy alta resolución, con píxeles de
aproximadamente diez micrones y un factor de llenado superior al 90%. Dado que
estas pantallas actúan por reflexión es posible, en principio, obtener una alta
modulación de fase sobre el haz incidente ya que el mismo efectúa un recorrido de ida
y vuelta dentro del cristal líquido.
205
Óptica de Fourier Claudio Iemmi
Polarizador
Cubierta de vidrio
Electrodo transparente
Lámina de alineado
Cristal líquido
Recubrimiento
CMOS reflectivo
Montura sobre
PCB
A pesar de estas características tan favorables, se encontró que existe un cierto grado
de depolarización de la luz y que este efecto puede llevar a una degradación en la
aplicación de un LCoS como modulador espacial de luz. Un estudio detallado mostró
que el problema está originado en la forma en la que se envía la señal eléctrica. En
estos dispositivos los distintos niveles de gris (o fase) se logran mediante pulsos que
no cambian su voltaje (como lo hacían en las LCTVs por transmisión) sino su serie
temporal. Esto es, se aplican más o menos pulsos para lograra el mismo efecto que
aplicando más o menos voltaje. Esto se traduce en que lo que aparentemente es un
nivel de gris o fase constante, en realidad está fluctuando en el tiempo. Una muestra
de este fenómeno puede verse en las siguientes figuras, donde se grafica la respuesta
temporal obtenidas con una pantalla tipo LCoS marca Philips modelo X97C3A0
controlada por el dispositivo LC-R2500 provisto por Holoeye [A. Lizana, I. Moreno, A.
Marquez, C. Iemmi, E. Fernandez, J. Campos, M. Yzuel, Time fluctuations of the
phase modulation in a liquid crystal on silicon display: Characterization and effects in
diffractive optics, Opt. Exp. 16, 16711 ( 2008)]. Esta pantalla, de 2.46 cm de diagonal,
es de cristal líquido nemático girado 45 grados, posee una resolución de 1024 X 768
pixeles con un factor de llenado del 93% y brinda 256 niveles de gris.
206
Óptica de Fourier Claudio Iemmi
-120
-180
-240
-300
GL=120 GL=211 GL=255
-360
0 5 10 15 20 25 30 35 40
(d) Time (ms)
VIII. HOLOGRAFÍA
INTRODUCCIÓN
En 1948 Dennis Gabor propuso un nuevo método formador de imágenes para el que
no hacían falta lentes. El denominó al mismo reconstrucción del frente de ondas. En
este proceso se deseaba registrar tanto la amplitud como la fase de las ondas
luminosas provenientes del objeto. Como es sabido los medios de registro son
sensible solamente a la intensidad de la onda luminosa por lo que resulta imposible
registrar una fase absoluta. Lo que si puede realizarse es convertir la información de
fase en variaciones de intensidad y así registrar estas últimas. Una técnica standard
que realiza esta tarea es la interferometría, en donde las variaciones de fase
introducidas por un objeto con respecto a un frente de ondas de referencia pueden ser
registradas como variaciones de intensidad. Para lograr esto, obviamente, el objeto
debe ser iluminado en forma coherente. Imaginemos por ejemplo un interferómetro de
Michelson en el que uno de los espejos es reemplazado por el objeto a testear y el
otro provee el haz de referencia
207
Óptica de Fourier Claudio Iemmi
Gabor demostró así que la amplitud y fase de la luz proveniente de un objeto podía
registrarse en una placa fotográfica (aunque esta sólo responde a cambios de
intensidad) y luego reconstruirse convenientemente. Llamó a tal registro holograma
(del griego holo:todo; grama:escritura).
Gabor publicó los primeros trabajos relacionados con el tema en los años 1949 y 1951
y estaban orientados a mejorar la resolución de los microscopios electrónicos.
[D.Gabor, Microscopy by Reconstructed Wavefronts, Proc. Roy. Soc. A197, 454,
(1949) ; D.Gabor, Microscopy by Reconstructed Wavefronts II, Proc. Phys. Soc. B64,
449, (1951)].
Debido a que en ese momento no se disponía de una fuente luminosa altamente
coherente el tema no fue tomado con interés. Por otra parte, como veremos, resultaba
muy difícil observar una imagen así obtenida.
Fue recién a comienzos de los años ´60 que, junto con el surgimiento del láser, dos
investigadores de la Universidad de Michigan, Emmeth Leith y Juris Upatnieks,
encontraron la forma de hacer del holograma algo realmente asombroso. La calidad y
realismo de las imágenes obtenidas pronto llamó la atención de cientos de
investigadores que se volcaron al tema. Surgió así una variedad enorme de
aplicaciones de la técnica (interferometría holográfica, filtrado espacial, reconocimiento
de caracteres, etc) que hicieron que su propiedad más efectista, esto es, la visión 3D
del objeto, fuese una de las menos importantes.
Gabor recibió el premio Nobel de física en 1971 y Leith y Upatnieks una mención
especial.
Analicemos algunas de las características del holograma de un objeto:
- Para su registro se utiliza, entre otros medios, material fotográfico cuya principal
diferencia con el utilizado habitualmente en forma comercial es su alta resolución.
208
Óptica de Fourier Claudio Iemmi
Sin embargo en una fotografía sólo se registra la intensidad proveniente del objeto
mientras que en un holograma registramos intensidad y fase relativa, en realidad
es la fotografía de un diagrama de interferencia.
- La holografía no necesita de lentes formadoras de imágenes por lo tanto no tiene
problemas de enfoque ni profundidad de campo.
- Un holograma no requiere pasarse de positivo a negativo pues sólo sería un
cambio de fase en que el ojo no registra.
- Si registramos el holograma de un objeto difusor y luego lo dividimos, a partir de
cada trozo es posible reconstruir la totalidad del objeto. Esta afirmación requiere
una aclaración: debemos imaginar que es como observar el objeto por una
ventana, si la misma es más chica lo seguiremos viendo pero se reduce el campo
de visión, también veremos que se pierde nitidez. La posibilidad de retener toda la
información en cada una de las partes se debe a que el objeto difusor emite luz en
todas las direcciones, luego cada punto del mismo ilumina aproximadamente a
toda la placa fotográfica, es decir sobre la misma queda registrado el diagrama de
interferencia del haz de referencia con la luz proveniente del objeto.
Onda objeto
Objeto
209
Óptica de Fourier Claudio Iemmi
Onda de referencia
Objeto
Onda objeto
Placa holográfica
I x, y O r O r O r 2 Oo ro cos o r
* 2 2
t x, y tb ´ I o x, y I b
Ib r
2
Tomemos, sólo para simplificar la expresión, que , luego reemplazando la
t x, y tb ´ O O r * O*r
2
Ahora bien, una vez que esta información está almacenada en la placa, debemos
i
analizar cómo extraerla. Para ello utilicemos un haz de reconstrucción R Ro e R
210
Óptica de Fourier Claudio Iemmi
Onda de reconstrucción
Observador
Holograma revelado
E x, y t x, y R tb ' O
2
R 'r O R 'r O R
* *
E1 E2 E3
R r , entonces
E2 ' r *O R ' r O
2
E3 ' r O* R ' r r O* E2 O
R=r
holograma
Obtenemos que E2 a menos de una constante provee el mismo campo que la onda
R r * , entonces
211
Óptica de Fourier Claudio Iemmi
E2 ' r *O R ' r *r *O R = r*
E3 ' r O* R ' r O*
2
E3 O*
holograma
R r*
R=r
holograma holograma
Sin embargo vemos que en este caso la curvatura de ambos haces no se cancelan
sobre el holograma.
Grafiquemos entonces el caso más general de dos fuentes conjugadas S y S’ que
emiten ondas esféricas
S S’
R = r*
holograma
212
Óptica de Fourier Claudio Iemmi
Vemos que sólo en el caso en el que la placa holográfica estuviese ubicada sobre la
línea de puntos que une las dos fuentes sucede lo anteriormente descripto, de manera
de cancelar sus curvaturas sobre el medio de registro. En caso contrario, si la placa
está ubicada tal como hemos visto con anterioridad, entonces r* debe converger
hacia S’, tal como se esquematiza en la figura.
DISPOSICIÓN DE GABOR
to xo , yo to to xo , yo ; to to
Donde to es la transmitancia promedio y to son los apartamientos de dicho
promedio.
xo, yo x, y
O r
objeto zo placa
Debemos recordar que en la época que Gabor diseñó esta experiencia no existían
fuentes altamente coherentes como los láseres, de hecho él utilizó una lámpara de
mercurio filtrada. De esta forma debió minimizar la diferencia de caminos a fin de que
el haz objeto y el de referencia pudiesen interferir. Así iluminó con un haz plano un
objeto cuya transmisión provocaba que una alta proporción de la luz incidente no se
I x, y r O x, y r O x, y r *O x, y r O* x, y
2 2 2
213
Óptica de Fourier Claudio Iemmi
r
2
placa y que tomamos, como lo hicimos anteriormente, I b . Además como r es
una onda plana que incide normalmente entonces r = r* . Luego podemos escribir la
transmisión de la misma como
t x, y tb ´ O O r O*r
2
Para la reconstrucción del holograma se emplea nuevamente una onda plana con
incidencia normal y amplitud R, con lo cual el campo transmitido viene dado por
Rt x, y R tb ´ O
2
´R r O x, y ´R r O x, y
*
El primer término es una onda plana que pasa por la placa casi sin dispersarse. El
zo zo
R
observador
214
Óptica de Fourier Claudio Iemmi
prisma
y
S
z
placa
lente zo
objeto
La luz emergente de la fuente puntual S es colimada por la lente. Una porción del
frente de ondas plano pasa a través del objeto con transmitancia to xo , yo . Otra
porción del frente de ondas pasa por el prisma que lo desvía formando un ángulo
con respecto a la normal a la placa. Así sobre el material de registro inciden dos ondas
coherentes, la que pasó por el objeto y la que deflectó el prisma con dirección
2
k 0, sin , cos . La distribución de campo sobre la placa será
2
E x, y ro exp i sin y O x, y
donde elegimos sobre la placa el plano z = 0. La distribución de intensidad será
215
Óptica de Fourier Claudio Iemmi
2 2
I x, y r O x, y ro exp i sin y O x, y ro exp i sin y O* x, y
2 2
Nuevamente vamos a suponer que trabajamos en una zona del material de registro
donde vale que la transmisión del campo es lineal con la intensidad de registro y
Ib r
2
tomamos . Así el factor de transmisión de la placa será
2 2
t x, y tb ' O x, y ro exp i sin y O x, y ro exp i sin y O* x, y
2
E x, y R t x, y Ro tb ' O x, y
2
' R r exp i 2 sin y O x, y
o o
2
' Ro ro exp i sin y O* x, y
El primer término es directamente transmitido y no aporta ninguna información ya que
R
-
zo zo
216
Óptica de Fourier Claudio Iemmi
Vemos que en este caso las imágenes están separadas angularmente y por lo tanto se
las puede observar en forma independiente. Por otra parte, dado que el fondo
coherente no se superpone con la imagen de interés, estos haces no interfieren y en
consecuencia no importa si el holograma es un positivo o un negativo. Otro detalle a
tener en cuenta es que la frecuencia espacial de la portadora, la cual está determinada
por el ángulo con que incidió el haz de referencia, debe elegirse de forma tal que sea
mucho mayor que las frecuencias espaciales propias del objeto, pues así quedan bien
proporcional a Ro ' O x, y .
2
217
Óptica de Fourier Claudio Iemmi
HOLOGRAMAS PLANOS
objeto x
O
d
do
y z
r
placa
218
Óptica de Fourier Claudio Iemmi
coordenadas xi , yi , zi .
y z
R
placa
Holograma de Fresnel
láser
filtro
espacial
espejo
placa
objeto
219
Óptica de Fourier Claudio Iemmi
Cada punto objeto es emisor de ondas esféricas por lo que la perturbación total que
produce en el plano del holograma viene dada por
ˆ , nˆ exp ik1d dx dy
O x, y, z
i
1
o o o
O x , y , z cos d d
o o
1
con d xo x yo y zo z .
2 2 2 2
exp ik1l
El haz de referencia sobre el holograma será r x, y, z ro donde
l
1
l xr x yr y zr z
2 2 2 2
Para simplificar un poco los cálculos los haremos en dos dimensiones y luego
extenderemos el resultado a tres.
En la aproximación de Fresnel y recordando que la placa se halla en Z = 0 tenemos
que la onda objeto viene dada por
ro xr2 x 2 x xr
r exp ik1 zr
zr 2 z r 2 z r z r
Ro xR2 x 2 x xR
R exp ik2 zR
zR 2 z R 2 z R zR
E x, y t x, y R tb ' O
2
R 'r O R 'r O R
* *
E1 E2 E3
220
Óptica de Fourier Claudio Iemmi
Para calcularlos debemos tener en cuenta que cuando iluminamos el holograma con el
haz de referencia R, es el mismo holograma el que actúa como una abertura
difractora. Así tenemos que
Escribamos esta ecuación de una manera más simplificada, para ello llamamos
k1 k k k xr x x x
a 2 2 1 ; b R o i
2 zr 2 zR 2 zi 2 zo 1 zr 2 zR 1 zo 2 zi
con lo cual podemos escribir
exp ik1 zo
ik1 xo2
E2 xi Cte O xo exp 2 zo exp ia x exp i 2 b x dx dxo
2
zo
O xo xo2
expresión exp ik1 zo es una onda esférica con centro en xo , zo
zo 2 zo
exp ik1 zo
ik1 xo2
E2 xi Cte o 2 zo
2
O x exp exp ia x exp i 2 b x dx dxo
zo
221
Óptica de Fourier Claudio Iemmi
La integral sobre dxo nos daría el campo que produce todo el objeto y la integral sobre
dx es la que contempla el efecto del holograma. Hagamos una analogía con lo visto
para sistemas formadores de imágenes y recordemos la expresión
Ei xi Eo xo h xi , xo dxo
En ese caso para que la imagen fuese lo más parecida posible al objeto, la función
respuesta al impulso del sistema debía aproximarse a una función delta de dirac.
En este caso si queremos que la onda difractada por el holograma sea lo más
parecida posible a la onda objeto, entonces debemos pedir que la integral entre llaves
tienda a una función delta. Luego eso es cierto si se cumple que a = 0. En ese caso
1
Pedir que a = 0 implica que zi . Por otra parte tenemos que
1 2 1 1
zR 1 zo zr
exp ik1 zo
ik1 xo2
E2 xi Cte O
o 2 zo b dxo
x exp
zo
Esta integral vale 0 siempre que b 0. Entonces obtenemos una onda difractada sólo
x x x
cuando b = 0 lo que implica que xi 2 o r R zi
1 zo zr zR
La ecuación hallada para zi es una condición de focalización análoga a la ecuación de
1 1 1
las lentes .
d o di f
Extendiendo los resultados hallados a tres dimensiones tenemos que
x x x
xi 2 o r R zi
1 zo zr z R
y y y
yi 2 o r R zi
1 zo zr z R
1
1 1 1
zi 2
z R 1 zo zr
Estas serán las coordenadas de la imagen primaria P xi , yi , zi .
222
Óptica de Fourier Claudio Iemmi
conjugada C xi , yi , zi
x x x
xi 2 o r R zi
1 zo zr z R
y y y
yi 2 o r R zi
1 zo zr zR
1
1 1 1
zi 2
z R 1 zo zr
Vemos, a partir de estas ecuaciones, que cuál imagen es real y cuál es virtual
depende de zi y en consecuencia de zo , zr y zR. Si zi > 0 la imagen es real, si zi < 0
1
M
z z donde el signo + es para P y el – para C
1 o 1 o
zr 2 zR
223
Óptica de Fourier Claudio Iemmi
Si 1 2 ; zr zR M 1 para P.
Si 1 2 ; zr zR M 1 para C.
Si 1 2 ; zr zR M 1 para P y C.
i
El aumento angular está dado por la relación M A . En la aproximación paraxial
o
xi x
podemos tomar tan i i ; tan o o o con lo cual podemos escribir
zi zo
xi zo
MA 2
xo zi 1
Vemos de esta relación porqué Gabor estaba interesado en el método holográfico para
aumentar la resolución de un microscopio electrónico. La relación entre la longitud de
onda de un electrón y la longitud de onda del visible es aproximadamente 100000 X.
Holograma de Fraunhofer
xo2 yo2
comparación de su distancia a la placa zo estamos en condiciones de
difracción de Fraunhofer. Una situación típica donde esto sucede es, por ejemplo, en
el holograma de un grupo de pequeñas partículas tales como las encontradas en una
suspensión o las producidas por un aerosol.
La ventaja de tener un holograma de una distribución de este tipo en vez de una
fotografía es que se puede ir estudiando por distintos planos a zo = cte.
Como sabemos la figura de difracción de Fraunhofer registrada en el holograma
corresponde a la transformada de Fourier del objeto. Existe otro tipo de holograma en
el que también se registra la transformada de Fourier en vez del objeto mismo y es el
que veremos a continuación.
Holograma de Fourier
Se define como holograma tipo Fourier a aquel en el que se registra sobre la placa
holográfica, la interferencia de dos ondas tales que corresponden a las transformadas
de Fourier del objeto y de la referencia, respectivamente. Normalmente esto implica
que el objeto sea bidimensional o de espesor limitado. Este tipo de hologramas, como
224
Óptica de Fourier Claudio Iemmi
xo x
referencia
a
r
objeto
placa
f f
2
O x, y to xo , yo ; r x, y ro exp i sin r x
fx f y
I x, y O x, y r x, y O x, y r x, y
*
2
2
to xo , yo ro2 ro to xo , yo exp i sin r x
fx f y fx f y
2
ro *
to xo , yo exp i sin r x
fx f y
225
Óptica de Fourier Claudio Iemmi
xi
placa
f f
2
E x, y Ro t x, y Ro tb ' to xo , yo
fx f y
2
' Ro ro to xo , yo exp i sin r x
fx f y
2
' Ro ro * to xo , yo exp i sin r x
fx f y
xi yi Rotb xi , yi
1
xi yi ' Ro f f to xo , yo . f f to xo , yo ' Ro xi yi f f to xo , yo
1 * 1
x y x y x y
x1 y * to xo , yo ' Ro to xi , yi to xi , yi
i i
fx f y
2
xi yi ' Ro ro f f to xo , yo exp i sin r x
1
x y
2
' Ro ro x1 y to xo , yo x1 y exp i sin r x
i i
fx f y i i
226
Óptica de Fourier Claudio Iemmi
a
Si tenemos en cuenta que sin r r y la definición de frecuencia espacial
f
tenemos que
2 2 a
exp i sin rx exp i x exp i 2 f x a
f
x1 y exp i 2 f x a xi a, yi
i i
Con lo cual
2
xi yi ' Ro ro f f to xo , yo exp i sin r x
1
x y
' Ro ro to xi , yi xi a, yi ' Ro ro to xi a, yi
2
xi yi ' Ro ro f f to xo , yo exp i sin r x ' Ro ro to* xi a, yi
1 *
x y
E xi , yi A xi , yi Bto xi , yi to xi , yi C to xi a, yi Dto* xi a, yi
con A, B, C, D constantes.
Las imágenes de interés, como siempre, son la primaria y la conjugada representadas
por el tercer y cuarto término respectivamente. Ambas imágenes son reales y están
sobre el mismo plano ya que el objeto y el haz de referencia también lo están.
xo primaria
a
yo yi
a
conjugada
objeto en el
plano de entrada xi
imágenes en el
plano de salida
227
Óptica de Fourier Claudio Iemmi
Luego veremos un dispositivo similar que permite sintetizar filtros para emplear en el
procesado de señales. Cabe destacar que acá estamos viendo los distintos tipos de
hologramas, qué es lo que se registra y cuál es la respuesta pero no entramos en el
detalle de su síntesis. En el caso particular de la generación de filtros aplicados al
procesado de señales, si bien se usan hologramas tipo Fourier, los mismos no se
sintetizan ópticamente sino que son generados por computadora. Es decir los
conceptos físicos involucrados son los mismos pero el proceso es numérico, luego el
holograma así calculado se implementa en un medio material que permite utilizarlo en
un procesador óptico. Luego veremos cómo hacerlo.
Los hologramas que hemos visto hasta ahora necesitan que tanto para su registro
como para su reconstrucción se emplee una fuente coherente. En el caso del registro,
esto es evidente, dado que se busca que la onda objeto interfiera con la utilizada como
referencia. En el proceso de reconstrucción también debe emplearse una fuente
puntual cuasimonocromática ya que el holograma, al actuar como una red de
difracción, genera una imagen para cada longitud de onda con un tamaño y posición
distintos.
La necesidad de utilizar iluminación coherente para su reconstrucción es una
característica que muchas veces limita el empleo de hologramas para distintos
propósitos. Por ejemplo, todos los hologramas comerciales no podrían existir si este
requerimiento fuese siempre necesario.
En el año 1969 Benton ideó un tipo de holograma que permitía evitar este requisito
[S.A.Benton Hologram reconstruction with extended light sources. JOSA 59, 1545
(1969)]. En un principio se llamaron white light transmisión holograms y luego fueron
conocidos como rainbow ya que la imagen que produce aparece coloreada como el
arco iris por la descomposición de la luz empleada para su iluminación.
Veremos entonces el procedimiento para realizar este tipo de hologramas.
En una primera etapa se registra un holograma tipo Fresnel, tal como los que hemos
estudiado. Llamaremos a este holograma H1, cuya síntesis se esquematiza en la
siguiente figura
228
Óptica de Fourier Claudio Iemmi
objeto
O1
H1
referencia
r1
holograma 1
utiliza un haz convergente tal que R1 r1* de forma tal que permita obtener una
*
imagen real proporcional a O . Por otra parte, y empleando la propiedad que poseen
los hologramas de objetos difusores, de poder reconstruir una imagen a partir de una
porción del mismo, se pone H1 en contacto con una ranura S tal como se muestra en
la figura
R1= r1*
O 1*
H1
H2
r2
S
t2 x, y tb ' O1* S r2* O1* S r2 O1 S *
2
Luego si se efectúa una reconstrucción de H2 con un haz R 2 r2* , obtenemos
229
Óptica de Fourier Claudio Iemmi
E R 2t2 x, y r2* tb ' O1* S
2
r r ' O S ' r
* *
2 2
*
1 2
2
O
1 S*
Vemos que en el último término aparece la onda objeto original O1 . Basta seguir la
R2 = r2*
S*
H2
O1
Así, si la reconstrucción se efectúa con una sola longitud de onda, el observador debe
“espiar “a través de la ranura S* para ver la imagen virtual O1 . Dado que la ranura es
muy fina, al mirar a través de ella sólo veremos una porción del objeto. Si en vez de
iluminar H2 con una luz monocromática se emplea una fuente de luz blanca, se verá un
continuo de ranuras, cada una con una posición caracterizada por la longitud de onda
con la que se ve coloreada y de acuerdo a esto se verá una porción del objeto O1 de
ese color.
S*
R2 = r2*
H2
O1
230
Óptica de Fourier Claudio Iemmi
Existen algunas variantes de esta técnica. Una de ellas es la que permite realizar el
holograma rainbow en un solo paso. Conceptualmente es el mismo fenómeno el que
se utiliza, sólo cambia la disposición experimental de registro.
Para ello se utiliza una lente convergente ubicada de forma tal que produzca una
imagen virtual de la ranura y una real del objeto
f
f
r
O’
S’ H
S
Holograma de fase
En los cálculos que hemos realizado para los distintos tipos de hologramas siempre
supusimos que la transmisión de los mismos era descripta por una función real
231
Óptica de Fourier Claudio Iemmi
amplitud según lo que hemos visto en materiales de registro. Esto es, si en general el
factor de transmisión de un elemento óptico puede escribirse como
dado por R x .
suponer que existe una relación lineal tal que x I x . Luego si la intensidad
sobre la placa es
232
Óptica de Fourier Claudio Iemmi
cos a cos J o a 2 1 J 2 n a cos 2n
n
n 1
sin a cos 2 1 J 2 n 1 a cos 2n 1
n
n 0
t x K J o a 2 1 J 2 n a cos 2n i 2 1 J 2 n 1 a cos 2n 1
n n
n 1 n 0
P ro exp i Oo2 ro2 J1 2Oo ro exp io
2
233
Óptica de Fourier Claudio Iemmi
P ro2 exp i Oo2 ro2 Oo exp io
2
Vemos que a menos de una amplitud constante y una fase global sin importancia se
reconstruye la onda objeto.
- Holograma de amplitud t x, y t x, y
Vamos a suponer que se registra el diagrama de interferencia más simple, esto es,
x
una red cosenoidal. Luego la transmisión será t x A B cos 2 donde
d
A es la transmitancia promedio, B es la modulación y d es el período. Dado que
1 1 x 1 1 x 1 x
t x cos 2 exp i 2 exp i 2
2 2 d 2 4 d 4 d
x x
x A B cos 2 de manera que t x exp iA exp iB cos 2 d
d
El factor de fase constante puede dejarse de lado. La otra exponencial, tal como
234
Óptica de Fourier Claudio Iemmi
0.6
0.4
0.2
B
/2
J12 B 0.339 o sea ronda el 34%. De este resultado vemos porqué interesan los
medios de fase como materiales de registro.
Existen varios factores que pueden limitar la resolución de una imagen holográfica, por
ejemplo: el tamaño y ancho de banda de las fuentes empleadas para el registro y la
reconstrucción, la resolución y tamaño del medio de registro, las aberraciones en el
frente de onda introducidas por la disposición geométrica empleada, speckle, etc.
Haremos un estudio muy superficial de algunos de estos inconvenientes,
analizándolos separadamente, aunque es evidente que en la práctica se pueden
presentar en forma conjunta.
235
Óptica de Fourier Claudio Iemmi
xo , xr xF
r
a
O
H
f f
I xF r xF O xF r xF O* xF r * xF O xF
2 2
Una vez registrado el holograma la reconstrucción se realiza con una fuente cuya
extensión es distinta de cero.
xR xF
R
H
f f f f xi
1
Se debe prestar atención a que hemos elegido una notación en la que la función y su
transformada de Fourier se denominan con la misma letra, se diferencian una de otra a partir
de la variable de la que dependen indicando en que plano se está trabajando.
236
Óptica de Fourier Claudio Iemmi
1 i 2
R xF R xR R x exp f x x dxR
i f
xF R R F
f
1
P xi xi r xF O xF R xF xi
1 * 1 *
r xr O xo R xR
i f
xF xF xF
f f f
i 2 x
P xi R x O x r x exp f x
R o
*
r
F
R xo xr xi dxF dxr dxR dxo
i 2 xF
P xi O x exp x a x dxF dxo
o
f
o i
O x x a x dx
o o i o O a xi
Supongamos que la extensión del holograma es 2LH. En ese caso para la imagen
primaria se obtiene
237
Óptica de Fourier Claudio Iemmi
xF i 2 xF
P xi O xo rect exp f xo a xi dxF dxo
2 LH
2L
O xo sinc H xo a xi dxo
f
argumento xo a xi por xo a xi ].
Así cada punto de la imagen correspondiente a un punto objeto no será
exactamente un punto sino la transformada de la forma del holograma, por ejemplo
en este caso un sinc2. El ancho característico de esta función (el ancho de su
f
máximo principal) vendrá dado por xi , luego habrá frecuencias espaciales
LH
en el objeto que no serán reproducidas en la imagen dado que existirá una
1 L
frecuencia de corte H (cabe aclarar que acá estamos suponiendo que la
xi f
lente no es un factor limitante y que la falta de resolución se debe solamente al
tamaño finito del holograma). Recordemos que con las redes de difracción sucedía
lo mismo, el tamaño de los órdenes de difracción dependía del la zona iluminada
de la red.
Por ejemplo si en este caso el holograma es una placa de 5 cm de lado, f = 20 cm
y = 632.8 nm, entonces la frecuencia de corte es aproximadamente 400 l/mm.
iii) Supongamos ahora que el tamaño del holograma es tal que no introduce
difracción y que el objeto es puntual, pero el tamaño de la fuente empleada como
referencia y el correspondiente al de reconstrucción vienen descriptos por las
x x
funciones r xr rect r y R xR rect R respectivamente.
2 Lr 2 LR
Tomemos que tanto la fuente usada como referencia así como la empleada para la
238
Óptica de Fourier Claudio Iemmi
xr xR i 2 xF
P xi rect 2 L rect xo A exp xR xo xr xi dxF dxr dxR dxo
r 2 LR f
xr xR i 2 xF
rect 2 L rect exp xR A xr xi dxF dxr dxR
r 2 LR f
xr xR
2 Lr
rect rect xr xR A xi dxr dxR
2 LR
xR A xi xR
2 Lr 2 LR dxR
rect rect
Dado que las dos funciones son reales, esto es una correlación, además al ser
pares coincide con la convolución. Así cada punto imagen se hallará ensanchado y
vendrá descripto por una función triangulo
L r+ L R
1 1
La frecuencia de corte en este caso será
xi 2 Lr LR
-(L r+ L R)
xi
siguiente
xo , xr xF
r
r
H
f f
239
Óptica de Fourier Claudio Iemmi
1 1
onda empleadas está en el rango 1 1 1 y la dispersión en
2 2
sin r
frecuencias espaciales viene dada por f H 1 .
12
Para reconstruir el holograma empleamos una onda plana incidente con un ángulo
2 2
R r y longitud de onda 2 tal que 2 2 2 , donde 2 es la
2 2
longitud de onda media de reconstrucción y 2 es el ancho de banda.
sin i sin i max
sin i min
donde
2 sin r sin r
sin i 2 max f H max 2
2
1
2 1 2
max
1
2 sin r sin r
sin i 2 min f H min 2
2
1
2 1 2
min
1
Luego
2 1 2
sin i
1 1
sin r
2
f sin i .
240
Óptica de Fourier Claudio Iemmi
El tema referente a las aberraciones escapa del tratamiento de este curso ya que es
muy compleja su formulación. Un estudio bastante riguroso acerca de aberraciones en
hologramas se encuentra en dos trabajos tradicionales [R.Meier, JOSA 56, 219 (1966);
R.Meier, JOSA 57, 895 (1967)].
Sólo diremos que para tener en cuenta las aberraciones deben considerarse términos
superiores en el desarrollo en serie que efectuamos para las distancias, es decir,
habría que trabajar con expresiones del tipo
r10 z
x
2
1 y12
x
2
0 y02 xx y y x x y y
4 4
1 0 1 0 1 3 0 1 3 0 ...........
2z 2z z z 8z 8z
241
Óptica de Fourier Claudio Iemmi
De acuerdo a la siguiente figura, asumamos que los vectores de onda k de los haces
que interfieren están en el plano x-z.
X
(x,z)
2 n
ko k r k n donde es la longitud de onda en el vacío y n el índice de
refracción del medio, entonces en dicho punto los campos serán:
O Oo exp i ko . Oo exp i o ; r ro exp i kr . ro exp i r
con
x xˆ z zˆ ; ko kn sin o xˆ coso zˆ ; kr kn sin r xˆ cosr zˆ
luego
242
Óptica de Fourier Claudio Iemmi
Así la diferencia de fase entre las dos ondas viene dada por
Tendremos máximos cuando esa diferencia de fase sea igual a 2m . Por lo tanto
resolviendo para x obtenemos una ecuación que define la posición de las franjas, la
cual viene dada por
cos o cos r m n
x z
sin r sin o sin r sin o
Vemos que consisten en franjas rectas cuya pendiente es el término entre corchetes y
m n
su ordenada al origen. En realidad para obtener exactamente la posición
sin r sin o
de los máximos deberíamos incluir el desfasaje de las dos ondas en el origen pero,
como veremos, una translación de las franjas no afecta los cálculos subsiguientes.
El espaciado de las franjas en la dirección x se puede deducir analizando cuál es la
variación en x cuando m cambia en una unidad. Así el período vendrá dado por
n
x d
sin r sin o
y la frecuencia espacial será
1 sin r sin o
f
d n
además que el ángulo que forman las franjas con el eje z puede deducirse de la
2sin o r sin r o
cos o cos r 2 2 tan o r o r
tan 2 2
sin r sin o o r r o
2 cos sin
2 2
243
Óptica de Fourier Claudio Iemmi
1 sin r sin o 2 r o 2
f cos o r sin 2 sin 2 cos
d n n 2 n
2
do d cos ; fo sin
n 2
X do
difractado en el punto xi , yi por una pantalla ubicada a una distancia z del origen.
Hasta ahora cuando calculábamos la difracción de Fraunhofer producida por una
abertura, considerábamos que la misma estaba ubicada en el origen de coordenadas.
Ahora debemos encontrar una expresión más general donde se tenga en cuenta que
la pantalla difractora es encuentre a una distancia z del origen.
244
Óptica de Fourier Claudio Iemmi
S L
So
i Lo
R Z
S 2 xR x zR z ; So2 xR2 z R2
2 2
Dado que las distancias So y Lo son grandes en comparación con x,z de la pantalla
difractora, podemos escribir en forma aproximada
x2 z 2 xR x zR z x2 z 2 xi x zi z
S So ....... ; L Lo .......
2 So 2 So So So 2 Lo 2 Lo Lo Lo
245
Óptica de Fourier Claudio Iemmi
A
E xi , zi exp ik So Lo .
i So Lo
xR x z R z xi x zi z x 2 z 2 x 2 z 2
. exp ik dx
o
S S o Lo Lo 2 S o 2 Lo
x x z z x x z z
E xi , zi Cte. exp ik R R i i dx
So So Lo Lo
Además, en esta aproximación, puede suponerse que el ángulo que forma So con el
eje z es aproximadamente igual al que forma S con dicho eje. Análogamente el ángulo
que forma Lo con el eje z es aproximadamente igual al que forma L. Luego, de
acuerdo a la convención de signos anteriormente tomada y teniendo en cuenta la
siguiente figura, podemos escribir
So
Lo
R+
i
R Z
zR zi
cos R cos R cos cos R ; cos i
So Lo
xR xi
sin R sin R cos sin R ; sin i
So Lo
Así el campo difractado en una dirección sin i ,cosi , debido a una pantalla
difractora plana (de espesor dz) ubicada a una distancia z del origen, viene dado por
246
Óptica de Fourier Claudio Iemmi
d E sin i , cos i Cte. exp ik sin i sin R x cos i cos R z dx dz
2
donde t x, z es la función transmisión de la red y knR , con nR la longitud
nR
de onda del haz de reconstrucción en el medio de registro.
X
H
(x,z)
-T/2 T/2
Z
-H
dz
en función de los resultados hallados anteriormente. Por un lado habíamos visto que
247
Óptica de Fourier Claudio Iemmi
Además
cos o cos r sin
tan cos r cos o sin r sin o
sin r sin o cos
sin
r x, z o x, z kn x sin r sin o z sin r sin o
cos
kn sin r sin o
x cos z sin
cos
sin r sin o 2
f sin cos f o cos
n n 2
resultando así
r x, z o x, z 2 fo x cos z sin
Ahora que tenemos esta ecuación veamos cuales son las condiciones de difracción de
Bragg. Para ello reemplacemos la expresión hallada para la función transmisión en la
ecuación que expresa el campo difractado.
248
Óptica de Fourier Claudio Iemmi
t
H 2
E sin i , cos i Cte. b ' Oo2 exp iknR sin i sin R x .
H T
2
obtenemos
E sin i , cos i Cte. tb ' Oo2 T 2 H sinc HknR sin i sin R .
T
.sinc knR cos i cos R
2
Cte. ' Oo ro T 2 H sinc H knR sin i knR sin R 2 f o cos .
T
.sinc knR cos i knR cos R 2 f o sin
2
Cte. ' Oo ro T 2 H sinc H knR sin i knR sin R 2 f o cos .
T
.sinc knR cos i knR cos R 2 f o sin
2
Esta es la solución a primer orden para la difracción por una red tridimensional.
Analicemos esta ecuación:
El primer término tiene un máximo cuando
sin i sin R
R i
cos i cos R
Esto corresponde al orden 0, es decir el haz incidente de reconstrucción es transmitido
por el holograma directamente sin desviarse
El segundo término es máximo cuando el argumento de ambas funciones sinc es cero,
esto es
2 f o
sin i sin R cos
knR
2 f o
cos i cos R sin
knR
249
Óptica de Fourier Claudio Iemmi
Vemos que reobtenemos la primera de las ecuaciones, vale decir que una de las
condiciones para que el segundo término del campo difractado por la red
tridimensional sea máximo, es justamente que se cumpla la ecuación de la red.
Para analizar la otra condición estudiemos cómo se reflejaría la onda incidente sobre
un plano coincidente con una de las franjas de interferencia. Sea io el ángulo de
incidencia y ro el de reflexión
cos io cos ro cos R cos sin R sin cos i cos sin i sin
sin
cos R cos i sin i sin R
cos
pero habíamos visto que
2 f o 2 f o
sin i sin R cos cos R cos i sin
knR knR
la dirección sin i , cosi cuando se cumplen las condiciones de Bragg, esto es,
250
Óptica de Fourier Claudio Iemmi
plano que, al igual que las franjas de interferencia, forma un ángulo con el eje Z.
2 f o 2 f o
sin i sin R cos ; cos i cos R sin
knR knR
y el segundo cuando
2 f o 2 f o
sin i sin R cos ; cos i cos R sin
knR knR
Veamos cómo podemos poner estas expresiones en función de los haces de registro.
nR
sin i sin R nR f o cos nR f sin r sin o
n
además
nR
cos i cos R nR f o sin nR f tan cosr coso
n
del haz objeto pero a través del mismo paso, esto es, se obtiene el objeto conjugado.
251
Óptica de Fourier Claudio Iemmi
O*
r
R= r*
O
ambos haces tienen igual dirección, esto es, cos r cos R ; sin r sin R ,
entonces sin i sin o ; cos i coso . Luego el haz imagen coincide con el haz
O
R= r
O
252
Óptica de Fourier Claudio Iemmi
Sensibilidad a la orientación
Para determinar el rango angular dentro del cual el holograma puede reconstruirse
P sin i , cos i Cte. ' Oo ro T 2 H sinc H knR sin i knR sin R 2 f o cos .
T
.sinc knR cos i knR cos R 2 f o sin
2
Llamemos I o Cte. ' Oo ro T 2 H
2
a la intensidad difractada cuando se cumplen las
condiciones de Bragg, esto es, cuando las dos funciones sinc valen 1.
Sea entonces
T
I I o sinc2 knR cos i knR cos R 2 f o sin
2
la intensidad difractada cuando se cumple la ecuación de la red pero no la de reflexión
en los planos de Bragg. La intensidad de la onda difractada valdrá cero cuando el
argumento de la función sinc sea igual a y esto nos indicará cuál es la rotación
angular tolerada por el holograma para poder reconstruir el frente de ondas. Vamos a
escribir dicho argumento en función de los haces de registro, para ello tengamos en
cuenta que
2
fo sin ; cos o cos r 2sin sin
n 2 2
2
luego 2 f o sin cos o cos r con lo cual el argumento toma la forma
n
253
Óptica de Fourier Claudio Iemmi
T 2 2
cos i cos R cos r cos o
2 nR n
Vamos a tomar que la longitud de onda de registro es igual a la de reconstrucción, así
la expresión anterior queda
T
cosi cos R cosr coso
n
Ahora bien, un cambio en la dirección del haz de reconstrucción cos R dará lugar
sin i sin R 0
Ahora bien, tenemos que
cos R
cos i i cos R R 0 i R
cos i
n cos R
sin i i sin R R sin i R sin R R
T cos i
254
Óptica de Fourier Claudio Iemmi
Con lo cual
1
n cos R
R sin R sin i
T cos i
n cos o n cos o
R
T sin r cos o sin o cos r T sin
Consideremos ahora el efecto que causa iluminar al holograma con una longitud de
onda diversa a la usada para el registro. Nuevamente estudiemos el término
correspondiente a la imagen primaria.
P sin i , cos i Cte. ' Oo ro T 2 H sinc H knR sin i knR sin R 2 f o cos .
T
.sinc knR cos i knR cos R 2 f o sin
2
255
Óptica de Fourier Claudio Iemmi
2
fo sin ; cos o cos r 2sin sin ; 2cos sin sin r sin o
n 2 2 2
máximo para sin i ,cosi sin o ,coso . Queremos ver ahora qué cambio en la
longitud de onda hace que este haz difractado se anule.
Sabemos que un cambio en la longitud de onda de reconstrucción
2 2
Ahora bien, por un lado tenemos que kn kn n y si
n n2
256
Óptica de Fourier Claudio Iemmi
n2 cos o n cos o
2
n
T 1 cos r cos o sin o sin r
T 1 cos
Tal como hemos dicho con anterioridad la mayoría de los hologramas deberían
considerarse con un espesor no despreciable. Sin embargo el hecho de que una sola
imagen se forme para una dada dirección de reconstrucción, se presenta solamente
en los hologramas verdaderamente gruesos. Así la mayoría se encuentra en una
situación intermedia.
Un holograma de volumen, ya sea de transmisión o reflexión, alcanza una eficiencia
de difracción teórica del 100% y experimentalmente (por ejemplo en gelatinas
dicromatadas) del 90%.
Hemos visto que hologramas de volumen de transmisión poseen una alta sensibilidad
a la orientación y al cambio de longitud de onda
n cos o n2 cos o
R
T sin
n
T 1 cos
257
Óptica de Fourier Claudio Iemmi
cambios de orientación. Esto hace que estos hologramas sean muy útiles ya que
pueden reconstruirse con luz blanca y desde casi cualquier dirección.
Los hologramas de reflexión fueron inventados por el ruso Denisyuk y sus primeras
investigaciones se encuentran en el trabajo [Y.Denisyuk, Soviet Phys. Doklady, 7, 543,
(1962)].
De acuerdo con los análisis que realizamos previamente veamos cómo son las
disposiciones geométricas para registrar y reconstruir hologramas de volumen.
Estudiemos primero el caso de los de transmisión
REGISTRO
planos de Bragg
Holograma
r
RECONSTRUCCIÓN
reflexión en los
planos de Bragg
O Imagen
primaria virtual
Holograma
R=r
RECONSTRUCCIÓN
reflexión en los
planos de Bragg
R = r* O* Imagen
conjugada real
Holograma
258
Óptica de Fourier Claudio Iemmi
REGISTRO
planos de Bragg
O
r
Holograma
RECONSTRUCCIÓN
reflexión en los
planos de Bragg
R=r O Imagen
primaria virtual
Holograma
Dado que este último tipo de hologramas es poco sensible a la orientación es posible
ver la imagen de un color para un ángulo de incidencia y de otro color (cuya longitud
de onda cumpla la nueva condición de Bragg) para otra orientación. Obviamente
suponemos que iluminamos con luz blanca. Cabe destacar que en el proceso
fotográfico pueden producirse dilataciones o contracciones que hacen que cambie el
espaciado de los planos de Bragg lo que se manifiesta en un cambio de la longitud de
onda reflejada con respecto a la de registro.
259
Óptica de Fourier Claudio Iemmi
INTERFEROMETRÍA HOLOGRÁFICA
Vamos a analizar el proceso planteando las ecuaciones en una dimensión para que
sean más sencillas. En una primera etapa se registra un holograma fuera de eje del
260
Óptica de Fourier Claudio Iemmi
i i kxsinr
referencia r x ro e r ro e , luego la transmisión del holograma
t x tb ' Oo2 x ro 2
' roOo x exp i o x kx sin r exp i o x kx sin r
El holograma una vez revelado, se repone exactamente en la misma posición que se
registró, ó se revela in situ y posteriormente se ilumina con la onda de reconstrucción
sólo que ahora O ' x es la onda correspondiente al objeto en su nuevo estado. Esto
i ' x
es O ' x Oo x e o , donde la amplitud es básicamente la misma que la
anterior, ya que la perturbación es mínima, y sólo se considera la variación en la fase.
' Oo2 x ro exp i o x o' x kx sin r ' ro2Oo x exp io x
De todos estos términos sólo nos interesan los que no contienen la información de la
fase del haz de referencia. Tales términos describen a ondas difractadas en
aproximadamente la misma dirección. Ellos son
261
Óptica de Fourier Claudio Iemmi
donde vemos que tanto los dos primeros términos como el coeficiente de la función
coseno del tercer término son magnitudes que no depende de la distribución de fase y
y final de fase.
Esta técnica tiene la ventaja de permitir ver en tiempo real una deformación ya que se
compara instante a instante la onda objeto inicial, registrada en el holograma, con la
onda objeto deformada que lo ilumina. Los inconvenientes residen en la reubicación
correcta de la placa y en el cuidado que se debe tener para que la emulsión fotográfica
no se deforme en el revelado y secado. Para esta técnica resultan ideales los
materiales de registro fototermoplásticos ya que, como vimos, se revelan in situ.
Luego
262
Óptica de Fourier Claudio Iemmi
' O1 x O2 x r * x O1 x O2 x r x
*
En este caso nuevamente se obtiene una imagen del objeto modulada por franjas de
interferencia que denotan las deformaciones sufridas por el mismo.
Otra forma de interpretar esta técnica es como un batido de las franjas registradas.
Cada holograma consiste en un registro de franjas de interferencia que varían muy
poco de uno a otro. El registro total consistirá en un diagrama de franjas que en
algunas zonas se refuerzan y en otras se borronean formando batidos o figuras de
moiré. Cuando reconstruimos el holograma, las zonas borroneadas no contribuirán con
luz difractada mientras que sí lo harán aquellas zonas donde las franjas se refuerzan.
La técnica de doble exposición es muy apropiada para el registro de fenómenos que
varían rápidamente en el tiempo (ondas de choque, flujo de fluidos, etc), se utiliza para
ello una fuente luminosa tal como un láser pulsado. Un ejemplo de este tipo de
interferogramas se muestra a continuación.
263
Óptica de Fourier Claudio Iemmi
264
Óptica de Fourier Claudio Iemmi
O F xI I
f f f f f yI
Recordemos que en estos casos la función de entrada sobre el plano objeto era
xF yF
Eo C1 , . Sobre dicho plano se ubicaba un filtro cuya transmisión
f f
x y
era tF xF , yF C2 H F , F donde C1 y C2 son constantes complejas. Así la
f f
distribución de campo detrás del filtro era proporcional a .H . Finalmente la lente LA
265
Óptica de Fourier Claudio Iemmi
x y 1 1 b 1 b
H F , F 1 exp i 2 yF exp i 2 yF
f f 2 2 f 2 f
b
La frecuencia espacial de dicha red es y tiene en cuenta la posición de la red
f
con respecto al eje z. Esto es, si no incluimos estamos asumiendo que para
yF 0 tenemos un máximo de transmisión, pero nosotros deseamos tener la libertad
de desplazar la red en el plano en el que se halla de manera que debemos incluir esa
variable.
Con este tipo de filtro un objeto ubicado en el plano de entrada será difractado en tres
imágenes sobre el plano de salida. Así dos objetos que no se superponen en el plano
o serán difractados en seis imágenes sobre i, algunas de las cuales pueden
superponerse.
Supongamos entonces que sobre o tenemos una distribución de campo dada por:
Eo xo , yo O1 xo , yo b ' O2 xo , yo b '
Luego
x f , y f
F F xF yF
,
exp i 2
b ' yF
xF y F
,
exp i 2
b ' yF
f
f f
1 2
f f f
donde las exponenciales aparecen debido al desplazamiento del objeto fuera del eje
óptico en b’. Si elegimos b’= b, después del filtro el campo es:
xF y F b yF 1 xF y F b yF
.H 12 1
,
exp i 2
2 2
,
exp i 2
f f f f f f
1 xF y F 1 xF y F b
, exp i , exp i 4 yF
f f f f f
1 2
4 4
1 xF y F b 1 xF y F
, exp i 4 yF , exp i
f f f
1 2
4 4 f f
266
Óptica de Fourier Claudio Iemmi
De todos estos términos sólo nos interesan los que se superponen en el origen
1 1 1
O1 xi , yi exp i O2 xi , yi exp i O1 xi , yi O2 xi , yi exp i 2 exp i
4 4 4
obtenemos la suma de los objetos. Si se desplaza la red de forma tal que se
2
obtiene la resta
O F I
LT LA
f f f f
267
Óptica de Fourier Claudio Iemmi
1
1 b 1 b
H1 1 exp i 2 yF exp i 2 yF
2
2 f 2 f
1
1 b 1 b
H2 1 exp
i 2 y F exp
i 2 y F
2
2 f 2 f
derivar O xo , yo . Luego cada red sobre el plano final producirá tres imágenes de
este objeto. Las correspondiente al orden 0 de difracción no sufren cambios, pero
analicemos por ejemplo que sucede en un orden difractado. Ahí tendremos una
1 1
Ei xi , yi ........... O xi , yi b exp i O xi , yi b ........
4 4
268
Óptica de Fourier Claudio Iemmi
La utilización de hologramas como filtros espaciales fue introducida por Vander Lugt,
quien propuso una nueva técnica para la síntesis de filtros complejos [A.Vander Lugt
Signal Detection by Complex Spatial Filtering, IEEE Trans. Inform. Theory, IT-10, 2
(1964)].
La principal característica de los filtros así sintetizados es que permiten registrar la
amplitud y fase mediante métodos interferométricos. Originalmente para la síntesis de
este tipo de filtros se utilizaban sistemas como los esquematizados en la figura.
Notemos la analogía con un holograma fuera de eje tipo Fourier.
Lc
prisma
S
yF
xF
r
placa
señal f f
LT
xo xF
referencia LT
señal
placa
f f
En este caso se toma como objeto la respuesta al impulso h deseada para el filtro. La
lente LT proporciona sobre el plano de registro una distribución de campo
269
Óptica de Fourier Claudio Iemmi
1 x y
H F , F proporcional a la transformada de Fourier de la señal de interés,
f f f
sin r
además sobre dicho plano incide el haz de referencia r ro exp i 2 yF .
Luego la distribución de intensidad sobre el film será
2
1 x y r x y sin r
I xF , yF ro2 2 2 H F , F o H F , F exp i 2 yF
f f f f f f
ro x y sin r
H * F , F exp i 2 yF
f f f
x y x y x y
H F , F H o F , F exp i H F , F queda registrada en el tercer
f f f f f f
término. Como supusimos siempre, se realiza una exposición tal que se trabaja en el
rango lineal del film de modo que la transmisión del filtro resulta
x y xF y F
t F , F I , .
f f f f
Veamos ahora como procesa una señal. El filtro holográfico de Vander Lugt se ubica,
como todo filtro en un procesador óptico coherente, en el plano transformado.
Si recordamos lo hecho para hologramas tipo Fourier, la placa se ubicaba en el foco
de una lente antitransformadora y se la iluminaba con un haz de reconstrucción plano.
Ahora en vez de un haz plano, sobre el filtro llega la transformada de Fourier de la
salida del filtro tendremos
x y 1 sin r
x f , y f r f 1f H o2
2
ro
t F , F . F F o
H exp i 2 yF
f f f 2 f 2
3 3
sin r
ro
2 f 2
H *
exp
i 2
yF
270
Óptica de Fourier Claudio Iemmi
1
Ei xi , yi ro2 Eo xi , yi 2 i
h x , yi h* xi , yi Eo xi , yi
f 2
ro
h xi , yi Eo xi , yi xi , yi f sin r
f
r
o h* xi , yi Eo xi , yi xi , yi f sin r
f
El primero y segundo término no son de interés y dan una distribución de luz centrada
en el origen de coordenadas. Analicemos el tercer y cuarto término.
h xi , yi Eo xi , yi xi , yi f sin r
h x x , y f sin
i o i r yo Eo xo , yo dxo dyo
h* xi , yi Eo xi , yi xi , yi f sin r
E x , y h x xi , yo yi f sin r dxo dyo
*
o o o o
271
Óptica de Fourier Claudio Iemmi
1 ro2 Eo xi , yi WE
1
2 i
2 h x , yi h* xi , yi Eo xi , yi 2Wh WE
f
2
ro
3 h x , y Eo xi , yi xi , yi f sin r Wh WE
f i i
ro
4 h* xi , yi Eo xi , yi xi , yi f sin r Wh WE
f
Luego esto nos da una idea de cuál es el ángulo que debemos elegir para el haz de
referencia r de modo tal que los términos no se superpongan.
convolución Wh + WE
f.sinr
2Wh + WE
WE
xi
f.sinr
correlación Wh + WE
yi
Se dice que un filtro está adaptado (matched filter) a una señal particular s xo , yo si
x y xF y F
H F , F *
f
transferencia resulta , . Analicemos la acción de este
f f f
filtro en un correlador tipo Vander Lugt, para ello centrémonos en el tercer término
mencionado en el punto anterior. El mismo da
s* xi , yi Eo xi , yi xi , yi f sin r
E x , y s x xi , yo yi f sin r dxo dyo
*
o o o o
272
Óptica de Fourier Claudio Iemmi
x f , y f
F F xF yF
, . En consecuencia el campo emergente del filtro viene
f f
dado por
x y xF y F ro2 sin r
.
1 ro
t F , F , 2
2
exp i 2 yF
f f f f f f3 f2
o o
3 2
ro sin r
exp i 2 yF
f2
2
2
Dado que tanto ro como o son reales, el tercer término representa una onda plana,
*
pues la curvatura de fue cancelada por la de . Luego tenemos que esta onda
plana, emergente del filtro, es antitransformada por la segunda lente del procesador
dando origen a un punto brillante en el plano focal de la misma. Si en cambio la señal
de entrada no coincide con la que fue sintetizado el filtro, las curvaturas no se
compensarán y dicho término dará origen a una distribución de luz tanto menos
intensa cuanto menor sea la correlación entre el objeto de entrada y la señal a
detectar. Otra forma de ver esto es la siguiente: el filtro fue sintetizado como un
holograma a partir de la interferencia del haz correspondiente a la señal y un haz plano
de referencia. En general se utiliza un haz de reconstrucción igual al de referencia con
el propósito de reobtener la onda objeto, sin embargo si iluminamos el holograma con
un haz de reconstrucción igual al haz objeto lo que reobtendremos es el haz de
referencia, esto es una onda plana.
273
Óptica de Fourier Claudio Iemmi
s x , y s x , y dx dy j 1...26
*
j o o 3 o o o o
integral. Por último cada señal s j xo , yo aportará distinta cantidad de luz así que
para independizarnos de esto debemos normalizar la expresión por
2
s x , y
2
3 o o dxo dyo
s3 xo , yo dxo dyo
2
s x , y
2
3 o o dxo dyo
2
s x , y s x , y dx dy
*
j o o 3 o o o o
j3
s x , y
2
j o o dxo dyo
2
Luego
274
Óptica de Fourier Claudio Iemmi
2
s x , y s x , y dx dy
*
j o o 3 o o o o
s3 xo , yo dxo dyo
2
s x , y
2
j o o dxo dyo
transferencia del filtro es H H o exp i H . J.Horner y P.Gianino [Appl. Opt. 23, 812
275
Óptica de Fourier Claudio Iemmi
burdos, en tanto que una cantidad considerablemente menor irá a las frecuencias más
altas que son las que contienen la información de los detalles más finos del objeto. Si
se construye un filtro en el cual eliminamos la información de la amplitud, estamos
dándole el mismo peso a todas las frecuencias de forma que logramos un aumento de
los detalles más finos. Obviamente este filtro también debe realizarse de forma digital
y luego transferirse a un medio que permita utilizarlo en el correlador.
A continuación se muestra una escena, su transformada de Fourier en módulo y fase y
la obtención de la antitransformada a partir de la información sólo de amplitud y sólo
de fase
Módulo TF Fase TF
Escena
276
Óptica de Fourier Claudio Iemmi
Filtro clásico
El correlador por transformada conjunta (JTC por Joint Transform Correlator), ideado
por Weaver y Goodman [C.Weaver, J.Goodman, A technique for optically concolving
two functions, Appl. Opt. 5, 1248 (1966)], ha sido diseñado originalmente para llevar a
cabo la operación de convolución entre dos funciones y por ello resulta de interés para
emplearlo en operaciones de filtrado.
La ventaja fundamental que presenta este dispositivo respecto del correlador 4f, se
debe a que su operación está basada en un proceso de registro holográfico, y luego
en un proceso de lectura empleando un haz plano, evitando así los inconvenientes de
alineación. La diferencia fundamental entre ambos dispositivos, radica en que tanto la
señal que se desea procesar, como la respuesta al impulso deseada, están presentes
simultáneamente durante el proceso de registro, por lo que basta emplear un haz de
lectura a fin de obtener la señal se salida deseada. En la figura se presenta un
esquema del proceso de registro en este dispositivo. Un haz plano de longitud de onda
277
Óptica de Fourier Claudio Iemmi
ilumina el plano de entrada del correlador, donde se representan, separadas por una
Lc LT
yF
S xF
f f Plano
transformado
conjunto
La señal de entrada está representada entonces por
d d
Eo xo , yo s xo , yo h xo , yo
2 2
Una lente convergente LT de distancia focal f es empleada para obtener en el plano
transformado conjunto JTP (Joint Transform Plane), ubicado a distancia focal de la
misma, la interferencia de las transformadas de Fourier de ambas señales,. Dicha
transformada viene representada por
xF yF d 1 xF yF d
x , yF
1
f , exp i yF
f f f f
H , exp i y
f
f f
F F
1 xF y F
2 2
xF y F
I xF , y F 2 2 H
f f f
, ,
f f
x y x y d
* F , F H F , F exp i 2 yF
f f f f f
x y x y d
H * F , F F , F exp i 2 yF
f f f f f
278
Óptica de Fourier Claudio Iemmi
Esta distribución debe ser registrada en algún medio sensible a la luz situado en el
plano transformado conjunto del sistema. Dicho registro constituye lo que
denominaremos transformada conjunta del correlador. Suponemos que se efectúa un
registro que posee una transmitancia lineal con la intensidad que llegó al mismo. El
proceso de lectura de la transformada conjunta se esquematiza en la figura
Lc LT
yF
S
xF
xi
yi
Plano f f
transformado
conjunto
De modo que si se la ilumina con un haz colimado de amplitud unitaria, tal como se
esquematiza en la figura, sobre el plano final tendremos
1
Ei xi , yi h xi , yi h* xi , yi s xi , yi s * xi , yi
f
h xi , yi s* xi , yi xi , yi d
h* xi , yi s xi , yi xi , yi d
Por lo tanto, la señal de salida está constituida en primer lugar, por los dos primeros
términos que dan respectivamente la autocorrelación de la escena y la referencia,
ambas distribuciones centradas en el origen y que no resultan de interés en este
análisis. Obtenemos además las correlaciones cruzadas
h xi , yi s* xi , yi xi , yi d
hx , y s x xi , yo yi d dxo dyo
*
o o o
279
Óptica de Fourier Claudio Iemmi
h* xi , yi s xi , yi xi , yi d
sx , y h x xi , yo yi d dxo dyo
*
o o o
Una ubicada en las coordenadas (0, -d) y la otra en (0, d). Así, una vez definida la
función transferencia asociada a un proceso particular, sólo es necesario calcular la
respuesta al impulso correspondiente e invertirla y conjugarla para emplearla como
señal de referencia en el correlador. En particular, para aplicaciones de
reconocimiento, donde se requiere una señal de salida dada por la correlación entre la
escena y la referencia, basta con representar ambas señales en el plano de entrada
del correlador. Notemos que en esa situación, las dos señales de correlación entre la
escena y la referencia que se obtienen, son una la reflexión especular de la otra.
Este correlador puede interpretarse como una experiencia de Young generalizada en
la cual en vez de simples aberturas iguales tenemos objetos. Si estos objetos (alguno
de la escena y el de referencia) son iguales entonces en el plano transformado los
espectros se superpondrán e interferirán dando lugar a franjas. Luego al iluminar estas
franjas, las mismas se comportarán como una red de difracción y darán lugar a
órdenes (picos de correlación). Si los objetos no son semejantes los espectros no
coinciden y por lo tanto al no interferir no se registrarán franjas.
Estos correladores poseen la ventaja de no utilizar filtros, cuyo posicionado siempre es
crítico, y además son aptos para procesado en tiempo real ya que la entrada puede
representarse en una pantalla de cristal líquido, el espectro se puede captar con una
CCD y luego representarse en una pantalla de cristal líquido para finalmente captar la
correlación con una CCD.
Cabe destacar que cuando este tipo de correlador es utilizado para el reconocimiento
de formas su poder de discriminación es equivalente al del filtro adaptado. Sin
embargo existen técnicas de procesado en la señal de entrada y en el plano de
registro que permiten llevar a cabo un filtrado equivalente al que se realiza con un filtro
POF.
MEMORIAS HOLOGRÁFICAS
El almacenado óptico de datos ha sido un tema de gran interés durante las últimas
décadas. Desde la aparición comercial de CDs y DVDs su uso creció
exponencialmente así como la necesidad de aumentar su capacidad de
almacenamiento y su velocidad de acceso. Los primeros CDs fueron evolucionando a
280
Óptica de Fourier Claudio Iemmi
Láser SLM
BS Lente
Espejo
Scanner Espejo
Haz de
referencia Señal
Cristal
CCD
281
Óptica de Fourier Claudio Iemmi
O haz objeto
R haz de reconstrucción
r haz de referencia
O
R O´ BS
r Holograma
r Holograma
r Holograma
de volumen de volumen de volumen
O
R R
O
O´ O´
282
Óptica de Fourier Claudio Iemmi
r'
r'
O O O
r r r
r' O
r r O
desplazamiento de un haz multiplexado angular +
esférico de referencia multiplexado espacial
283
Óptica de Fourier Claudio Iemmi
284
Óptica de Fourier Claudio Iemmi
por lo que se hace necesario también analizar cómo se transforma Fourier una señal
digital.
DIGITALIZACIÓN
Para que una imagen del tipo analógico pueda pasar de su estado continuo a una
representación numérica en la computadora deben realizarse ciertos pasos
intermedios de transformación.
En primer lugar la imagen debe ser registrada por algún dispositivo electrónico o
generada numéricamente y luego debe almacenarse en la memoria de una
computadora. La resolución espacial de esta imagen dependerá del número de píxeles
que contenga, es decir de cómo se realizó el muestreo de la señal. En la siguiente
figura se muestra un ejemplo de una imagen digitalizada usando distintos tamaños de
píxel.
285
Óptica de Fourier Claudio Iemmi
El nivel de cuantización dependerá del número de bits con el que trabaje el conversor
analógico digital. En cuanto al número de píxeles que deben escogerse para que la
pérdida de información se reduzca al mínimo o sea inexistente, existe un teorema que
contempla este punto en base al estudio del proceso de discretización de una función
continua. Recibe el nombre de teorema del muestreo o de Shannon y es el tema a
tratar a continuación.
Sea f(x,y) una función continua cuya transformada de Fourier es finita. El objetivo es
representar esta función a través de un arreglo de puntos.
Muestrearla de forma regular a un intervalo determinado involucra el hecho de aceptar
de la función solo aquellos valores que coinciden con múltiplos de este intervalo.
Matemáticamente, esta operación equivale a multiplicarla por la función peine de
Dirac:
p( x , y ) x nT y mQ
n m
286
Óptica de Fourier Claudio Iemmi
f M ( x, y) f ( x, y) x nT y mQ
n m
FM f,f
x y F fx , f y P fx , f y
1 n m
donde P( f x , f y ) fx fy
TQ n m T Q
corresponde a la
1 n m
FM ( f x , f y )
TQ n m
F ( fx , f y ) fx f y
T Q
1 n m
F fx , f y
TQ n m T Q
287
Óptica de Fourier Claudio Iemmi
Espectro de la
función f(x,y)
uo
fx
vo
fy
Espectro de la
función p(x,y)
fx
fy
Espectro de la
función fM(x,y)
fx
fy
Analizando el espectro de la función muestreada se puede ver que debe existir alguna
relación entre el tamaño de F( fx, fy ) y las distancias 1/T y 1/Q entre los distintos
espectros como para que estos no se superpongan. Si uo y vo son las máximas
frecuencias de F( fx, fy ) en las direcciones fx y fy respectivamente, las situaciones
planteadas son las siguientes:
Si los intervalos de muestreo son tal que 1/T 2uo, 1/Q 2vo entonces los
espectros no se superpondrán, dando lo que se muestra en la figura anterior.
288
Óptica de Fourier Claudio Iemmi
Si por el contrario, los intervalos son tal que 1/T < 2uo, 1/Q < 2vo entonces los
espectros tendrán porciones superpuestas y se formará lo que se llama aliasing1,
haciendo imposible la reconstrucción de la señal. Esta situación está representada
en forma unidimensional en la siguiente figura
Espectro de la función
Espectro de la función
muestreada
1
El término aliasing se refiere a un falso traslado. Es decir, si se considera uno solo de los
espectros, la porción del mismo que queda afuera del rango de frecuencias determinado por
los intervalos de muestreo se dice que es “falsamente trasladada” dentro de ese rango. La
expresión “falsamente” viene a cuenta de que en realidad estas frecuencias espurias presentes
en el espectro provienen de los componentes vecinos.
289
Óptica de Fourier Claudio Iemmi
repetida con período 1/T y multiplicada por 1/T como indica el teorema del muestreo.
En la siguiente figura se grafican las cuatro funciones.
transformada de
función analógica s(x) Fourier S(f)
(a) (b)
(c) (d)
Si ahora a la señal sM(x) se la repite con período mT se obtendrá una función b(x)
dada por
b x s x x kmT
k
M
1 k
B f
mT k
S M f f
mT
290
Óptica de Fourier Claudio Iemmi
función b(x) obtenida
de replicar sM(x)
transformada de
Fourier B(f)
s M ( x) snT x nT
n
291
Óptica de Fourier Claudio Iemmi
1 k
B f
mT k
S M f f
mT
se obtiene
1 k
B f snT exp i 2 nTf f
mT k n mT
1 n k
mT k n
snT exp i 2 k f
m
mT
Comparando esta expresión con la anteriormente hallada para B(f) en término de los
Bk se obtiene que
1 n
Bk
mT n
s(nT ) exp i 2 k
m
n
snT T
k ( en periodo 1/ T )
Bk exp i 2 k
m
Estas dos últimas ecuaciones expresan el hecho de que s(nT) (sk en la figura) y Bk son
transformada y antitransformada discreta respectivamente.
Una vez establecida la forma en que una imagen es expresada como un arreglo
numérico y luego llevada al plano de Fourier de manera también digital, el paso
siguiente sería generar el holograma computacionalmente. Existen varias técnicas
para esto último, algunas de las cuales veremos a continuación.
292
Óptica de Fourier Claudio Iemmi
Los trabajos en los que se describe este método son [B.Brown, A.Lohmann Complex
spatial filtering with binary mask Appl. Opt. 5, 967 (1966)], [A.Lohmann, D.Paris, Binary
Fraunhofer Holograms, generated by computer, Appl. Opt. 6, 1739 (1967)]
En este método tenemos una malla de M x N píxeles (número de puntos en los que
fue muestreado el objeto) de los cuales debe codificarse la amplitud y la fase. Dicha
codificación se realiza mediante aberturas rectangulares cuyo tamaño, dentro del
píxel, determina la amplitud y su posición la fase.
Y
x
y
n.y A.y
m.x X
El píxel (m, n) tiene, como todos, tamaño x, y y su posición es (mx, ny); la
cantidad de luz que deja pasar depende del alto de la abertura Ay (el ancho se deja
293
Óptica de Fourier Claudio Iemmi
Este método se halla detallado en las publicaciones [W.Lee Sampled Fourier transform
hologram generated by computer, Appl. Opt. 9, 639 (1970)], [C.Burckhardt, A
simplification of Lee´s method of generating holograms by computer, Appl. Opt. 9,
1949 (1970)].
El mismo se basa en que un número complejo puede descomponerse como
2 2
A exp i A1 exp i Ao A1 exp i
3 3
con
2 2
exp i
*
exp i
A1 3 3
3
*
Ao
3
2 2
exp i
*
exp i
A1 3 3
3
294
Óptica de Fourier Claudio Iemmi
295
Óptica de Fourier Claudio Iemmi
Una vez generados estos hologramas deben transferirse a un medio físico que permita
utilizar el filtro en un procesador óptico. Hay medios más o menos sofisticados que
permiten llevar a cabo este propósito. Una de las formas más elementales en la que
esto puede hacerse es realizar una impresión del holograma en papel y luego realizar
una fotorreducción. Otro, un poco más elaborado, es realizar directamente la
impresión utilizando una impresora de artes gráficas de alta resolución
(aproximadamente 3000 dpi) que proporciona el filtro sobre película fotográfica.
Como dato curioso, existe una fotografía en donde se ve a Lohmann y sus
colaboradores, pintando con pincel un holograma sobre una pared para después
reducirlo fotográficamente. Recordemos que cuando Lohmann propuso su método
para sintetizar CGHs no existían impresoras y prácticamente no había computadoras.
Constituye todo un desafío estudiar y proponer la utilización de un elemento para el
que aún no están desarrollados los métodos de implementación.
Independientemente de cuál sea el método de transferencia estos hologramas binarios
presentan algunos inconvenientes ya que si se desea superponer hologramas para
hacer un filtro compuesto cada vez hay más zonas oscuras que bloquean las claras.
Por otra parte al ser las franjas representadas por rectángulos, es muy difícil elegir
arbitrariamente la frecuencia y dirección de la portadora.
El advenimiento de nuevos métodos de transferencia permitió la representación de
CGHs directamente en niveles de gris [C. Iemmi, S. Ledesma, J. Campos, M.
Villarreal, Gray-level computer generated hologram filters for multiple-object
correlation, Appl. Opt. 39, 1233 (2000)]. Esto hace que sean verdaderos
interferogramas que pueden sumarse y que permiten elegir arbitrariamente, a través
de la portadora, donde irá a parar la correlación. Esta técnica permite implementar el
procesado multicanal y es aplicable ya sea utilizando material fotográfico (mediante
máquina de copiado AGFA Alto) o directamente en pantallas de cristal líquido. A
continuación se ve una porción ampliada de un holograma generado en niveles de gris
296
Óptica de Fourier Claudio Iemmi
HOLOGRAFÍA DIGITAL
297
Óptica de Fourier Claudio Iemmi
dicho haz sin la presencia del orden cero de difracción, ni del haz conjugado. Esto es,
el holograma sería lo que se conoce como un kinoform.
Existe una técnica de interferometría, llamada por corrimiento de fase [K. Creath
Phase-measurement Interferometry Techniques, E. Wolf, Progress in Optics XXVI
Cap. 5, (1988)], que permite obtener la distribución de amplitud y fase del haz de
interés. En este método se utiliza un interferómetro por división de amplitud, tipo
Michelson, en uno de cuyos brazos hay un espejo que proporciona el haz de
referencia y en el otro se halla el objeto que proporciona la onda a la cuál se le quiere
medir la distribución de fase. La imagen de interferencia es adquirida por una cámara
CCD. A diferencia de la técnicas interferométricas usuales en este método la fase del
haz de referencia es desplazada, en forma controlada, con respecto a la fase del haz a
testear, y la intensidad resultante es registrada para diferentes valores de corrimiento
relativo. Algunas técnicas cambian la fase escalonadamente entre mediciones y otras
integran la intensidad a medida que la fase es corrida. La cantidad de corrimientos
también varía.
Se necesita un mínimo de tres mediciones, ya que son tres las incógnitas a develar,
surgidas de la ecuación para la intensidad en el patrón de interferencia:
I I 0 [1 cos( )]
I1 x, y I 0 x, y 1 cos x, y
1
2
I 2 x, y I 0 x, y 1 cos x, y I 0 x, y 1 sin x, y
I3 x, y I 0 x, y 1 cos x, y I 0 x, y 1 cos x, y
3
2
I 4 x, y I 0 x, y 1 0 cos x, y I 0 x, y 1 0 sin x, y
298
Óptica de Fourier Claudio Iemmi
I 4 ( x, y ) I 2 x, y I1 ( x, y ) I 3 ( x, y )
2 2
( x, y )
2I0
I 4 ( x, y ) I 2 x, y
( x, y ) arctan
I1 ( x, y ) I 3 ( x, y )
299
Óptica de Fourier Claudio Iemmi
el plano del holograma, viene dada por U ( x, y) exp i ( x, y) , luego es posible
reconstruir el objeto, numéricamente, aplicando dicha integral de propagación hasta el
plano objeto. La misma resulta:
U ,
iA
d
exp i
d
( x, y )
i 2
2 2 U ( x, y) exp
d
x y dx dy
2 2
cumple la condición d 3
4
x 2 y 2 .
2
Para una dimensión
muestreo en cada eje de coordenadas sobre el plano de la CCD, luego las variables
continuas y deben reemplazarse por r y s, donde r y s son números enteros.
300
Óptica de Fourier Claudio Iemmi
2 2 2 2
U r , s exp i
d
r s .
kr ls
. k 0 l 0 U (k , l ) exp i
N 1 N 1
d
k 2 x 2 l 2 y 2 exp i 2
N N
301