FR 2944935 A1 20101029 FR 0902068 A 20090428 PROCEDE D'ESTIMATION DU DEBIT ET DE LA DISTORSION DE DONNEES D'IMAGES CODEES A POSTERIORI DE L'ENCODAGE La présente invention a pour objet un procédé d'estimation du débit et de la distorsion d'un flux de données d'images compressées, ladite estimation étant effectuée a posteriori de l'étape de compression. Le champ d'application de l'invention est celui de la transmission d'images à travers un canal de communication à débit contraint et soumis à des erreurs, par exemple un canal sans fil de type hertzien. Le débit limité disponible pour la transmission nécessite une compression du flux de données d'images, à cet effet un des standards de compression d'image couramment utilisé et pour lequel l'invention s'applique avantageusement est le standard JPEG 2000 (Joint Photographic Expert Group) développé par l'Organisation Internationale de normalisation connue sous l'abréviation anglo-saxonne ISO. L'invention s'applique également dans le cadre de réseaux de communications hétérogènes comportant de nombreux noeuds de communication intermédiaires pour lesquels il est important de pouvoir déterminer les caractéristiques du flux d'images transmis sans avoir à effectuer une étape de décompression trop pénalisante en terme de complexité de calculs et de délai introduit dans la transmission. La connaissance des caractéristiques telles que le débit ou la distorsion d'un flux d'images compressées sont utiles lorsqu'on cherche, par exemple, à les modifier au cours d'une étape de transcodage. Les utilisations possibles de l'invention concernent tous les systèmes mettant en oeuvre des transmissions de flux d'images compressées à travers des canaux de propagation sujets à erreurs, par exemple les systèmes d'imagerie à la demande, de diffusion, de réseaux de capteurs sans fils ou encore de vidéo surveillance. Dans un contexte de transmission de données d'images compressées où le débit disponible et les perturbations impactant le canal sont variables, il est nécessaire de pouvoir adapter à tout moment le débit utile d'un flux de données compressées ainsi que sa résistance aux erreurs afin de s'adapter aux variations des contraintes dues à l'environnement. A cet effet, il est donc important de pouvoir estimer les caractéristiques d'un flux d'images compressées telles que le débit et la distorsion de compression, a posteriori de l'étape de compression et sans effectuer de décodage dudit flux. Une autre caractéristique importante à estimer est la distorsion de l'image décodée résultant de la perte ou de la corruption lors de sa transmission d'un paquet contenu dans le flux compressé. Cette information permet de mesurer l'influence sur le rendu visuel des erreurs de propagation qui impactent le canal de transmission et d'en déduire, notamment, une stratégie appropriée de protection du flux compressé par ajout de redondance. La problématique de l'estimation de débit et de distorsion d'un flux d'images compressées est notamment traitée dans le contexte de l'allocation débit/distorsion qui vise à déterminer les meilleurs paramètres de quantification à appliquer par un mécanisme de compression d'images afin d'obtenir une distorsion sur l'image décodée la plus faible possible pour un débit donné. Le document [1] propose plusieurs modèles d'allocation de débit et/ou de distorsion pour le codage d'images et de vidéos, cependant il ne traite pas le problème de l'estimation de la distorsion sur des flux images déjà compressés et sans connaissance a priori de l'image encodée. Des méthodes de modélisation de la sensibilité aux erreurs d'un flux vidéo existent telle que celle divulguée dans le brevet français [2], mais ce dernier se limite au cas d'application au standard vidéo H.264. Le procédé mis en oeuvre dans ce brevet suppose, de plus, une interaction avec le codeur source et/ou la connaissance, du point de vue du récepteur, de la séquence vidéo d'origine transmise. Afin de compenser, notamment, les limitations de l'art antérieur précédemment introduites, l'invention propose une solution permettant d'estimer les paramètres de débit et de distorsion d'une séquence d'images préalablement compressées et ce sans nécessiter de décompression ni de connaissance a priori de l'image encodée. II repose notamment sur l'utilisation d'un modèle de prédiction de la distorsion d'une image décodée lorsqu'un paquet compris dans le flux transmis est perdu. A cet effet, l'invention a pour objet un procédé d'estimation du débit et de la distorsion d'au moins une image compressée à l'aide d'au moins une opération de transformée en ondelettes délivrant au moins une sous-bande b contenant des coefficients, ladite sous-bande étant subdivisée en au moins un paquet de coefficients, ledit paquet comportant en plus un entête de paquet et étant transmis à travers un canal de propagation, ledit procédé étant caractérisé en ce qu'il comprend au moins les étapes suivantes: o une étape d'estimation de la distorsion D liée à la compression 15 des coefficients dudit paquet, o une étape d'estimation de la distorsion D 'b liée à la corruption ou à la perte dudit paquet lors de sa transmission, o une étape d'estimation du débit Lb de compression de ladite sous-bande b, 20 o une étape d'estimation de la sensibilité D b dudit paquet en fonction des distorsions D et D 'b et du débit Lb de ladite sous-bande b, o une étape d'estimation de la distorsion de ladite image liée à la décompression dudit paquet à partir de sa sensibilité D 25 Selon une variante de réalisation de l'invention, ledit paquet est associé à une couche i de qualité de la compression et l'étape d'estimation de la sensibilité est réalisée à l'aide de la relation suivante : ~~yc k l k ( i-1 ~-c Db ù n(1ù P e r r b ) D k +E P e r r b n\1 P e r r b J . LJbi-1 + Perrb • DTb b , i=1 i=2 j=1 J Ib où PerrQi b est la probabilité que le paquet d'indice c; correspondant à la couche de qualité i de la sous-bande b soit erroné, est sa distorsion de compression, Db est la distorsion de corruption de la sous-bande b, et T sont les nombres de coefficients respectivement contenus dans ledit paquet c et ladite sous-bande b. Selon une variante de réalisation de l'invention, l'estimée de la distorsion liée à la compression est obtenue à l'aide de la relation suivante: Ob ù Ob .Atr Db 12 12 où Ab est un pas de quantification lié à la compression et Ni: est un pas de quantification lié à une étape de troncature des coefficients résultant de la compression. Selon une variante de réalisation de l'invention, le pas de quantification 15 L\trb est obtenu à partir de paramètres caractéristiques de ladite sous-bande b tels que le nombre de bits de quantification des coefficients de ladite sous-bande b, le nombre de plans de bits à zéro, et le nombre de passes de codage appliquées aux plans de bits desdits coefficients, les valeurs desdits paramètres caractéristiques étant déterminées à partir desdits entêtes dudit 20 flux de données compressées. Selon une variante de réalisation de l'invention, l'estimée de la distorsion Db liée à la corruption de la sous-bande b est obtenue à l'aide de la relation suivante :10 Db = Db .22L" où Db est la distorsion liée à la compression de la sous-bande b Selon une variante de réalisation de l'invention, l'estimation du débit Lb de la sous-bande b est obtenue en divisant la somme des longueurs exprimées en octets de tous les paquets inclus dans la sous-bande b par le nombre de coefficients contenus dans la sous-bande b. Selon une variante de réalisation de l'invention, le procédé comprend en plus les sous étapes suivantes: o une étape d'estimation de l'erreur quadratique moyenne de ladite image compressée à partir de la sensibilité D b de chaque paquet dudit flux compressé. o une étape d'estimation du débit total L de compression de ladite image obtenu en faisant la somme des débits Lb de chaque sous-bande b pondérés par le rapport entre le nombre de coefficients de la sous-bande b et la taille de l'image en pixels. Selon une variante de réalisation de l'invention, le mécanisme de compression d'images mis en oeuvre respecte le standard JPEG 2000 et ledit paquet est un code-bloc. Selon une variante de réalisation de l'invention, le procédé comprend 20 en plus une étape de transcodage au cours de laquelle lesdites estimées de débit L b et de sensibilité D b sont utilisées pour supprimer, protéger ou adapter le débit d'au moins un desdits paquets du flux compressé. Selon une variante de réalisation de l'invention, ladite étape de transcodage respecte le standard JPWL. L'invention a également pour objet une utilisation du procédé décrit précédemment caractérisée en ce que les distorsions D et D 'b sont utilisées afin de déterminer une mesure de la qualité de ladite image. 25 D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés qui représentent: La figure 1, une illustration de la syntaxe d'un flux de code JPEG 5 2000, La figure 2, une représentation du résultat d'une transformée en ondelettes sur deux niveaux de résolution, La figure 3, une représentation de la notion de plan de bits dans un code-bloc, 10 La figure 4, un exemple de système adapté à exécuter les étapes selon l'invention. Le standard JPEG 2000 (Joint Picture Expert Group 2000), tel que décrit notamment dans [3] définit un mécanisme de compression d'images 15 avec ou sans perte utilisant une transformation en ondelettes. La partie 11 de ce standard référencée sous l'abréviation anglo- saxonne JPWL (JPEG 2000 Wireless) et décrite dans le document [4] définit une syntaxe et des outils permettant de protéger les données encodées au format JPEG 2000 en amont de leur transmission à travers un canal bruité 20 hertzien. Cette phase de protection des données peut s'effectuer ultérieurement à la phase de compression dans une étape connue de l'Homme du métier sous le terme de transcodage. Un encodeur respectant le standard JPWL se décompose principalement en trois modules : un encodeur au standard JPEG 2000 qui délivre un flux compressé, un module 25 de génération de la sensibilité aux erreurs dudit flux compressé et un module de protection qui permet d'insérer de la redondance dans le flux compressé à partir des informations de sensibilité aux erreurs. Ces deux derniers modules sont indépendants du codeur source JPEG 2000 et peuvent s'appliquer directement sur un flux déjà compressé lors d'une étape de transcodage. 30 Cette particularité présente l'avantage de permettre l'insertion de protection aux erreurs à un endroit quelconque de la chaine de transmission et non nécessairement immédiatement après l'étape de compression. Dans la suite de la description, le terme code-bloc est utilisé en référence à un ensemble de coefficients obtenus par transformée en ondelettes et appartenant à une sous-bande donnée. Ce terme est introduit dans la définition du standard [3] JPEG 2000. L'utilisation du terme sensibilité en référence à un paquet de données d'image compressée ou à un code-bloc contenant des données d'image compressée désigne la distorsion sur l'image décodée résultant de la 1 o décompression de ce paquet ou code-bloc. Le terme � sensibilité � est utilisé car il s'agit ici de mesurer l'influence d'erreurs impactant les données contenues dans ledit paquet ou code-bloc sur le rendu visuel de l'image transmise. Le standard JPWL définit des moyens d'insérer dans le flux 15 compressé des informations de sensibilité aux erreurs et des informations de redondance. Il décrit notamment la structure et l'encapsulation des données ajoutées a posteriori de l'encodage. L'information de sensibilité aux erreurs est utilisée en particulier dans un mécanisme de protection inégale aux erreurs qui vise à répartir la protection par ajout de redondance de façon 20 inégale sur les différents paquets qui composent le flux compressé de façon à prendre en compte l'impact différent des erreurs de propagation sur chacun de ces paquets. Par exemple, les entêtes de paquets contenant des informations de signalisation ou des paramètres de décompression indispensables au décodage du flux ont intrinsèquement une sensibilité aux 25 erreurs plus accrue que certaines parties des données d'images compressées. Afin de mettre en oeuvre le type de mécanisme de protection inégale aux erreurs introduit précédemment, il est nécessaire de définir des valeurs de sensibilité pertinentes affectées aux différents paquets qui composent un flux compressé ce que ne décrit pas le standard JPWL. La 30 présente invention propose donc à cet effet une méthode d'estimation du débit et de la distorsion d'un flux d'images compressées permettant de différencier les paquets composant un flux de code JPEG 2000 en fonction de leur sensibilité aux erreurs. La figure 1 schématise la structure d'un flux de données au format résultant de la compression JPEG 2000 d'une image. Ce flux est notamment composé d'un entête principal 101 suivi par un ou plusieurs entêtes de tuiles 102,103,104 qui délimitent les données associées à une tuile qui est une zone rectangulaire de l'image considérée définie en tant que telle par le standard JPEG 2000. Un entête de tuile 101 est suivi de plusieurs paquets 1 o 106 qui comprennent chacun un entête de paquet 107 et des données entropiques 108 associées à une zone de l'image et à une sous bande. Un paquet 107 peut contenir plusieurs codes-bloc. Un flux de code JPEG 2000 se termine toujours par un marqueur de fin de flux 105 connu sous le terme anglo-saxon � End of Codestream � EOC. 15 La figure 2 schématise le principe connu de l'Homme du métier de la décomposition d'une image en sous bandes suite à une transformée en ondelettes, technique connue dans le domaine de la compression de données numériques. Une image 200 de résolution donnée est soumise à 20 une première étape de décomposition en ondelettes aboutissant aux sous- bandes 201a et 201b. Cette première étape consiste à appliquer, selon l'axe vertical 200b de l'image 200, respectivement un filtrage passe bas permettant d'obtenir une première sous-bande 1 L 201 a puis un filtrage passe haut permettant d'obtenir une seconde sous-bande 1H 201b. Dans une 25 seconde étape les mêmes opérations de filtrage passe bas et filtrage passe haut sont appliquées suivant l'axe horizontal 200a afin d'obtenir au total quatre sous-bandes 1 LL,1 HL,1 LH et 1HH. Ces quatre sous-bandes correspondent respectivement à l'image d'origine avec une résolution divisée par deux 202a, et trois images de détails 202b,202c,202d qui donnent les 30 hautes fréquences de l'image originale. L'application du mécanisme précédemment décrit sur l'image 202a de résolution divisée par deux permet d'obtenir, in fine, une image 204 de résolution divisée par quatre par rapport à l'image 200 d'origine. Dans la suite de la description, une sous-bande désignera le résultat d'une étape de filtrage passe haut ou passe bas, selon une transformée en ondelettes, appliquée à l'image d'origine 200 ou à une autre sous-bande 201 a,201 b,202a,202b,202c,202d. Les coefficients d'une sous bande désignent les coefficients obtenus par une opération de transformée en ondelettes. Dans un procédé de l'art antérieur de compression d'une image selon le standard JPEG 2000, les coefficients de chaque sous-bande obtenue sont ensuite quantifiés et divisés en plusieurs codes-blocs qui contiennent les coefficients issus de la décomposition en ondelettes. Les codes-blocs correspondent à l'unité de base du codage entropique opéré par le standard JPEG 2000 et sont codés indépendamment. Les coefficients de ces code- blocs sont quantifiés et le résultat de cette quantification peut être représenté de manière tridimensionnelle, afin de faire apparaître la notion de plans de bits. Un paquet JPEG2000 peut contenir plusieurs code-blocs. Dans l'exemple qui est décrit par la suite, le procédé d'estimation du débit/distorsion selon l'invention est appliqué à un code-bloc mais il peut être généralisé à tout paquet contenant des données compressées associées à une sous-bande obtenue par une opération de transformée en ondelettes d'une image. La figure 3 illustre une représentation en trois dimensions de la notion de plan de bits. Un ensemble de coefficients, par exemple des coefficients d'une sous bande, sont représentés dans un repère à trois axes. Le premier axe 305 représente la dynamique d'un coefficient, c'est-à-dire la représentation binaire du coefficient quantifié. Les deux autres axes 306,307 correspondent aux positions à l'intérieur du code-bloc du coefficient. Un coefficient quantifié comporte un nombre donné de bits que l'on range classiquement du bit le moins significatif connu par l'acronyme anglo-saxon LSB (Least Signifiant Bit) au bit le plus significatif ou MSB (Most Signifiant Bit). Un plan de bit 301,302,303,304 correspond aux bits de même rang de l'ensemble des coefficients d'un code-bloc. Par exemple, le plan MSB 301 contient les bits les plus significatifs de tous les coefficients pour un code-bloc considéré. De façon similaire le plan LSB 304 contient les bits les moins significatifs. Un coefficient 308 est représenté par l'ensemble de ses bits appartenant chacun à un plan de bit distinct. Dans une étape de compression selon le standard JPEG 2000 les coefficients quantifiés sont ensuite encodés, notamment par un codage arithmétique réalisé en trois passes et appliqué à chaque plan de bits en 1 o commençant par le plan du bit le plus significatif. Le procédé selon l'invention consiste, notamment, à estimer le débit et la distorsion de tout ou partie d'un flux d'images compressées, par exemple selon le standard JPEG 2000, à partir d'informations extraites des différents 15 entêtes de ce flux. En particulier, l'invention s'applique à un ou plusieurs codes-bloc ou à tout ensemble regroupant une pluralité de coefficients issus d'une transformée en ondelettes. La distorsion d'une image ou partie d'image transmise au travers d'un canal sans fil puis décompressée résulte de deux phénomènes distincts. 20 Tout d'abord il existe une distorsion D inhérente à l'opération de compression qui impacte chaque code-bloc même si celui-ci est correctement reçu à l'issue de sa transmission sur un canal sans fil donné. Cette première distorsion est essentiellement due aux opérations de quantification qui ont lieu lors de l'étape de compression et qui induisent des pertes d'information. 25 Ensuite, il existe une distorsion D' liée à la perte ou a la corruption d'un paquet ou d'un code bloc lors de la transmission sur un canal soumis à erreurs de propagation. A partir des deux termes de distorsion D et D' introduits, il est possible de déterminer la sensibilité intrinsèque d'un paquet de données compressées à l'aide de la relation suivante : 30 D=(1-Prr)•D+Prr•D' (1) avec Perr la probabilité que le paquet transmis sur un canal de transmission ne soit pas correctement reçu. La sensibilité d'un paquet désigne la contribution dudit paquet à la distorsion engendrée sur l'image décodée. Le procédé selon l'invention permet notamment de déterminer, pour chaque paquet, les deux termes de distorsion D et D' utiles au calcul de la sensibilité intrinsèque dudit paquet. Le document [1] introduit un modèle asymptotique d'évaluation de la distorsion de compression d'une sous-bande en fonction du pas de quantification utilisé pour quantifier les coefficients de cette sous-bande obtenus par une étape de transformée en ondelettes. Sous une hypothèse de haute résolution, il est connu que la densité de probabilité des coefficients d'une transformée en ondelettes peut être approchée par une distribution Gaussienne généralisée de moyenne nulle pour chaque sous-bande. Les coefficients d'une transformée en ondelettes sont supposés indépendants et identiquement distribués dans chaque sous-bande. Les hypothèses de stationnarité et d'ergodicité sont également retenues de façon à pouvoir estimer les caractéristiques statistiques des différentes sous-bandes simplement. La distorsion D due à la compression peut être reliée au pas de quantification A à l'aide de la relation suivante : D=12(2) Afin de déterminer la distorsion liée à la quantification scalaire et à l'étape de troncature opérée lors du codage par plan de bits des coefficients de transformée en ondelettes, la relation (2) est appliquée à un code-bloc. L'étape de troncature a pour but notamment de diminuer le débit en supprimant les plans de bits les moins significatifs. Dans le cadre d'une compression respectant le standard JPEG 2000 cette étape est mise en oeuvre par un mécanisme d'allocation en débit/distorsion. Elle est assimilable à une quantification scalaire car la troncature de N plans de bits est équivalente à une division par 2N. On obtient ainsi le pas de quantification équivalent Ob associé à un code-bloc c d'une sous bande b à l'aide de la relation suivante : A'b ùOb.LVYb (3) avec Ab le pas de quantification appliqué à la sous-bande b et Atib , le pas de quantification lié à l'étape de troncature. Le pas de quantification appliqué à une sous-bande b est défini dans le standard JPEG 2000 par la relation suivante : i + Pb l 2") (4) où RB est la dynamique nominale de la sous-bande b , RB = RI + log 2 \Yb) • (5) Les valeurs des différents paramètres R,, Eb, Pb et yb utiles à la détermination du pas de quantification de la sous-bande b peuvent être extraites des 15 entêtes des paquets du flux compressé au format JPEG 2000. R, est le nombre de bits utilisés pour représenter les composantes des pixels de l'image source et peut être déterminé par le marqueur SIZ(� Image and Tile Size Marker �) présent dans l'en-tête principal 101 et qui définit, selon le standard JPEG 2000, la taille d'une tuile et de l'image associée. 20 Eb est l'exposant du pas de quantification de la sous-bande b qui peut être déterminé à partir du marqueur QCC (� Quantization Component Marker �) inclus dans l'entête principal 101 et dans l'entête de tuile 102,103,104. Pb est la mantisse du pas de quantification de la sous-bande b et est spécifiée dans le marqueur QCD (� Quantization Default Marker �) présent 25 dans l'entête principal 101 et dans l'entête de tuile 102,103,104. yb est le gain de la sous-bande b qui est défini dans le standard JPEG 2000 en fonction du type de sous-bande. 10 A ù 2RBùeb b Afin de déterminer le second terme Nrb nécessaire au calcul du pas de quantification associé à un code-bloc, il faut connaître le nombre de plans de bits tronqués lors de l'étape d'allocation en débit et distorsion selon le standard JPEG 2000. Cette information peut être obtenue à partir du nombre de plans de bits conservés qui peut être calculé grâce aux données stockées dans les en-têtes des paquets. Soient nb bitplanesbc le nombre de plans de bits conservés par l'algorithme d'allocation en débit/distorsion mis en oeuvre dans le cadre d'une compression respectant le standard JPEG 2000 et nb bitplanes trb° le nombre de plans de bits tronqués dans le code-bloc 1 o d'indice c de la sous-bande d'indice b. Chaque plan de bits est codé en trois passes, excepté le premier plan de bits qui est codé en une passe et le bit de signe qui est codé lors de la première passe dite de signifiance. Le dernier plan de bits est supposé conservé si le résultat d'une au moins des passes de codage appliquée à ce plan de bits est conservé. Enfin le nombre de 15 plans de bits à zéro qui permet de déterminer le bit le plus significatif est codé séparément. Le nombre de plans de bits conservés nb bitplanesb' est donc donné par la relation suivante: nb_bitplanesb` = IMSBbc + nb _passesb -1 + 2 (6), 3 20 où IMSBb° est le nombre de plans de bits à zéro situés avant le plan du bit le plus significatif. La valeur de ce nombre de plans de bits à zéro peut être extraite de l'entête du paquet contenant le code-bloc d'indice c de la sous-bande b, nb passesb° est le nombre de passes de codage conservées par l'algorithme 25 d'allocation en débit/distorsion qui peut être également extrait de l'entête du paquet. Le nombre de plans de bits tronqués correspond à la différence entre le nombre de bits Mb restants après l'étape de quantification scalaire et le nombre de bits conservés lors de l'étape d'allocation en débit/distorsion nb_bitplanesbc. nb _ bitplanes _ trbc = Mb û nb _ bitplanesb° , (7) Le nombre de bits Mb peut être obtenu à partir des marqueurs QCC et QCD 5 présents dans l'entête principale 101 et les entêtes de tuiles 102,103,104,105. A partir de la relation (7) on peut en déduire le pas de quantification lié à la troncature : Ac nb_bitplanes tri At-r' = 2 (8) 10 A partir des relations (2) à (8) introduites précédemment, on peut déterminer la distorsion Db° induite par l'étape de compression d'un codeur JPEG 2000 pour un code-bloc d'indice c appartenant à la sous-bande d'indice b en fonction d'informations extraites des marqueurs présents dans les entêtes du 15 paquet qui contient ce code-bloc. Les relations suivantes résument les dépendances, représentées par l'intermédiaire des fonctions f et g, entre la distorsion, les deux termes de quantification et les informations extraites des entêtes : Dc Ob Ob.Otrb b 12 12 Ab = f (RB,eb,,ub,gainb) (9) Atrb = g(Mb,IMSBb,nb_ passes La deuxième étape du procédé selon l'invention consiste à estimer la 20 distorsion liée à l'impact de la perte d'un paquet JPEG 2000 lors de sa transmission. Cette distorsion dépend en partie de la notion de couche de qualité associée à un code-bloc. Une couche de qualité d'indice i correspond 25 au plan de bits de rang i tel qu'explicité précédemment à l'appui de la figure 3. A titre d'exemple, le plan des bits les plus significatifs correspond à une couche de qualité de base alors que tout plan de bits supplémentaire permet d'introduire une couche de qualité supplémentaire encore appelée couche de qualité de rehaussement ou de raffinement. Deux cas de figure peuvent se présenter. Tout d'abord, si le paquet perdu 5 est associé à une couche de qualité de rehaussement, alors la distorsion Db'' résultante de la perte de ce paquet est égale à la distorsion de compression D, d'une couche de qualité inférieure transmise correctement. Si c; est l'indice du code-bloc correspondant à la ième couche de qualité d'une sous bande d'indice b, on peut expliciter ce résultat à l'aide 1 o de la relation suivante : D'bC~ = Dbr ' pour i >1 (10) Si, au contraire, le paquet perdu est associé à une couche de qualité de base, alors la distorsion résultante est équivalente à l'énergie du code- 15 bloc de la sous-bande considérée, c'est-à-dire à la distorsion de compression induite par un pas de quantification de taille infinie. Pour estimer, dans ce dernier cas, la distorsion induite lors de la corruption ou la perte d'un code-bloc, on pose l'hypothèse que les coefficients issus de la transformée en ondelettes d'une sous-bande b 20 suivent une loi statistique Gaussienne généralisée de variance ab et de paramètre ab. II est connu de l'Homme du métier, notamment à travers le document [1] qu'une relation existe entre le débit Lb de la sous-bande b et la distorsion Db due à la compression dans chaque sous-bande b : Db = F(ab, Qb) X 2-2Lb (11) 25 avec F(ab, ab) une fonction dépendant de la statistique d'une répartition gaussienne généralisée associée à une sous-bande b. La distorsion totale Db de la sous-bande b est obtenue en moyennant les distorsions de compression D°b de chaque code-bloc de la sous-bande b estimées précédemment à l'aide des relations (9) et compte tenu du nombre Tb de coefficients de la sous-bande b : Ncbb -1 E Dbc Db= c 0 Tb (12), où Ncbb est le nombre de code-blocs contenus dans la sous-bande b. La relation (11) nécessite également de déterminer le débit Lb de la sous-bande b qui est par ailleurs également nécessaire afin de modéliser la courbe de débit/distorsion de la sous-bande b. Le débit Lb est égal à la somme des débits Lb' de chaque code-bloc appartenant à la sous-bande b et est obtenu à l'aide de la relation suivante : Ncbb -1 Lb = LbC = c=0 Ncbb -1 1 L' c c=0 Tb (13) où L'b° est la longueur en octets du code-bloc d'indice c qui peut être extraite de l'entête du paquet contenant ledit code-bloc. Le débit Lb est exprimé en nombre de bits par pixels. La distorsion liée à la perte de la sous-bande b équivaut à la distorsion 15 résultante lorsque la sous-bande n'est pas transmise c'est-à-dire lorsque son débit Lb vaut 0, ainsi on peut écrire : Db = F(ab,6b) = Db.224 (14) A partir des relations (11) à (14) on peut donc déterminer F(ab, ab) puis D'b. Enfin la distorsion Db résultante de la perte d'un code-bloc d'indice c1 20 correspondant à la première couche de qualité, c'est-à-dire la couche de qualité de base, s'écrit : Dbcl = Db Tb (15) Tb où Tb° est le nombre de coefficients contenus dans le code-bloc d'indice c de la sous-bande b. Les relations (10) et (15) obtenues précédemment permettent donc de 5 déterminer la distorsion résultante de la perte d'un code-bloc quel que soit la couche de qualité considérée. Si on considère que la compression selon le standard JPEG 2000 a engendré k couches de qualité pour un même code-bloc c d'une sous-bande 10 b, alors la distorsion totale dudit code-bloc résultant de l'encodage puis de la transmission à travers un canal bruité peut s'écrire à l'aide de la relation suivante obtenue en regard de la relation (1) : C Ci Ck Db ù~(iùPerrb) Db Db'c� (16) k ( i 1 Ci 1ù1 (1 Ci Perrb (i ù Perrb ) i=1 i=1 j=1 où Perf1 b est la probabilité que le code-bloc d'indice c correspondant à la ième 15 couche de qualité de la sous-bande b soit erroné et 11 désigne un produit. La relation (16) peut s'écrire encore de la façon suivante : r~C k r~C k ci Perrb ' 1 Db ). 11(1ù vb ù J(i ù p vbk + E Perrb . Ubi Perrb l b 17 + • ( ) i=1 i=2 j=1 T ~ b La relation (17) permet donc de déterminer la sensibilité intrinsèque de 20 chaque code-bloc JPEG 2000 à partir des informations contenues dans les en-têtes présents dans le flux compressé. Cette information de sensibilité peut être utilisée de façons diverses. Notamment, cette information permet de mettre en oeuvre un mécanisme de protection différenciée aux erreurs par ajout de redondance selon le standard JPWL ou tout autre mécanisme 25 permettant d'allouer un code correcteur à chaque zone du flux compressé dont on a préalablement déterminé une estimation de sa sensibilité et de son débit à l'aide du procédé selon l'invention. On peut également choisir de ne pas transmettre un code-bloc qui présente une sensibilité trop importante en cherchant à minimiser la distorsion globale d'une image en regard d'une contrainte de débit total. A partir des relations précédemment établies, on peut déduire pour une image entière le débit total Ltot et l'erreur quadratique moyenne totale MSEtot en tenant compte des gains de reconstruction des filtres de synthèse. Nsb-1 MSEtot = abDb (18) b=0 où ab est le gain du filtre de reconstruction de la sous-bande b et Nsb est le 10 nombre de sous-bandes dans une image. Le débit total de l'image peut également être estimé à partir de la relation suivante : Nsb-1 Ltot = ab Lb (19) b=0 avec ab le rapport entre le nombre de coefficients de la sous bande b et la 15 taille de l'image en pixels. Dans le cas d'une transmission avec ajout de redondance, il faut tenir compte de la redondance ajoutée dans le calcul du débit total ainsi que de l'influence du code correcteur associé à cette redondance dans la détermination de la probabilité d'erreur d'un code-bloc. 20 De façon générale, l'invention permet d'adapter le flux de données compressées en fonction de la nature et de l'état du canal de transmission. La figure 4 schématise un exemple de système adapté à exécuter les étapes du procédé selon l'invention. 25 Un serveur 401 de diffusion d'images compressées, par exemple un serveur JPEG 2000, transmet un flux d'images compressées à travers un réseau 402 lequel est relié à une station de base 403 comprenant un module 404 de transcodage du flux reçu et des couches d'accès radio 405. Le module de transcodage 404 est adapté à mettre en oeuvre le procédé d'estimation de débit/distorsion selon l'invention à partir du flux compressé reçu. Ce module peut également recevoir depuis les couches d'accès radio 405 des informations sur l'état du canal ou sur le débit disponible permettant d'utiliser l'estimation de débit/distorsion selon l'invention afin d'adapter le flux compressé ou de lui appliquer un mécanisme de protection contre les erreurs. La station de base 403 est reliée à une antenne d'émission 406 apte à transmettre le flux d'images compressées à travers un canal sans fil 407 à 1 o destination d'une station mobile 409 qui reçoit ledit flux par l'intermédiaire d'une antenne réceptrice 408. La mobile 409 comprend notamment des couches d'accès radio 411 ainsi qu'un décodeur image 410 apte à décompresser le flux reçu et produire une ou plusieurs images 412. 15 Le procédé selon l'invention présente notamment l'avantage d'estimer le débit et la distorsion d'un flux d'images sans opérer de décompression de ce flux. Les estimées de débit et de distorsion selon l'invention permettent d'adapter le flux de données en fonction des conditions de transmission et notamment d'appliquer une stratégie de protection aux erreurs qui prend en 20 compte la sensibilité intrinsèque de chaque zone du flux compressé. L'invention a encore comme avantage de permettre d'estimer a posteriori de la compression l'impact de la perte d'un paquet de données d'images. Le procédé selon l'invention peut s'appliquer directement au niveau du codeur source ou être intégré dans un module de transcodage ou bien 25 encore à tout autre niveau d'une couche protocolaire tel que la couche transport ou réseau. Cet avantage permet d'adapter le flux de données à n'importe quel emplacement de la chaîne de transmission globale. Enfin l'invention peut être intégrée dans un mécanisme dit de Qualité de service connu sous l'abréviation anglo-saxonne QoS (Quality Of Service ). En 30 particulier, l'invention permet d'estimer une mesure de la qualité image en fonction des distorsions D b et D 'b pondérée par un facteur de type psycho-visuel. Cette mesure, ou métrique, de la qualité psycho-visuelle d'une image peut être utilisée dans un mécanisme de gestion de la qualité de service puisqu'elle donne accès, à chaque point ou noeud de la chaîne de transmission d'une image, à une estimation de l'impact dû aux étapes précédentes en termes de dégradations de la qualité psycho-visuelle de l'image. Références [1] Christophe Parisot, � Allocations basées modèles et transformée en ondelettes au fil de l'eau pour le codage d'images et de vidéos �, Thèse de l'université de Nice ûSophia Antipolis 20 janvier 2003 [2] Cyril Bergeron, Catherine Lamy-Bergot, � Procédé permettant de déterminer des paramètres de compression et de protection pour la transmission de données multimédia sur un canal sans fils � Brevet FR2903272 [3] International standard ISO/IEC 15444-1 ITU-T recommendation T.800 (08/2002) [4] Information technology ù JPEG 2000 image coding system ù Part 11: 20 Wireless, ISO/IEC FCD 15444-11, 03/2005 REVENDICATIONS1. Procédé d'estimation du débit et de la distorsion d'au moins une image (200) compressée à l'aide d'au moins une opération de transformée en ondelettes délivrant au moins une sous-bande b (202a,202b,202c,202d) contenant des coefficients, ladite sous-bande étant subdivisée en au moins un paquet (106) de coefficients, ledit paquet (106) comportant en 1 o plus un entête (101,102,103,104,107) de paquet et étant transmis à travers un canal de propagation, ledit procédé étant caractérisé en ce qu'il comprend au moins les étapes suivantes: o une étape d'estimation de la distorsion D b liée à la compression des coefficients dudit paquet (106), 15 o une étape d'estimation de la distorsion D 'b liée à la corruption ou à la perte dudit paquet (106) lors de sa transmission, o une étape d'estimation du débit Lb de compression de ladite sous-bande b, o une étape d'estimation de la sensibilité D b dudit paquet (106) en fonction des distorsions D b et D 'b et du débit Lb de ladite sous-bande b, o une étape d'estimation de la distorsion de ladite image (200) liée à la décompression dudit paquet (106) à partir de sa sensibilité D b 2. Procédé selon la revendication 1 caractérisé en ce que ledit paquet (106) est associé à une couche i de qualité de la compression et que l'étape d'estimation de la sensibilité est réalisée à l'aide de la relation suivante : 20 25c, Db +Prrb'T Tb b k k ( fl(i u = rrb ci E ci l k + errb ù P r rb l e b i=1 i=2 j=1 ~ où PeCf 1 b est la probabilité que le paquet d'indice c; correspondant à la couche de qualité i de la sous-bande b soit erroné, est sa distorsion de compression, Db est la distorsion de corruption de la sous-bande b, b c et Tb sont les nombres de coefficients respectivement contenus dans ledit paquet c et ladite sous-bande b. 3. Procédé selon la revendication 2 caractérisé en ce que l'estimée de la distorsion Deb' liée à la compression est obtenue à l'aide de la relation suivante: D1 = Ob _ ob.otrb b 12 12 où Ab est un pas de quantification lié à la compression et Itrb est un pas de quantification lié à une étape de troncature des coefficients résultant de la compression. 4. Procédé selon la revendication 3 caractérisé en ce que le pas de quantification Ltrb est obtenu à partir de paramètres caractéristiques de ladite sous-bande b tels que le nombre de bits de quantification des coefficients de ladite sous-bande b, le nombre de plans de bits à zéro, et le nombre de passes de codage appliquées aux plans de bits desdits coefficients, les valeurs desdits paramètres caractéristiques étant déterminés à partir desdits entêtes (101,102,103,104,107) dudit flux de données compressées. 5. Procédé selon la revendication 2 caractérisé en ce que l'estimée de la distorsion Db liée à la corruption de la sous-bande b est obtenue à l'aide de la relation suivante : D' - D 22' b ù b' où Db est la distorsion liée à la compression de la sous-bande b 6. Procédé selon l'une des revendications précédentes caractérisé en ce que l'estimation du débit Lb de la sous-bande b est obtenu en divisant la somme des longueurs exprimées en octets de tous les paquets inclus 1 o dans la sous-bande b par le nombre de coefficients contenus dans la sous-bande b. 7. Procédé selon l'une des revendications précédentes caractérisé en ce qu'il comprend en plus les sous étapes suivantes: 15 o une étape d'estimation de l'erreur quadratique moyenne de ladite image (200) compressée à partir de la sensibilité D b de chaque paquet (106) dudit flux compressé. o une étape d'estimation du débit total L de compression de ladite image (200) obtenu en faisant la somme des débits Lb de 20 chaque sous-bande b (202a,202b,202c,202d) pondérés par le rapport entre le nombre de coefficients de la sous-bande b et la taille de l'image (200) en pixels. 8. Procédé selon l'une des revendications précédentes caractérisé en ce 25 que le mécanisme de compression d'images mis en oeuvre respecte le standard JPEG 2000 et ledit paquet (106) est un code-bloc. 9. Procédé selon l'une des revendications précédentes caractérisé en ce qu'il comprend en plus une étape de transcodage (404) au cours de laquelle lesdites estimées de débit L b et de sensibilité D b sont utilisées pour supprimer, protéger ou adapter le débit d'au moins un 5 desdits paquets (106) du flux compressé. 10. Procédé selon la revendication 9 caractérisé en ce que ladite étape de transcodage (404) respecte le standard JPWL. io 11. Utilisation du procédé selon l'une des revendications précédentes caractérisée en ce que les distorsions D b et D sont utilisées afin de déterminer une mesure de la qualité de ladite image (200).