Brevet FR3013490A1 — publié le 20150522

Procédé et dispositif de compression d'image numérique La présente invention concerne un procédé de compression d'image numérique représentée par au moins une matrice bidimensionnelle d'échantillons d'image. Elle concerne également un dispositif de compression associé. L'invention se situe dans le domaine de la compression d'images numériques, et plus particulièrement dans le domaine de la compression d'images avec un objectif de qualité à atteindre. Les images numériques obtenues par divers capteurs sont en général représentées par des matrices d'échantillons d'image ou pixels, chaque échantillon d'image ayant une intensité de dynamique donnée. De nos jours, les images numériques ont de très grandes tailles, qu'il s'agisse des images dans le domaine de la photographie utilisant des appareils de haute qualité, et dans le domaine de la prise d'images par télédétection, ou dans le domaine de l'imagerie médicale. Afin de stocker et transmettre de telles images, il est connu d'appliquer une compression d'images numériques. Diverses techniques de compression existent. D'une manière générale, les méthodes de compression comprennent l'application de trois traitements à une image numérique source, qui sont les suivants : décorrélation, quantification, codage. Dans plusieurs standards de compression d'image numérique, la décorrélation est obtenue par application d'une transformation, appliquée aux pixels pour obtenir des coefficients transformés. La norme de codage JPEG (« Joint Photographic Expert Group ») utilise une transformation en cosinus discrète (DOT) appliquée sur des blocs d'image de taille 8x8, suivie d'une quantification par blocs et d'un codage entropique. La quantification consiste à appliquer une matrice de quantification aux coefficients d'un bloc, et il est possible de définir un niveau de qualité de compression, un facteur multiplicatif étant alors appliqué à la quantification de tous les blocs de l'image numérique source. Cependant, un tel niveau de qualité résulte en une compression plus ou moins forte selon les blocs de l'image numérique source, et ne permet pas de contrôler finement la qualité de l'image à coder. La norme de codage JPEG2000 a été mise au point après JPEG. Elle utilise une transformation en ondelettes discrètes (DWT) de l'image numérique source. La quantification et le codage sont effectués en une seule étape, grâce à un codage par couches de qualité, obtenues par un codage par plans de bits des coefficients transformés. Afin de limiter la quantité de données d'information relatives aux paramètres de codage, la norme JPEG2000 utilise des unités de codage qui sont appelées « code- blocs », typiquement de taille au moins égale à 32x32 pixels. Ainsi, JPEG2000 ne permet pas un contrôle fin de la qualité de codage. Dans le domaine du traitement des images obtenues par satellite, la norme de codage IDC (« Image Data Compression ») définie par le CCSDS (« Consultative Committee for Space Data Systems ») utilise également une transformation en ondelettes discrète, suivie d'un regroupement spatial des coefficients transformés obtenus en blocs de taille prédéterminée, typiquement de taille 8x8, et un codage par segments regroupant un nombre prédéterminé de blocs consécutifs. Chaque bloc de coefficients transformés comprend un coefficient dit basse fréquence, représentatif de la moyenne de l'intensité des pixels du bloc, et une pluralité de coefficients dit fréquentiels, correspondant à trois niveaux de décomposition en ondelettes pour un bloc de taille 8x8. Les segments sont codés indépendamment les uns des autres, les paramètres nécessaires pour la compression étant encodés en en-tête de chaque segment. Les coefficients fréquentiels des blocs d'un segment sont codés par plans de bits. Comme dans le cas de JPEG2000, l'ensemble des données binaires de codage associées à un segment, appelé également train binaire associé au segment, est organisé hiérarchiquement et peut être tronqué en tout point, pour obtenir un volume de données plus faible mais à qualité dégradée. La norme de codage CCSDS IDC propose deux modes, permettant un contrôle de la qualité et/ou du volume de données pour chaque segment. La qualité est contrôlée notamment par le nombre de plans de bits codés pour le segment. De ce fait, la qualité effective de l'image reconstruite dépend des caractéristiques de l'image numérique source. De ce fait, en pratique, il est particulièrement difficile de trouver des paramètres permettant d'assurer une qualité souhaitée pour un ensemble d'images. En effet, à paramètres constants, on peut obtenir de fortes différences en termes de rapport signal à bruit entre deux images numériques sources différentes. Il existe donc un besoin d'améliorer les méthodes de compression existantes en permettant un pilotage fin du niveau de qualité local, tout en garantissant un débit minimum pour un niveau de qualité local cible donné. A cet effet, l'invention propose, selon un premier aspect, un procédé de compression d'image numérique représentée par au moins une matrice bidimensionnelle d'échantillons d'image, chaque échantillon ayant une valeur d'intensité associée, comportant des étapes de : - application d'une transformation afin d'obtenir une représentation de l'image comportant un ensemble de blocs de coefficients, chaque bloc de coefficients comportant un coefficient dit coefficient basse fréquence correspondant à une valeur moyenne de l'intensité des échantillons d'image dudit bloc et une pluralité de coefficients dits coefficients fréquentiels ; - quantification et codage des blocs de coefficients, afin d'obtenir un ensemble de données numériques représentatives de l'image numérique compressée. Ce procédé de compression comprend, après l'obtention de blocs de coefficients, les étapes de : -pour au moins un bloc traité sélectionné parmi les blocs de coefficients : - calcul d'une première valeur représentative dudit bloc traité à partir des coefficients fréquentiels dudit bloc traité, - obtention d'une deuxième valeur représentative dudit bloc traité calculée en fonction d'un modèle de bruit d'acquisition d'image associé audit bloc traité - quantification et codage dudit bloc traité en fonction desdites première et deuxième valeurs associées audit bloc traité. Avantageusement, le procédé de compression de l'invention permet d'appliquer un codage différencié pour des blocs traités en fonction d'une première valeur calculée à partir des coefficients fréquentiels d'un bloc traité et d'une deuxième valeur représentative d'un modèle de bruit associé au bloc traité. Avantageusement, un pilotage fin de la qualité est obtenu, grâce à la sélection d'une cible de qualité par bloc. Le procédé de compression selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prises isolément ou en combinaison, suivant toutes les combinaisons techniquement possibles. Le procédé de comporte, pour un bloc traité, avant quantification et codage, une étape de détermination d'une cible de qualité en fonction desdites première et deuxième valeurs associées audit bloc traité. Selon une caractéristique, pour un bloc traité, ladite cible de qualité est une valeur représentative de l'erreur quadratique moyenne entre les coefficients du bloc traité avant quantification et les coefficients du bloc traité après quantification. Le procédé comporte en outre, pour chaque bloc traité, une étape de détermination d'un niveau de quantification permettant de respecter la cible de qualité. Selon un mode de réalisation, le procédé comprend, pour chaque bloc traité, une comparaison de ladite première valeur à un premier niveau seuil calculé à partie de ladite deuxième valeur, et lorsque ladite première valeur est inférieure audit premier niveau seuil, une étape de mise à zéro des coefficients fréquentiels dudit bloc traité. Lorsque ladite première valeur est supérieure audit premier niveau seuil et inférieure à un deuxième niveau seuil, ledit bloc traité est considéré comme étant un bloc d'exception pour lequel ladite cible de qualité est sélectionnée pour minimiser l'erreur quadratique moyenne entre les coefficients du bloc traité avant quantification et les coefficients du bloc traité après quantification. Selon une caractéristique, la première valeur représentative d'un bloc traité est représentative de la variance des coefficients fréquentiels du bloc traité, et la deuxième valeur représentative d'un bloc traité est fonction d'une variance d'un modèle de bruit fonction de la valeur du coefficient basse fréquence du bloc traité. Les blocs traités sont sélectionnés par application d'un masque binaire. Le procédé comporte, pour chaque bloc non sélectionné comme bloc traité par application dudit masque, une étape de mise à zéro des coefficients fréquentiels dudit bloc. Le procédé comporte une étape de codage d'une indication relative à un niveau de quantification à appliquer pour chaque bloc traité, une dite indication de codage étant associée à l'ensemble de données numériques représentatives de l'image numérique compressée. Selon une variante, l'étape de quantification et de codage est mise en oeuvre par un codage par plans de bits représentatifs des coefficients, et ladite indication est une indication relative au nombre de plans de bits à coder. Selon une variante, l'étape de quantification est mise en oeuvre par application d'une matrice de quantification par bloc, et ladite indication est une indication relative à un coefficient de pondération à appliquer à une matrice de quantification prédéterminée. Les blocs de coefficients sont obtenus par application d'une transformation en ondelettes discrète et regroupement spatial des coefficients de transformation obtenus. Les blocs de coefficients sont obtenus par application d'une transformation en cosinus discrète appliquée sur des blocs de taille prédéterminée. Selon un deuxième aspect, l'invention concerne un dispositif de compression d'image numérique représentée par au moins une matrice bidimensionnelle d'échantillons d'image, chaque échantillon ayant une valeur d'intensité associée, comportant des modules aptes à réaliser : -une application d'une transformation afin d'obtenir une représentation de l'image comportant un ensemble de blocs de coefficients, chaque bloc de coefficients comportant un coefficient dit coefficient basse fréquence correspondant à une valeur moyenne de l'intensité des échantillons d'image dudit bloc et une pluralité de coefficients dits coefficients fréquentiels ; - une quantification et un codage des blocs de coefficients, afin d'obtenir un ensemble de données numériques représentatives de l'image numérique compressée. Ce dispositif comprend, des modules, aptes à réaliser, après l'obtention de blocs de coefficients, pour au moins un bloc traité sélectionné parmi les blocs de coefficients : - un calcul d'une première valeur représentative dudit bloc traité à partir des coefficients fréquentiels dudit bloc traité, - une obtention d'une deuxième valeur représentative dudit bloc traité calculée en fonction d'un modèle de bruit associé audit bloc traité - une quantification et un codage dudit bloc traité en fonction desdites première et deuxième valeurs associées audit bloc traité. Le dispositif de compression comprend des modules aptes à mettre en oeuvre les étapes du procédé de compression d'images numériques brièvement décrit ci-dessus. Ce dispositif de compression apporte des avantages analogues aux avantages du procédé de compression brièvement décrit ci-dessus. D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : -la figure 1 est un synoptique d'un système de compression d'images mettant en oeuvre un procédé de compression selon l'invention ; - la figure 2 représente schématiquement un ensemble de sous-bandes obtenues par application d'une transformation en ondelettes discrète à une image numérique ; - la figure 3 représente schématiquement un bloc de coefficients transformés obtenus par une transformation en ondelettes ; -la figure 4 est un schéma représentant les blocs fonctionnels d'un dispositif programmable apte à mettre en oeuvre l'invention ; - la figure 5 est un synoptique d'une étape de détermination d'une cible de qualité par bloc selon un mode de réalisation de l'invention ; - la figure 6 est un tableau illustrant l'association de cibles de qualité et une valeur représentative des coefficients fréquentiels de blocs à coder ; - la figure 7 illustre schématiquement un ensemble de plans de bits d'un bloc ; - la figure 8 est un synoptique d'une étape de quantification et codage utilisant une cible de qualité selon un premier mode de réalisation de l'invention ; - la figure 9 est un synoptique d'un procédé de compression mettant en oeuvre la détermination d'une cible de qualité selon un deuxième mode de réalisation de l'invention. La figure 1 illustre un système 1 de compression/décompression d'images. Le système est composé d'une première partie 2 effectuant la compression d'images numériques, d'une deuxième partie de stockage/transmission 3 et d'une troisième partie de décompression 4. Les parties 2 et 4 sont respectivement mises en oeuvre par des dispositifs tels que brièvement décrits ci-après en référence à la figure 4, comportant une ou plusieurs unités de traitement, par exemple des processeurs, et aptes à effectuer des calculs. La deuxième partie de stockage et transmission 3 est réalisée par un stockage de données dans une mémoire d'un appareil programmable ou d'un dispositif électronique et/ou par un réseau de communications apte à effectuer une transmission de données numériques. En variante, les procédés de l'invention sont mis en oeuvre par des dispositifs électroniques de type circuits logiques programmables, tels des cartes électroniques à base de FPGA ou ASIC. Une image numérique à compresser I est présentée en entrée de la première partie 2 de compression. L'image numérique I est composée d'une ou plusieurs matrices d'échantillons d'image ou pixels, chaque échantillon d'image ayant une valeur d'intensité associée. Par la suite on considère pour simplifier l'explication qu'une image numérique à compresser est composée d'une matrice de dimensions MxN de pixels, chaque pixel ayant une intensité initialement codée sur P bits. Par exemple, P=8. L'image numérique I est fournie à un premier module 6 de transformation de l'image numérique à compresser en blocs de coefficients transformés, de taille prédéterminée. Les blocs de coefficients sont agencés spatialement. Selon un mode de réalisation, une transformation en ondelettes discrète (DWT) est appliquée sur l'ensemble de la matrice d'image, sur 3 niveaux de décomposition, comme illustré schématiquement à la figure 2. De manière connue, l'application d'une DWT s'effectue d'abord sur la matrice d'image complète, permettant d'obtenir quatre sous-bandes, chaque sous-bande comportant des coefficients transformés correspondant à l'application d'un filtrage passe-bas L et d'un filtrage passe-haut H, d'abord selon une première direction de la matrice, puis selon une deuxième direction, perpendiculaire à la première direction. L'application d'un premier niveau de décomposition à la matrice I permet d'obtenir les sous-bandes de coefficients ni, LH1, HL, et HH,. La sous-bande LL, comprend une version passe-bas de l'image initiale I à un niveau de résolution inférieur au niveau de résolution de l'image initiale I. Si l'image initiale est composée de M lignes et N colonnes, chaque sous-bande ni, LH1, HL1, HH, de premier niveau de décomposition comporte M/2 lignes et N/2 colonnes. L'application des filtres respectifs L et H est itérée sur la sous-bande LI-1, permettant d'obtenir les sous-bandes LL2, LH2, HL2 et HH2, de deuxième niveau de décomposition et de taille M/4 lignes et N/4 colonnes. La sous-bande LL2 est à nouveau traitée par filtrage pour obtenir les sous-bandes de troisième niveau de décomposition LL3, LH3, HL3 et HH3. Ainsi, après application d'une transformation en ondelettes sur trois niveaux de décomposition, comme illustré à la figure 2, l'image numérique est représentée par un ensemble de sous-bandes de coefficients transformés. A un niveau de décomposition p, la taille des sous-bandes est de MI 2P x NI2P . Les coefficients transformés des sous-bandes sont regroupés en blocs de coefficients localisées spatialement de taille 2P x2 pour p niveaux de décomposition, donc des blocs de taille 8x8 pour une décomposition à 3 niveaux de décomposition. Ainsi, comme illustré aux figures 2 et 3, un bloc Bk de coefficients transformés est obtenu par regroupement spatial de coefficients de la transformation DWT correspondant à une même position spatiale de l'image initiale I. Chaque bloc de coefficients Bk comprend un coefficient basse fréquence DC, issue de la sous-bande LL3 dans l'exemple, et plusieurs coefficients fréquentiels AC issus des sous-bandes LH3, HL3, HH3, LH2, HI-2, HH2, LH1, HL1, HH, comme illustré. En variante, selon un mode de réalisation alternatif, l'image I est divisée en blocs de pixels de taille 8x8 et une transformation en cosinus discrète (DCT) est appliquée sur chaque bloc de pixels, permettant d'obtenir des blocs de taille 8x8 localisés spatialement de coefficients transformés. De manière générale, toute transformation de décorrélation, permettant d'obtenir des blocs de coefficients de taille mxn est applicable. De retour à la figure 1, suite à la transformation, une représentation de l'image par blocs Bk de coefficients transformés incluant un coefficient basse fréquence C0,0, également appelé coefficient DC, et des coefficients fréquentiels est obtenue. Les blocs de coefficients sont ensuite traités par un module 8 de détermination de cible de qualité par bloc de l'image. Ce module 8 de détermination de cible de qualité par bloc est spécifique à l'invention. Il prend en entrée des paramètres P1 relatifs à un modèle de bruit associé à l'image traitée et des paramètres P2 relatifs à la cible de qualité. Optionnellement, un ou plusieurs masques binaires MK permettant de distinguer entre blocs d'intérêt et blocs de non-intérêt est également fourni. Par exemple, si l'image numérique I à compresser est une image de télédétection, acquise par un capteur embarqué à bord d'un satellite, un masque distinguant des blocs appartenant à des zones de nuages est fourni, les blocs de nuages étant des blocs sans intérêt pour la mission assignée au satellite car sans contenu d'information. D'une manière plus générale, pour une application cible donnée, dans un contexte d'utilisation donné, il est possible de fournir de tels masques MK permettant de définir l'intérêt de chaque bloc de manière à pouvoir adapter la qualité image cible associée. Par la suite, seuls les blocs d'intérêt sont traités pour établir une cible de qualité associée. Les blocs de non-intérêt sont simplement codés par leur coefficient basse fréquence. La détermination d'une cible de qualité par bloc d'intérêt, appelé également bloc traité, sera détaillée ci-après en référence à la figure 5. Dans un mode de réalisation de l'invention, le module 8 de détermination d'une cible de qualité par bloc traité fournit, pour chaque bloc traité Bk, une valeur cible d'erreur quadratique moyenne MSEk (Mean Square Error) de codage pour le bloc Bk. Chaque coefficient du bloc Bk est codé par un coefficient quantifié C, * L'erreur quadratique moyenne de codage est donnée par : MSEk =1 - ) 2 (Eq 1) mxn (i,.1)#(0,0) pour un bloc Bk de taille mxn. Ces valeurs MSEk ainsi que les blocs de coefficients transformés sont fournis en entrée d'un module 10 de quantification/codage. Pour chaque bloc traité Bk, le module 10 de quantification/codage détermine statistiquement ou par mesure de l'erreur quadratique moyenne (MSE), un niveau de quantification associé et réalise le codage entropique en appliquant le niveau de quantification déterminé. Selon un premier mode de réalisation, le module 10 réalise la quantification et le codage par un codage des coefficients d'un bloc traité par plans de bits. Selon un deuxième mode de réalisation, la quantification est réalisée par application d'une matrice de quantification par bloc, et le codage est un codage entropique des coefficients quantifiés. En sortie du module 10, un ensemble de données binaires ou train binaire représentatif des coefficients quantifiés codés 12 ainsi qu'un ensemble d'indications 14 relatives aux paramètres de quantification/codage sont obtenus. Dans un mode de réalisation, des groupes de blocs ou segments d'images sont codés indépendamment. Un ensemble d'informations 14 et un ensemble de données binaires 12 sont associés à chaque segment pour former un train binaire de données compressées associées au segment. L'ensemble de ces trains binaires par segment forme un ensemble de données numériques compressées F représentatives de l'image I. En variante, un seul ensemble d'information 14 est associé à l'ensemble des blocs de l'image I. L'ensemble de données numériques compressées F est soit stocké, par exemple sous forme d'un fichier, soit transmis par l'intermédiaire d'un réseau de communications, par une partie de stockage/transmission 3 du système 1. L'ensemble de données numériques compressées est reçu par la troisième partie de décompression 4. Un module de décodage et déquantification 16 reçoit, pour chaque segment de l'image, les données numériques compressées F comportant les données binaires 12 et les ensembles d'information 14 associées, et décode des coefficients transformés associés aux blocs Bk de l'image I. Des coefficients quantifiés C,: pour chacun des blocs Bk de l'image représentée par les données numériques compressées F sont obtenus en sortie du module 16. Ensuite une transformation inverse par blocs de coefficients est appliquée par le module 18. Si la transformation appliquée par le module de transformation 6 en blocs de coefficients est une DWT comme expliqué ci-dessus, les coefficients des blocs Bk sont redistribués par sous-bandes, et une transformation DWT inverse est appliquée par le module 18 d'application d'une transformation inverse, en fonction du nombre de niveaux de décomposition. Si la transformation appliquée par le module de transformation 6 est une DOT par blocs, une DOT inverse par blocs est appliquée par le module 18 d'application d'une transformation inverse. En sortie du module 18, une image reconstruite I* est obtenue, les valeurs d'intensité des échantillons de l'image I* étant obtenues par application de la transformation inverse comme expliqué ci-dessus. Les procédés de l'invention sont mis en oeuvre par un dispositif électronique, par exemple un dispositif programmable de type ordinateur, station de travail, dont les principaux blocs fonctionnels sont illustrés à la figure 4. Un dispositif électronique 20 apte à mettre en oeuvre les procédés de l'invention comprend une unité centrale de traitement 22, par exemple un CPU, apte à exécuter des opérations préprogrammées ou des instructions de programme informatique lorsque le dispositif 20 est mis sous tension. Dans un mode de réalisation, une unité centrale de traitement multi-processeurs est utilisée, permettant d'effectuer des calculs parallèles. Le dispositif 20 comporte également des moyens de stockage d'informations 24, par exemple des registres, aptes à stocker des instructions de code exécutable permettant la mise en oeuvre de programmes comportant des instructions de code aptes à mettre en oeuvre les procédés selon l'invention. Le dispositif 20 comporte des moyens de commande 26 permettant de mettre à jour des paramètres et de recevoir des commandes d'un opérateur. Lorsque le dispositif programmable 20 est un dispositif embarqué, les moyens de commande 26 comprennent un dispositif de télécommunication permettant de recevoir des commandes et des valeurs de paramètres à distance. Alternativement et de manière optionnelle, les moyens de commande 26 sont des moyens de saisie de commandes d'un opérateur, par exemple un clavier. De manière optionnelle, le dispositif programmable 20 comprend un écran 28 et un moyen supplémentaire de pointage 30, tel une souris. Les divers blocs fonctionnels du dispositif 20 décrits ci-dessus sont connectés via un bus de communication 32. En variante, les procédés de l'invention sont mis en oeuvre par des dispositifs électroniques de type circuits logiques programmables, tels des cartes électroniques à base de FPGA ou ASIC. La figure 5 illustre un mode de réalisation d'une étape de détermination de cible de qualité par blocs de coefficients transformés d'une image numérique I, mise en oeuvre par des modules du processeur 22 d'un dispositif 20 apte à mettre en oeuvre l'invention ou par une carte électronique développée pour la mise en oeuvre de l'invention. La détermination d'une cible de qualité est mise en oeuvre pour chaque bloc de coefficients transformés Bk issus de l'image I. Un bloc Bk est fourni en entrée d'une étape optionnelle 34 de sélection de blocs d'intérêt. Si un masque MK binaire définissant des blocs d'intérêt et des blocs de non- intérêt est fourni, l'étape 34 consiste à appliquer le masque MK pour déterminer si le bloc Bk est un bloc d'intérêt. En cas de réponse négative, l'étape 34 est suivie d'une étape 35 de détermination d'au moins un paramètre de codage pour le bloc Bk. Lors de l'étape 35, le bloc Bk est modifié en un bloc ak à coder dans lequel tous les coefficients fréquentiels AC sont mis à zéro. Avantageusement, le débit de codage associé aux blocs de non-intérêt est ainsi minimisé, une information minimale relative à ces blocs étant néanmoins codée par codage du coefficient basse-fréquence du bloc. En variante, une erreur quadratique moyenne MSEk cible est associée au bloc Bk à l'étape 35. Si le bloc Bk est un bloc d'intérêt, l'étape 34 est suivie d'une étape 36 de calcul d'une première valeur représentative des coefficients fréquentiels du bloc Bk. Dans le mode de réalisation préféré, la variance des coefficients fréquentiels du bloc est calculée : VI (Bk) = 1 (mxn -1) ,,kEB, Le calcul de la première valeur prend en compte tous les coefficients C, du bloc Bk sauf le coefficient basse fréquence C0,0, C étant la valeur moyenne des coefficients fréquentiels du bloc Bk. Alternativement, l'écart-type des coefficients fréquentiels du bloc Bk ou une autre valeur représentative de la statistique des coefficients, comme un moment d'ordre supérieur, est calculée. Lors d'une étape 38, une deuxième valeur V2, représentative d'un modèle de bruit d'acquisition d'image associé au bloc Bk, est calculée. Dans le mode de réalisation préféré, un modèle de bruit paramétré est utilisé. La deuxième valeur V2 est fonction de ce modèle de bruit et de la valeur du coefficient basse fréquence C0,0, représentatif de la moyenne de l'intensité des pixels du bloc Bk, selon la formule suivante : V2 (Bk ) = a2 + fi - C0,0 La deuxième valeur V2 est la variance du bruit local, et les paramètres a et 13 sont fournis dans l'ensemble des paramètres P1 en entrée du procédé. Ces paramètres a et 13 sont des constantes prédéterminées. Par exemple, les constantes a et 13 sont représentatives du capteur d'acquisition de l'image I à compresser et sont déterminées par étalonnage du capteur et mémorisées. Pour le satellite PLEIADES-1A, les valeurs a et 13 sont égales respectivement à 2.267 et 0.0393. Alternativement, pour un bruit de détection non « blanc », le calcul du bruit peut prendre en compte des modèles différents pour chaque sous-bande ou coefficient de transformé. Le bruit est donc obtenu avec la formule suivante : V2 (Bk ) at,j 2 ± )6, - Co 0 1,JEB (1,.1)(0,0) Il est à noter que les étapes de calcul 36, 38 peuvent être effectuées en parallèle ou leur ordre peut être inversé. Selon l'invention, la quantification et le codage du bloc traité sont effectuées en fonction des première valeur V1 et deuxième valeur V2 calculées pour le bloc. Lors d'une étape de détermination du type de bloc 40, la première valeur Vl(Bk) est comparée à la deuxième valeur calculée, pondérée par un coefficient K. Le coefficient K est une constante. Ce coefficient est un paramètre fourni par l'utilisateur, et permet de déterminer un niveau seuil permettant de départager les blocs contenant de l'information utile des blocs de bruit. Dans un mode de réalisation, K est choisi de préférence entre 0,5 et 4, par exemple K=1. Si la première valeur Vl(Bk) est inférieure au premier niveau seuil K - V2 (Bk ) , il en est conclu que le bloc traité Bk est un bloc contenant du bruit. Dans ce cas, l'étape 40 est suivie d'une étape 42, analogue à l'étape 35 décrite ci-dessus, dans laquelle le bloc Bk est modifié en un bloc B' k à coder dans lequel tous les coefficients fréquentiels AC sont mis à zéro. En effet, il est inutile d'allouer un débit de codage important à un bloc qui ne comporte que du bruit. Si la première valeur Vl(Bk) est supérieure au premier niveau de seuil K - V 2(B k) , il en est conclu que le bloc Bk contient de l'information utile. L'étape 40 est suivie d'une étape 44 de comparaison de la première valeur Vl(Bk) à un deuxième niveau de seuil prédéterminé SE. Si la valeur Vl(Bk) est inférieure à ce deuxième niveau de seuil, il en est conclu que le niveau de signal du bloc Bk, estimé par la valeur Vl(Bk) est faible, légèrement plus élevé qu'un niveau de signal de type bruit. Par exemple, SE est égal compris entre 16 et 32. Ces blocs sont des blocs dits d'exception, comportant une information tout en ayant un faible niveau de signal. Ce type de blocs est fortement compressé par les méthodes de compression classiques, même pour des faibles taux de compression. Pour certaines applications, notamment dans le domaine de la télédétection, il a été noté qu'il est important de préserver de tels blocs. De ce fait un traitement particulier leur est appliqué. Dans le mode de réalisation préféré, la cible de qualité pour ces blocs est mise à un niveau élevé, ce qui est réalisé par exemple par la mise à 0, lors de l'étape 46, de l'erreur quadratique moyenne cible MSEk pour le bloc Bk. Le modèle de bruit présenté ci-dessus, peut être utilisé, lors de l'étape optionnelle 47, pour réaliser un débruitage par seuillage ou pondération des coefficients pour ce type de bloc, afin de garantir un codage plus efficace de l'information utile. Dans ce cas, la MSEk = 0 doit être garantie entre le bloc débruité obtenu en sortie de l'étape de débruitage 47 et le bloc après codage. Si la valeur V, (Bk) est supérieure au seuil SE, le bloc Bk est considéré comme bloc régulier, pour lequel une cible de qualité est fixée lors de l'étape 48. Lors de cette étape 48, une erreur quadratique moyenne de codage cible MSEk est fixée pour le bloc Bk. De préférence, la valeur MSEk est fixée en fonction de la première valeur Vl(Bk) égale à la variance des coefficients fréquentiels dans ce mode de réalisation. Des valeurs d'erreur quadratique moyenne cible MSE cible en fonction des valeurs de variance ou d'écart-type sont tabulées et mémorisées afin d'être fournies parmi les paramètres P2 en entrée du procédé. Par exemple, ces valeurs sont stockées dans des tables, comme par exemple la table montrée en figure 6 reliant des valeurs de variance de coefficients fréquentiels de bloc à des valeurs de cible de qualité MSE cible. Ainsi, par exemple, en utilisant la table T de la figure 6, si le bloc Bk a une première valeur V1 associée strictement inférieure à 64, l'erreur quadratique moyenne cible est MSEk=2. En variante, plusieurs tables reliant la variance des coefficients fréquentiels à une valeur cible de qualité sont mémorisées et utilisées par exemple en fonction de l'appartenance du bloc traité à un type de région. Pour cela, un classement préalable en régions doit être effectué. Une fois la détermination de cible de qualité effectuée pour un bloc Bk, elle est mémorisée afin d'être fournie au module 10 de quantification et de codage, pour un codage du bloc respectant cette cible de qualité. En variante, le modèle de bruit présenté ci-dessus, peut également être utilisé, lors d'une étape optionnelle 49, pour réaliser un débruitage par seuillage ou pondération des coefficients du bloc Bk, afin de garantir un codage plus efficace de l'information utile. Dans ce cas, la valeur d'erreur quadratique moyenne cible MSEk associée au bloc Bk doit être garantie entre le bloc débruité obtenu après l'étape de débruitage 49 et le bloc après codage. Dans le cas particulier où les coefficients fréquentiels AC sont mis à 0, le bloc modifié B'k est directement fourni pour quantification et codage, sans indication de cible de qualité particulière associée. Les figures 7 à 8 illustrent la mise en oeuvre des cibles de qualité et paramètres de codage déterminés dans un premier mode de réalisation, dans lequel les blocs de coefficients sont codés par plans de bits. Ce mode de réalisation est en particulier compatible avec le standard de compression CCSDS IDC. La figure 7 illustre schématiquement le principe d'un codage par plan de bits. Pour un bloc Bk donné, un bloc 4x4 étant illustré à la figure 7, chaque coefficient C, a une valeur associée représentée sur un nombre B de bits en représentation binaire, rangés dans l'ordre entre le bit 0P-1 de poids fort ou MSB pour « most significant bit » et le bit 0,° de poids faible, ou LSB pour « least significant bit ». Ainsi, il est possible de représenter un bloc par les plans de bits de ses coefficients, un plan de bits regroupant l'ensemble des bits de même poids de tous les coefficients d'un bloc. Le nombre maximal de plans de bits nécessaire pour représenter tous les coefficients fréquentiels d'un bloc, c'est-à-dire le nombre de bits nécessaire pour coder le plus grand coefficient fréquentiel d'un bloc de coefficients, est appelé profondeur ou « bitdepth ». Dans l'exemple de la figure 7, le bloc Bk est représenté par B=5 plans de bits, notés BP0 à BP4, donc sa profondeur est égale à B=5. Le codage par plans de bits consiste à coder successivement les plans de bits représentatifs des coefficients d'un bloc ou de plusieurs blocs consécutifs d'une image, en partant du plan MSB, comme indiqué par le sens de la flèche de la figure 7, ce qui permet d'obtenir un ensemble binaire de données ou train binaire hiérarchiquement emboîté, l'information la plus importante se situant en tête du train binaire. La quantification est réalisée par l'arrêt du codage à un certain nombre de plans de bits pour un bloc Bk. Par exemple, le codage des plans de bits BP4 et BP3 peut s'avérer suffisant pour obtenir une valeur d'erreur quadratique cible. Selon le standard de compression CCSDS IDC, des segments regroupant plusieurs blocs sont codés indépendamment. Pour chaque segment, les indications de codage sont codées dans l'en-tête du segment. Le coefficient basse-fréquence ou coefficient DC de chaque bloc est codé par codage prédictif, séparément des coefficients fréquentiels ou coefficients AC des blocs du segment. De préférence, un codage de Rice est appliqué aux différences entre coefficients DC quantifiés successifs des blocs de coefficients du segment. Ensuite, les profondeurs BitDepthAC Block, associées au codage des coefficients fréquentiels des blocs B, successifs de coefficients du segment sont codées par un codage entropique, par exemple le codage de Rice. Ensuite, pour chaque bloc, les plans de bits de tous les coefficients fréquentiels AC successifs, en partant du plan de bits de poids fort, sont codés, par étapes de codage ou « stages » en anglais, les étapes de codage permettant une hiérarchisation des données codées par plan de bits. Un nombre maximal de 4 étapes de codage pour chaque plan de bits est prévu par le standard CCSDS IDC. Ainsi, un train binaire à couches de qualité est obtenu pour un segment de blocs de coefficients. Selon l'invention, le codage pour chaque bloc de coefficients d'un segment prend en compte la cible de qualité ou la pré-quantification par mise à 0 des coefficients fréquentiels mises en oeuvre dans l'étape de détermination de cible de qualité par blocs de coefficients décrite ci-dessus. La figure 8 illustre les principales étapes d'un procédé de quantification et codage des coefficients fréquentiels prenant en compte les cibles de qualité déterminées et mémorisées par bloc, selon un premier mode de réalisation de l'invention. Pour chaque bloc traité, un niveau de quantification par bloc traité est déterminé en fonction de la cible de qualité préalablement calculée pour le bloc. Pour un bloc Bk d'un segment d'image à compresser donné, une première étape 50 d'initialisation des paramètres de codage est mise en oeuvre. Lors de cette étape d'initialisation 50, la profondeur B du bloc Bk à coder est calculée, un indice b du plan de bit courant BPb à coder est mis à B-1, et un indice s représentatif de l'étape de codage du plan de bit b est mis à 0. Ensuite, lors d'une étape de quantification 52, les coefficients fréquentiels du bloc Bk, quantifiés selon les paramètres b et s courants, sont calculés par la fonction notée quant(Bk,b,$). L'erreur quadratique moyenne de quantification courante, notée mse, est également calculée lors de l'étape de quantification 52. Lors d'une étape de comparaison 54, l'erreur quadratique moyenne mse du bloc calculée est comparée à l'erreur quadratique moyenne cible MSEk pour le bloc Bk, préalablement déterminée et mémorisée. Si l'erreur quadratique moyenne de quantification mse courante, avec les paramètres de codage courant, est inférieure à l'erreur quadratique moyenne cible MSEk, l'étape 54 est suivie d'une étape de mémorisation d'une indication de codage d'indice b de plan de bits à coder et d'étape de codage s associée pour le bloc courant Bk. L'indication de codage par bloc est associée aux données binaires codées et transmises avec les données compressées pour utilisation ultérieure lors du décodage, comme expliqué ci-dessus. Si l'erreur quadratique moyenne de quantification courante, avec les paramètres de codage courant, est supérieure à l'erreur quadratique moyenne cible, l'étape 54 est suivie d'une étape 58 de comparaison de l'indice de l'étape de codage à 3. Si la valeur de l'indice s a atteint la valeur 3, alors l'étape 58 est suivie d'une étape 60 de mise à jour des paramètres. L'indice de plan de bits à coder est décrémenté de 1 et l'indice s est de nouveau mis à 0. Sinon, l'indice s est incrémenté de 1 lors de l'étape 62 d'incrémentation d'étape de codage. Ainsi, quatre étapes de codage par plan de bits sont autorisées, ces étapes de codage étant appliquées conformément au standard de compression CCSDS IDC. Les étapes respectives 60 et 62 sont suivies de l'étape de quantification 52 précédemment décrite. Les étapes 50 à 62 sont appliquées pour chaque bloc Bk à coder de chaque segment. Pour les blocs pour lesquels les coefficients fréquentiels AC ont tous été mis à 0, seul le coefficient basse-fréquence est quantifié et codé, le niveau de quantification étant déterminé de manière classique pour tout le segment. Dans un mode de réalisation, un fichier réunissant l'ensemble d'indications de codage pour l'ensemble des blocs d'un segment est créé pour chaque segment codé et stocké/transmis avec les données codées du segment, par exemple de la même manière que la séquence BitDepthAC Block,. En variante, deux simplifications peuvent être introduites dans l'étape de quantification/codage visant à limiter la complexité calculatoire. Selon une première variante, le calcul de l'erreur quadratique moyenne mse est réalisé une seule fois par plan de bits, pour une même étape « s» de codage. De cette manière les étapes 58 et 62 peuvent être supprimées. Selon une deuxième variante, nécessitant encore moins de ressources calculatoires, les valeurs de « b» et « s» sont associées à des niveaux de qualité image cible en fonction de la variance du bloc. Ces valeurs sont déterminées statistiquement pour une application ou capteur donnés. Ainsi, les « b» et « s » nécessaires à atteindre (statistiquement) la qualité image cible pour le bloc peuvent être appliqués directement dès l'étape 50. Les étapes 52, 54, 58, 60 et 62 deviennent alors inutiles. La figure 9 illustre un deuxième mode de réalisation de l'invention, compatible avec le standard de compression JPEG. L'image numérique à compresser I est fournie en entrée. Une première étape 80 d'obtention de blocs 8x8 de l'image numérique I et d'application d'une DCT par blocs 8x8 est mise en oeuvre. Des blocs de coefficients transformés sont obtenus. Ensuite, une étape de détermination 82 de cible de qualité, selon l'invention, pour l'ensemble des blocs de coefficients est appliquée. Les cibles de qualité par blocs sont mémorisées. Ensuite, pour chaque bloc Bk de coefficients de l'image, une quantification des coefficients par application d'une matrice de quantification yQ(u,v) est appliquée à l'étape de quantification 86. La matrice de quantification Q(u,v) est une matrice prédéterminée mémorisée, et y est un coefficient de pondération appliqué à tous les coefficients de la matrice Q(u,v). Dans un mode de réalisation tous les coefficients de la matrice Q(u,v) sont multipliés par y. Lors de l'initialisation, y=1. Un bloc de coefficients quantifiés, noté qBk, est obtenu. L'erreur quadratique moyenne mse entre le bloc de coefficients quantifié qBk et le bloc Bk est calculée à l'étape de calcul d'erreur de quantification 88. L'erreur quadratique moyenne mse calculée est comparée à l'erreur quadratique moyenne cible lors de l'étape 90. Si l'erreur quadratique moyenne mse est inférieure à la cible de qualité MSEk, la valeur du coefficient de pondération courante y est mémorisée comme indication de codage lors d'une étape 92. Si l'erreur quadratique moyenne mse est supérieure à la cible de qualité MSEk pour le bloc Bk courant, le coefficient de pondération y est diminué lors d'une étape 94 de mise à jour du coefficient de pondération. Par exemple, y est divisé par 2. En variante, une nouvelle valeur de y est choisie dans un ensemble de valeurs prédéterminées mémorisé. Les étapes 86 à 94 sont mises en oeuvre pour chaque bloc, et permettent de déterminer un niveau de quantification en fonction de la cible de qualité MSEk à respecter. L'indication de codage pour chaque bloc est mémorisé et indiquée en association avec les données compressées, par exemple dans l'en-tête d'un fichier de données compressées. Les étapes décrites ci-dessus sont suivies d'une étape de codage d'évènement plage-longueur classique appliquée sur les coefficients fréquentiels de chaque bloc quantifié. En variante, plusieurs matrices de quantification Q*(u,v) sont préalablement mémorisées, et une matrice parmi les matrices mémorisées est choisie pour la quantification. Une indication relative à la matrice choisie parmi l'ensemble des matrices mémorisées est alors codée. Selon une autre variante, permettant de simplifier les calculs, les valeurs de y ou les matrices à appliquer peuvent également être déterminées statistiquement en fonction de la variance du bloc et de la qualité image cible. Elles peuvent donc être appliquées dès l'étape 82. Les étapes 88, 90 et 94 deviennent alors inutiles Le procédé de l'invention a été décrit ci-dessus en particulier dans son adaptation avec le standard de compression CCSDS IDC et avec le standard de compression JPEG. Il est cependant clair que le procédé selon l'invention s'applique en général avec tout autre procédé de compression, notamment avec le standard de compression JPEG2000. En particulier, une taille quelconque des blocs traités peut être mises en oeuvre avec le procédé de l'invention. REVENDICATIONS1.- Procédé de compression d'image numérique représentée par au moins une matrice bidimensionnelle d'échantillons d'image, chaque échantillon ayant une valeur d'intensité associée, comportant des étapes de : - application d'une transformation afin d'obtenir une représentation de l'image comportant un ensemble de blocs de coefficients, chaque bloc de coefficients comportant un coefficient dit coefficient basse fréquence correspondant à une valeur moyenne de l'intensité des échantillons d'image dudit bloc et une pluralité de coefficients dits coefficients fréquentiels ; - quantification et codage des blocs de coefficients, afin d'obtenir un ensemble de données numériques représentatives de l'image numérique compressée, caractérisé en ce qu'il comprend, après l'obtention de blocs de coefficients, les étapes de : -pour au moins un bloc traité sélectionné parmi les blocs de coefficients : - calcul (36) d'une première valeur représentative dudit bloc traité à partir des coefficients fréquentiels dudit bloc traité, - obtention (38) d'une deuxième valeur représentative dudit bloc traité calculée en fonction d'un modèle de bruit d'acquisition d'image associé audit bloc traité, - quantification et codage dudit bloc traité en fonction desdites première et deuxième valeurs associées audit bloc traité. 2.- Procédé de compression selon la revendication 1, caractérisé en ce qu'il comporte, pour un bloc traité, avant quantification et codage, une étape de détermination (46, 48) d'une cible de qualité en fonction desdites première et deuxième valeurs associées audit bloc traité. 3.- Procédé de compression selon la revendication 2, caractérisé en ce que, pour un bloc traité, ladite cible de qualité est une valeur représentative de l'erreur quadratique moyenne entre les coefficients du bloc traité avant quantification et les coefficients du bloc traité après quantification. 4.- Procédé de compression selon l'une quelconque des revendications 2 à 3, caractérisé en ce qu'il comporte en outre, pour chaque bloc traité, une étape dedétermination (56) d'un niveau de quantification permettant de respecter la cible de qualité. 5.- Procédé de compression selon l'une quelconques des revendications 2 à 4, caractérisé en ce qu'il comprend, pour chaque bloc traité, une comparaison (40) de ladite première valeur à un premier niveau seuil calculé à partir de ladite deuxième valeur, et lorsque ladite première valeur est inférieure audit premier niveau seuil, une étape de mise à zéro (46) des coefficients fréquentiels dudit bloc traité. 6.- Procédé de compression selon la revendication 5, caractérisé en ce que, lorsque ladite première valeur est supérieure audit premier niveau seuil et inférieure à un deuxième niveau seuil, ledit bloc traité est considéré comme étant un bloc d'exception pour lequel ladite cible de qualité est sélectionnée pour minimiser l'erreur quadratique moyenne entre les coefficients du bloc traité avant quantification et les coefficients du bloc traité après quantification. 7.- Procédé de compression selon l'une quelconque des revendications 1 à 6, caractérisé en ce que ladite première valeur représentative d'un bloc traité est représentative de la variance des coefficients fréquentiels du bloc traité, et en ce que ladite deuxième valeur représentative d'un bloc traité est fonction d'une variance d'un modèle de bruit fonction de la valeur du coefficient basse fréquence du bloc traité. 8.- Procédé de compression selon l'une quelconque des revendications 1 à 7, caractérisé en ce que lesdits blocs traités sont sélectionnés par application (34) d'un masque binaire. 9.- Procédé de compression selon la revendication 8, caractérisé en ce qu'il comporte, pour chaque bloc non sélectionné comme bloc traité par application dudit masque, une étape de mise à zéro (35) des coefficients fréquentiels dudit bloc. 10.- Procédé de compression selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape de codage (56, 92) d'une indication relative à un niveau de quantification à appliquer pour chaque bloc traité, une dite indication de codage étant associée à l'ensemble de données numériques représentatives de l'image numérique compressée. 11.- Procédé de compression selon la revendication 10, dans lequel l'étape de quantification et de codage est mise en oeuvre par un codage par plans de bits représentatifs des coefficients, caractérisé en ce que ladite indication est une indication relative au nombre de plans de bits à coder. 12.- Procédé de compression selon l'une quelconque des revendications précédentes, caractérisé en ce que les blocs de coefficients sont obtenus par application d'une transformation en ondelettes discrète et regroupement spatial des coefficients de transformation obtenus. 13.- Procédé de compression selon la revendication 10, dans lequel l'étape de quantification est mise en oeuvre par application d'une matrice de quantification par bloc, caractérisé en ce que ladite indication est une indication relative à un coefficient de pondération à appliquer à une matrice de quantification prédéterminée. 14.- Procédé de compression selon à la revendication 13, caractérisé en ce que les blocs de coefficients sont obtenus par application d'une transformation en cosinus discrète appliquée sur des blocs de taille prédéterminée. 20 15.- Dispositif de compression d'image numérique représentée par au moins une matrice bidimensionnelle d'échantillons d'image, chaque échantillon ayant une valeur d'intensité associée, comportant : -un module (6) d'application d'une transformation afin d'obtenir une représentation de l'image comportant un ensemble de blocs de coefficients, chaque bloc de coefficients 25 comportant un coefficient dit coefficient basse fréquence correspondant à une valeur moyenne de l'intensité des échantillons d'image dudit bloc et une pluralité de coefficients dits coefficients fréquentiels ; - un module (10) de quantification et de codage des blocs de coefficients, afin d'obtenir un ensemble de données numériques représentatives de l'image numérique 30 compressée, caractérisé en ce qu'il comprend, des modules, aptes à réaliser, après l'obtention de blocs de coefficients, pour au moins un bloc traité sélectionné parmi les blocs de coefficients : - un calcul (8) d'une première valeur représentative dudit bloc traité à partir 35 des coefficients fréquentiels dudit bloc traité, 10 15- une obtention (8) d'une deuxième valeur représentative dudit bloc traité calculée en fonction d'un modèle de bruit associé audit bloc traité - une quantification et un codage (10) dudit bloc traité en fonction desdites première et deuxième valeurs associées audit bloc traité.5