Chapter 2
Chapter 2
Optique de Fourier
Principe
z
= z
2π
k= k2x + k 2y + k 2z = λ
1
pz = – p 2x – p 2y
λ2
U(x,y,z) = A exp[i2π(px x + py y + pz z)]
Angles:
sinθ x = kx/k = λ px
sinθ y = ky/k = λ py
Exemple (voir figure suivante): objet périodique, période Λ x = 1/px
⇒ déflexion d'une onde plane sinθ x = λ px = λ/Λ x
x k
illumination θx
objet périodique
Λ x = 1/px
Spectre angulaire
illumination
y
U(x,y)
∞ ∞
^ (p ,p ) =
U x y ∫ ∫ U(x,y) exp[–i2π (px x + py y)] dx dy
–∞ –∞
∞ ∞
U(x,y,z) = ∫ ∫ ^ (p ,p ) exp[i2π (p x + p y)] exp[ik z] dp dp
U x y x y z x y
–∞ –∞
1
avec kz = 2π – p 2x – p 2y
λ2
1
Condition pour les ondes qui se propagent: p2x + p2y <
λ2
1
Ondes évanescentes: p2x + p2y >
λ2
kz = iα ⇒ U
^ (p ,p ) exp[i2π (p x + p y)] exp[-α z]
x y x y
x0 x
Système
optique
y0 y
^ (p ,p ) = h
^ (p ,p ) U
^ (p ,p )
U x y x y 0 x y (multiplication)
∞ ∞
U(x,y) = ∫ ∫ h^ (px,p y) U
^ (p ,p ) exp[i2π (p x + p y)] dp dp
0 x y x y x y
–∞ –∞
1
h^ (px,p y) = exp i2π z – p 2 – p2
y
λ2 x
1 1 z exp[ikr]
h(x,y) = – 2π ik – r r
r
r2 = x2 + y2 + z2
1 z exp[ikr]
Si r >> λ ⇒ Rayleigh-Sommerfeld: h(x,y) = iλ r r
où z/r est le facteur d'obliquité.
∞ ∞
U(x,y) = ∫ ∫ U 0(x0,y0) h(x – x0, y – y0) dx0 dy0 ⇒
–∞ –∞
1 ∞ ∞ exp[ikr10] z
U(x,y;z) = iλ ∫ ∫ U 0(x0,y 0)
–∞ –∞ r10 r10 dx 0 dy 0
Approximation de Fresnel:
1 z
z – p 2x – p 2y =λ 1 – λ2 p 2x – λ2 p 2y
λ2
z λ2 p2x λ2 p2y λ4 p4x λ4 p4y
= λ (1 – 2 – 2 + 8 + 8 ........
| Fresnel |
Fresnel:
∞ ∞
1 ik
U(x,y;z) = iλz exp[ikz] ∫ ∫ U 0(x0,y 0) exp2z[(x – x0)2 + (y – y 0)2] dx 0 dy 0
–∞ –∞
1 ik
U(x,y;z) = iλz expikz + 2z (x2 + y2)
∞ ∞
ik
∫ ∫ U 0(x0,y 0) exp2z (x20 + y 20) exp[–i2π (px x 0 + py y 0)] dx0 dy 0
–∞ –∞
1 ik ik
U(x,y;z) = iλz expikz + 2z (x2 + y2) TF U 0(x0,y 0) exp2z (x20 + y 20)
Approximation de Fraunhofer
π
z >> λ (x20 + y20)
ik
exp2z (x20 + y 20) → 1
1 ik
U(x,y;z) = iλz expikz + 2z (x2 + y2)
∞ ∞
∫ ∫ U 0(x0,y0) exp[–i2π (px x0 + py y0)] dx0 dy0
–∞ –∞
1 ik
U(x,y;z) = iλz expikz + 2z (x2 + y2) TF{U 0(x0,y 0)}
théorème de translation:
∫
TF{U(x – x 0)} = dx U(x – x 0) e– i2π x px
^ (p )2
TF{U(x – x 0)}2 = TF{U(x)}2 = U invariant
x
théorème de similitude
TF{U(ax)} ∫
= U(ax) e– i2π x px dx
⌠ x' 1 1 ^ px
= U(x') exp[– i2π p ) dx' =
⌡ a x a a U( a )
fonction de Dirac:
π x2
δ(x) = lim (1/a) exp – a2
a→0
∫ δ(x) dx= 1
π (x2 + y2)
δ(x, y) = lim (1/a2) exp –
a→0 a2
∫ δ(x, y) dx dy = 1
^ (p )
U x
Gauss, Dirac:
a δ(p x )
U(x)
1 1
1/e 1/e
x px
1 a
0 0 1
π π π
U(x) = e – π x2 ^ (p ) = e– π px2
U x
U(x) = 1 ^ (p ) = δ(p )
U x x
(théorème de similitude)
^ (p )
U(x) U x
2
1
1
1/2 1/2
x px
0
0 g=1/p0 –p 0 0 +p0
U(x) = 1 + cos(2π p0 x) = ^ (p ) =
U x
1 1 1
1 + 2 (ei2π p0 x + e–i2π p0 x ) 2 δ(px+p0) + δ(px) + 2 δ(px–p0)
(théorème de translation)
Couples de transformée de Fourier
U(x) ^ (p )
U x
a
1
x px
–a/2 0 a/2
–2/a 0 1/a 3/a
U(x,y)=1 x
D
p
2
1.22 p 2 = p2
x +p y
0
D
U1 U2
f (distance focale)
r2 + f 2 = (f + δ) 2 = f 2 + 2fδ + δ 2
1 1
approximation quadratique (δ << f): δ ≅ 2f r2 = 2f (x2 + y2)
k
U 1 = eikz → U 2 = eik(z – δ) = eikz exp–i 2f (x2 + y2)
k
U 2(x,y) = tL(x,y) U 1(x,y) tL(x,y) = P(x,y) exp –i 2f (x2 + y2)
pour le suivant:
x1 xA x2
αx = λ/Λ = λ px
U1
z x2 = f α x = λ f px
αx
Λ U2 sinα ≈ tanα ≈ α
d1 f
amplitudes
x 1, y1 complexes
xA , yA x2, y2
U1 → U2
U1 UA U2
z
d1 = 0
d1 d2 = f xi = (x i, yi, 0)
k
U A(xA) = U1(xA) tL(xA) = U1(xA) P(xA) exp–i 2f (xA2)
1 ik
U(x2;z) = iλz expikz + 2z (x22)
k ik
x TF U 1(xA) P(xA) exp–i 2f (xA2) exp2z (xA2)
avec z = f
1 ik
U(x2;f) = iλf expikf + 2f (x22) TF {U 1(xA) P(xA) }
Si le diamètre de la pupille est plus grand que l'objet on obtient pour l'intensité:
exp[ikf] k ^ x2
i (x22) U
1 λf
U 2(x2) = exp 2f
iλf
1 ^ x22
* U
(λf)2 1λf
I(x2) = U2(x2) U 2 (x2) =
x
Fréquences spatiales: p = λf2
x 1, y1 xA , yA x2, y2 amplitudes
complexes
U1 UL UA U2
z U1 → U2
xi = (x i, yi, 0)
d1 d2 = f
Admettons que P(xA) = 1 partout (diamètre de la lentille est plus grand que le
faisceau).
^ (p) = U
^ (p) h
^ (p)
UL 1 h^ (p): fonction de transfert
^ (p) exp[ikd ] exp –iπλd p 2 x
=U1 1 { 1 } avec p = λf2
Selon a):
exp[ikf] k x
i x 22 U^ 2
U 2(x2) = exp 2f L λf
iλf
k x x2 2
^ 2 exp
exp[ik(f + d1)]
i x 22 U
= exp 2f 1 λf –iπλd 1 λf
iλf
exp[ik(f+d1)] k d1 ∞ x x
U 2(x2) = exp i 1 – x 2 ∫ U 1(x1) exp –i2π 2 1 d 2x1
2
iλf 2f f –∞ λf
1 ^ x22
I(x2) = U2(x2) U *2 (x2) = U
(λf)2 1λf
exp[ik2f] ^ x2
pour d1 = f U 2(x2) = U1λf
iλf
x 1, y1 x2, y2
U1 U2
z
d1 d2 = f
∞ d x x
U 2(x2) = h0(x2) ∫ U 1(x1) P(x1 + f1 x 2) exp–i2π 2λf 1 d2x1
–∞
eik(f+d1) k d
avec h0(x2) = exp i 1 – 1 x22
iλf 2f f
c) objet derrière la lentille
illumination
x1, y1 x 2, y2
U1 UA
U2
z
d1
d2 = f
x 2
Illumination avec une onde sphérique: US (x1) = exp–ik 2d1
1
f
U A(x1) = US (x1) U 1(x1) P( d x 1)
1
1 ik
U(x,y;z) = iλz expikz + 2z (x22)
x 2 f ik
x TF exp–ik 2d1 U1(x1) P( d x 1) exp2z (x12)
1 1
aved z = d1:
∞
exp[ikd ] x 2 ⌠ x x
U 2(x2) = iλd 1 expik 2d2 U (x ) P( f x ) exp –i2π 2 1 d 2x
1 1 d1 1 λd1 1
1 1 ⌡
–∞
Si tout l'objet est illuminé on obtient:
1 ^ x2 2
I(x2) = U2(x2) U *2 (x2) =
(λd1)2 1λd1
U
filtrage optique (cohérent):
x1 x2 = λ f p x3
^
F(p)
^ ^
U1 U1 U3 U3
f f f f
plan de Fourier x1
∫
^ (p) = ^F(p)·U
U ^ (p) ⇒ U 3(x) = F(x') U1(x – x') d2x' = F(x) ° U 1(x)
3 1
multiplication convolution
plan de Fourier plan image
x1 xP x2 = xI image géométrique:
x I = – M x1 = M xO
agrandissement:
z
M = d I/dO
Q loi de conjugaison:
dO d2 = d I 1 1 1
+
xO dI dO = f
k
2
exp[ikdO]
devant la lentille: UP (xP ) = exp i (x – xQ
)
iλdO 2dO P
k
après la lentille: U'P (xP ) = UP (xP ) P(xP ) exp–i 2f (xP 2)
exp[ikdI] ∞ k
2
U Q(xI,x Q) = iλd ∫ U' (x
P P ) exp
i 2d (xI – xP
)
d 2xP
I –∞ I
⇒
eik(d O+dI) xQ2 xI2 ⌠ xI xQ
U Q(xI,x Q)=– λ 2d d expik + P(x ) exp –ikx – d d 2xP
2dO 2dI P
P d
I O
O I ⌡
MxQ2 x – M x
eik(d O+dI) ^ I Q
U Q(xI,x Q) = – λ 2d d expik P λd
O I 2f I
z
Q
f dI
dO
xO
xI – M xQ
réponse impulsionnelle h: U Q(xI,xQ) ∝ ^P λd ∝ h(xI – M xQ)
I
⌠ x
1 xP ·x I 2 1 ^ I
P(xP ) exp–i2π d xP =
A P λdI λdI
h(xI) = P
A P λdI λdI
⌡
AP MxO2
γ = λd exp ik dI + dO +
O 2f
représentation spectrale: (convolution → multiplication)
^ (p) = ^h(p) · O
U ^ (Mp)
I filtrage
⇒ ^I = h^ O
^ ^h O
*
^
^I (p) = H
^ (p) · ^I (Mp)
représentation spectrale: I O filtrage
fonction de transfert incohérente: autocorrélation de la
fonction pupillaire P(xP )
fonction
pupillaire
xP
-D/2 D/2
^h(p) ^
H(p)
fonction de transfert 1
fonction de transfert incohérente
cohérente
P P
sinθ / λ 2 sinθ / λ
D / (2λz) D / (λz)
Remarques
incohérent cohérent
^ (0.0) = 1 ^h(p) libre
H
^ (–p) = H
H ^ *(p)
^ (p)| ≤ H
|H ^ (0.0)
Exemples
a) O(x) = cos(2π p0 x)
1 1
b) O(x) = |cos(2π p0 x)| ≈ 2 + 2 cos(4π p 0 x) + ...
1
a) O(x) = cos(2π p0 x) = 2 {exp(i2πp0x) + exp(–i2πp0x) }
^ (p ) = 1 δ(p +p ) + 1 δ(p –p )
U x x 0 x 0
2 2
1 1 1 1
b) O(x) = 2 + 2 cos(4π p 0 x) = 2 + 4 {exp(i4πp0x) + exp(–i4πp0x) }
^ (p ) = 1 δ(p +2p ) + 1 δ(p ) + 1 δ(p –2p )
U x x 0 x x 0
4 2 4
cohérent incohérent
^ ^ ^ O ^
O IO = O * (autocorrélation)
1/2 1/2
1/2
1/4 1/4
px px
^ ^
h H
px px
-p 0 +p -2p0 0 +2p0
0 0
^I = Oh
^^ Oh^^ ^ ^ ^
* I = H IO
1/2 1/2
1/4 1/4
<1/4 <1/4
px px
y y
P P
Forme de la pupille Autocorrélation
(normalement ronde)
xP xP
λ dI p y
aire totale λ dI p x
aire commune
∫∫ dx P dy P
' '
ˆ (x P ) ∝ ∩( xP )
H
∫∫ dx P dyP
' '
∩(0,0)
∫∫
∩( xP )
{[ ( ) (
exp ik W x 'P − W x 'P − x P )]}dx'Pdy'P
ˆ (x P ) ∝
H
∫∫ dx P dy P
' '
∩(0,0)
avec xP = λ dI p
xP " x P " "
∫∫ + − x P − dx P dy P
"
exp ik W P
x W
2 2
ˆ (x P ) ∝ ∩( xP )
H
∫∫ dx Pdy P
" "
∩(0,0)
xP
avec x"P = x'P − et x P = λd I p
2
^ (p)| ^
|H avec aberrations ≤ |H(p)|sans aberrations
preuve
substitution:
x x
( )
X x"P = exp ikW x"P + P
2
et ( )
Y x"P = exp − ikW x "P − P
2
X2 = Y2 =1
on trouve:
2 2
∫∫ XYdx P dy P ∫∫ dx P dy P
" " " "
2 ∩( x P ) ∩( x P ) 2
ˆ ( xP )
H = ≤ = H
ˆ (x P )
∫∫ dx P dy P ∫∫ dx P dy P
avec aberration " " " " sans aberration
∩(0,0) ∩(0,0)
^
IHI
1
-p p __
D
λz
Y électronique
x
illumination
objective à dédecteur
réseau avec
tester fente d'analyse
fréquence variable
Spatiale:
∞ ∞
^ (p
U x,p y) = ∫ ∫ U(x,y) exp[-i2π (px x + py y)] dx dy
–∞ –∞
∞ ∞
U(x,y) = ∫ ∫ ^ (p ,p ) exp[i2π (p x + p y)] dp dp
U x y x y x y
–∞ –∞
^ (p ,p )
TF{U(x,y)} = U x y
1) Théorème de linéarité
TF{a U1(x,y) + b U 2(x,y)} = a TF{U 1(x,y)} + b TF{U2(x,y)}
2) Théorème de similitude
1 ^ px py
TF{U 1(ax,by)} = |ab| Ua ‚ b
3) Théorème de translation
^ (p ,p ) exp[-i2π (p a + p b)]
TF{U 1(x – a, y – b)} = U x y x y
4) Théorème de Parseval
∞ ∞ ∞ ∞
∫ ∫ |U(x,y)|2 dx dy =∫ ∫ ^ (p ,p )|2 dp dp
|U x y x y
–∞ –∞ –∞ –∞
5) Théorème de convolution
∞ ∞
TF{ ∫ ∫ ^ (p ,p )
U 1(x0,y 0) U 2(x – x0, y – y0) dx0 dy 0 = U1 x y
–∞ –∞
^ (p ,p )
U2 x y
5) Théorème d'autocorrélation
∞ ∞
TF{ ∫ ∫ ^ (p ,p )|2
U(x 0,y 0) U *(x0 – x, y0 – y) dx0 dy 0 = |U x y
–∞ –∞
Convolution
∞
g( x) = f ( x) o h( x) = ∫ f ( x' )h( x − x' )dx'
−∞
∞
TF {g ( x ) } = gˆ ( p ) = ∫ g ( x ) exp[ −2 πixp] dx
–∞
∞ ∞
gˆ ( p ) = ∫ exp[ −2 πixp ] ∫ f ( x' ) h ( x − x ' )dx' dx
–∞ −∞
∞ ∞
gˆ ( p ) = ∫ ∫ h ( x − x ' ) exp[ −2 πixp ] dx f ( x ' ) dx'
– ∞ −∞
substitution w = x – x'
∞ ∞
gˆ ( p ) = ∫ ∫ h ( w ) exp[ −2 πiwp] exp[ −2 πix'p ] dw f ( x' ) dx'
– ∞ −∞
∞ ∞
gˆ ( p ) = ∫ f ( x ' ) exp[ −2 πix'p ] dx' ∫ h ( w ) exp[ −2 πiwp] dw
–∞ −∞
gˆ ( p ) = ˆf ( p ) h
ˆ (p )
g ( x ) = f ( x ) o h( x ) convolution
gˆ ( p ) = ˆf ( p ) h
ˆ (p )
−∞
∞
TF {g ( x ) } = gˆ ( p ) = ∫ g ( x ) exp[ −2 πixp] dx
–∞
∞ ∞
gˆ ( p ) = ∫ exp[ −2 πixp] ∫ f ( x' ) h ( x ' − x )dx' dx
*
–∞ −∞
∞ ∞
g ( p ) = ∫ ∫ h* ( x ' − x ) exp[ −2 πixp ] dx f ( x' ) dx'
ˆ
– ∞ −∞
substitution w = x' – x
∞ ∞
gˆ ( p ) = ∫ ∫ h* ( w ) exp[ +2πiwp] exp[ −2 πix'p ]dw f ( x ' ) dx'
– ∞ −∞
∞ ∞
gˆ ( p ) = ∫ f ( x ' ) exp[ −2 πix'p ] dx' ∫ h ( w ) exp[ +2πiwp] dw
*
–∞ −∞
∞ ∞
∫ f ( x ' ) exp[ −2 πix'p ] dx' ∫ [ h ( w ) exp[ −2 πiwp]] dw
*
gˆ ( p ) =
–∞ −∞
gˆ ( p ) = ˆf ( p ) h
ˆ * ( p)
g ( x ) = f ( x )∗ h ( x ) corrélation
gˆ ( p ) = ˆf ( p ) h
ˆ * ( p)
g ( x ) = f ( x )∗ f ( x ) autocorrélation
2
gˆ ( p ) = ˆf ( p ) ˆf * ( p ) = ˆf ( p )