Diffraction de la lumière

Cours 5

1
1 Introduction

2
 La diffraction est un phénomène propre
aux ondes (Huygens)
Vision très grossière du phénomène:

Écran opaque
a vaut quelques λ
a >> λ

a θ ~λ/a
a

Onde plane Onde plane Onde après le trou 3

progressive progressive
 Ce phénomène s’observe effectivement
avec la lumière

Faisceau lumineux Ce que prévoit Ce que prévoit la

incident= Cylindre l’optique géométrique théorie ondulatoire
et la théorie
corpusculaire

C’est ce qui est observé

si le trou est très petit
4
 Interprétation électromagnétique
L’onde incidente Si(M,t) excite les atomes de l’écran qui génèrent à leur tour un
champ Se(M,t)

A) Lorsque l’écran est « plein », le champ de l’autre côté de l’écran est tel
que: S(M,t) = Se(M,t)+Si(M,t) = 0. L’écran est opaque.

B) Lorsqu’on fait un trou dans l’écran, on retire des « sources » qui

généraient avant un champ St(M,t) (petite contribution au champ Se(M,t))

Après l’écran on a donc:

S’(M,t) = Se(M,t) + Si(M,t)-St(M,t) = -St(M,t)

Si le trou est tout petit, la portion retirée rayonnait avant un champ d’onde
quasi-sphérique.

Donc S’(M,t) = -St(M,t) est pratiquement un champ d’onde sphérique.

5
On retrouve ainsi ce que nous avions supposé sur les trous d’Young.
 Interprétation quantique





p = ℏk = ℏ(k x ex + k z ez )
Impulsion d’un photon, quantum du champ EM:



= ℏk (sin θex + cos θez )

Le trou dans l’écran permet de dire que la position du photon au niveau de

l’écran était –a/2 < x < a/2
L’incertitude sur la position est ∆x ~ a

Inégalité de Heisenberg:

Plus ∆x (et donc a) est petit, plus ∆px est grand et plus l’incertitude sur ө est
grande
plus la divergence de l’onde sortante sera grande.

6
2 Principe de Huygens-Fresnel

7
 Interprétation « Mécanique »
Une onde se propage par ébranlement de proche en proche des points
du milieu matériel (ou ici du champ).

Un point atteint par l’onde se met à vibrer et rayonne un champ d’onde

sphérique. La vibration est proportionnelle à la perturbation incidente

Chaque point de Σ émet une

petite onde sphérique

Onde incidente… Et après Σ?

La superposition de ces
ondelettes sphériques donne
8
une onde plane…
 Huygens-Fresnel: un « principe »?

On peut montrer (avec beaucoup de mathématiques…) que ce

« principe » découle en fait des équations de Maxwell.

Les équations de Maxwell sont des équations aux dérivées partielles,

qui décrivent localement le comportement du champ.

On peut cependant les traduire sous une forme intégrale. Un exemple

célèbre est le théorème de Gauss de l’électrostatique. Ce dernier fait
intervenir l’intégrale d’un champ sur une surface.

9
 Principe de Huygens-Fresnel (Formulation Moderne)
y Y
M(X,Y)
P(x,y)
x X
R(x,y)

Plan d’arrivée (z)

Plan de départ (z=0) Champ = S(X,Y,t)
Champ = s(x,y,t) =S(X,Y)e -iωt
=s(x,y)e -iωt

Intégrale de Kirchhoff
ikPM
e
S ( X , Y ) = ∫∫ s ( x, y ) K dxdy
PM
Champ dans le plan de départ
Champ en M d’une
K: Facteur de Kirchhoff, dépendant de
ondelette sphérique émise
10
l’inclinaison de PM par rapport à (Oz)
par P
 Interprétation
Si on connaît le champ dans un plan, on peut le déterminer
dans un autre plan situé à la distance D>>λ

Le champ en M est la somme d’une infinité de champs

élémentaires émis par les points P du plan de départ.

Ces champs élémentaires sont des ondelettes sphériques

d’amplitudes complexes initiales proportionnelles à celles
du champ en P.

11
3-Approximation de Fresnel
et de Fraunhofer

12
 Approximation de Fresnel
• On suppose que z >> X,Y,x,y > λ ALORS

Au dénominateur, PM~R~z

• En revanche, dans le terme de phase on a:

PM = z 2 + ( X − x) 2 + (Y − y ) 2
PM = z 2 + X 2 + Y 2 + x 2 + y 2 − 2 Xx − 2Yy
PM = R + x 2 + y 2 − 2 Xx − 2Yy
 ( x 2 + y 2 ) ( xX + yY ) 
PM ≈ R1 + 2
− 2

 2R R 
( x 2 + y 2 ) ( xX + yY )
PM ≈ R + −
2R R
13
 Approximation de Fresnel (Suite)
i
De plus K≈ (admis)
λ
L’intégrale de Huygens Fresnel devient alors…
ikR k ( x2 + y2 ) k ( Xx +Yy )
i e i −i
S ( X ,Y ) =
λ z ∫∫ s( x, y)e 2R
e R
dxdy
R~ z dans les conditions de Fresnel

Intégrale de Fresnel

ikR k ( x2 + y2 ) k ( Xx +Yy )
i e i −i
S ( X ,Y ) =
λ z ∫∫ s( x, y)e 2z
e R
dxdy
14
« terme ou noyau de Fresnel »
 Approximation de Fraunhofer
Situation dans laquelle le terme de Fresnel devient proche de 1, c’est à dire
La phase de k(x2+y2)/2z << π (1/100ème par exemple)

x2 y2
, << 1 Condition très drastique !!!
λ z λz

Alors l’intégrale de Fresnel devient:

k ( Xx +Yy )
i e ikR −i
S ( X ,Y ) =
λ z ∫∫ s( x, y)e R
dxdy

Intégrale de Fraunhofer
x2 y 2
, << 1
λz λ z 15
 Diffraction de Fraunhofer comme diffraction à « l’infini »
Si R ou z tendent vers l’infini, le noyau de Fresnel tend vers 1. L’intégrale de
Fraunhofer garde son sens si on fait apparaître les angles θx et θy.
M
y
X
s(x,y) R=OM
θy
O
z

X sin θ x Intégrale de Fraunhofer

On pose u= =
λR λ
Y sin θ y
v= =
λR λ
On obtient i e ikR −i 2π ( xu + yv )
S (u , v) =
λ z ∫∫ s ( x , y ) e dxdy
16
 Importance de la diffraction de Fraunhofer (I)
Intérêt des coordonnées u,v:
1) Ne fait plus apparaître z et X mais un angle
permet de décrire le
comportement « à l’infini » car coordonnées indépendante de la distance.

2) La diffraction de Fraunhofer est souvent qualifiée de diffraction à l’infini.

On regarde la répartition de l’amplitude du champ diffracté
ANGULAIREMENT

3) Permet de faire apparaitre l’intégrale de Fraunhofer comme la transformée

de Fourier du champ s(x,y)

i eikR X sin θ x
−i 2π ( xu + yv ) u= =
S ( X , Y ) = S (u , v) =
λ z ∫∫ s( x, y)e dxdy λR λ
Y sin θ y
i eikR v= =
S (u, v) = TF ( s) λR λ
λ z
17
 Approximation de Fraunhofer (III)

θy
z
s(x,y)

18
 Diffraction de Fraunhofer en
pratique

sin θ X
θX X
y Y u= ≈
≈
λ λ λf '
sin θY θY Y
v= ≈ ≈
λ λ λf '
S(u,v)
θy θy
z
s(x,y)
f’

d
Lentille mince, « parfaite »
(diamètre infini,
stigmatique…)

19
Résultat admis (pour l’instant)
Bilan

20
Calcul d’une figure de diffraction dans
l’approximation de Fraunhofer

21
 Un problème concret…
Plan juste après le
Onde incidente
diaphragme
CONNUE Amplitude,
phase,
INTENSITE
A TROUVER !!!

Diaphragme Écran
de forme extrêmement
quelconque loin ou au
MAIS foyer d’une
CONNUE lentille
22
 Résolution
1) On détermine s(x,y) le champ juste après le diaphragme

s(x,y) = t(x,y) si(x,y)

Où si(x,y) est le champ de l’onde incidente sur le plan du diaphragme,

Et t(x,y) est la transmission en amplitude du diaphragme.

Exemples: si t(x,y)=0 , le diaphragme est opaque en x,y.

si t(x,y)=1, le diaphragme est transparent en x,y
si |t(x,y)|<1 le diaphragme est partiellement transparent en x,y
si t(x,y) est complexe, le diaphragme déphase l’onde incidente
(fine plaque de verre..)

2) On calcule le champ diffracté sur l’écran en intégrant l’intégrale de

Fraunhofer
i e ikR −i 2π ( xu + yv )
S (u , v) =
λ z ∫∫ s ( x , y ) e dxdy
Avec X sin θ x Y sin θ y
u= = v= =
λR λ λR λ 23
 Utilisation des transformées de Fourier
On peut intégrer directement l’intégrale dans quelques cas simples.
Pour des cas plus complexes (et plus intéressants !) cela n’est pas possible.
Heureusement, l’arsenal mathématique des transformées de Fourier va nous aider.

i e ikR −i 2π ( xu + yv )
S (u , v) =
λ z ∫∫ s ( x , y ) e dxdy

i e ikR
S (u , v) = TF ( s )
λ z
I (u, v) = S (u, v) ∝ TF ( s )
2 2

24
 Rappel: propriétés des TF
+∞
Définition (1D) F (u ) = TF ( f ) = ∫ f ( x)e −i 2πux dx
−∞
+∞
i 2πux
f ( x) = TF −1 ( F ) = ∫ F (u ) e du
−∞
Propriétés

TF −1 (TF ( f )) = f
Linéarité TF (af ( x) + bg ( x)) = aF (u ) + bG (u )
translation TF ( f ( x − a )) = e −i 2πau F (u )
x
Similitude TF ( f ( )) = aF (au )
a
TF d’une TF TF (TF ( f ( x))) = f (− x) 25
 TF de fonctions de base
Fonction rectangle (« rect(x) »)

1 TF (rect ( x / a )) = a sin(πua) / πua = a sin c(πua)

-a/2 a/2

Fonction gaussienne TF (exp(−πx 2 )) = exp(−πu 2 )

26
 TF à deux dimensions
F (u , v) = ∫∫ f ( x, y )e − i 2π ( xu + yv ) dxdy

Cas particulier important fonction f à variables séparées:

:
Alors: f ( x, y ) = h( x ) g ( y )
F (u , v) = H (u )G (v)

27
Exemple 1: diffraction par une fente éclairée
par une onde plane
Y

si ( x, y ) = s0 eikz = s0 ( z = 0)
b X
t ( x, y ) = rect ( x / a)rect ( y / b)
s ( x, y ) = t ( x, y ) si ( x, y ) = s0
a

S (u, v) ∝ TF (rect ( x / a)rect ( y / b))

S (u, v) ∝ s0 ab sin c(πau ) sin c(πbv)
I (u, v) = S (u, v)
2

= I 0 a 2b 2 sin c 2 (πau ) sin c 2 (πbv)

Intensité maximale si u=0, v=0.

Intensité s’annule une première fois pour u=1/a et/ou v=1/b 28
 Exemple 1: diffraction par une fente éclairée
par une onde plane
Intensité maximale si u=0, v=0.
Intensité s’annule une première fois pour u=1/a et/ou v=1/b
Ceci correspond à des angles :
λ λ
sin θ x = sin θ y =
a b
Si les angles sont petits, la divergence du faisceau vaut:

2λ 2λ
∆θ x ≈ ∆θ y ≈
a b

29
TF 2D de deux fentes (cf TD)
Calcul effectué avec Scilab

30
TF 2D d’un « trou circulaire »

r = x2 + y2
r circ( x, a) = 0 r > a
circ( x, a) = 1 r < a

2a 2 J1 (2πaρ )
TF (circ( x, a )) = ρ = u 2 + v2
2πaρ
J1 : Fonction de Bessel

31
 TF 2D d’un « trou circulaire »
2 J1 (2πaρ )
TF (circ( x, a )) = ρ = u 2 + v2
2πaρ

J1 : Fonction de Bessel
La fonction de Bessel s’annule en plusieurs
points. Le premier zéro correspond à:

Position du premier zéro : ρ=1.22/2a

Angle correspondant au premier zéro:

2
TF (circ( x, a )) Θ= =1.22λ/2a

32
 Application: le réseau de diffraction (voir TD)
Succession de N fentes identiques parallèles entre elles et séparées de la même
distance d
x
Direction d’observation

O θx z

diaphragme

d
Y
Fente 1
Onde incidente

… b … X

33
a
 Application: le réseau de diffraction
(voir TD)
• Succession de N fentes identiques parallèles entre elles
et séparées de la même distance d
• Formellement: interférences à N ondes
• Peut être traité totalement à l’aide de la diffraction

Pics d’intensité (Interférences constructives) dans les directions telles que:

sin θ X − sin i = pλ / d = PNλ

Nombre de fentes par
unité de longueur
Direction incidente
34
 Application: le réseau de diffraction
(voir TD)


La direction du pic d’intensité dépend de la longueur d’onde

Un réseau disperse la lumière

Ce phénomène est à l’origine de l’irisation des disques CD

Spectroscopes à réseaux

manipulation d’impulsions lasers courtes

Les spectroscopes à réseau sont utilisés pour détecter les exoplanètes par
effet Doppler.

35
Filtrage spatial en lumière
cohérente

36
 Montage conceptuel
f’ f’ f’ f’

Onde plane

Objet diffractant

Plan
d’observation
Montage Montage
classique de la réciproque de la
diffraction de diffraction de
Fraunhofer Fraunhofer

Ici on observe la Figure de

diffraction de Fraunhofer Plan « de Fraunhofer » ou « de Fourier »

37
 Point de vue géométrique
f’ f’ f’ f’

Onde plane

Objet diffractant

Plan
d’observation

L’optique géométrique prévoit que le plan d’observation est le

plan image du plan de sortie de l’objet diffractant.

L’image à la même taille que l’objet, mais apparaîtra renversée.

38
 Point de vue ondulatoire
f’ f’ f’ f’

Onde plane Plan

d’observation

Objet diffractant

Fraunhofer Fraunhofer
s(x,y)=si(x,y)t(x,y)
S (u , v) ∝ TF ( s ( x, y ))
s' (x, y) ∝TF(TF(s(x, y)))
Le champ dans le plan
d’observation est bien similaire à s' (x, y) ∝ s(−x,−y)
39
celui du plan objet et renversé!
 Illustration sur un exemple simple
Objet diffractant de transmission sinusoidale et limité par une fente
Transmission de l’objet modulée par une fréquence spatiale 1/Λ

t ( x) = 1 rect ( x / a )(1 + cos(2πx ))

2
Λ
t ( x) = 1 rect ( x / a ) + 1 2 rect ( x / a ) exp(i 2πx )) + 1 2 rect ( x / a ) exp(−i 2πx ))
2
Λ Λ
Amplitude du champ diffracté dans le plan de Fourier (diffraction de Fraunhofer)

S (u ) = TF ( si t ( x))
S (u ) = 1 s (sin c(πau ) + sin c(πa (u − 1 / Λ)) + sin c(πa (u + 1 / Λ )) )
2 0

Un pic = une fréquence spatiale

40
Chaque pic est convolué par un sinus cardinal (effet de rect(x/a))
 Illustration sur un exemple simple
Objet diffractant de transmission sinusoidale et limité par une fente

t(x)

S(u)

Un pic = une fréquence spatiale41

Chaque pic est convolué par un sinus cardinal (dimension finie de l’objet)
 Généralisation: signification du signal dans
le plan de Fourier

− i 2π ( xu + yv )
S (u , v) ∝ ∫∫ s ( x, y )e dxdy
Dans cette intégrale de Fourier, u et v ont le caractère de fréquences spatiales

S(u,v) donne le contenu en fréquences spatiales du « signal » s(x,y).

Plus s(x,y) contient de détails à petite échelle, plus il est composé de signaux
hautes fréquences.

S(u,v) aura alors une valeur importantes pour u et v loin de zéro.

42
 Principe du filtrage spatial

S(u,v) donne le contenu en fréquences spatiales du « signal » s(x,y).

On va modifier le signal dans le plan de Fourier afin de modifier l’image

Il suffit pour cela de placer un diaphragme tF(u,v) dans le plan de Fourier

S’(u,v)= tF(u,v) S(u,v)

Comme dans un filtre en électronique, on peut ainsi selectionner/éliminer

certaines fréquences.

Si tF(u,v) est faible pour les valeurs de u et v grandes, on parle de filtre passe
bas…etc..;

On peut alors éliminer ou renforcer certains détails de l’image

43
Exemple 1D
Amplitude du champ au
niveau de l’objet

Diffraction de Fraunhofer
sans filtrage

Diffraction de Fraunhofer
après filtrage par une fente

Intensité dans le plan image

=> La modulation sinusoidale
a disparue !!! Il reste les
fréquences spatiales basses 44
 Exemple 2D
1) Objet initial 4) Image finale

2) Diffraction dans le plan de Fourier

3) Filtrage par une fente
Basses fréquences

Hautes fréquences

45
 Strioscopie
Certains objets sont transparents et se contentent de déphaser la lumière (cellules…)
t ( x , y ) = e iφ ( x , y )
L’intensité de l’image d’un tel objet est uniforme
on ne voit aucun détail !

Supposons que φ(x,y)<<1 alors: t ( x, y ) = e iφ ( x , y ) ≈ 1 + iφ ( x, y )

~
Et dans le plan de Fourier on a: TF (t ) ≈ δ (u , v) + iφ (u , v)

Si on place un petit diaphragme près du centre, on élimine la fonction δ !

Le champ dans le plan image est alors : S ( x , y ) ≈ iφ ( x , y )

Et l’intensité dans le plan image: I ( x , y ) ≈ φ 2 ( x, y )

Les variations de phase (épaisseur, indice…) deviennent visibles !!!

Mais on ne « voit » pas le signe de φ(x,y)…
46
 Strioscopie
f’ f’ f’ f’

Onde plane

Objet diffractant
transparent
Plan
d’observation

Plan de Fourier

Ondes de choc produites par une balle de fusil se

déplaçant dans l’air observées par strioscopie
47
 Contraste de phase
Analogue à la strioscopie mais cette fois, on déphase de π/2 la fonction δ
au lieu de la supprimer

t ( x , y ) = e i φ ( x , y ) ≈ 1 + iφ ( x , y )
~
Contraste
TF (t ) ≈ δ (u , v) + iφ (u , v)
~
c(u , v)TF (t ) ≈ iδ (u , v) + iφ (u , v)
de phase

Champ dans le
plan image
S ( x, y ) ≈ i (1 + φ ( x, y ))
Intensité dans le I ( x, y ) ≈ 1 + 2φ ( x, y ) + φ 2 ( x, y ) ≈ 1 + 2φ ( x, y )
plan image

On est cette fois sensible au signe de la phase

48
 Concrètement…
microscopie à contraste de phase

L’opération de filtrage
Microscopie classique Microscopie à contraste de
spatial se fait dans l’objectif
de microscope phase
49
Formation des images et
diffraction

50
 Notre objectif
Peut-on généraliser la notion de diffraction à des situation où la lumière passe par
un système optique ? (lentilles…)

Peut-on reformuler la notion d’image de l’optique géométrique dans le contexte

des ondes?

Intérêt pratique: les optiques doivent se comporter comme des diaphragmes à

cause de leur taille finie. Conséquences sur les performances d’un système
d’imagerie ?

Champ connu dans ce plan

lentille
Champ ici ?

« L’objet » 51
 Lentille « parfaite »
Lentille infiniment grande et fonctionnant dans les conditions de Gauss.
Elle transforme une onde incidente plane en une onde sphérique convergeant
vers le foyer F’

F’

f’
La lentille est équivalente à un « diaphragme » purement déphasant qui transforme
une onde plane en une onde sphérique convergent en f’. La transmission de la
lentille est formellement égale à:

− ik ( x 2 + y 2 ) / 2 f '
t L ( x, y ) = e 52
 Lentille de diamètre fini
En réalité, la lentille a un certain diamètre d’entrée (une certaine pupille)

F’

f’

On peut modéliser cette pupille par un diaphragme de transmission tO(x,y)

Exemple:
tO(x,y)=circ(r,a)

53
 Champ produit par un point source image

F’

p’
Point source
f’ Plan d’observation
=> Éclaire la lentille
avec une onde
sphérique
Pour calculer le champ dans le plan d’observation, on va utiliser le principe de Huygens
Fresnel:

1) on calcule le champ juste après la lentille

2) on intégre l’intégrale de Fresnel pour obtenir le champ dans le plan d’observation.
54
Ici, fraunhofer n’est a priori pas valable ! (diffraction à distance finie)
 Champ produit par un point source objet
ikR k ( x2 + y2 ) k ( Xx +Yy )
i e i −i
Intégrale de Fresnel S ( X ,Y ) =
λ z ∫∫ s( x, y)e 2z
e R
dxdy

ikR k ( x2 + y2 )
i e i
−i ( 2πux + 2πvy )
S ( X ,Y ) = ∫∫ s( x, y)e
2 p'
Ici z=p’ e dxdy
λ p'
On fait par ailleurs apparaître la partie « Fraunhofer »
s ( x , y ) = t ( x , y ) si ( x , y )
( x2 + y2 )
s0 ikr s0 ikp + ik
si ( x , y ) = e ≈ e 2p
Onde sphérique incidente
On calcule p p
s(x,y) t ( x , y ) = t O ( x , y )t L ( x , y )
( x2 + y2 )
− ik
t ( x, y ) = tO ( x, y )e 2f'

ik ( p + p ') k ( x2 + y2 )  1 1 1 
i  + − 
i s0 e  p p ' f '  − i ( 2 πux + 2 πvy )
S ( X ,Y ) = ∫∫ t
2
On obtient… ( x , y )e e dxdy
λ p' p O

55
Très compliqué…sauf si…
 Champ produit par un point source objet

1 1 1
Très compliqué…sauf si… + − =0
p p' f '

i s0 e ik ( p + p ') −i ( 2πux + 2πvy )

S ( X ,Y ) =
λ p' p ∫∫ O
t ( x , y ) e dxdy

S ( X , Y ) ∝ ∫∫ tO ( x, y )e −i ( 2πux + 2πvy ) dxdy

Dans ces conditions: la figure de diffraction produite dans le plan

d’observation est alors simplement la figure de diffraction de Fraunhofer de
la pupille de la lentille tO(x,y)

56
 Champ produit par un point source objet
A quoi correspond le plan d’observation tel que:
1 1 1
+ − =0
p p' f'
Ou encore avec des 1 1 1
notations algébriques: − =
p' p f'

C’est la relation de conjugaison de l’optique géométrique!!!

Le plan d’observation est dans ces conditions le plan image!

A cause de la diffraction , l’image d’un point est une tache lumineuse

Dans le plan image, la figure de diffraction d’un point objet est la figure de
diffraction de Fraunhofer de la pupille de la lentille tO(x,y).

Plus la pupille est petite, plus la tache sera grande!

57
 Conclusion: limite de diffraction

La résolution d’un système d’imagerie est liée à la taille de la tache de

diffraction produite par un point source.

Si deux points objets sont trop proches, leurs taches de diffraction

peuvent se chevaucher
on ne peut pas les discerner sur l’image.

La taille de la tache lumineuse produite par un point est :

λp'
δ '∝ Où a est le diamètre de l’optique
a
Pour pouvoir résoudre deux points objets, il faut que la distance entre eux
soit supérieure à

p λp λ
δ= δ '∝ ∝
p' a N .A
Ouverture numérique du
système d’imagerie
58
 Généralisation à toute une image:
Notion de réponse impulsionnelle
En anglais: PSF ou Point Spread Function

Point
source
objet
champ
diffracté
H(X,Y)

Plan objet Plan image

(x,y) (X,Y)
Système optique Σ 59
centré, aplanétique
 Rappel mathématique I: Distribution de Dirac
1
δ ( x) = lim rect ( x / ε )
Propriétés ε − >0 ε
+∞
1 = ∫ δ (t )dt
−∞

+∞
f ( x) = ∫ f (t )δ (t − x)dt
−∞

Transformée de Fourier TF (1) = δ (u )

TF (ei 2πax ) = δ (u − a )
TF (δ ( x)) = 1
TF (δ ( x − a )) = e −i 2πua 60
 Rappel mathématique 2:
Le produit de convolution
+∞

Définition h( x ) = ∫ f ( x' ) g ( x − x' )dx' = f * g = g * f

−∞

Propriété TF (h) = TF ( f )TF ( g )

Remarque utile:
+∞

−∞
∫ f ( x' )δ ( x − x' )dx' = f ( x) = f * δ

61
 Notion de réponse impulsionnelle
Quel est le champ dans le plan image généré par un point source objet ?

Système Σ
Point objet en (x=0,y=0) Tache de diffraction
dont le champ est H(X,Y)

Si le système est aplanétique:

Système Σ
Point objet en (x=x0,y=y0) Tache de diffraction
dont le champ est
H(X-γx0,Y- γ y0)
avec γ le grandissement
Modélisation du champ du point objet:

Champ de taille infiniment petite et donc d’intensité infiniment

grande (sa puissance reste finie)

fonction de Dirac : δ(x-x0,y-y0)
62
 Intérêt de la réponse impulsionnelle

Si on connait le champ dans le plan image (« la tache de diffraction »)

produit par un point source, on peut le connaitre pour un objet plus
complexe formé par une grande quantité de points sources.

Formellement on a:
Produit de convolution

G ( X , Y ) = ∫∫ s ( x, y ) H ( X − γx, Y − γy )dxdy

Réponse impulsionnelle du système optique

Champ dans le plan Champ dans le plan objet

image

63
 Quelques exemples (I)
Propagation libre sur un distance d
C’est Fresnel !!!
( X − x ) 2 + (Y − y ) 2
ikd + ik
2z
e
G ( X , Y ) ∝ ∫∫ s ( x, y ) dxdy
d

Champ « objet » Réponse impulsionnelle

Champ « image » c’est une onde sphérique !

Écrit ainsi, l’intégrale de Fresnel est un produit

de convolution !

Remarque:
( X − x ) 2 + (Y − y ) 2 X 2 +Y 2
ikd + ik ikd + ik
2z 2z
e e
∫∫ δ ( x, y) d
dxdy =
d
Le champ émis par nos points sources se comporte bien comme une onde sphérique
64
dans un autre plan!!!
 Quelques exemples (II)
Point source à l’infini

Champ d’un Diffraction:

f’
point objet à Diffraction de Fraunhofer
l’infini
de la pupille:
onde plane H(X,Y) proportionnel à
TF(t)

Pupille de la lentille :
modélisée par une
transmission t(x,y)

65
 Quelques exemples (III)
Point source à distance finie
Champ émis p p’ Diffraction:
par le point Diffraction de
objet
Fresnel du champ
Onde après
sphérique lentille+pupille

Plan objet Plan image

Pupille de la lentille : modélisée par
une transmission t(x,y)

Champ après lentille=t(x,y)*champ d’une onde convergente au

point image
Ce dernier compense le terme de Fresnel !

A nouveau H(X,Y) proportionnel à TF(t) (diffraction de 66

Fraunhofer)
 Fonction de transfert du système
(En anglais : Modulation Transfert Function)

On peut réécrire le formalisme précédent en utilisant la transformée de Fourier et

le fait que la TF d’un produit de convolution de deux fonctions est le produit des
leurs TF. Ainsi:

G( X , Y ) = s * H
TF (G ) = TF ( s )TF ( H )
« Réponse fréquentielle de
l’optique »
fonction de
Contenu en fréquences transfert du système
spatiales de l’image Contenu en fréquences
spatiales de l’objet

67
 Fonction de transfert du système
Plus la PSF d’un système optique est « fine »…

 plus le système est résolvant (on peut distinguer deux points objets très
proches)
Plus la MTF est large
Plus on peut distinguer des fréquences spatiales élevées dans l’image

Si on connait la PSF on connait la MTF et vice-versa…

68
 Fréquences spatiales ? Rappel…
Voici un objet dont on veut faire l’image…

Modulation de fréquence spatiale

élevée, dans la direction y

Modulation dont la fréquence spatiale est basse, dans la direction x

69
 Fonction de transfert: interprétation
Indique quelles modulations seront discernables ou pas (« floutées »)
Critère objectif: perte de contraste de la modulation

Plan objet Plan image

Fonction de transfert : équivalent 2D du diagramme de Bode des filtres

électroniques !

70
 Exemple numérique
1) Objet 3) Pupille de
Sa transmission l’optique
est modulée par formant
deux fréquences
l’image
spatiales (une
basse et une
haute)

2) TF (Objet)
-> pics éloignés du
centre
4) PSF de
=fréquence
spatiale haute l’optique
-> pic proche du =TF(pupille)
centre = =fonction
fréquences
spatiales d’Airy
basses

5) Produit de 6) Intensité
TF(Objet) et de dans le plan
la MTF de image La
l’optique modulation haute
Les fréquences fréquence n’est plus
hautes sont visible => résolution
atténuées dégradée 71par la
pupille
Formation des images en
lumière incohérente

72
 Notion de réponse impulsionnelle
en lumière incohérente

Point
source
champ
objet 1
diffracté 2
champ
Point diffracté 1
source H(X,Y)
objet 2

Plan objet Plan image

(x,y) (X,Y)
Système optique Σ 73
centré, aplanétique
 Notion de réponse impulsionnelle
en lumière incohérente
Les champs générés par 2 objets mutuellements incohérents se superposent EN
INTENSITE !!!

K ( X ,Y ) = H ( X ,Y )
2
Soit Réponse impulsionnelle incohérente

I ' ( X , Y ) = ∫∫ I ( x, y ) K ( X − γx, Y − γy )dxdy

Intensité dans le plan image

Intensité dans le plan objet

74