Frédéric Legrand Licence Creative Commons 1

Filtrage numérique par convolution

1. Introduction
Un signal numérique est une suite de nombres, obtenue par exemple par échantillon-
nage d’un signal analogique. Une fréquence d’échantillonnage est généralement associée
au signal numérique.
Le filtrage numérique est un traitement effectué sur cette suite de nombres, pour
transformer le signal. Ce calcul peut être fait sur un ordinateur à partir de signaux
stockés en mémoire. Dans les circuits d’électronique numérique, il est fait en temps réel
par des microprocesseurs spécialisés (DSP Digital Signal Processor), que l’on rencontre
par exemple dans les téléphones portables pour le traitement du son.
Nous allons aborder ici un type particulier de filtrage, le filtrage par convolution. Les
filtres obtenus sont à réponse impulsionnelle finie (filtre RIF).
Les images numérisées sont aussi des signaux numériques, mais à deux dimensions.
Une image numérisée est donc une matrice de nombres. Le filtrage par convolution se
pratique couramment sur les images numériques, par exemple pour rendre les images
légèrement floues afin de réduire le bruit, ou bien au contraire pour accentuer les détails.
Nous verrons comment se fait le filtrage par convolution sur une image.

2. Filtrage des signaux

2.a. Convolution
On considère un signal numérique donné par la suite xn , où n est un indice, compris
en 0 et N − 1.
Lorsque le signal est stocké en mémoire, le nombre d’échantillons N est fixé dès le
départ. Dans un système de traitement numérique temps réel, un nouvel échantillon doit
être traité à intervalle de temps régulier (la période d’échantillonnage). Dans ce cas, le
nombre d’échantillons N augmente au cours du temps.
Pour comprendre le filtrage par convolution, nous commençons par le cas du filtrage
temps réel. La figure suivante en montre le principe. Les indices sont notés entre crochets.
Frédéric Legrand Licence Creative Commons 2

Te t[N-2] t[N-1] t[N]

temps

x[N-5] x[N-4] x[N-3] x[N-2] x[N-1] x[N] Entrée

h[2] h[1] h[0] Noyau

y[N-6] y[N-5] y[N-4] y[N-3] y[N-2] y[N-1] y[N] Sortie

Seuls les 6 derniers échantillons sont représentés. La convolution se fait avec un noyau,
qui dans cet exemple comporte trois termes h[0], h[1], h[2]. Le signal filtré est noté yn .
Pour obtenir le dernier échantillon du signal filtré, on effectue une combinaison linéaire
des trois derniers échantillons du signal :

yN = h0 xN + h1 xN −1 + h2 xN −2 (1)
Dans un système temps réel, le dernier échantillon de la sortie est décalé d’une durée égale
à la période d’échantillonnage par rapport au dernier échantillon de l’entrée. En effet, le
calcul de cet échantillon prend un certain temps (qui doit être inférieur à la période). Sur
cet exemple, le nombre d’opérations à effectuer est très réduit (3 multiplications et la
somme de trois nombres), mais les cadences d’échantillonnage peuvent être très élevées,
de l’ordre de plusieurs gigahertz.
Plus généralement, si M est le nombre de termes du noyau, l’échantillon n de la sortie
se calcule en faisant une combinaison linéaire des M derniers échantillons de l’entrée :
M
X −1
yn = hi xn−i (2)
i=0

Cette relation définit la convolution discrète, et permet donc de réaliser un filtrage par
convolution. L’effet du filtre dépend évidemment des coefficients choisis pour le noyau.

2.b. Réponse impulsionnelle

Une impulsion est un signal numérique défini par :

x0 = 1 (3)
xn = 0 (n 6= 0) (4)

La sortie est alors :

Frédéric Legrand Licence Creative Commons 3

yn = 0 (n < 0) (5)
y0 = h0 x0 = h0 (6)
y1 = h1 x0 = h1 (7)
··· (8)
yM −1 = hM −1 x0 = hM −1 (9)
yn = 0 (n ≥ M ) (10)

Les coefficients h0 , h1 ...hN −1 (le noyau du filtre) constituent donc la réponse impulsion-
nelle du filtre.
Un filtre fonctionnant par convolution est dit à réponse impulsionnelle finie (RIF)
car le nombre de termes de la réponse impulsionnelle est fini.
Les filtres numériques RIF sont toujours stables : si l’entrée reste bornée, la sortie
l’est aussi.

2.c. Exemples
Filtre moyenneur
Le filtre moyenneur effectue la moyenne arithmétique de M échantillons consécutifs
du signal d’entrée. Par exemple, pour M = 3, le noyau est :
1 1 1
h=( , , ) (11)
3 3 3
Le filtre moyenneur est un filtre passe-bas, qui réduit les hautes fréquences.

Filtre passe-bas en sinus cardinal

Le filtre moyenneur a l’inconvénient d’être très peu sélectif : si l’on augmente la taille
du noyau, on réduit encore plus les hautes fréquences mais on altère aussi les basses
fréquences.
Le filtre passe-bas en sinus cardinal peut être très sélectif (bien plus que les filtres
électroniques analogiques). Pour le définir, il faut tout d’abord calculer le rapport de la
fréquence de coupure sur la fréquence d’échantillonnage (qui doit être inférieur à 1/2) :
fc
a= (12)
fe
On appelle filtre idéal un filtre qui laisse passer sans atténuation toutes les fréquences
inférieures à la fréquence de coupure, et qui annule les fréquences supérieures à la coupure.
Malheureusement, la réponse impulsionnelle d’un filtre passe-bas idéal est infinie. Elle
est donnée par :

sin(2π ka)
g(k) = 2a (13)
2π ka
Cette réponse impulsionnelle est rendue finie par troncature. Soit P l’indice de tronca-
ture. La réponse impulsionnelle finie comporte alors M = 2P + 1 termes. Elle est définie
par :

hk = g(k − P ) (14)
Frédéric Legrand Licence Creative Commons 4

l’indice k variant de 0 à 2P .
L’indice de troncature devra être d’autant plus grand que a est faible. Il faut donc
éviter une fréquence d’échantillonnage trop grande par rapport à la fréquence de coupure
(le sur-échantillonnage n’est pas toujours souhaitable en filtrage numérique). Un critère
simple est P a > 1. Pour un coefficient a donné, la sélectivité du filtre augmente avec P .
À la limite P → ∞, on obtient un filtre idéal.
L’inconvénient d’un indice P grand est la perte des M = 2P +1 échantillons au début
du signal.

Filtre dérivateur
La dérivation d’un signal numérique repose sur la formule de Taylor :

f (t + ) = f (t) + f 0 (t) + O(2 ) (15)

d’où l’on déduit :

f (t + ) − f (t)
f 0 (t) = + O() (16)

En terme de signal numérique, cela donne :
xn − xn−1
yn = (17)
Te
où Te est la période d’échantillonnage. Cette évaluation de la dérivée est appelée différence
finie. La réponse impulsionnelle est donc (on omet la période d’échantillonnage) :

h = (1, −1) (18)

Il y une autre manière d’évaluer la dérivée. En considérant le développement suivant :

f (t − ) = f (t) − f 0 (t) + O(2 ) (19)

La différence avec le premier donne :

f (t + ) − f (t − )
f 0 (t) = + O() (20)
2
ce qui donne la différence finie suivante :
xn − xn−2
yn = (21)
2Te
dont le noyau est :

h = (1, 0, −1) (22)

Dans ce cas, on considère généralement que l’échantillon de sortie yn est simultané avec
l’échantillon xn−1 . Cette différence finie est appelée différence finie centrée.
Frédéric Legrand Licence Creative Commons 5

3. Filtrage des images

3.a. Convolution
Les images numériques sont des signaux numériques à deux dimensions, qui s’appa-
rentent aux signaux de taille fixe discutés plus haut. Les échantillons sont appelés les
pixels.
Le filtrage par convolution se fait sur une image en niveaux de gris. Les images en
couleur sont composées de trois couches, par exemple rouge, vert et bleu. Chaque couche
peut être vue comme une image en niveau de gris. Il faudra donc opérer sur chaque
couche séparément pour filtrer une image en couleur (pas nécessairement avec le même
filtre).
On considère une image numérique en niveaux de gris, dont la valeur du pixel (i, j)
est notée Xi,j . On adopte la convention qui consiste à placer l’indice de colonne (indice
i) en premier, car cet indice correspond à l’axe x sur l’image. La réponse impulsionnelle
(ou noyau) dans le cas d’un signal numérique à deux dimensions est elle-même un signal
à deux dimensions, c’est-à-dire une matrice, dont les coefficients seront notés hm,n . Pour
simplifier, on se limite au cas des réponses impulsionnelles carrées, avec un nombre de
termes impair sur chaque dimension M = 2P + 1. Par exemple, pour P = 1 :
 
h00 h10 h20
H =  h01 h11 h21  (23)
h02 h12 h22
L’image filtrée Y s’obtient en effectuant le produit de convolution entre X et H, ce qui
donne pour un pixel (i, j) de l’image finale :

Yi,j = h00 Xi−1,j−1 +h10 Xi,j−1 +h20 Xi+1,j−1

+h01 Xi−1,j +h11 Xi,j +h21 Xi+1,j (24)
+h01 Xi−1,j+1 +h12 Xi,j+1 +h22 Xi+1,j+1
Le pixel Yi,j est donc obtenu en faisant une combinaison linéaire (ou moyenne pondérée)
du pixel Xi,j de l’image initiale et de ses 8 proches voisins. La relation ci-dessus se
généralise aux noyaux de convolution 5x5, 7x7, etc.
Frédéric Legrand Licence Creative Commons 6

Ligne 0

Ligne 1

Ligne 2

Ligne 3

i,j+1

i-1,j i,j i+1,j

i,j-1

h00 h10 h20

h01 h11 h12

h02 h12 h22

Remarque : il faudra faire attention au fait que les images sont stockées dans des tableaux,
dont le premier indice se réfère aux lignes horizontales. Le premier indice est donc l’indice
j. On remarque aussi que les lignes des images sont numérotées en partant du haut. Ainsi
la ligne 0 se trouve en haut de l’image. La correspondance entre les indices définis ci-
dessus et ceux du tableau est donc

Xi,j → image[Ny − 1 − j][i] (25)

où Ny est le nombre de lignes de l’image.
On voit que la convolution ne peut s’appliquer sur les bords de l’image, plus préci-
sément sur les P rangées horizontales et verticales des bords (coloriées en bleu sur la
figure). On évite généralement de réduire la taille des images lors du filtrage. Il faut donc
prévoir un traitement spécial pour ces rangées. La solution la plus simple consiste à ne
pas modifier ces pixels.
Bien que les valeurs des pixels soient souvent des entiers (par exemple codés sur 8
bits), les calculs de filtrage sont faits avec des nombres à virgule flottante. Cela évite des
erreurs d’arrondi grossières lorsqu’on applique plusieurs filtrages successifs, ou d’autres
types de traitement. Les valeurs de l’image finale sont converties en entier uniquement
au moment de l’affichage, ou du stockage dans un fichier.
Frédéric Legrand Licence Creative Commons 7

3.b. Exemples
Filtre moyenneur
Voici la réponse impulsionnelle d’un filtre moyenneur 3x3 :
 1 1 1 
9 9 9
1 1 1
H= 9 9 9
 (26)
1 1 1
9 9 9
Comme dans le cas des signaux à une dimension, le filtre moyenneur est un filtre passe-
bas, qui atténue les fins détails de l’image, ce qui peut être utile pour réduire le bruit.
Ce filtre ne modifie pas la plage des valeurs des niveaux de gris. Par exemple pour une
image codée en 8 bits (par couche), les niveaux sont compris entre 0 et 255, et restent
dans cette plage après application du filtre.

Filtres dérivateurs
La dérivation est une opération courante en traitement d’image, car elle permet de
faire ressortir les bords des formes. La dérivation par rapport à la direction horizontale
(X) est :

Hx = (−1, 0, 1) (27)
Il s’agit d’une différence finie centrée car la dérivée en un pixel est évaluée en utilisant
le pixel situé à sa droite et le pixel à gauche.
La dérivation par rapport à la direction verticale (Y) est :
 
1
Hy =  0  (28)
−1
La dérivation a l’inconvénient d’amplifier le bruit. Pour réduire cet effet, on associe un
filtrage passe-bas à la dérivation. C’est ce que fait le filtre de Sobel, par exemple pour la
direction X :
 
−1 0 1
Sx =  −2 0 2  (29)
−1 0 1
L’effet de filtrage passe-bas est obtenu ici par un moyennage pondéré dans la direction
Y.
Pour en savoir plus sur les filtres dérivateurs et l’application à la détection de bords,
voir Détection des bords.
Pour voir l’effet de plusieurs filtres sur des images, voir cette simulation.
Les filtres dérivateurs ne conservent pas l’intervalle [0, 255] des images 8 bits. En
effet, certaines dérivées peuvent être négatives. Il faut en tenir compte lorsqu’on veut
afficher l’image de sortie. La fonction [Link] convertit les valeurs
de manière à les faire entrer dans la plage [0, 255]. Une autre solution consiste à prendre
la valeur absolue pour afficher la dérivée. C’est ce qui est fait dans cette simulation.