IUT Chby, Dpt SRC, Guyet T.

2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

En pratique …
On retrouve la transformée de Fourier dans
Transformée de Fourier et – (Re-)Echantillonnage
• Scanner un document!
Images • Ré-échantillonner une image (Et donc pour toute
transformation géométrique!)
– Le traitement d’image (Filtres)
• Débruiter une image = Passe-Bas
• Augmenter le contraste = Passe-Haut
– Format JPEG : A la base du fonctionnement

1 2

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

L’idée de Fourier

Tout signal périodique peut être décrit

Rappels de TF du signal comme une combinaison linéaire (CL) de
signaux élémentaires
Cosinus
Sinus
Caractérisées par une fréquence

3 4

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Un exemple … Un exemple …
cos(w1t) cos(w2t)
1 1
S(t) 0.5 0.5
4
0 0

2 -0.5 -0.5

-1 -1
0 5 10 15 20 25 30 35 0 5 10 15 20 25 30 35
0

-2 C’est une
combinaison linéaire
-4
0 5 10 15 20 25 30 35 S(t) = 2. cos(2πw2t) + cos(2πw1t) de deux signaux
4 élémentaires
2

-2

-4
0 5 10 15 20 25 30 35

5 6

1
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Un exemple … La transformé de Fourier

• Spectre en amplitude (ne pas oublier qu’il y a • Comment trouver la combinaison linéaire
aussi une phase!) ?
S(t) = 2. cos(2πw2t) + cos(2πw1t)
• Comment trouver quels cosinus et sinus
4=2² interviennent et quel est leur poids ?

1=1² Grâce à la transformé de Fourier !

w1 w2

7 8

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Pour les signaux finis Pour les signaux finis

• Tous les signaux ne sont pas • On fait comme si on recopiait le signal
nécessairement fini
– Infini Il devient :
– Périodique • Infini
• Périodique
• Et même aucun en pratique …
• On peut alors calculer sa TF !
Comment définir une TF ?

9 10

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Ondes élémentaire dans une

image
• Caractérisée par
Le cas des images en niveaux – Sa direction (Orthogonale à l’onde)
de gris – Sa longueur d’onde

11 12

2
 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Ondes élémentaire dans une Ondes élémentaire dans une

image image
• Caractérisée par • Caractérisée par
– Sa direction (Orthogonale à l’onde) – Sa direction (Orthogonale à l’onde)
– Sa longueur d’onde – Sa longueur d’onde

13 14

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Aide à la visualisation d’images

Images et fréquences
artificielles : Représentation 3D
• De même que pour les signaux on peut
10

20
retrouver des cosinus et des sinus dans
30 les images
40

50 10

20
60
30
70 40

80 50

60
90
70

100 80
10 20 30 40 50 60 70 80 90 100 90

100
10 20 30 40 50 60 70 80 90 100

15 16

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Image et fréquences Interprétation d’une FFT 2D

• Basse fréquence : Dégradé lent • Les hautes fréquences correspondent à des
changements d’intensité rapides
10

20

30

40

50

60

70

80

90

100
10 20 30 40 50 60 70 80 90 100

17 18

3
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Les extrémités en image Mise en garde

Front : Nécessite des hautes fréquences

Onde = Onde spatiale

Fréquence = Fréquence spatiale
Aplat : Uniquement une composante continue
Ceci n’a RIEN à voir avec les couleurs !!!!

19 20

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Exemple sur une image

Construction d’une
représentation des ondes
contenues dans une image
La Transformée de Fourier (TF)
d’Images

21 22
Image « Barbara »

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Transformé de Fourier d’une image Illustration de l’expression

• Toute image finie peut être décrite comme
une CL de sinus et cosinus = ( π (ω +ω )) + ( π (ω +ω ))
= =
• Expression complexe
πω πω Programme MatLab correspondant :
= = =
= = Icos=zeros(100);
for(i=1:100)
= (
π ω +ω ) for(j=1:100)
= =
Icos(i,j)=128+127*cos(2*pi*(10*i+20*j)/128);
end
= ( π (ω +ω )) + ( π (ω +ω )) end
= =
• Le rapport entre Wi et Wj donne l’orientation de l’onde
• La DFT permet de calculer les ci,j
23
• La valeur de + donne la pulsation d’onde associée 24

4
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Transformé de Fourier d’une image Représentation d’une DFT 2D

• Calcul des ci,j , pour tout i,j dans INn et INm • Pour un signal, le spectre est représenté par
une courbe
− π ω +ω • Le spectre d’une image peut être représenté par
=
+ + = =
une nouvelle ‘image’ (une matrice) (ci,j)
= = – Cette ‘image’ comprend les valeurs complexes des
+ coefficients de Fourier.
+
– La taille de la transformée est la même que celle de
− π + l’image
= + +
• De même, on parle de domaine fréquentiel,
+ + = = par opposition au domaine spatial (de l’image)

25 26

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Représentation d’une DFT 2D Importance des deux composantes

• Image complexe : On ne peut pas l’afficher
50

100

150

θ
=ρ
200

250

300

350

400

450

• Décomposition du spectre en images affichables 500

50 100 150 200 250 300 350 400 450 500

– Images d’amplitude 50

• Image constituée des (ρi,j) (i.e. On attribut cette valeur comme

100

150

intensité de gris)
200

250

– et de phase
300

350

• Image constituée des (θi,j)

400

450

500
50 100 150 200 250 300 350 400 450 500

27 28

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Image d’amplitude Vue en 3D

• Image de la même taille que
l’image du domaine spatial !
• S(i,j) = ρi-n/2,j-m/2 50 50

– S(i,j) est le coefficient

100 100
150

correspondant aux pulsations

150
200
200
suivant x et y: 250

250
300

− − 350 300
= = 400
+ + 450
350

• Plus un point est blanc, plus 500

50 100 150 200 250 300 350 400 450 500
400

sa valeur numérique est

450

500
importante 50 100 150 200 250 300 350 400 450 500

Plus il intervient dans l’image

29 30

5
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Transformée inverse Récapitulatif

• Transformée de Fourier inverse
– Opération inverse de la TF • Toute image numérique peut être décrite
Permet de reconstruire une image à partir d’un
spectre complexe comme un mélange d’ondes élémentaires
– Une image et son spectre sont équivalente
On peut passer de l’un à l’autre sans perdre d’informations
!
• Attention : il faut les deux composantes (amplitudes
et phases) pour reconstruire correctement le signal

31 32

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Récapitulatif Récapitulatif
Deux représentations pour un même objet/ • A quoi ça peut servir ?
Deux vues sur un même objet
– le spectre peut être
• Interprété : Le spectre donne des informations sur
l’image d’un caractère différent
• Modifié : En modifiant le spectre on va pouvoir
modifier les image
50 50

• Etudié : Les résultats théoriques sur les TFs seront

100 100

150 150

200 200

250 250

300 300

transposable aux images (e.g. critère de Nyquist)

350 350

400 400

450 450

500 500
50 100 150 200 250 300 350 400 450 500 50 100 150 200 250 300 350 400 450 500

Image Amplitude Image Phase

33 34

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Interprétation d’une FFT 2D

• Trois zones
importantes
Propriétés et interprétation des – Hautes fréquences : loin
spectres de Fourier d’une image du centre de l’image
– Basses fréquences :
proche du centre de
l’image
– Composante continue
(moyenne de l’image) sur
le pixel du centre

35 36

6
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Exemple de DTF Exemple de DFT

37 38

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Exemple de DFT Exemple de DFT

39 40

IUT Chby, Dpt SRC, Guyet T. 2006-

2006-2007 IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Rotation d’images Retour sur Jasmine

5
! "# ! $ 10

15

20

25

30

35

40

45

50

55

5 10 15 20 25 30 35 40 45 50 55

10

15

20

25

30

35

40

45

50

55

5 10 15 20 25 30 35 40 45 50 55

41 42

7
IUT Chby, Dpt SRC, Guyet T. 2006-
2006-2007

Interprétation

50

100

150

200

250

300

350

400

450

500
50 100 150 200 250 300 350 400 450 500

43