Tweet
Tweet

Slide 1

Slide 1 text

Master Mathematical Analysis and Applications Course M1 - S1 Image processing - Restoration - Week 8-10 Mohammed Hachama [email protected] http://hachama.github.io/home Saad Dahlab University - Blida 1 -February 2021-

Slide 2

Slide 2 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (2/21) M. Hachama ([email protected])

Slide 3

Slide 3 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (3/21) M. Hachama ([email protected])

Slide 4

Slide 4 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Diff´ erentes perturbations d´ egradant les images • ´ etalonnage de l’appareil, • mauvaises conditions d’acquisition : nuages pour les photos a´ eriennes,manque de luminosit´ e etc ... • vieillissement des supports (photographies ou films anciens), • artefacts num´ eriques lors de reconstructions ` a partir des donn´ ees physiques (comme en tomographie par exemple). Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 5

Slide 5 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Estimation du bruit • Rapport signal sur bruit (SNR : Signal to Noise Ratio) d’un signal x entach´ e d’un bruit b sur un intervalle de temps fini I (en d´ ecibels dB) : SNR = 20 · log10 x b • Le SNR peut ˆ etre obtenu (exp´ erimentalement) grˆ ace ` a l’´ etalonnage des instruments de mesure. Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 6

Slide 6 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Bruit additif gaussien (centr´ e) d’´ ecart-type σ : ub = u + bσ • Bruit additif uniforme : ub = u + b • Bruit non lin´ eaire poivre et sel : ub = f (u) Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 7

Slide 7 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Bruit multiplicatif (donc proportionnel ` a l’intensit´ e lumineuse) : ub = u · (1 + bσ ) Exp. Le chatoiement (speckle) : petites taches rapidement fluctuantes qui donnent un aspect granuleux. Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 8

Slide 8 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 9

Slide 9 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Densit´ e gaussienne : P(z) = 1 σ √ 2π e−(z−µ)2 2σ2 • Bruit thermique, appareils, ... • Densit´ e de Rayleigh : P(x) = x σ2 exp −x2 2σ2 • Image de profondeur (depth), ... Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 10

Slide 10 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Densit´ e exponentielle : f (x) = λe−λx ; si x 0 0 ; si x < 0 • Laser, ... • Poivre et sel : P(w) = pw , P(b) = pb , P(x) = 0 sinon • Uniforme : P(x) = cst Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 11

Slide 11 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Bruit p´ eriodique • Des pics anormaux dans le domaine de Fourier Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 12

Slide 12 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Bruit p´ eriodique • Des pics anormaux dans le domaine de Fourier Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 13

Slide 13 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Description du bruit • Flou : convolution ` a effet fortement r´ egularisant ub = u ∗ gσ Image processing (week 8-10) -Restoration- (4/21) M. Hachama ([email protected])

Slide 14

Slide 14 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Restauration d’images • Am´ eliorer la qualit´ e visuelle des images • M´ ethodes objectives • Filtrage passe-bas/Filtrage spatial • Filtres non lin´ eaires : m´ edian, ... • Filtre de rehaussement • M´ ethodes subjectives • Mod´ eliser la d´ egradation de l’images • M´ ethodes en fonction de la d´ egradation Image processing (week 8-10) -Restoration- (5/21) M. Hachama ([email protected])

Slide 15

Slide 15 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Restauration d’images • Mod` ele g´ en´ erique : I0 = I ∗ h + n • Cas particulier (image bruit´ ee) : I0 = I + n Image processing (week 8-10) -Restoration- (5/21) M. Hachama ([email protected])

Slide 16

Slide 16 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (6/21) M. Hachama ([email protected])

Slide 17

Slide 17 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre spatial : utilisation de masque. • Filtre passe-bas : diminue le bruit mais att´ enue les d´ etails de l’image (flou plus prononc´ e). • Filtre passe-haut : accentue les contours et les d´ etails, mais amplifie le bruit. • Filtre passe-bande : ´ elimine certaines fr´ equences ind´ esirables. Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 18

Slide 18 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Propri´ et´ es des filtres • Filtre s´ eparable Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 19

Slide 19 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Strat´ egies pour le lissage d’images • Les images ne sont pas lisses lorsque les pixels voisins sont ”tr` es” diff´ erents. • Lissage = rendre les pixels adjacents similaires • Exemple : Filtre de moyenne • Chaque pixel rec ¸oit la moyenne de ses voisins Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 20

Slide 20 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre de moyenne Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 21

Slide 21 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre gaussien Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 22

Slide 22 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre gaussien Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 23

Slide 23 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre m´ edian • Combien de points bruit´ es peut-il filtrer ? Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 24

Slide 24 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre m´ edian Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 25

Slide 25 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre m´ edian Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 26

Slide 26 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre m´ edian Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 27

Slide 27 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre m´ edian Image processing (week 8-10) -Restoration- (7/21) M. Hachama ([email protected])

Slide 28

Slide 28 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • moyenne vs gaussien Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 29

Slide 29 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • moyenne vs gaussien Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 30

Slide 30 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • moyenne vs gaussien Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 31

Slide 31 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Filtre de moyenne Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 32

Slide 32 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Filtre gaussien Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 33

Slide 33 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Filtre de moyenne Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 34

Slide 34 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Filtre gaussien Image processing (week 8-10) -Restoration- (8/21) M. Hachama ([email protected])

Slide 35

Slide 35 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Filtre gaussien Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 36

Slide 36 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • ´ Eviter de moyenner sur les contours Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 37

Slide 37 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral • Poids adaptatifs Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 38

Slide 38 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 39

Slide 39 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 40

Slide 40 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 41

Slide 41 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 42

Slide 42 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 43

Slide 43 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 44

Slide 44 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 45

Slide 45 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 46

Slide 46 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 47

Slide 47 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 48

Slide 48 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 49

Slide 49 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 50

Slide 50 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 51

Slide 51 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral : Effet escalier Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 52

Slide 52 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral : Suppression de la texture Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 53

Slide 53 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral : Effet cartoon Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 54

Slide 54 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Filtre bilat´ eral : Effet cartoon Image processing (week 8-10) -Restoration- (9/21) M. Hachama ([email protected])

Slide 55

Slide 55 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local • Comme le filtre bilat´ eral : moyenner des pixels ”similaires” • Bilat´ eral : pixels voisins avec des intensit´ es proches • NL-Means : pixels ayant des voisinages similaires Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 56

Slide 56 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 57

Slide 57 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 58

Slide 58 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 59

Slide 59 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 60

Slide 60 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 61

Slide 61 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 62

Slide 62 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 63

Slide 63 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 64

Slide 64 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 65

Slide 65 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 66

Slide 66 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 67

Slide 67 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre non-local Image processing (week 8-10) -Restoration- (10/21) M. Hachama ([email protected])

Slide 68

Slide 68 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener • On connait I0 (image observ´ ee) et σ (variance du bruit). • Pour chaque pixel, on calcule la moyenne µL et la variance σL locales, calcul´ ees sur une fenˆ etre W • Reconstruction : J = µL + σ2 L σ2 L +σ2 (I0 − µL) • σ = 0 : J = I0 . Pas de bruit. • σL = 0 : J = µL . Pas de bruit. • σL >>= σ : J = I0 . Pas de r´ egularisation car pr´ esence d’un bord. • σL ≈ σ : J = 1 2 (I + µL ). Tirer la valeur du pixel vers la moyenne de ses voisins dans les ”r´ egions normales”. • Estimation de σ ` a partir de r´ egion homog` ene de l’image. Image processing (week 8-10) -Restoration- (11/21) M. Hachama ([email protected])

Slide 69

Slide 69 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener • Exemple Image processing (week 8-10) -Restoration- (11/21) M. Hachama ([email protected])

Slide 70

Slide 70 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (12/21) M. Hachama ([email protected])

Slide 71

Slide 71 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage • Types de filtre fr´ equentiel/spatial • Filtre fr´ equentiel F(g) = HF(f ) Image processing (week 8-10) -Restoration- (13/21) M. Hachama ([email protected])

Slide 72

Slide 72 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre passe-bas • Passe-bas id´ eal Image processing (week 8-10) -Restoration- (14/21) M. Hachama ([email protected])

Slide 73

Slide 73 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre passe-bas • Passe-bas id´ eal Image processing (week 8-10) -Restoration- (14/21) M. Hachama ([email protected])

Slide 74

Slide 74 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre passe-bas • Passe-bas id´ eal Image processing (week 8-10) -Restoration- (14/21) M. Hachama ([email protected])

Slide 75

Slide 75 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre passe-bas • Passe-bas de Butterworth Image processing (week 8-10) -Restoration- (14/21) M. Hachama ([email protected])

Slide 76

Slide 76 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique • Input Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 77

Slide 77 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique • Transform´ ee de Fourier Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 78

Slide 78 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique • Transform´ ee de Fourier Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 79

Slide 79 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique • Filtre Notch Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 80

Slide 80 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique • Filtrage fr´ equentiel Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 81

Slide 81 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtrage d’un bruit p´ eriodique Image processing (week 8-10) -Restoration- (15/21) M. Hachama ([email protected])

Slide 82

Slide 82 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (16/21) M. Hachama ([email protected])

Slide 83

Slide 83 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Application sur des images binaires et en niveaux de gris • Modification des morphologies des objets • Nettoyer le r´ esultat de la segmentation • Remplir les trous, ´ eliminer le bruit • Lisser le r´ esultat de la segmentation • Caract´ erisation • un ´ el´ ement structurant • des transformations : ´ erosion, dilatation, ouverture (´ erosion-dilatation), fermeture (dilatation-´ erosion) Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 84

Slide 84 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Soit une image binaire avec fond=0 et objet=1 • On glisse l’´ el´ ement structurant ”glisse” sur les pixels • Autres : ´ el´ ements structurants pas forc´ ement sym´ etriques • Application dune transformation • ´ Erosion : Si l’un des pixels du masque est fond (valeur 0), alors le pixel central devient fond (i.e., min) • Dilatation : Si l’un des pixels du masque fait partie de l’ objet (valeur> 0), alors le pixel central devient objet (i.e., max) • Autres : ´ Erosion puis dilatation (Ouverture) ; Dilatation puis ´ erosion (Fermeture) Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 85

Slide 85 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • ´ Erosion Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 86

Slide 86 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Dilatation Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 87

Slide 87 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Nombre d’objets ? Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 88

Slide 88 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Exemple de Dilatation Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 89

Slide 89 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • ´ Erosion-Dilatation Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 90

Slide 90 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Ouverture-Fermeture Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 91

Slide 91 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Exemple Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 92

Slide 92 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Influence de l’´ el´ ement structurant Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 93

Slide 93 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Exemple Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 94

Slide 94 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Gradient morphologique Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 95

Slide 95 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtres morphologiques • Gradient morphologique Image processing (week 8-10) -Restoration- (17/21) M. Hachama ([email protected])

Slide 96

Slide 96 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Plan 1. Bruit 2. Filtrage spatial 3. Filtrage fr´ equentiel 4. Filtrage morphologiques 5. D´ econvolution Image processing (week 8-10) -Restoration- (18/21) M. Hachama ([email protected])

Slide 97

Slide 97 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Mod´ elisation de la d´ egradation : I0 = I ∗ h + η • Filtre inverse • Formulation f = h ∗ u , f = hu • Calcul de 1/h. C’est le plus simple et peut donner de bons r´ esultats pour une image non bruit´ ee. Toutefois il n’est pas toujours possible d’inverser h car elle peut s’annuler. Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 98

Slide 98 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Filtre inverse : Image originale Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 99

Slide 99 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Filtre inverse : Image Floue Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 100

Slide 100 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Filtre inverse : D´ econvolution Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 101

Slide 101 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Filtre inverse : Image floue + Bruit Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 102

Slide 102 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre inverse • Filtre inverse : D´ econvolution Image processing (week 8-10) -Restoration- (19/21) M. Hachama ([email protected])

Slide 103

Slide 103 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Algorithme de Van Cittert • Formulation f = h ∗ u • R´ esolution : Posons ˆ g = 1 − ˆ h de sorte que formellement on obtient ˆ u = ˆ f 1 − ˆ g = +∞ k=0 ˆ gk ˆ f Si on pose uo = f et ˆ un = n k=0 ˆ gk ˆ f pour tout n = 1, on obtient ˆ un+1 = ˆ f + ˆ gˆ un = ˆ f + (1 − ˆ h)ˆ un , ou de mani` ere ´ equivalente un+1 = f + un − h ∗ un Image processing (week 8-10) -Restoration- (20/21) M. Hachama ([email protected])

Slide 104

Slide 104 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Algorithme de Van Cittert Image processing (week 8-10) -Restoration- (20/21) M. Hachama ([email protected])

Slide 105

Slide 105 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Algorithme de Van Cittert • Choix a priori du filtre h. • Ne prend pas en compte le bruit. Si f est une image flout´ ee et bruit´ ee :f = h ∗ u + b, o` u u est un bruit blanc gaussien. Le passage ` a un filtre inverse donne ˆ u = ˆ f ˆ h − ˆ b ˆ h . Image processing (week 8-10) -Restoration- (20/21) M. Hachama ([email protected])

Slide 106

Slide 106 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener • Formulation variationnelle de type Tychonov min v∈L2(R) f − h ∗ u 2 L2 + q ∗ u 2 L2 min v∈L2(R) ˆ f − ˆ h v 2 2 + ˆ qv 2 2 • D´ erivation : ˆ h∗(ˆ h ˆ u − ˆ f ) + |ˆ q|2ˆ u = 0 Et donc : ˆ u = ˆ h∗ˆ f |ˆ h|2 + |ˆ q|2 Image processing (week 8-10) -Restoration- (21/21) M. Hachama ([email protected])

Slide 107

Slide 107 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener • Choix de q : Id´ ealement, il faudrait choisir ˆ q(w) = |ˆ b(w)| |ˆ u(w)| car, dans ce cas, le deuxi` eme terme p´ enalise le bruit. En effet, : q ∗ u 2 = ˆ qˆ u 2 = ˆ b 2 = b 2 . • Apr` es calcul : W = ¯ ˆ h |ˆ h|2 + |ˆ q|2 , Dans le cas o` u b est un bruit blanc gaussien d’´ ecart-type σ : |ˆ q|2 = σ2 < |ˆ h| >2 < |ˆ f | >2 −σ2 . Image processing (week 8-10) -Restoration- (21/21) M. Hachama ([email protected])

Slide 108

Slide 108 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener Image processing (week 8-10) -Restoration- (21/21) M. Hachama ([email protected])

Slide 109

Slide 109 text

Bruit Filtrage spatial Filtrage fr´ equentiel Filtrage morphologiques D´ econvolution Filtre de Wiener Image processing (week 8-10) -Restoration- (21/21) M. Hachama ([email protected])