0% ont trouvé ce document utile (0 vote)
13 vues28 pages

Chapter 2

Le document traite de l'optique de Fourier, en se concentrant sur la propagation de la lumière dans l'espace libre et les transformations de Fourier associées. Il décrit les principes fondamentaux, les ondes planes, le spectre angulaire, ainsi que les approximations de Fresnel et de Fraunhofer. De plus, il aborde les effets des lentilles minces sur les fronts d'onde et la transmission d'amplitude complexe.

Transféré par

h.a.dresser93
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
13 vues28 pages

Chapter 2

Le document traite de l'optique de Fourier, en se concentrant sur la propagation de la lumière dans l'espace libre et les transformations de Fourier associées. Il décrit les principes fondamentaux, les ondes planes, le spectre angulaire, ainsi que les approximations de Fresnel et de Fraunhofer. De plus, il aborde les effets des lentilles minces sur les fronts d'onde et la transmission d'amplitude complexe.

Transféré par

h.a.dresser93
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

2.

Optique de Fourier

2.1 Propagation de la lumière dans l'espace libre

Principe

Une onde quelconque peut


être représentée par une
superposition d'ondes planes

z
= z

Onde planes ⇔ fréquences spatiales

Onde plane: U(x,y,z) = A exp[i(kx x + ky y + kz z)]


k= k2x + k 2y + k 2z = λ

Fréquences spatiales: px = kx/2π et py = ky/2π

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

Construction avec la Sphère d'Ewald (voir chapitre 4)

Système linéaire, fonction de transfert

x0 x

Système
optique

y0 y

Système général: U 0(x0,y 0) ⇒ U(x,y)


Propagation libre: U 0(x0,y 0;z=0) ⇒ U(x,y) = U 0(x0,y 0;z=d)

Réponse impulsionnelle h(x,y):


Objet: δ(0,0) ⇒ Image: h(x,y)

Pour la propagation libre on peut écrire:


∞ ∞
U(x,y) = ∫ ∫ U 0(x0,y0) h(x – x0, y – y0) dx0 dy0 (convolution)
–∞ –∞

^ (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
–∞ –∞

h(x,y): réponse impulsionnelle h^ (px,p y): fonction de transfert


Propagation libre:
Distance z = d (voir spectre angulaire)

 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

r10 = (x – x0)2 + (y – y 0)2 + z2

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 |

h^ (px,p y) = exp[ikz] exp[–iπλz (p2x + p 2y)]


1  x2 + y2 
h(x,y) = iλz exp[ikz] expik 2z 
 
∞ ∞
U(x,y) = ∫ ∫ U 0(x0,y0) h(x – x0, y – y0) dx0 dy0 ⇒
–∞ –∞

Fresnel:
∞ ∞
1  ik 
U(x,y;z) = iλz exp[ikz] ∫ ∫ U 0(x0,y 0) exp2z[(x – x0)2 + (y – y 0)2] dx 0 dy 0
–∞ –∞  

L'approximation de Fresnel peut être reécrit comme une transformation de


Fourier.
x y
Avec px = λz et py = λz on trouve:

1  ik 
U(x,y;z) = iλz expikz + 2z (x2 + y2)
 
∞ ∞
 ik 
∫ ∫ U 0(x0,y 0) exp2z (x20 + y 20) exp[–i2π (px x 0 + py y 0)] dx0 dy 0
 
–∞ –∞

1  ik    ik  
U(x,y;z) = iλz expikz + 2z (x2 + y2) TF U 0(x0,y 0) exp2z (x20 + y 20)
    

Approximation de Fraunhofer

π
z >> λ (x20 + y20)

 ik 
exp2z (x20 + y 20) → 1
 

1  ik 
U(x,y;z) = iλz expikz + 2z (x2 + y2)
 
∞ ∞
∫ ∫ U 0(x0,y0) exp[–i2π (px x0 + py y0)] dx0 dy0
–∞ –∞
1  ik 
U(x,y;z) = iλz expikz + 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

= ∫ dx' U(x') e– i2π (x' + x0)px = U


^ (p ) e– i2π x0 p x
x

^ (p )2
TF{U(x – x 0)}2 = TF{U(x)}2 = U invariant
x

intensité de la diffraction Fraunhofer: invariante à la translation

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

source ponctuelle: fonction de Dirac bidimensionnelle

 π (x2 + y2) 
δ(x, y) = lim (1/a2) exp –
a→0  a2

 ∫ δ(x, y) dx dy = 1

∫ U(x, y) δ(x – x0, y – y0) dx dy = U(x0, y0)

∫ e– i2π(px x+py y) dxdy = δ(px,py), ∫ e+ i2π(px x+py y) dpx dpy = δ(x, y)


U(x,y) =1 ⇒ ^ (p ,py) = δ(p ,p )
U x x y
Couples de transformée de Fourier

^ (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) = e – π x2/a2 ^ (p ) = a e– π px2 a2


U x

U(x) = 1 ^ (p ) = δ(p )
U x x

(théorème de similitude)

réseau (de transmission) sinusoïdal:

^ (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

fente (fonction rectangulaire):

U(x) ^ (p )
U x
a
1

x px

–a/2 0 a/2
–2/a 0 1/a 3/a

|x|< a/2: U(x) = 1 ^ (p ) = a sin(a π px)


U
|x|> a/2; U(x) = 0 x (a π p ) x

U(x) = rect(x/a) ^ (p ) = a sinc(a p )


U x x
(théorème de similitude)
sin(π p)
sinc(p) = (π p) (définition)

pupille circulaire: fonction Airy:


^
y U (p ) 2
U(x,.y)=0 πD
4

U(x,y)=1 x
D
p

2
1.22 p 2 = p2
x +p y
0
D

x2 + y2 < D/2 : U(x, y) = 1


^ (p , p ) = π D2 2 J1(π D p)
U x y 4 (π D p)
x2 + y2 > D/2 : U(x, y) = 0
2 J1(π p)
airy(p) = (π p)
2.2 Transformation de Fourier par une lentille
Lentille mince

lentille mince idéale: transformation d'un front d'onde:


δ
onde sphérique
f
F (foyer) r F
z f z

U1 U2
f (distance focale)
r2 + f 2 = (f + δ) 2 = f 2 + 2fδ + δ 2

1 1
approximation quadratique (δ << f): δ ≅ 2f r2 = 2f (x2 + y2)

changement de phase (front d'onde) par la lentille:

 k 
U 1 = eikz → U 2 = eik(z – δ) = eikz exp–i 2f (x2 + y2)
 

transmission d'amplitude complexe de la lentille


(élément mince)

 k 
U 2(x,y) = tL(x,y) U 1(x,y) tL(x,y) = P(x,y) exp –i 2f (x2 + y2)
 

fonction pupillaire: P(x,y) = |P|(x,y) eiρ(x,y)

|P|(x,y): transmission de l'ouverture, ρ(x,y): aberrations (de phase)

pour le suivant:

lentille sans aberrations (approximation paraxiale), ouverture idéale:

à l'intérieure de l'ouverture: P(x,y) = 1, ailleurs: P(x,y) = 0


transformation de Fourier (par lentille)
principe: 1) diffraction Fraunhofer: fréquences spatiales → directions
2) lentille: directions → points dans le plan focal

x1 xA x2
αx = λ/Λ = λ px
U1
z x2 = f α x = λ f px
αx
Λ U2 sinα ≈ tanα ≈ α
d1 f

a) objet devant la lentille (d1 = 0)

amplitudes
x 1, y1 complexes
xA , yA x2, y2
U1 → U2

U1 UA U2
z
d1 = 0

Objet (x1,y 1) = (x A,y A)

d1 d2 = f xi = (x i, yi, 0)

xi2 = xi2 + yi2

 k 
U A(xA) = U1(xA) tL(xA) = U1(xA) P(xA) exp–i 2f (xA2)
 

voir page 2-5, diffraction Fresnel reécrite comme une transformation de


Fourier:

1  ik 
U(x2;z) = iλz expikz + 2z (x22)
 
  k   ik  
x TF U 1(xA) P(xA) exp–i 2f (xA2) exp2z (xA2)
    
avec z = f

1  ik 
U(x2;f) = iλf expikf + 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  ^ x22
* U 
(λf)2  1λf
I(x2) = U2(x2) U 2 (x2) =

x
Fréquences spatiales: p = λf2

b) objet devant la lentille (d 1)

x 1, y1 xA , yA x2, y2 amplitudes
complexes

U1 UL UA U2
z U1 → U2

xi = (x i, yi, 0)
d1 d2 = f

xi2 = xi2 + yi2

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 


= exp  i  1 –  x 2  ^  2
U
iλf  2f  f  2  1λ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  ^ x22
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  

limitations par les dimensions finies de l'ouverture de la lentille (vignetage)

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 expikz + 2z (x22)
 

  x 2  f  ik  
x TF exp–ik 2d1  U1(x1) P( d x 1) exp2z (x12)
  1  1  

aved z = d1:


exp[ikd ]  x 2  ⌠  x x
U 2(x2) = iλd 1 expik 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

objet: U1(x) filtre: ^F(p) image: U3(x)


Microscopie à contraste de phase (champs noir) ϕ(x) << 1
(F. Zernike, 1935, Prix Nobel 1953)

Problème: visualiser un objet de phase, par exemple pour détecter des


inhomogenités dans un bloc de verre.
a) filtre d'amplitude:
objet
U 1(x) = U 0 eiϕ(x) ≅ U0 [1 + iϕ(x)], I1(x) = U 02
inhomogène
x2 ^ (p) = U [δ(p) + iϕ
U ^ (p)],
1 0
^F(p)·U
^ (p) = i U ϕ
^
1 0 (p)
^
F(p) U 2(x) = i U0 ϕ(x), I2(x) = U 02 ϕ 2(x)

b) filtre de phase: ^F(p):= i δ(p)


objet
U 1(x) ≅ U0 [1 + iϕ(x)], I1(x) = U 02
inhomogène
x2 ^ (p) = U [δ(p) + iϕ
U ^ (p)],
1 0
^F(p)·U
^ (p) = i U [δ(p) + ϕ
^ (p)]
0
^
F(p) U 2(x) = i U0 [1 + ϕ(x)], I2(x) ≅ U02 [1 + 2ϕ(x)]
2.3 Formation d'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

diffraction Fresnel pour une source ponctuelle:

Object: UQ(–x1) = UQ(xO) = δ(xO – xQ)

 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 expik  +    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 expik  P  λd 
O I  2f   I 
 

élargissement de l'image géométrique d'une source ponctuelle par la diffraction


^P(x /λd ) sur l'ouverture de la lentille P(x )
I I P
x1 xP xI
P(x P ) ^
P h(x I – Mx Q )

z
Q

f dI
dO
xO

pupille P(xP ) circulaire ⇒ h(xI) ∝ ^P(xI/λdI) = airy(D|xI|/λdI)


(diamètre D) résolution: h(xI) = 0 ⇒ |xI| = 1.22 λdI/D

Formation d'image cohérente, superposition cohérente:

superposition cohérente au plan image: U I(xI) = ∫ U 1(xQ) UQ(xI,xQ) d2xQ

 
xI – M xQ
réponse impulsionnelle h: U Q(xI,xQ) ∝ ^P  λd  ∝ h(xI – M xQ)
 I 

U I(xI) = ∫ O(xO)h(xI – MxO) d2xO


image cohérente convolution

fonction pupillaire: réponse impulsionnelle:


P(xP ), ∫ |P(xP )|2 d 2xP = A P h(xI), ∫ |h(xI)|2 d2xI = 1

⌠  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  

fonction objet: O(xO) = γ U1(–x1) = γ U1(xO)

AP   MxO2  
γ = λd exp ik dI + dO +  
O   2f  
représentation spectrale: (convolution → multiplication)

spectre de l'objet: spectre de l'image:


^ (p) = ∫ O(x ) exp{–i2π p·x } d 2x
O ^ (p) = ∫ U (x ) exp{–i2π p·x } d 2x
U
O O O I I I I I

^ (p) = ^h(p) · O
U ^ (Mp)
I filtrage

U=hoO⇒ I = (h o O) (h o O)* ⇒ ^I = TF(h o O) * TF(h o O)

⇒ ^I = h^ O
^ ^h O
*
^

fonction de transfert (cohérente):

^ (x /λd )} = λdI P(– x )


^h(p) = ( A λd )–1 TF{P xP
P I I I
AP P p = λd
I

Formation d'image incohérente, superposition incohérente:

superposition incohérente au plan image, addition des intensités:

II(xI) = ∫ |O(x O)|2 |h(x I – M xO)| 2 d 2xO

II(xI) = ∫ IO(xO)H(xI – M xO) d2xO


image incohérente convolution

réponse impulsionnelle incohérente: ^ (x /λd )|2


H(xI) = |h(xI)|2 ∝ |P I I

H(xI) est proportionnel au spectre de puissance de la pupille

objet incohérent (intensité): IO(xO) = |O(xO)|2

^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 )

^ (p) = TF{|h(x )|2} ∝ TF{P


H I
^P^ *} ∝
∫ P(x'P) P*(x'P – xP) d2x'P
^ (p) est proportionnel à l'autocorrélation de la pupille
H

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

OTF: Optical Transfer Function = H^ (p)

MTF: Modulation Transfer Function ^ (p)|


= |H

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

a) L'illumination cohérente est meilleure


b) L'illumination incohérente est meilleure
Exemple a)

cohérent incohérent

^ ^ ^ O ^
O IO = O * (autocorrélation)
1/2 1/2
1/2
1/4 1/4

px px

-p0 0 +p0 -2p 0 0 +2p0

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

-2p 0 0 +2p0 -2p 0 +2p


0 0
Interprétation géométrique de la fonction de transfer

y y
P P
Forme de la pupille Autocorrélation
(normalement ronde)

xP xP

λ dI p y

aire totale λ dI p x

aire commune

Système sans aberrations: ^ (p) ∝ aire commune


H aire totale

∫∫ dx P dy P
' '

ˆ (x P ) ∝ ∩( xP )
H
∫∫ dx P dyP
' '

∩(0,0)

Systèmes avec aberrations

système avec aberrations: P(xP ) = |P(xP )| exp{ikW(x P )}


système sans aberrations: W(xP ) = 0
W: l'erreur du chemin optique

∫∫
∩( 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

exemples d'aberrations (monochromatique):

défocalisation W(x,y) = a(x2 + y2)


astigmatisme b x2 + c y2
coma dx(x2 + y2)
aberration sphérique e(x2 + y2)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

avec l'inégalité de Schwarz


2
∫∫ XYdx P dy P ≤ ∫∫ X 2 dx "P dy "P ∫∫ Y
" " 2
dx "P dy "P

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)

Mesure de la fonction de transfert (contraste)

Réseau en intensité: I = 1 + cos(2πpx)


⇒ I = 1 + m cos(2πpx) avec m < 1

^
IHI
1

-p p __
D
λz

^ (p)| = Imax – Imin = contraste


MTF = |H Imax + I min

Y électronique
x

illumination
objective à dédecteur
réseau avec
tester fente d'analyse
fréquence variable

G. Schröder, Technische Optik, Vogel Buchverlag Würzburg, 1990, p. 196


Exemples calculés de la MTF

aberrations sphériques W(r) = ar4 , r=0...1


(A) a = 0, (B) a= λ/4, (C) a= λ/2, (D) a= λ

système sans aberrations avec fonction pupillaire en forme d'anneau


(A) s 0/sm = 0, (B) s 0/sm = 0.25, (C) s 0/sm = 0.5, (D) s0/sm = 0.75
Transformation de Fourier

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| Ua ‚ 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 )

Valable pour TF et TF-1


Corrélation

g ( x ) = f ( x )∗ 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
*

–∞ −∞
∞ ∞
∫ 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 )

Valable pour TF et TF-1

Exemple: Objet O Intensité I = O O*


TF{ I } = TF{ O } * TF{ O } = autocorrélation

Vous aimerez peut-être aussi