FR 2840087 A1 20031128 FR 0206208 A 20020522 DEPOLLUTION D'IMAGE NUMERIQUE BRUITEE ET COMPRIMEE. DOMAINE TECHNIQUE GENERAL. La présente invention concerne un procédé de dépollution d'image numérique comprimée et bruitée. Plus précisément, elle concerne un procédé permettant de diminuer jusqu'à une valeur choisie la pollution des images numériques comprimées. ETAT DE L'ART. Lors de l'acquisition d'images par des instruments d'acquisition et lors de la compression des images acquises, un certain nombre de bruits s'additionnent au signal correspondant à l'image. En effet, d'une part les instruments d'acquisition ne sont jamais parfaits. Ils introduisent des signaux parasites que l'on appelle bruit instrumental. D'autre part, une compression d'une image consiste en l'attribution d'une valeur uniforme du signal sur un intervalle donné. Elle introduit donc une erreur sur le signal. On la qualifie d'uniforme et on l'appelle bruit de compression ou de quantification. On cherche à diminuer le bruit instrumental et le bruit de quantification dans les images comprimées. Les procédés de débruitage des images comprimées de l'état de l'art permettent le débruitage du bruit instrumental. Il existe une multitude de tels procédés de débruitage des bruits instrumentaux. Certains d'entre eux utilisent la décomposition du signal comprimé bruité dans une base choisie de l'espace de fréquentiel de compression. Cette base peut être une base DCT ( Discrete Cosinus Transform selon la terminologie anglo-saxonne) ou dans une base d'ondelettes. Les techniques précédentes présentent cependant des inconvénients. En effet, elles ne permettent pas d'effectuer le débruitage du bruit de compression en même temps que le débruitage instrumental. PRESENTATION DE L'INVENTION. L'invention propose de pallier ces inconvénients. Un des buts de l'invention est de proposer un procédé qui permet de diminuer jusqu'à une valeur choisie la pollution des images numériques comprimées. La dépollution de l'image élimine le bruit de compression et le bruit instrumental par un même procédé. Un autre but de l'invention est de proposer un procédé de détermination de la meilleure base d'expression de l'image numérique pour cette dépollution. Un autre but de l'invention est de proposer une méthode de suréchantillonnage pour éviter la distorsion par repliement des images dépolluées. A cet effet, l'invention propose un procédé de dépollution d'image numérique bruitée et comprimée caractérisé en ce qu'on effectue le déflouage de ladite image en effectuant une déconvolution du signal et en effectuant le débruitage du bruit de compression et du bruit instrumental de l'image en exprimant le signal numérique de l'image dans une base de paquets d'ondelettes sur laquelle ces deux bruits sont gaussiens. L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leur combinaison techniquement possible : - la détermination de la base dans laquelle les bruits de compression et instrumental sont gaussiens est effectuée durant un apprentissage des caractéristiques de la chaîne instrumental associée à l'image et à sa compression, l'apprentissage comportant les étapes selon lesquelles : - on construit un tableau de coefficients, ces coefficients étant fonction de l'effet de la déconvolution de déflouage sur un signal de bruit moyen exprimé dans une base de Walsh généralisée ; - on construit une base de décomposition la moins décomposée en comparant la valeur des coefficients de l'effet de cette déconvolution à une valeur seuil de finesse ; - on construit une base d'ondelettes permettant le meilleur débruitage du bruit moyen décomposé dans la base de Walsh généralisée ; - on construit la meilleure base en effectuant la fusion des deux bases précédentes, le coefficient de décomposition le plus fin de chaque base étant conservé pour la construction de cette meilleure base. - pour construire le tableau de coefficients fonction de l'effet de la déconvolution, on effectue les opérations suivantes : - on décompresse un signal de bruit uniforme moyen ; - on ajoute audit bruit décompressé un bruit instrumental moyen pour obtenir un signal de bruit total ; - on effectue d'une part une déconvolution du signal de bruit de compression uniforme moyen avant de l'exprimer dans une base de Walsh généralisée ; - on effectue d'autre part une décomposition du signal de bruit total non déconvolué dans la même base de Walsh généralisée ; - on construit les deux tableaux des normes L2des différents coefficients trouvés ; - on construit un tableau résultant à partir du rapport des normes des coefficients. - pour construire la base de décomposition la moins décomposée, on effectue les opérations suivantes : - on calcule une valeur de finesse en décomposant l'image du bruit déconvoluée en base de Walsh d'une profondeur d'analyse immédiatement supérieure à l'analyse déjà effectuée pour la construction du tableau de l'effet de la déconvolution ; - on compare à la valeur de finesse choisie, à travers tous les niveaux de décomposition en base de Walsh de tous les paquets, d'abord en ligne puis en colonne, le rapport en norme L2des coefficients du tableau de la décomposition du bruit total déconvolué pour chaque paire de paquets frères associés ; - si le rapport desdits coefficients est inférieur à la valeur de la finesse, on recompose ces paquets frères en un seul paquet ; - on effectue le procédé de recomposition avec des valeurs de la finesse de plus en plus importantes jusqu'à quasi éliminer les paquets décomposés au dernier niveau choisi, cette valeur de la finesse étant ensuite utilisée au niveau immédiatement inférieur. - l'étape de dépollution comporte une étape de détection du module des coefficients de la décomposition de l'image dans la meilleure base par rapport à une valeur d'un seuil, suivie d'une étape de traitement de ces coefficients ; - on calcule la valeur du seuil, cette valeur étant fonction de la valeur de l'écart type du bruit en tout point de l'image, d'un paramètre de seuillage, de la valeur du filtre de déconvolution, de la valeur en niveau de gris de l'image et de la valeur locale de l'intervalle de quantification ; - les petits coefficients par rapport au seuil sont multipliés par un facteur de débruitage et les forts coefficients sont déconvolués ; - on effectue en outre une étape de suréchantillonnage des paquets avant le débruitage et/ou la déconvolution afin d'éviter le phénomène de distorsion de repliement ; - le suréchantillonnage comporte les étapes selon lesquelles, sur chaque paquet : - on décale d'une ligne, d'une colonne et d'une ligne et d'une colonne, les paquets au niveau immédiatement inférieur de la meilleure base ; - on décompose les quatre paquets ainsi trouvés ; - on les débruite au niveau immédiatement supérieur après décomposition, une déconvolution pouvant en outre à cette étape de débruitage ; - on recompose les paquets ainsi obtenus ; - on moyenne les divers paquets après recalage. PRESENTATION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit qui est purement illustrative et non limitative et qui doit être lue en regard des dessins annexés sur lesquels : - la figure 1 représente schématiquement une base de Walsh généralisée ; - la figure 2 représente schématiquement une base de paquets d'ondelettes, - la figure 3 représente schématiquement un exemple de base de décomposition la moins décomposée trouvée pendant l'étape d'apprentissage ; et - la figure 4 représente schématiquement un exemple de meilleure base issue de la fusion des bases des figures 2 et 3. DESCRIPTION DETAILLEE DE L'INVENTION. Le procédé de débruitage selon l'invention utilise trois informations en entrée. La première entrée est l'image image polluée (bruitée et compressée). La deuxième entrée est le doublet (FTM, bruit instrumental) relatif à l'image bruitée et floue. Le premier élément du doublet est la fonction FTM (Fonction de Transfert de Modulation). C'est une fonction de l'espace fréquentiel. Elle représente l'aptitude du système instrumental à transmettre les différentes fréquences. Elle est la fonction de flou caractérisant les pertes de transfert instrumentales. Le deuxième élément du doublet est caractérisé par son écart type local. Le bruit de compression est non stationnaire car son écart type dépend localement du facteur F de quantification de la compression. La troisième entrée est l'ensemble des facteurs de quantification relatifs à la compression dans l'espace transformé. L'espace transformé peut être un espace DCT (Discrete Cosinus Transform), un espace d'ondelettes, ou un espace d'une autre transformation de compression. De plus, les images traitées peuvent être définies sur une grille d'échantillonnage en géométrie carrée ou quinconce (entrelacement de deux sous-images d'échantillonnage carré décalées d'un demi pas en ligne et en colonne). Le procédé utilise l'ensemble de ces données pour fournir in fine une image déflouée et débruitée dont on peut diminuer à volonté le niveau de pollution uniforme sur les zones uniformes. On appelle pollution uniforme une pollution du type bruit blanc. On rappelle que la compression concentre l'information contenue dans une image sur un petit nombre de coefficients par transformation de l'image dans un espace où la représentation de l'image est compacte, c'est- à-dire décorrélée. La décorrélation fait référence au fait que l'expression de l'image sur une base de transformation permet d'obtenir des coefficients dont peu ont une grosse valeur, porteuse d'information. La plupart des coefficients de la décomposition ont un faible coefficient, et peuvent donc être éliminés par la compression. Lors de la compression, tout coefficient appartenant à un intervalle de taille F est représenté par un même point appartenant à cet intervalle. Ceci peut être formalisé par l'ajout d'un bruit uniforme de quantification local d'intervalle F(l',p') dans l'espace transformé des coefficients. (l',p') représente la position dans l'espace transformé. D'autre part, l'écart type du bruit instrumental dépend de la position (l, p) dans l'image. Son expression est la suivante. 2(l,p) = A + B x L(l, p) où A et B sont des constantes connues pour les images considérées, et L(l,p) est le niveau de gris de l'image au point (l,p). Il est constitué d'un processus de Poisson pour la partie B#L(l,p) et d'un processus gaussien pour la partie A. Le procédé de dépollution est fondé sur le résultat suivant : le bruit de compression ou de quantification, qui est uniforme lorsqu'il est exprimé dans une base DCT, ondelettes ou autre de l'espace transformé, devient quasi Gaussien s'il est exprimé dans l'espace image, ou une autre base d'ondelettes par exemple, que celle choisie pour la compression multirésolution. De plus, le bruit de compression est décorrélé du bruit instrumental. Le procédé de dépollution des images comprimées selon l'invention est globalement formalisé par la transformation des doublets caractérisant respectivement l'image brute et l'image dépolluée et déflouée caractérisée par la transformation : (FTM, bruit -blanc(l, p)) -> (P, k x bruit -blahc(l, p)) . La déconvolution ou déflouage (passage de FTM à P) permet d'obtenir une image déflouée. L'image déflouée possède une fonction de transfert rehaussée P, le bruit reste blanc sur les zones uniformes et son écart type local est contrôlé par le coefficient de débruitage k. On remarque que l'écart-type du bruit dépend de la position sur l'image. Par contre, la déconvolution est globale et elle ne dépend pas de l'image. Le déflouage de l'image consiste en une opération de déconvolution de l'image par filtrage. Cette opération augmente particulièrement les Hautes Fréquences (HF) qui avaient été en priorités affaissées par la FTM instrumentale. L'information HF contenue dans l'image est rehaussée, mais également la partie HF du bruit blanc qui se transforme en un bruit coloré. Le procédé selon l'invention traite par reblanchissement du bruit des images déconvoluées. Le procédé selon l'invention prend en compte les deux espaces (espace image et espace transformé). On considère toujours les deux paramètres espace-fréquence. En effet, le bruit dépend d'une part de la localisation spatiale. D'autre part, la décoloration du bruit - suite à la déconvolution - dépend de son spectre fréquentiel. Les ondelettes présentent une bonne représentation spatiale mais sont spectralement mal résolues dans les hautes fréquences. Elles ne conviennent donc pas pour la décomposition du bruit dès que le débruitage est couplé à la déconvolution. Il faut donc trouver une autre base que les ondelettes pour la représentation du bruit. Une bonne représentation conjointe (espace et fréquence) est réalisée par la transformation en paquets d'ondelettes. Elle consiste en une décomposition dyadique complète de l'image en alternant des convolutions par un filtre passe-bas et passe-haut séparable suivant les lignes et suivant les colonnes. Le procédé de décomposition selon l'invention est itératif. Le nombre d'itérations s'appelle le nombre d'échelles de décomposition choisi. Il correspond à une profondeur d'analyse. La décomposition en paquets est très souple. En effet, le nombre de bases possibles dans la décomposition en paquets est immense. Le nombre de paquets varie en croissance exponentielle avec le nombre d'échelles. De plus, chaque paquet à chaque échelle peut être ou ne pas être décomposé. La figure 1 représente schématiquement une analyse en paquets où tous les paquets sont décomposés jusqu'au niveau 4, c'est-à-dire en effectuant 4 itérations. Cette base particulière uniforme en décomposition (chaque paquet à été décomposé jusqu'au dernier niveau) est nommée base de Walsh généralisée. Chaque paquet de cette base peut être interprété comme représentant essentiellement une zone spectrale de l'espace de Fourier pourvue que l'ondelette choisie soit régulière. Cette ondelette régulière peut être du type spline, Daubechies bi-orthogonal avec beaucoup de moments nuls par exemple. Un des paquets représente la position. Ce paquet correspond à la seule décomposition de l'image à l'aide des passes-bas. Il est appelé paquet résumé de l'image. Il s'agit en fait d'une image sous résolue construite à partir d'un moyennage pondéré de l'image de départ. La figure 2 représente schématiquement la base d'ondelettes qui pour un niveau de décomposition donné, par exemple de 4, débruite le mieux un bruit blanc se sur-rajoutant à l'image. Cette base est générée par la seule décomposition à chaque niveau du paquet basse fréquence. Elle possède donc une très faible résolution spectrale pour les hautes fréquences. La base d'ondelettes ne permet pas le blanchissement du bruit introduit par la déconvolution. La meilleure base permettant le blanchissement du bruit introduit par la déconvolution est celle qui utilise des paquets d'ondelettes. Au final, la meilleure base pour la déconvolution (qui correspond au déflouage) puis le débruitage (qui correspond à la dépollution) doit répondre à deux problèmes complémentaires. La meilleure base est celle où chaque paquet d'une part épouse au mieux le spectre du filtre de déconvolution, et d'autre part est de niveau de décomposition minimal. La première caractéristique permet le blanchissement du bruit de déconvolution. La deuxième permet de rester proche de la décomposition en ondelettes qui permet de débruiter le bruit blanc se sur rajoutant à l'image. Ainsi le procédé de dépollution des images comprimées selon l'invention est fondé sur la décomposition en paquets d'ondelettes, pour prendre en compte le débruitage, et un seuillage des petits coefficients de l'image sur ladite base, pour prendre en compte la déconvolution. La prise en compte de l'échantillonnage quinconce peut être réalisée directement par une mise à 0 des coefficients de Fourier du filtre de déconvolution en dehors de la cellule réciproque quinconce de l'échantillonnage. PRINCIPE DU PROCEDE. La dépollution des images se décompose en deux étapes. La première étape est constituée de l'étape d'apprentissage sur une image de bruit moyen généré par la chaîne image et la seconde étape est l'opération de dépollution de l'image proprement dite. APPRENTISSAGE. L'étape d'apprentissage a deux objectifs principaux. Le premier objectif est le calcul de la meilleure base de décomposition de l'image comprimée. Cette meilleure base de décomposition comporte des paquets d'ondelettes. Le deuxième objectif est le calcul des rapports de décomposition en chaque paquet d'ondelettes. L'apprentissage prend en compte l'ensemble des bruits indépendants de l'image. Ils peuvent donc être modélisés une fois pour toutes pour chaque instrument. ETAPES PRELIMINAIRES. L'apprentissage s'effectue sur une image quelconque de bruit moyen. Dans le cas général, on génère un bruit blanc uniforme, qui correspond à une quantification moyenne, dans la base de représentation des coefficients de la compression. On décompresse ce bruit. On lui ajoute un bruit instrumental moyen. Bien entendu, si l'image est constituée de sous images l'opération précédente est appliquée à chacune d'entre elles et les sous images de bruit sont entrelacées pour constituer l'image de bruit final. La somme du bruit uniforme décompressé et du bruit instrumental moyen génère ainsi un écart type du bruit total Gaussien (ou quasi Gaussien) d'expression : #@2=k (A+B##+#2/12). L'étape d'apprentissage permet de connaître les caractéristiques des appareils de la chaîne d'acquisition d'image, représentées par la FTM. Grâce à la connaissance de la FTM, on calcule un filtre de déconvolution PIFTM où P est une fonction de flou objectif. Typiquement, la fonction P peut être la fonction sphéroïdale aplatie en anglais prolate - de support spectral égal à la cellule réciproque de l'échantillonnage de l'image. Selon le procédé de l'invention, on peut générer une image du bruit total issu de la somme de taille 1024x1024. Cette taille permet d'effectuer des statistiques sur les paquets jusqu'à la dernière échelle de décomposition. Le nombre d'échelle ne dépasse pas 5. Par conséquent, la taille d'un paquet au niveau 5 est de 32x32 (25x25). Cette image de bruit total est ensuite déconvoluée par le filtre de déconvolution. On réalise ensuite une décomposition simultanée de l'image du bruit total et de l'image du bruit total déconvoluée. La décomposition des deux images s'effectue dans la base de Walsh. La décomposition n'excède pas un niveau 5. Les normes L2de chaque paquet pour chaque image de bruit sont alors calculées. On appelle Tab~bruit~L2~Walsh le tableau des normes L2des paquets de l'image de bruit, et Tab~bruit~deconv~L2-Walsh celui correspondant à l'image de bruit déconvoluée. On construit ensuite un autre tableau appelé Tab deconv Walsh dans lequel on place le rapport des normes L2des paquets correspondants. En chaque case du tableau, on a donc le rapport : norme image de bruit déconvoluée norme image de bruit pour chaque paquet correspondant. Par ailleurs, on calcule le paquet résumé de l'image. On normalise ensuite le tableau Tab~deconv~Walsh par la norme L2 du paquet résumé de même niveau. Les paquets représentent essentiellement une zone fréquentielle précise de l'espace spectral. Par conséquent, le tableau Tab~deconv~Walsh représente l'effet du filtre de déconvolution dans la décomposition dyadique en base de Walsh. On constate que pour une base de Walsh de profondeur 4, une analyse de Fourier, qui correspond à la meilleure déconvolution possible, et les résultats donnés par le tableau Tab~deconv~Walsh, donnent des fonctions très proches. Ce résultat justifie la possibilité d'effectuer la déconvolution dans les paquets d'ondelettes. CALCUL DE LA MEILLEURE BASE. Le principe de calcul de la meilleure base se décompose en deux temps. On recherche d'abord la base de déconvolution la moins décomposée. Dans un deuxième temps, on effectue une fusion avec la meilleure base de débruitage : la base d'ondelettes. Le premier temps consiste à se donner une profondeur d'analyse en échelle (typiquement 4) et un coefficient de finesse Deconv~finesse d'analyse du filtre de déconvolution. Ce premier temps est nécessaire du fait que la fonction de déconvolution ne doit pas varier de façon trop brutale entre différents paquets successifs. Lorsque la variation d'un paquet à l'autre est trop brutale, la décomposition en paquets doit être plus fine. La valeur de Deconv finesse est calculée tout d'abord en décomposant l'image du bruit déconvoluée en base de Walsh d'une profondeur d'analyse immédiatement supérieure à l'analyse déjà effectuée. On effectue ainsi typiquement une décomposition de Walsh au niveau 5 si l'on s'intéresse au niveau 4. Par ailleurs, les paquets d'ondelettes sont reliés dans l'analyse ligne et colonne par paires. Ils sont appelés paquets frères. Les étapes du premier temps sont les suivantes. On examine à travers tous les niveaux de décomposition et tous les paquets, d'abord en ligne puis en colonne, et en partant de la base de Walsh, si deux paquets frères sont tels que leur rapport - en norme L2 - de leur coefficient du tableau Tab~deconv~Walsh ou l'inverse de leur rapport 1/Tab~deconv~Walsh est inférieur à Deconv finesse. On compare donc la valeur du rapport : Tab~deconv~Walsh~premier paquet Tab~deconv~Walsh~deuxième paquet ou l'inverse de ce rapport à la valeur du coefficient de finesse que l'on s'est fixée. Si le rapport est inférieur au coefficient de finesse, on recompose ces paquets frères en un seul paquet. Ce paquet résultant peut être rectangulaire, avec un niveau d'analyse différent en ligne et en colonne. Ensuite, on effectue le procédé de recomposition avec des valeurs de Deconv finesse de plus en plus importantes jusqu'à quasi éliminer les paquets décomposés au dernier niveau (typiquement 5). Cette dernière valeur de Deconv finesse est choisie pour l'analyse en paquet au niveau inférieur (typiquement 4). La figure 3 représente schématiquement la base de déconvolution la moins décomposée trouvée lors du premier temps. Le second temps de la détermination de la meilleure base consiste en une union entre la base de déconvolution et la base d'ondelettes. Dans notre exemple, la base d'ondelettes à considérer est représentée à la figure 2. On conserve in fine la décomposition la plus fine entre ces deux bases dans la décomposition dyadique. On obtient donc la base de décomposition représentée à la figure 4. On a donc atteint un premier objectif, à savoir le calcul de la meilleure base de décomposition. On rappelle que l'apprentissage doit également remplir un autre objectif, à savoir le calcul des rapports de décomposition en chaque paquet d'ondelettes. CALCUL DES RAPPORTS DE DECOMPOSITION. On calcule, comme on l'a fait pour la base de Walsh, le tableau des rapports des normes L2des paquets correspondants dans la meilleure base trouvée. On effectue alors les mêmes opérations que précédemment. On génère un bruit blanc uniforme, qui correspond à une quantification moyenne, dans la base de représentation des coefficients de la compression. On décompresse ce bruit. On ajoute un bruit instrumental moyen à l'image du bruit uniforme décompressée. On garde d'une part cette image. On déconvolue d'autre part cette image de bruit par le filtre P/FTM. On construit ensuite les deux tableaux des normes L2des coefficients après avoir décomposé les images des deux bruits dans la meilleure base trouvée figure 4. On appelle Tab bruit deconv L2 mb et Tab~bruit~deconv~L2~mb ces deux tableaux. On appelle Tab~deconv~mb le tableau des rapports correspondants. Ce tableau interviendra dans le débruitage de l'image. Après l'étape d'apprentissage, on effectue la dépollution proprement dite de l'image. DEPOLLUTION DE L'IMAGE. On rappelle que la dépollution de l'image est constituée d'une étape de détection du module des coefficients de la décomposition de l'image dans la meilleure base par rapport à une valeur d'un seuil, suivie d'une étape de traitement de ces coefficients. CALCUL DU SEUIL - PRESENTATION DES VARIABLES. Premièrement, le calcul du seuil nécessite le calcul de la valeur de l'écart type du bruit en tout point de l'image. On a donc besoin de la connaissance locale des facteurs de quantification et du niveau de bruit de l'image. Deuxièmement, la dépollution sera essentiellement gouvernée par le paramètre de seuillage S. Il est choisi par l'utilisateur. Il a une valeur comprise entre 0 et 3. 5. La valeur 3. 5 correspond à un fort débruitage, tandis que la valeur 0 correspond à une déconvolution sans débruitage. Par ailleurs, la valeur du seuil des petits coefficients après déconvolution par le filtre spectral de rapport PIFTM dépend de la représentation spectrale correspondante du paquet. Elle est donc quasi proportionnelle à la valeur spectrale correspondante du filtre PIFTM. La valeur du seuil des petits coefficients dépend également de la valeur en niveau de gris de l'image et de la valeur locale de l'intervalle de quantification. Ce caractère local peut être facilement appréhendé car le bruit est quasi Gaussien et la valeur du seuil est proportionnelle à l'écart type du bruit, calculable en tout point (l,p) dès que l'on connaît la valeur F du facteur de quantification et de L, le niveau de gris en ce point. Ainsi, le seuil est proportionnel à A+B#L (l,p)+F2(l,p)/12 A+B##+#2/12 CALCUL DU RAPPORT DE L'ECART TYPE LOCAL. On décompose l'image comprimée en paquets dans la meilleure base. On calcule par ailleurs les résumés des tableaux des facteurs de quantification de la compression à tous les niveaux de décomposition. On calcule par ailleurs, par le même procédé de décomposition, le résumé de l'image. Grâce à ce résumé de l'image et à celui des facteurs de quantification, il est alors possible de calculer le rapport de l'écart type local : k~sig(l', p')= A+BxL(l',p')+F2(l',p')l12 ~ ~ ~SIg A+B##+#2/12 où (l', p') est la position dans les 2 résumés (image et tableau des facteurs de quantification) correspondante au point à dépolluer dans les paquets. DETECTION DU MODULE DES COEFFICIENTS ET TRAITEMENT. Une fois que les petits coefficients par rapport au seuil sont calculés, trois types de traitement sont proposés. La déconvolution est effectuée en Fourier avant ou après le débruitage. La déconvolution ou le débruitage peuvent également être effectués simultanément dans la décomposition en paquets. Le premier type de traitement est un débruitage dans les paquets puis une déconvolution dans Fourier. La détection des petits coefficients est effectuée par le calcul du seuil Seuil=S#k~sig#Tab~bruit~L2. Les petits coefficients (inférieurs à Seuil) sont comprimés dans le rapport : 1/ sig Tab - bruit - deconv - 1 (k -sig 7o& ~ M/y ~ Le second type de traitement est une déconvolution dans Fourier puis un débruitage dans les paquets. La détection des petits coefficients est effectuée par le calcul du seuil : Seuil = S x k-sig x Tab-bruit-deconv-L2 . Les petits coefficients (inférieurs à Seuil) sont comprimés dans le rapport : 1/(k~ sig x Tab - bruit ~ deconv - L2) ' Le troisième type de traitement est un débruitage et une déconvolution dans les paquets. La détection des petits coefficients est effectuée par le calcul du même seuil que précédemment : Seuil = S x k - sig x Tab - bruit - deconv ~ L2 . Le traitement se sépare ici en deux parties. Les petits coefficients sont comprimés dans le rapport : 1/ sig Tab - bruit - deconv - l (k -sig x Tab -bruit ~ deconv- L2) Les gros coefficients sont dilatés par multiplication par le tableau : Tab bruit deconv L2. On rappelle que la valeur du paramètre de seuillage est comprise entre 0 et 3.5. La valeur 0 correspond à la seule déconvolution, on n'applique alors pas de débruitage. La valeur 3. 5 correspond à un débruitage sur la quasi totalité des petits coefficients (ceci en conformité avec les propriétés d'une distribution quasi Gaussienne). Les valeurs intermédiaires du paramètre de seuillage S pondèrent le traitement entre ces deux extrêmes. Le coefficient de débruitage k sur les zones uniformes est en général maintenu égal à 1. PHENOMENE DE DISTORSION DE REPLIEMENT. On rappelle que la distorsion de repliement, ou aliasing selon la terminologie anglo-saxonne généralement employée, est une distorsion d'un signal, qui se produit lors de son échantillonnage à une fréquence d'échantillonnage inappropriée. La distorsion de repliement est due à un chevauchement des bandes latérales créées dans le spectre du signal échantillonné. Les paquets obtenus après le seuillage sont aliasés. Ils subissent un sous échantillonnage dans un rapport 2. Une amélioration pour éviter les artefacts provoqués par le phénomène d'aliasing lors du seuillage est d'effectuer un suréchantillonnage des paquets avant débruitage (et éventuellement déconvolution). Le seuillage est alors appliqué seulement au dernier niveau de décomposition de la meilleure base. PROCEDE DE SURECHANTILLONNAGE. La procédure consiste à décaler les paquets au niveau immédiatement inférieur de la meilleure base, les débruiter au niveau immédiatement supérieur après décomposition, recomposer et moyenner les divers paquets après recalage. Pour chaque paquet du niveau inférieur de profondeur d'analyse correspondant à la meilleure base, on décale le paquet ainsi trouvé d'une ligne, d'une colonne et d'une ligne et d'une colonne. On décompose les 4 paquets ainsi trouvés (si l'on tient compte également du paquet non décalé). On débruite (et déconvolue) ces 4 paquets. On les recompose et on les recale. Le résultat est constitué d'un moyennage de ces 4 paquets recomposés au niveau inférieur et débruités au niveau immédiatement supérieur. Ceci est appliqué sur chaque paquet. Il ne reste plus alors qu'à parachever la recomposition pour obtenir l'image dépolluée. Tous les développements qui précèdent s'appliquent avantageusement à une image comprimée dans un espace DCT. EXCEPTION DANS LE CAS D'UNE COMPRESSION EN ONDELETTES. Dans le cas où la compression est réalisée par ondelettes et non pas par DCT, la compression multirésolution et la restauration sont effectuées dans des bases d'ondelettes ou de paquets. Une solution proposée par l'invention est de choisir pour la compression la meilleure base définie plus haut. La compression s'effectue sur les coefficients dans une meilleure base de profondeur donnée (typiquement 3) avec une ondelette de faible nombre de moments nuls. Il s'agit en général d'ondelettes Daubechies biorthogonales. Avantageusement, le facteur de quantification choisi en compression multirésolution tient compte du filtre de déconvolution. Chaque paquet ou sous bande est quantifié par un facteur inversement proportionnel à la valeur du filtre de déconvolution relatif à la fréquence moyenne représenté par le paquet. Ainsi, le bruitage par compression restera-t-il blanc après déconvolution. Cette partie ne dépend pas de l'image. Comme précédemment une seconde partie dépend de l'image par l'intermédiaire d'un tableau de facteurs de quantification. Le débruitage-déconvolution s'effectue, comme nous l'avons précédemment indiqué, avec une ondelette spline ou une ondelette de Daubechies possédant un nombre important de moments nuls. Le procédé de dépollution consiste à effectuer le débruitage (et la déconvolution) par changement de base en décomposant chaque paquet compressé d'une profondeur de 1 ou 2 niveaux supplémentaires avec une ondelette spline ou une Daubechies possédant un nombre important de moments nuls. APPRENTISSAGE. La partie apprentissage sera effectuée sur les bruits indépendants de l'image. On a ainsi un bruit instrumental et un bruit de quantification des paquets (sous bandes) par bruitage dans la base de compression (typiquement niveau 3) avec un bruit uniforme d'intervalle égal à la valeur de la quantification correspondante à chaque paquet. Soit C(vl,vp) le bruit global non blanc correspondant au bruit de quantification des paquets (sous bandes). Alors, le bruit d'apprentissage est quasi Gaussien et la variance égale à : A+B##+C (vl,vp). Les normes L2des bruits déconvolués et non déconvolués sont alors calculées comme précédemment. Cette formulation du bruit permet d'effectuer la détection des petits coefficients en déterminant le seuil pour chaque paquet et chaque point. Ce bruit est ensuite analysé par décomposition d'une profondeur de 1 ou 2 niveaux supplémentaires (base de débruitage et de déconvolution). Le seuillage (et la déconvolution) de l'image est effectué dans cette nouvelle base. DEPOLLUTION. La dépollution de l'image est réalisée en tenant compte du tableau des facteurs de quantifications F propres à l'image par calcul du résumé du tableau des F. Pour la détection nous utiliserons le facteur : k~sig{l\p\vhv +xZ(/',')+C( )+(/'')/12 k sig(I', p,,, VI, vp) A + B x L(1', p') C(vl, v,) + F' p') /.)=J#######=####-####### . k~Slg(l,p A+BxL - +C(V,vp)+F2 - l12 et pour la dépollution : - . , , A+BxL(l',p')+C(V,vp)+F2(l',p')112 k~Sig(l',p',vl,vp)= A+BxL(l',p')+C(vl,v ~ ~ p~2 )+F2(l',p')l12 # A+B##+#+#2/12 Ceci permet d'atteindre le bruit blanc #A+B##+#+#2/12. Cette opération est effectuée dans la base de débruitage déconvolution et comme dans la partie précédente avec décalage. TRAITEMENT. Ici aussi trois types de traitement sont possibles. Le premier type de traitement est comme précédemment un débruitage dans les paquets puis une déconvolution dans Fourier. La détection des petits coefficients est effectuée par le calcul du seuil : Seuil = S x k -sig(l', p', v" vn ) x Tab -bruit -L2 . Les petits coefficients sont comprimés dans le rapport : 1/ (k~sig#Tab~bruit~deconv~L2). Le second type de traitement est une déconvolution dans Fourier puis un débruitage dans les paquets. La détection des petits coefficients est effectuée par le calcul du seuil : Seuil = S x k -sig(l', p', v" v ) x Tab-bruit-deconv- L2 . Les petits coefficients sont comprimés dans le rapport 1/~ 1/(k~sig x Tab~bruit~deconv- L2) Le troisième type de traitement est un débruitage et une déconvolution dans les paquets. La détection des petits coefficients est effectuée par le calcul du même seuil que précédemment : Seuil =S#k~sig(l',p',vl,vp)#Tab~bruit~L2 Le traitement se sépare ici en deux parties. Les petits coefficients sont comprimés dans le rapport : 1/ (k~sig) alors que les gros coefficients sont dilatés par multiplication par le tableau : Tab~bruit~deconv~L2. REVENDICATIONS. 1. Procédé de dépollution d'image numérique bruitée et comprimée caractérisé en ce qu'on effectue le déflouage de ladite image en effectuant une déconvolution du signal et en effectuant le débruitage du bruit de compression et du bruit instrumental de l'image en exprimant le signal numérique de l'image dans une base de paquets d'ondelettes sur laquelle ces deux bruits sont gaussiens. 2. Procédé selon la revendication 1, caractérisé en ce que la détermination de la base dans laquelle les bruits de compression et instrumental sont gaussiens est effectuée durant un apprentissage des caractéristiques de la chaîne instrumental associée à l'image et à sa compression, l'apprentissage comportant les étapes selon lesquelles : - on construit un tableau de coefficients, ces coefficients étant fonction de l'effet de la déconvolution de déflouage sur un signal de bruit moyen exprimé dans une base de Walsh généralisée ; - on construit une base de décomposition la moins décomposée en comparant la valeur des coefficients de l'effet de cette déconvolution à une valeur seuil de finesse ; - on construit une base d'ondelettes permettant le meilleur débruitage du bruit moyen décomposé dans la base de Walsh généralisée ; - on construit la meilleure base en effectuant la fusion des deux bases précédentes, le coefficient de décomposition le plus fin de chaque base étant conservé pour la construction de cette meilleure base. 3. Procédé selon la revendication 2, caractérisé en ce que pour construire le tableau de coefficients fonction de l'effet de la déconvolution, on effectue les opérations suivantes : - on décompresse un signal de bruit uniforme moyen ; - on ajoute audit bruit décompressé un bruit instrumental moyen pour obtenir un signal de bruit total ; - on effectue d'une part une déconvolution du signal de bruit de compression uniforme moyen avant de l'exprimer dans une base de Walsh généralisée ; - on effectue d'autre part une décomposition du signal de bruit total non déconvolué dans la même base de Walsh généralisée ; - on construit les deux tableaux des normes L2des différents coefficients trouvés ; - on construit un tableau résultant à partir du rapport des normes des coefficients. 4. Procédé selon les revendications 2 et 3, caractérisé en ce que pour construire la base de décomposition la moins décomposée, on effectue les opérations suivantes : - on calcule une valeur de finesse en décomposant l'image du bruit déconvoluée en base de Walsh d'une profondeur d'analyse immédiatement supérieure à l'analyse déjà effectuée pour la construction du tableau de l'effet de la déconvolution ; - on compare à la valeur de finesse choisie, à travers tous les niveaux de décomposition en base de Walsh de tous les paquets, d'abord en ligne puis en colonne, le rapport en norme L2des coefficients du tableau de la décomposition du bruit total déconvolué pour chaque paire de paquets frères associés ; - si le rapport desdits coefficients est inférieur à la valeur de la finesse, on recompose ces paquets frères en un seul paquet ; - on effectue le procédé de recomposition avec des valeurs de la finesse de plus en plus importantes jusqu'à quasi éliminer les paquets décomposés au dernier niveau choisi cette valeur de la finesse étant ensuite utilisée au niveau immédiatement inférieur. 5. Procédé selon la revendication 1, caractérisé en ce que l'étape de dépollution comporte une étape de détection du module des coefficients de la décomposition de l'image dans la meilleure base par rapport à une valeur d'un seuil, suivie d'une étape de traitement de ces coefficients. 6. Procédé selon la revendication 5, caractérisé en ce qu'on calcule la valeur du seuil, cette valeur étant fonction de la valeur de l'écart type du bruit en tout point de l'image, d'un paramètre de seuillage, de la valeur du filtre de déconvolution, de la valeur en niveau de gris de l'image et de la valeur locale de l'intervalle de quantification. 7. Procédé selon la revendication 6, caractérisé en ce qu'on calcule la valeur de l'écart type en fonction de la connaissance locale des facteurs de quantification et du niveau de bruit de l'image. 8. Procédé selon la revendication 5, caractérisé en ce que les petits coefficients par rapport au seuil sont multipliés par un facteur de débruitage, et les gros coefficients sont déconvolués. 9. Procédé selon la revendication 5, caractérisé en ce qu'on effectue en outre une étape de suréchantillonnage des paquets avant le débruitage et/ou la déconvolution afin d'éviter le phénomène de distorsion de repliement. 10. Procédé selon la revendication 9, caractérisé en ce que le suréchantillonnage comporte les étapes selon lesquelles, sur chaque paquet : - on décale d'une ligne, d'une colonne, et d'une ligne et d'une colonne, les paquets au niveau immédiatement inférieur de la meilleure base ; - on décompose les quatre paquets ainsi trouvés ; - on les débruite au niveau immédiatement supérieur après décomposition, une déconvolution pouvant en outre à cette étape de débruitage ; - on recompose les paquets ainsi obtenus ; - on moyenne les divers paquets après recalage.