FR 2736741 A1 19970117 FR 9508632 A 19950711 Le domaine de l'invention est celui des images numériques, que celles-ci soient générées directement par un système d'acquisition numérique ou bien numérisées après leur acquisition. Plus précisément, l'invention concerne un procédé de traitement d'une séquence source d'images numériques ayant été endommagées (c'est-à-dire présentant un bruit élevé), de façon à obtenir une séquence de sortie d'images numériques corrigées. D'une façon générale, le bruit qui affecte les images numériques se traduit par une dégradadation du contraste de ces images en fonction de l'intensité du bruit et par une perte d'informations dans les zones endommagées de ces images. Il existe plusieurs types de bruit, tels que notamment les bruits dus aux mouvement ou défocalisation du capteur rendant les images floues (motion/focus blur en anglo-saxon), les bruits sous forme de tâches (sperme noise en anglo-saxon) ou encore les bruits dus aux mesures (bruits additifs ou multiplicatifs). Il convient donc de compenser les effets du bruit, de façon à permettre la restitution d'une séquence d'images présentant une qualité suffisante. I1 existe plusieurs techniques connues de correction de séquences d'images endommagées, telles que notamment le filtrage inverse ou le calcul de la moyenne de la séquence d'images. La technique connue du filtrage inverse consiste à appliquer l'inverse de la fonction de transfert à l'image endommagée. Ceci suppose que le type de bruit (ou une approximation de celui-ci soit connu. Or cela n'est pas toujours le cas. Par ailleurs, cette technique connue du filtrage inverse est peu efficace lorsque le rapport bruit / signal est élevé. D'une façon générale, la seconde technique connue, qui consiste à calculer la moyenne de la séquence d'images, ne permet pas d'apporter une correction suffisante aux images endommagées. La technique connue du filtrage inverse présente quant à elle l'inconvénient majeur d'être spécifique à un type de bruit particulier et d'introduire des erreurs lorsque la séquence d'images contient des objets mobiles. Par conséquent, lorsqu'elle est mise en oeuvre pour corriger un type de bruit différent de celui pour lequel elle est prévue, elle engendre des pertes d'informations supplémentaires. Par ailleurs, dans la plupart de ces techniques connues, la correction de chaque image de la séquence est effectuée en plusieurs passes successives, chaque passe correspondant à un filtrage de correction dont les paramètres peuvent être modifiés par un opérateur (on parle dans ce cas de filtrage interactif). Le temps nécessaire à l'exécution de ces passes successives rend impossible une correction en temps réel de la séquence d'images numériques endommagées. L'invention a notamment pour objectif de pallier ces différents inconvénients de l'état de la technique. Plus précisément, l'un des objectifs de la présente invention est de fournir un procédé de traitement d'une séquence source d'images numériques endommagées, du type permettant d'obtenir une séquence de sortie d'images numériques corrigées, qui soit indépendant du type de bruit affectant les images de la séquence source. L'invention a également pour objectif de fournir un tel procédé qui permette un traitement en temps réel de la séquence source. Un autre objectif de l'invention est de fournir un tel procédé qui permette de restituer une séquence source d'images représentant un panorama (avec ou sans objets mobiles), ainsi que le contraste de ces images. Ces différents objectifs, ainsi que d'autres qui apparaîtront par la suite, sont atteints selon l'invention à l'aide d'un procédé de traitement d'une séquence source d'images numériques endommagées, du type permettant d'obtenir une séquence de sortie d'images numériques corrigées, chacune des images numériques brouilées ou corrigées étant décrite pixel par pixel, chacun desdits pixels étant caractérisé par un niveau d'amplitude parmi une pluralité de niveaux d'amplitude possibles, caractérisé en ce que ledit procédé comprend notamment, pour chaque image numérique endommagée de ladite séquence source, les étapes suivantes:: - calcul, pour l'ensemble des pixels de l'image endommagée, d'un premier paramètre, dit paramètre global, d'estimation de la correction à apporter à chacun des pixels de l'image endommagée, en fonction dudit ensemble des pixels de l'image endommagée et de l'ensemble des pixels d'une image de référence, ladite image de référence étant une image corrigée qui résulte de la correction de l'image précédant ladite image endommagée dans ladite séquence source; - calcul, pour chaque pixel donné de l'image endommagée, d'un second paramètre, dit paramètre local, d'estimation de la correction à apporter audit pixel donné, en fonction dudit pixel donné et d'autres pixels, dits pixels voisins, localisés dans un voisinage prédéterminé dudit pixel donné;; - calcul, pour chaque pixel donné de l'image endommagée, d'un troisième paramètre, dit paramètre temporel, d'estimation de la correction à apporter audit pixel donné, en fonction dudit pixel donné et du pixel de l'image de référence possédant les mêmes coordonnées spatiales que ledit pixel donné; - calcul, pour chaque pixel donné de l'image endommagée, d'un coefficient de correction, par combinaison desdits premier, second et troisième paramètres d'estimation associés audit pixel donné; - correction, selon une stratégie de correction prédéterminée, de chaque pixel donné de l'image endommagée de façon à obtenir le pixel correspondant de l'image corrigée de la séquence de sortie, en fonction du coefficient de correction associé audit pixel donné, du pixel donné de l'image endommagée et du pixel de l'image de référence possédant les mêmes coordonnées spatiales que ledit pixel donné. Ainsi, le principe général de l'invention consiste à traiter successivement chacune des images endommagées formant la séquence source. Pour chaque image endommagée, on corrige chacun des pixels en fonction d'un coefficient de correction qui lui est associé et qui résulte de la combinaison de trois paramètres distincts (respectivement global, local et temporel) d'estimation de la correction à apporter au pixel. Grâce à cette combinaison de trois paramètres fournissant chacun des indications distinctes sur la correction à apporter, le procédé de l'invention n'est pas interactif et ne nécessite qu'une seule passe de traitement pour chaque image endommagée. De plus, la plupart des calculs peuvent se présenter sous forme matricielle. Par conséquent, le procédé de l'invention est parfaitement adapté à une mise en oeuvre en temps réel. Par ailleurs, cette combinaison de trois paramètres rend le procédé de l'invention pratiquement indépendant du type de bruit affectant les images de la séquence source. Ceci permet de ne faire qu'un nombre limité d'hypothèses quant au type de bruit ayant affecté les images de la séquence source. On peut par exemple se limiter aux hypothèses selon lesquelles le bruit est distribué et intermittent. En d'autres termes, le procédé de l'invention est du type adaptatif et convient quel que soit le type de bruit à compenser. Le premier paramètre (ou paramètre global) est commun à tous les pixels d'une même image endommagée, et fournit une indication sur la qualité générale de cette image endommagée, par comparaison avec l'image corrigée précédente (ou image de référence). Avantageusement, ledit calcul, pour l'ensemble des pixels de l'image endommagée, du premier paramètre Pl s'écrit selon la formule suivante : P1 = K + fE (H1, H2) où: - K est une valeur de décalage prédéterminée; - H1 est un premier histogramme des niveaux d'amplitude des pixels de l'image endommagée; - H2 est un second histogramme des niveaux d'amplitude des pixels de l'image de référence; - fE est une fonction d'erreur prédéterminée permettant de calculer un écart entre deux fonctions. La fonction d'erreur prédéterminée est par exemple basée sur la méthode des moindres carrés. Elle peut également être basée sur la différence des variances des histogrammes de l'image endommagée et de l'image de référence. La valeur de décalage prédéterminée (ou offset) permet, si elle est choisie différente de zéro, d'éviter que le premier paramètre soit nul lorsque les deux histogrammes sont identiques (cas où fE (H1, H2) = 0). Le second paramètre (ou paramètre local) est spécifique à chacun des pixels d'une même image endommagée. I1 fournit une indication sur la présence d'une discontinuité spatiale, par comparaison avec les pixels voisins de la même image endommagée. En effet, le bruit se traduit par une discontinuité spatiale assez prononcée dans une même image. Par conséquent, la détection d'une telle discontinuité spatiale au niveau d'un pixel indique une probabilité assez forte pour que ce pixel soit bruité. De façon avantageuse, ledit calcul du second paramètre est effectué simultanément pour tous les pixels de l'image endommagée, et s'écrit selon la formule suivante: [Mp2] = (l/ pixel distinct de l'image endommagée; - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée; - F2 est un filtrage moyen ou médian ou passe-bas ou tout autre filtrage adapté au bruit traité; - a2 est un premier facteur de normalisation. On peut éventuellement utiliser une combinaison de filtreages moyen, médian et passe-bas, au lieu d'utiliser un seul filtrage, de façon à combiner les avantages des différents types de filtrage. Le troisième paramètre (ou paramètre temporel) est spécifique à chacun des pixels d'une même image endommagée. Il fournit une indication sur un changement de valeur d'amplitude, par comparaison avec le pixel de l'image corrigée précédente (image de référence) possédant les mêmes coordonnées spatiales. En effet, le bruit se traduit également par une discontinuité temporelle assez prononcée entre deux images succesives. Par conséquent, la détection d'une telle discontinuité temporelle au niveau d'un pixel indique une probabilité assez forte pour que ce pixel soit bruité. Il convient toutefois de s'assurer que cette discontinuité temporelle n'est pas due à un objet en mouvement entre les images successives. Préférentiellement, ledit calcul du troisième paramètre est effectué simultanément pour tous les pixels de l'image endommagée, et s'écrit selon la formule suivante: [Mp3] = (1/a3) . F3 { [1t] - [1t-i] } où: - [Mp3] est une matrice des troisièmes paramètres P3 associés chacun à un pixel distinct de l'image endommagée; - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée; - [1t-i] est une matrice des valeurs d'amplitude des pixels de l'image de référence; - F3 est un filtrage moyen ou médian ou passe-bas; - Ot3 est un second facteur de normalisation. Le filtrage permet de ne pas prendre en compte des changements de valeurs d'amplitude de pixel dus à des mouvements d'objet entre les images successives. De façon préférentielle, ledit calcul du coefficient de correction s'écrit selon la formule suivante: C = max [1, P1 . fc (P2, P3)] où: - P1, P2 et P3 sont lesdits premier, second et troisième paramètres normalisés à 1; - fc est une fonction de combinaison prédéterminée. Avantageusement, ladite stratégie de correction prédéterminée consiste, pour chaque pixel donné de l'image endommagée, à calculer une valeur d'amplitude de pixel corrigé It' (x, y ) selon la formule suivante: It'(x,y ) R1 (C). .It(x,y)+R2(C) 1t-I (x, y) où: - It (x, y ) est la valeur d'amplitude du pixel donné de l'image endommagée; - 1t-1 (x, y ) est la valeur d'amplitude du pixel de l'image de référence possédant les mêmes coordonnées spatiales (x, y) que le pixel donné; - C est le coefficient de correction normalisé à 1; - R1 et R2 sont des fonctions de répartition respectant les contraintes suivantes:R1 (C) + R2 (C) R1 (O) = 1, R2 est une fonction croissante telle que R2 (i) = 1. Ainsi, on peut distinguer les trois cas suivants: - si C = O, la valeur d'amplitude du pixel corrigé est égale à la valeur d'amplitude du pixel de l'image endommagée (pas de remplacement); - si C = 1, la valeur d'amplitude du pixel corrigé est égale à la valeur d'amplitude du pixel de l'image de référence (remplacement complet); - si 0 pondérée des valeurs d'amplitude du pixel de l'image endommagée et du pixel de l'image de référence (remplacement pondéré). Selon une variante avantageuse, ladite stratégie de correction prédéterminée consiste, pour chaque pixel donné de l'image endommagée, à calculer une valeur d'amplitude de pixel corrigé It" (x, y ) selon la formule suivante: I"(x,y) R1 (C). It(x,y)+R(C) eIt, (x,y)+R3(C) It,F(x,y) où:: - it (x, y ) est la valeur d'amplitude du pixel donné de l'image endommagée; - It l (x, y ) est la valeur d'amplitude du pixel de l'image de référence possédant les mêmes coordonnées spatiales (x, y) que le pixel donné; - 1t,F (x, y ) est la valeur d'amplitude du pixel d'une image filtrée possédant les mêmes coordonnées spatiales (x, y) que le pixel donné, ladite image filtrée étant obtenue par un filtrage moyen ou médian ou passe-bas ou tout autre filtrage adapté au bruit traité de l'image endommagée; - C est le coefficient de correction normalisé à i - R1, R2 et R3 sont des fonctions de répartition respectant les contraintes suivantes :R1 (C) + R2 (C) + R3 (C) S 1, R1 est une fonction décroissante telle que R1 (0) = 1, R2 est une fonction croissante telle que R2 (1) = 1. Préférentiellement, ledit procédé comprend, pour chaque pixel donné de l'image endommagée, les étapes supplémentaires suivantes: - calcul d'une probabilité d'erreur sur la valeur d'amplitude du pixel donné, en fonction de l'écart entre le nombre de pixel de l'image endommagée et le nombre de pixel de l'image de référence possédant la même valeur d'amplitude que ledit pixel donné, ladite probabilité d'erreur constituant un quatrième paramètre d'estimation de la correction à apporter audit pixel donné - pondération par ledit quatrième paramètre des premier et second paramètres ou du coefficient de correction associé audit pixel donné. Ces deux étapes supplementaires permettent d'affiner le calcul du coefficient de correction. Dans un mode de réalisation préférentiel de l'invention, ledit procédé comprend les étapes préalables suivantes: - calcul d'un premier histogramme des niveaux d'amplitude des pixels de l'image endommagée; - calcul d'un second histogramme des niveaux d'amplitude des pixels de l'image de référence; - calcul de l'écart entre les valeurs moyennes desdits premier et second histogrammes; - correction des valeurs d'amplitude des pixels de l'image endommagée en fonction dudit écart de valeur moyenne, de façon à équilibrer les valeurs d'amplitude des pixels de l'image endommagée avec celles des pixels de l'image de référence, et obtenir une image endommagée précorrigée qui est utilisée à la place de l'image endommagée dans toutes les autres étapes dudit procédé. Ces étapes préalables correspondent à la mise en oeuvre d'une correction de gain sur les valeurs d'amplitude de l'image endommagée en fonction des valeurs d'amplitude de l'image de référence. Avantageusement, le procédé de l'invention s'applique au traitement en temps réel d'une séquence source d'images numériques endommagées. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante d'un mode de réalisation préférentiel de l'invention, donné à titre d'exemple indicatif et non limitatif, et des dessins annexés, dans lesquels: - la figure 1 illustre l'application du procédé selon l'invention; - la figure 2 présente un mode de réalisation particulier du procédé selon l'invention; - la figure 3 illustre le fonctionnement du procédé de la figure 2; - les figures 4A et 4B présentent respectivement, de façon schématique, une image endommagée et une image de référence correspondante ; et - les figures 5, 6 et 7 illustrent chacune un mode de réalisation distinct du calcul du coefficient de correction d'un pixel en fonction des trois ou quatre paramètres associés à ce pixel. L'invention concerne donc un procédé de traitement d'images numériques. Plus précisément, comme présenté sur la figure 1, l'invention s'applique au traitement 1 d'une séquence source 2 d'images numériques endommagées, de façon à d'obtenir une séquence de sortie 3 d'images numériques corrigées. De façon connue, les images numériques (endommagées ou corrigées) sont composées d'une pluralité de pixels caractérisé chacun par un niveau d'amplitude parmi plusieurs possibles. Par niveau d'amplitude, on entend par exemple, pour une image en noir et blanc, le niveau de gris (généralement compris entre 0 et 255). Dans le cas d'une image couleur, chacune des trois couleurs de base peut être traitée séparément. Ceci revient à effectuer, sur chacune des trois images associée à une couleur distincte, un traitement du type noir et blanc. Selon une variante, le niveau d'amplitude est égal au niveau d'intensité de l'image couleur. La figure 2 présente un mode de réalisation particulier du procédé selon l'invention, et la figure 3 en illustre le fonctionnement. Dans un mode de réalisation particulier, le procédé selon l'invention comprend, pour chaque image numérique endommagée 31 de la séquence source 2, les étapes suivantes: - calcul (21) d'un premier paramètre P1, dit paramètre global, d'estimation de la correction à apporter à chacun des pixels de l'image endommagée 31; - calcul (22) d'un second paramètre P2, dit paramètre local, d'estimation de la correction à apporter à un pixel donné de l'image endommagée 31; - calcul (23) d'un troisième paramètre P3, dit paramètre temporel, d'estimation de la correction à apporter à un pixel donné de l'image endommagée 31; - calcul (24) d'un coefficient de correction C d'un pixel donné de l'image endommagée 31, par combinaison des premier, second et troisième paramètres P1, P2, P3 associés à ce pixel donné;; - correction (25) d'un pixel donné de l'image endommagée 31, en fonction du coefficient de correction C associé à ce pixel donné. L'étape 21 de calcul du premier paramètre P1 est commune à tous les pixels d'une même image endommagée 31. Par contre, les étapes suivantes 22, 23, 24 (calcul des second et troisième paramètres P2, P3, calcul du coefficient de correction C et correction de pixel) sont répétées pour chacun des pixels de l'image endommagée 31,jusqu'à ce que tous les pixels aient été traités. On détaille maintenant successivement chacune de ces étapes 21 à 25. L'étape 21 de calcul du premier paramètre P1, qui est commune à l'ensemble des pixels de l'image endommagée 31, est fonction de l'ensemble des pixels d'une image de référence 32. Cette image de référence 32 est l'image corrigée appartenant à la séquence de sortie 3, qui résulte de la correction de l'image précédant l'image endommagée 31 dans la séquence source 2. En d'autres termes, l'image de référence 32 d'une image endommagée 31 de rang n dans la séquence source 2 est le résultat du traitement de l'image endommagée de rang n - 1 dans cette même séquence source 2. Dans le présent mode de réalisation, le premier paramètre P1 est un scalaire tel que: P1 =K+fE(Hl,H2) où: - K est une valeur de décalage prédéterminée; - H1 est un premier histogramme des niveaux d'amplitude des pixels de l'image endommagée; - H2 est un second histogramme des niveaux d'amplitude des pixels de l'image de référence; - fE est une fonction d'erreur prédéterminée permettant de calculer un écart entre deux fonctions. La valeur de décalage K est par exemple égale à 1. La fonction d'erreur fE est par exemple la fonction des moindres carrés. D'une façon générale, fE (H1, H2) indique la distance entre les premier et second histogrammes H 1, H2. L'étape 22 de calcul du second paramètre P2, qui est réalisée pour chaque pixel donné de l'image endommagée 31, est fonction de ce pixel donné et de pixels voisins de ce pixel donné. La figure 4A présente de façon schématique l'image endommagée 31 et permet d'expliciter la notion de pixels voisins. Le pixel donné à traiter 41 possède ici les coordonnées spatiales (x, y). Les pixels voisins 44 de ce pixel donné 41 sont ceux compris dans une zone de voisinage 43 prédéterminée (par exemple une fenêtre de N pixels sur N pixels, centrée sur le pixel donné 41). Dans le présent mode de réalisation, le calcul 22 du second paramètre P2 est effectué simultanément pour tous les pixels de l'image endommagée (31), et s'écrit selon la formule suivante: [Mp2] = (l/a2) . ([1t] - F2 [1t] ) où: - [Mp2] est une matrice des seconds paramètres P2 (chacun étaassociés chacun à un pixel distinct de l'image endommagée 31; - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée 31; - F2 est un filtrage moyen ou médian ou passe-bas ou tout autre filtrage adapté au bruit traité; - a2 est un premier facteur de normalisation (par exemple égal à la moyenne ou à la médiane pour un filtrage F2 moyen ou médian respectivement). L'étape 23 de calcul du troisième paramètre P3, qui est réalisée pour chaque pixel donné 41 de l'image endommagée 31, est fonction de ce pixel donné 41 et du pixel de l'image de référence possédant les mêmes coordonnées spatiales que ce pixel donné 41. La figure 4B présente, de façon schématique, l'image de référence 32 correspondant à l'image endommagée 31 de la figure 4A. On y a représenté le pixel 42 possédant les mêmes coordonnées spatiales que le pixel donné 41 de l'image endommagée 31. Dans le présent mode de réalisation, le calcul 23 du troisième paramètre P3 est effectué simultanément pour tous les pixels de l'image endommagée 31, et s'écrit selon la formule suivante: [Mp3] = (l/a3) . F3 { [It] - [It l] } où: - [Mp3] est une matrice des troisièmes paramètres P3 associés chacun à un pixel distinct de l'image endommagée 31; - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée 31; - [1t-i] est une matrice des valeurs d'amplitude des pixels de l'image de référence 32; - F3 est un filtrage moyen ou médian ou passe-bas; - a3 est un second facteur de normalisation (par exemple égal à la moyenne ou à la médiane pour un filtrage F3 moyen ou médian respectivement). Dans un premier mode de réalisation, illustré sur la figure 5, le coefficient de correction C est tel que : C = max [1, P1 . fc (P2, P3)] où: - P1, P2 et P3 sont les premier, second et troisième paramètres normalisés à 1 ; - fc est une fonction de combinaison prédéterminée (par exemple, une moyenne pondérée). Les paramètres P2 et P3 sont combinés (51), par exemple par une sommation pondérée, et le résultat de cette combinaison 53 est multiplié (52) par le premier paramètre P1, de façon à obtenir le coefficient de correction C. Un second et un troisième modes de réalisation du calcul du coefficient de correction C sont maintenant explicités, en relation avec les figures 6 et 7. Ces deux variantes nécessitent, pour chaque pixel donné de l'image endommagée 31, le calcul d'un quatrième paramètre P4 d'estimation de la correction à apporter à ce pixel donné 41. Ce quatrième paramètre P4 est une probalité d'erreur sur la valeur d'amplitude du pixel donné 41, en fonction de l'écart E entre le nombre de pixel de l'image endommagée 31 et le nombre de pixel de l'image de référence possédant la même valeur d'amplitude que le pixel donné 41. L'écart E s'écrit par exemple: où: - H1 et H2 sont les histogrammes des niveaux d'amplitude des pixels respectivement de l'image endommagée 31 et de l'image de référence 32; - n est la valeur d'amplitude du pixel donné 41. Dans le second mode de réalisation du calcul du coefficient de correction C (cf fig.6), les second, troisième et quatrième paramètres P2, P3 et P4 sont combinés (61), par exemple par une sommation pondérée, et le résultat de cette combinaison 63 est multiplié (62) par le premier paramètre P1, de façon à obtenir le coefficient de correction C. Dans le troisième mode de réalisation du calcul du coefficient de correction C (cf fig.7), les paramètres P2 et P3 sont combinés (71), par exemple par une sommation pondérée, et le résultat de cette combinaison 73 est multiplié (72) par le premier paramètre P1. Le résultat de cette multiplication 74 est combiné à ou multiplié avec (75) le quatrième paramètre P4, de façon à obtenir le coefficient de correction C. L'étape 25 de correction d'un pixel donné 41 de l'image endommagée 31, en fonction du coefficient de correction C associé à ce pixel donné 41, est effectuée selon une stratégie de correction prédéterminée. L'image corrigée correspondant à l'image endommagée 31 est obtenue après que tous les pixels de cette image endommagée ont été corrigés. On rappelle ici que cette image corrigée servira d'image de référence pour le traitement de l'image endommagée suivante de la séquence source 2. Un premier mode de réalisation de la stratégie de correction prédéterminée consiste, pour chaque pixel donné 41 de l'image endommagée 31, à calculer une valeur d'amplitude de pixel corrigé It' (x, y ) selon la formule suivante: It (xty)=Rl(c) It(xsy)+R2(c) It-l (xsy) où: - It (x, y ) est la valeur d'amplitude du pixel donné 41 de l'image endommagée 31; - It l (x, y ) est la valeur d'amplitude du pixel 42 de l'image de référence 32 possédant les mêmes coordonnées spatiales (x, y) que le pixel donné 41; - C est le coefficient de correction normalisé à 1; - R1 et R2 sont des fonctions de répartition. Les fonctions de répartition R1 et R2 doivent respecter un certain nombre de contraintes, à savoirs: - Rl (C) + R2 (C) - R1 est une fonction décroissante telle que R1 (O) = 1, - R2 est une fonction croissante telle que R2 (1) = 1. Les deux paramètres R1 (C) et R2 (C) sont par exemple calculés lors d'une étape particulière 28 effectuée après le calcul du coefficient de correction C. La contrainte R1 (O) = 1 garantit que lorsque C =0, la valeur d'amplitude du pixel corrigé est égale à la valeur d'amplitude du pixel de l'image endommagée (pas de remplacement). De même, la contrainte R2 (1) = 1 garantit que lorsque C = 1, la valeur d'amplitude du pixel corrigé est égale à la valeur d'amplitude du pixel de l'image de référence (remplacement complet). Entre ces deux cas extrêmes, c'est-à-dire si 0 En d'autres termes, le pixel donné 41 de l'image endommagée 31 est remplacé par le pixel correspondant 42 de l'image de référence 32 proportionnellement au bruit qu'il est supposé présenter, ce dernier étant indiqué par le coefficient de correction C associé à ce pixel donné 41. Un second mode de réalisation de la stratégie de correction prédéterminée consiste, pour chaque pixel donné 41 de l'image endommagée 31, à calculer une valeur d'amplitude de pixel corrigé It" (x, y ) selon la formule suivante: It" (x, y) = R1 (C). I(x, y) + R2 (C). 1t-I (x, y) + R3 (C).It,(x, y) où: - It(x, y ) est la valeur d'amplitude du pixel donné 41 de l'image endommagée 31; - It l (x, y ) est la valeur d'amplitude du pixel 42 de l'image de référence 32 possédant les mêmes coordonnées spatiales (x, y) que le pixel donné 41; - It,F (x, y ) est la valeur d'amplitude du pixel d'une image filtrée possédant les mêmes coordonnées spatiales (x, y) que le pixel donné 41, l'image filtrée étant obtenue par un filtrage moyen ou médian ou passe-bas ou tout autre filtrage adapté au bruit traité de l'image endommagée 31; - C est le coefficient de correction normalisé à 1; - R1, R2 et R3 sont des fonctions de répartition. Les fonctions de répartition doivent respecter un certain nombre de contraintes, à savoir: - R1(C) + R2 (C)+ R3 (C)I 1, - R1 est une fonction décroissante telle que Ru (0) = 1, - R2 est une fonction croissante telle que R2 (1) = 1. Le paramètre R3 (C), qui constitue un paramètre supplémentaire par rapport au premier mode de réalisation présenté auparavant, est par exemple calculé lors d'une étape particulière 29 effectuée après le calcul du coefficient de correction C. Selon cette variante, l'étape de correction de chaque pixel donné de l'image endommagée 31 est réalisée non seulement en fonction des images source et référence, mais également en fonction d'une image source préfiltrée. Optionnellement, le procédé de l'invention peut comprendre une étape préalable 26 (cf fig. 2 et 3) de correction de gain, visant à équilibrer les valeurs d'amplitude des pixels de l'image endommagée 31 avec celles des pixels de l'image de référence 32. L'image endommagée précorrigée 33 ainsi obtenue est utilisée à la place de l'image endommagée 31 dans toutes les autres étapes 21 à 25, 25', 28 du procédé. Cette étape 26 de correction de gain consiste à: - calculer l'écart entre les valeurs moyennes des premier et second histogrammes H1, H2 des niveaux d'amplitude des pixels respectivement de l'image endommagée 31 et de l'image de référence 32, cet écart pouvant s'exprimer sous la forme d'un gain G; - corriger les valeurs d'amplitude des pixels de l'image endommagée 31 en fonction de cet écart de valeur moyenne, en multipliant la valeur d'amplitude de chaque pixel de l'image endommagée 31 par l'inverse du gain G, c'est-à-dire G-l. Le procédé de l'invention n'est pas dépendant du type de bruit affectant les images de la séquence source 2. I1 s'applique notamment, mais non exclusivement, au traitement en temps réel d'images numériques endommagées. REVENDICATIONS 1. Procédé de traitement d'une séquence source (2) d'images numériques endommagées, du type permettant d'obtenir une séquence de sortie (3) d'images numériques corrigées, chacune des images numériques brouilées ou corrigées étant décrite pixel par pixel, chacun desdits pixels étant caractérisé par un niveau d'amplitude parmi une pluralité de niveaux d'amplitude possibles, caractérisé en ce que ledit procédé comprend notamment, pour chaque image numérique endommagée (31) de ladite séquence source (2), les étapes suivantes:: - calcul (21), pour l'ensemble des pixels de l'image endommagée (31), d'un premier paramètre (P1), dit paramètre global, d'estimation de la correction à apporter à chacun des pixels de l'image endommagée, en fonction dudit ensemble des pixels de l'image endommagée (31) et de l'ensemble des pixels d'une image de référence (32), ladite image de référence (32) étant une image corrigée qui résulte de la correction de l'image précédant ladite image endommagée dans ladite séquence source (2); - calcul (22), pour chaque pixel donné (41) de l'image endommagée (31), d'un second paramètre (P2), dit paramètre local, d'estimation de la correction à apporter audit pixel donné, en fonction dudit pixel donné (41) et d'autres pixels, dits pixels voisins (44), localisés dans un voisinage prédéterminé (43) dudit pixel donné (41);; - calcul (23), pour chaque pixel donné (41) de l'image endommagée (31), d'un troisième paramètre (P3), dit paramètre temporel, d'estimation de la correction à apporter audit pixel donné (41), en fonction dudit pixel donné (41) et du pixel (42) de l'image de référence (32) possédant les mêmes coordonnées spatiales (x, y) que ledit pixel donné; - calcul (24), pour chaque pixel donné (41) de l'image endommagée (31), d'un coefficient de correction (C), par combinaison desdits premier, second et troisième paramètres (P1, P2, P3) d'estimation associés audit pixel donné (41);; - correction (25), selon une stratégie de correction prédéterminée, de chaque pixel donné (41) de l'image endommagée (31) de façon à obtenir le pixel correspondant de l'image corrigée de la séquence de sortie (3), en fonction du coefficient de correction (C) associé audit pixel donné (41), du pixel donné (41) de l'image endommagée (31) et du pixel de (42) L'image de référence (32) possédant les mêmes coordonnées spatiales (x, y) que ledit pixel donné (41). 2. Procédé selon la revendication 1, caractérisé en ce que ledit calcul (21), pour l'ensemble des pixels de l'image endommagée (31), du premier paramètre (P1) s'écrit selon la formule suivante: P1 = K + fE (H1, H2) où: - K est une valeur de décalage prédéterminée; - H1 est un premier histogramme des niveaux d'amplitude des pixels de l'image endommagée (31); - H2 est un second histogramme des niveaux d'amplitude des pixels de l'image de référence (32); - fE est une fonction d'erreur prédéterminée permettant de calculer un écart entre deux fonctions. 3. Procédé selon l'une quelconque des revendications 1 et 2, caractérisé en ce que ledit calcul (22) du second paramètre (P2) est effectué simultanément pour tous les pixels de l'image endommagée (31), et s'écrit selon la formule suivante: [Mp2] = (l/a2) . ( [1t] - F2 [1t]) où: - [Mp2] est une matrice des seconds paramètres P2 associés chacun à un pixel distinct de l'image endommagée (31); - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée (31); - F2 est un filtrage prédéterminé; - a2 est un premier facteur de normalisation. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que ledit calcul (23) du troisième paramètre (P3) est effectué simultanément pour tous les pixels de l'image endommagée (31), et s'écrit selon la formule suivante: [Mp3] = (1/a3) F3 F3 ( [1t] - [It-i] J où: - [Mp3] est une matrice des troisièmes paramètres P3 associés chacun à un pixel distinct de l'image endommagée (31); - [It] est une matrice des valeurs d'amplitude des pixels de l'image endommagée (31); - [1t-i] est une matrice des valeurs d'amplitude des pixels de l'image de référence (32); - F3 est un filtrage moyen ou médian ou passe-bas; - a3 est un second facteur de normalisation. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ledit calcul (24) du coefficient de correction (C) s'écrit selon la formule suivante: C = max [1, Pi . f, (P2, P3)1 où: - P1, P2 et P3 sont lesdits premier, second et troisième paramètres normalisés à 1; - fc est une fonction de combinaison prédéterminée. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que ladite stratégie de correction prédéterminée consiste, pour chaque pixel donné (41) de l'image endommagée (31), à calculer une valeur d'amplitude de pixel corrigé It' (x, y ) selon la formule suivante: I'(x,y) = R1 (C).It(x,y)+R2(C).It-1 (x,y) où: - It (x, y ) est la valeur d'amplitude du pixel donné (41) de l'image endommagée (31); - It l (x, y ) est la valeur d'amplitude du pixel (42) de l'image de référence (32) possédant les mêmes coordonnées spatiales (x, y) que le pixel donné (41); - C est le coefficient de correction normalisé à 1; - R1 et R2 sont des fonctions de répartition respectant les contraintes suivantes :R1 (C) + R2 (C) S 1, R1 est une fonction décroissante telle que R1 (0) = 1, R2 est une fonction croissante telle que R2 (1) = 1. 7. Procédé selon selon l'une quelconque des revendications 1 à 5, caractérisé en ce ladite stratégie de correction prédéterminée consiste, pour chaque pixel donné (41) de l'image endommagée (31), à calculer une valeur d'amplitude de pixel corrigé It" (x, y) selon la formule suivante: 4" (x, Y)= R1(C) . 4 (x, y )+ R2(C) est, (x, y) + R3 (C).Ity(x, y) où: - It (x, y ) est la valeur d'amplitude du pixel donné (41) de l'image endommagée (31); - It l (x, y ) est la valeur d'amplitude du pixel (42) de l'image de référence (32) possédant les mêmes coordonnées spatiales (x, y) que le pixel donné (41); - It,F (x, y ) est la valeur d'amplitude du pixel d'une image filtrée possédant les mêmes coordonnées spatiales (x, y) que le pixel donné (41), ladite image filtrée étant obtenue par un filtrage prédéterminé de l'image endommagée (31); - C est le coefficient de correction normalisé à 1; - R1, R2 et R3 sont des fonctions de répartition respectant les contraintes suivantes :R1 (C) + R2 (C) + R3 (C) telle que Ru (0) = 1, R2 est une fonction croissante telle que R2 (1) = 1. 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce qu'il comprend, pour chaque pixel donné (41) de l'image endommagée (31), les étapes supplémentaires suivantes: - calcul d'une probalité d'erreur sur la valeur d'amplitude du pixel donné (41), en fonction de l'écart entre le nombre de pixel de l'image endommagée (Hl(n)) et le nombre de pixel de l'image de référence (H2(n)) possédant la même valeur d'amplitude (n) que ledit pixel donné (41), ladite probabilité d'erreur constituant un quatrième paramètre (P4) d'estimation de la correction à apporter audit pixel donné - pondération (61; 71) par ledit quatrième paramètre (P4) des premier et second paramètres (P1, F2) ou du coefficient de correction (C) associés audit pixel donné (41) 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce qu'il comprend les étapes préalables suivantes: - calcul d'un premier histogramme (H1) des niveaux d'amplitude des pixels de l'image endommagée (31); - calcul d'un second histogramme (H2) des niveaux d'amplitude des pixels de l'image de référence (32); - calcul de l'écart entre les valeurs moyennes desdits premier et second histogrammes; ; - correction des valeurs d'amplitude des pixels de l'image endommagée (31) en fonction dudit écart de valeur moyenne, de façon à équilibrer les valeurs d'amplitude des pixels de l'image endommagée (31) avec celles des pixels de l'image de référence (32), et obtenir une image endommagée précorrigée (33) qui est utilisée à la place de l'image endommagée (31) dans toutes les autres étapes (21 à 25, 25', 28) dudit procédé. 10. Application du procédé selon l'une quelconque des revendications 1 à 9 au traitement en temps réel d'une séquence source d'images numériques endommagées.