Brevet FR2769453A1 — publié le 19990409

PROCEDE D'EVALUATION DE LA DEGRADATION D'UNE IMAGE VIDEO INTRODUITE PAR UN SYSTEME DE CODAGE ET/OU DE STOCKAGE ET/OU DE TRANSMISSION NUMERIQUE La présente invention a pour objet un procédé d'évaluation de la dégradation d'une image vidéo, introduite par un système de codage et/ou de stockage et/ou de transmission numérique, notamment un système mettant en oeuvre un codage et/ou un stockage et/ou une transmission de signal video numérique avec un bas débit. Les dégradations sont en général dues aux erreurs introduites par les algorithmes de réduction de débit et/ou par les liaisons défectueuses de transmission ou bien encore aux défauts des codeurs et des décodeurs. La numérisation des signaux vidéo a ouvert la possibilité de pouvoir copier, stocker ou transmettre l'information vidéo en maintenant une qualité constante. Cependant, la grande quantité d'information véhiculée par les images vidéo nécessite en pratique l'utilisation de méthodes de compression numérique pour réduire le débit binaire. Une méthode de compression très répandue dans le domaine de la vidéo est décrite dans la norme ISO/CEI 13918 MPEG2. Cet algorithme est dit du type "avec pertes", car l'image restituée après le décodaae n'est pas identique à l'origInal. Afin de maintenir une qualité acceprable pour le téléspectateur final, les algorithmes de réduction de débit tiennent compte des propriétés perceptuelles du système humain de vision. En dépit de ceci, les contraintes imposées de débit ou de largeur de bande disponible pour la transmission ainsi que le contenu du signal vidéo impliquent l'apparition de dégradations caractéristiques sur le signal après décodage. Ces dégradations introduites par la chaîne, par exemple MPEG2, influent directement sur la qualité perçue de l'image finale. L'évaluation automatique de la qualité des signaux audiovisuels a un large champ d'applications dans la chaîne de télévision numérique : production, distribution, et évaluation des performances des systèmes. Par ailleurs, les dispositifs existants ont été élaborés pour des tests en laboratoires et ne sont pas adaptés pour la télésurveillance des réseaux de distribution. Des séquences d'images codées MPEG et transmises à bas débit pour la diffusion de la télévision numérique ou pour d'autres applications multimédia, vont présenter un certain nombre de défauts ou de déformations par rapport à la séquence originale. En effet, une liste non exhaustive de dégradations visibles peut être dressée. Les plus perceptibles d'entre elles sont les erreurs granulaires, les déformations de contours, les pertes d'information, les contours "exotiques", les effets de bloc, etc. . Cependant, des erreurs de transmission peu importantes peuvent se traduire par des effets plus ou moins localisés sur l'image. En cas de perturbation importante, elles peuvent se traduire par des difficultés d'accès aux inrormations, par exemple coupures de service ou gels d'images pendant un temps plus ou moins long selon ia perturbatIon. L'envergure des erreurs dépend du niveau de pertinence et de la structure des données qu'elles affectent : mots de synchronisation, vecteurs mouvement, images codées avec ou sans prédiction ou images de base pour les prédictions. En plus des coupures ou gels d'images, les dégradations constatées se matérialisent par des blocs ou des macro-blocs erronés ou mal positionnés. Ceci a pour effet de propager les dégradations sur toute la séquence vidéo jusqu'à l'image codée sans prédiction, donc codée indépendamment des autres. Un procédé d'évaluation été proposé par le NTIA (National Telecommunications and Information Administration) dans l'article de A.A. WEBSTER et Collaborateurs, intitulé "An objective videoquality assesment system based on human perception" et publié en juin 1993 dans la revue SPIE vol.13, p.15-26. Ce procédé met en oeuvre une analyse des images dégradées et des images d'origine, après les avoir filtrées par deux opérateurs de SOBEL (matrices 3x3) vertical et horizontal. L'image filtrée est obtenue par convolution en faisant glisser les matrices de SOBEL horizontalement et verticalement, et les résultats obtenus représentent les gradients verticaux et horizontaux de l'image. En d'autres termes, l'image filtrée met en valeur les contours verticaux et horizontaux contenus dans l'image initiale (non filtrée). Une mesure basée sur cette information permet de mettre en évidence le changement de contenu entre l'entrée du système vidéo et sa sortie. Le procédé proposé par le NTIA met en oeuvre deux paramètres - d'une part, l'informatIon spatiale SI qui représente l'écart moyen mesuré sur les pixels de l'image filtrée par l'opérateur de SOBEL. Il s'agit ici ae déterminer l'écart moyen au niveau des contours de l'image filtrée, en considérant cue les contours sont importants pour la vision et qu'ils sont affectés par les différents traitements des systèmes numériques à réduction de débit - d'autre part, l'information temporelle TI qui représente l'écart-type de L'image de différence entre deux images successives, cet écart-type étant calculé à partir des différences entre les valeurs des mêmes pixels de deux trames successives. Le paramètre TI peut révéler un mouvement saccadé dû à un défaut du codeur. Le procédé proposé par le NTIA met en oeuvre un calcul comparatif des paramètres SI et TT sur le signai vidéo numérique, entre une image d'entrée et une mage de sortie d'un système. Ce procédé présente un certain nombre d'inconvénients. Le filtre de SOBEL conserve seulement certaines fréquences de l'image qui permettent de prendre en compte la perte de contour; ce qui fait que la perte de définition qui pourrait être prise en compte est celle qui se situe dans la gamme des fréquences conservées. En d'autres termes, la perte de définition ne peut être prise en compte que de manière partielle. D'autre part, le paramètre SI prend en compte des défauts qui tendent à se compenser. En effet, la perte d'information de l'image tend à faire diminuer le paramètre SI alors que les faux contours et les effets de bloc tendent au contraire à le faire augmenter, ce qui fait que le paramètre SI n'est significatif que si l'un ou l'autre des phénomènes est dominant. Enfin, la méthode de calcul des paramètres SI et TI, en utilisant un écart type calculé de manière globale sur l'image entière, réduit de manière drastique 1 impact sur ces paramètres des dégradatIons localsees. La présente invention a pour objet un procédé qui permet de remédier au moins en partIe aux défauts précites. Le procédé selon la présente invention exploite une transformée par blocs, par exemple la transformée en cosinus discrète utilisee notamment dans le standard MPEG, pour mettre en valeur des signatures caractéristiques des défauts identifiés. Cette approche originale permet non seulement de mesurer finement les erreurs introduites, mais elle permet également de tenir compte du contenu initial du signal vidéo et des algorithmes mis en oeuvre dans MPEG. Les transformations par blocs d'une image (transformée de 'ourler, transformée en cosinus discrète TCD, etc...) sont obtenues par l'opération [Fn,m] = [T] . [fn,m] . [T]t où f(x,y) désigne le bloc d'image à transformer et > (x,y) la matrice de la transformation. Une autre transformée par blocs est réalisée à partir de la transformée en ondelettes de l'image, en réorganisant les coefficients d'ondelettes de manière à obtenir des blocs transformés ayant la taille désirée et en particulier de même taille que les blocs obtenus par les méthodes précitées. Un tel mode de réorganisation est indiqué dans l'article de R. de QUEIROZ et Collaborateurs, intitulé "Wavelet transforms in a JPEG - like Image Coder", publié en avril 1997 dans la revue IEEE Transactions On Circuits and Systems for Video Technology, vol. 7 n"2 p. 419-424. L'idée de base de 1 invention est en particulier de réaliser des calculs sur les blocs selon lesquels l'image transmise a été codée, de manière à générer un paramètre significatif exempt de l'effet de bloc. L'invention concerne ainsi un procédé d'évaluation de la dégradation d'unie image vidéo codée par blocs de points images ou pixels, cette dégradation étant générée par un système de codage et/ou de stockage et/ou de transmission produlsant une image de sortie à partir d'une image d'entrée, caractérisé en ce qu'il comporte les étapes suivantes a) sélectionner une image d'entrée et déterminer une activité spatiale (SA) de l'image d'entrée dans une fenêtre d'analyse représentant au moins une partie de l'image qui présente un ensemble de dits blocs de points images ou pixels, cette détermination mettant en oeuvre les sous-étapes suivantes l) déterminer pour chaque bloc (n,m) du dit ensemble de blocs de pixels, les coefficients transformés Fn,m (i,j) par une transformée par blocs selon les dits blocs ll) déterminer à partir des coerilcients transformés Fn,m (i,j), l'activité spatiale bsa de chaque oloc dudit ensemble de blocs de pixels iii) déterminer à partir de l'activité spatiale de bsa de chaque bloc l'activité spatiale globale SA1 de l'ensemble des blocs de pixels constituant la fenêtre d'analyse b) sélectionner l'image de sortie correspondant à l'image d'entrée et déterminer ladite activité spatiale globale SA2 de l'image de sortie dans ladite fenêtre d'analyse, en mettant en oeuvre as sousétapes a i) à a iii) c) comparer l'activité spatiale globale (SA2) de l'image de sortie dans la fenêtre d'analyse et l'activité spatiale globale (SA1) de l'image d'entrée dans la fenêtre d'analyse. L'activité spatiale bsa d'un bloc (n,m) peut être déterminée par la formule suivante L'activité spatiale bsa peut être obtenue comme indiqué ci-dessus par combinaison des carrés des composantes de l'image à partir de la quasi-totalité des composantes de la TCD. Toute autre fonction de combinaison des composantes de a t ransformée utilisée est applicable pour caractériser le contenu de ~'limage. Une fonction plus générale est où k(i,j) est un coefficient constant de pondération de la composante i,j utilisée, et p une constante. Il y a plusieurs possibilités pour choisir les constantes k(i,j), er le choix de l'un ou l'autre se fait en fonction de l'application recherchée. En effet, la fonction k(i,j) est choisie suivant la caractéristique à mettre en exergue : tenir compte du système visuel ou extraire une partie de l'information pertinente. Les cas suivants sont proposés 1) le procédé numéro un est à utiliser lorsque l'on privilégie une bonne corrélation avec la perception humaine. Les valeurs des constantes k(i,j) sont initialisées en fonction de l'importance relatIve des coefficients de la transformée pour l'oeil humain, notamment sa sensibilité fréquentielle, afin de fournir un paramètre d'activité représentatif de ce qui est perçu. Par exemple, dans le cas de la TCD, on prend k(i,j) = 1/Q(i,j). Les Q(i,j) sont les composantes de la matrice de quantification exploitée pour la réduction de débit, indiquées dans le document suivant qui est extrait de la norme MPEG-2 : ISO/IEC CD 13818-2 : "Information technology - Generic coding of moving pictures and associated audio information - Part 2 : video", Dec 1993, p 45, 5 2) le second procédé est utilisé lorsque les coefficients par exemple TCD affectés par la dégradation sont identifiables. Les constantes a) les coefficients correspondant aux hautes fréquences spatiales sont souvent les plus affectés par la compression. Un exemple de matrice de pondération suivant le rang du coefficient par exemple TCD est donné dans le tableau ci dessous 00001111 00001111 O 0 o-o 1 1 1 1 00001111 11111111 11111111 11111111 11111111 b) certains coefficients de faible amplitude sont ramenés à zéro lors de l'étape de compression. Afin de choisir ces coefficients, une moyenne pondérée de chacun des coefficients sur la zone d'image analysée (de taille M.N. blocs) est effectuée ,où Q(i,j) est défini comme ci-dessus, en 1) Les coefficients dont les moyennes sont parmi les 48 plus faibles valeurs sont retenus. Pour ces derniers, Dans le cas du paramètre "activité spatiale", on fixe k(i,j)=1, et p=2 pour la suite du document relative à la description des figures. L'activité spatiale globale SA de l'ensemble des blocs peut être alors déterminée par la formule suivante (HxW) représentant le nombre de blocs de pixels dans la fenêtre d'analyse. Ladite comparaison (étape c) est avantageusement effectuée à l'aide du paramètre R défini de la manière suivante LR = gj[fi(SAl,SA2)] avec i, j#{1,2} avec f(t,y)=x-y ou f2(x,y)= x/y et g1(x)= 100H ou et par exemple : LR = long10 (SAe/SAs). Le procédé peut être alors caractérisé en ce qu'il met en oeuvre le calcul du paramètre LR sur une séquence de M images d'entrée et de M images de sorties correspondants, et en ce qu'il met en oeuvre le calcul d'un paramètre d'évaluation ou note de qualité MLR défini de la manière suivante MLR = Maximum (LR) Le procédé peut être avantageusement caractérisé en ce qu'il comporte la détermination d'un indicateur d'activité temporelle TA pour un groupe de M images, déterminé de la manière suivante rsA(U), pour u variant de o à M-1, désignant les M coefficients d'une transformée par blocs, par exemple une transformée en cosinus discrète appliquée à M échantillons temporels de l'actIvité spatiale SA. Le procédé peut être caractérisé en ce qu'il comporte, en vue de la détermination de l'effet de bloc, les étapes suivantes d) déterminer l'activité spatiale Adl de l'image d'entrée dans une fenêtre d'analyse décalée présentant des blocs de points image ou pixels qui sont décalés d'au moins un pixel dans la direction des lignes de l'image et/ou dans une direction perpendiculaire à la direction des lignes de l'image par rapport auxdits blocs de pixels de ladite fenêtre d'analyse, cette détermination mettant en oeuvre les sous-étapes suivantes i) appliquer ladite transformée par blocs sur chaque bloc décalé (n,m) dudit ensemble de blocs de pixels de la fenêtre d'analyse décalée pour déterminer les coeeficients transformés Fdn,m (i,j), ii) déterminer à partir des coefficients transformés Fdn,m (i,j) des blocs décalés, l'activité spatiale bsad de chaque bloc dudit ensemble de blocs décalés, iii) déterminer l'activité spatiale globale SAdi de l'ensemble des blocs décalés constituant la fenêtre d'analyse décalée, e) déterminer l'activité spatiale globale SAd2 de l'image de sortie dans ladite fenêtre d'analyse décalée en mettant en oeuvre pour l'image de sortie les sous-étapes d i) à d iii) f) comparer d'une part l'activité spatiale globale SAd2 de l'image de sortie dans la fenêtre d'analyse décalée et l'activité spatiale globale SA2 de l'image de sortie dans la fenêtre d'analyse, pour évaluer l'effet de bloc dans l'image de sortie et d'autre part l'activité spatiale SAd1 de l'image d'entrée dans la fenêtre d'analyse décalée et l'activité spatiale SA1 de l'image d'entrée dans la fenêtre d'analyse pour évaluer l'effet de bloc dans l'image d' entrée. L'effet de blocs se caractérise de deux manières différentes : - intrinsèquement, c'est-à-dire par une fonction représentative du contenu de l'image et de l'image seule. - de manière différenciée, c'est-à-dire qu'il est nécessaire de comparer deux valeurs d'une fonction représentative du contenu de l'image : la première calculée sur l'image de référence (image d'entrée) et la seconde sur l'image dégradée (image ce sortie). BM =gk[f@(f@(SAd2,Sd2),f@(SAd1,SA1))],avec i,j,k # {1,2} et on a f1(x,y) = x - y ou f2(x,y)=x/y g1(x)=100.|x| ou g2(x)= 100.|log(|x|)| En particulier, la formule est utilisée. L'étape f met avantageusement en oeuvre la détermination d'un indicateur BM d'effet de bloc, selon la formule suivante SAd2 / SA2 BM = - x 100 SAd1 / SA1 D'autres caractéristiques et avantages de l'invention apparaitront mieux à la lecture de la description qui va suivre, donnée à titre d'exemple non limitatif, en liaison avec les dessins ci-annexés, dans lesquels - la figure ' est une représentation schématique d'un dispositif pour la mise en oeuvre du procédé selon la présente invention, - la figure 2 représente une grille alignée pour le calcul du paramètre d'évaluation SA et une grille décalée pour le calcul du paramètre Sauf, - la figure 3 représente l'évolution du paramètre BM au cours d'une séquence filmée particulière, - et les figures 4a à 4c représentent respectivement l'évolution temporelle du paramètre SA (4a) et la TCD appliquée à M échantillons successifs de SA (4b et 4c). Selon la figure 1, des images d'entrée 1e sont codées par un codeur CCD et transmises par un système de transmission TRANS, par exemple un système de transmission par satellite, puis décodées à la réception par un décodeur DEC pour donner une image de sortie décodée s. Un circuit TCDE applique å titre d'exemple une transformée en cosinus discrète à l'image d'entrée Ie et un circuit PARE détermine les paramètres Pe de l'image d'entrée 1e qui seront utiles à l'évaluation à réaliser. De même, un circuit TCDS applique à titre d'exemple une transformée en cosinus discrète à l'image de sortie Is et un circuit PARS détermine les paramètres P5 de l'image de sortie Is qui seront utiles à l'évaluation à réaliser. Un circuit de comparaison COMP effectue des comparaisons entre les paramètres d'évaluation Pe et Ps er produit en sortie des résultats d'évaluation R. Les paramètres Pe sont transmis vers la sortie par le système TRANS ou indépendamment de celui-ci. Alternativement, il est possible de faire "remonter" les paramètres Ps vers l'entrée afin qu'ils y soient comparés avec les paramètres Pe. On remarquera également qu'il est possible d'utiliser les coefficients TCD disponibles au niveau du codeur CCD et/ou du décodeur DEC, mais alors l'évaluation ne tiendra pas compte des défauts du codeur CCD et/ou du décodeur DEC. La transtormee en cosinus discrète bicuirectlonnelle TCD fait correspondre à un bloc de N x N pixels (en général N = 8) de l'espace pixel un bloc de N x N coefficients du domaine transformé. Ces derniers expriment la représentation fréquentielle de l'amplitude de la luminance des pixels. La transformée TCD présente des coefficients E(u,v) dans lesquels - u,v = index de la fréquence spatiale horizontale et verricale, - x, = index de position horizontale et -,erticale du pixel dans le bloc, - f(x,y) représente la luminance du pixel au point (x,y) dans l'image, - F(u,v) représente le coefficient transformé au point (u,v) du domaine transformé, - avec u,v,x,y variant de O à N-l (N = 8). Les paramètres exploités pour représenter les dégradations contenues dans les images numériques codées notamment en MPEG2 utilisent cette transformation. Une information spatiale est fournie à partir du domaine transformé par le calcul de l'activité spatiale globale SA de l'image. Ce paramètre informe précisément sur les détails dans le signal image (It) indexé à l'instant t. (It) contient (HxW) blocs de pixels. La taille du bloc f,(x,y) est 8x8 pixels (voir fig.l). Dans un premier temps la transformée TCD par blocs est appliquée à l'image entièrement, ou partiellement (sur une zone d'intérêt ou fenêtre d'analyse) . Le premier facteur représentatif du contenu de l'image, l'activité spatiale, est obtenu à partir de l'image tranformée TCD (il). Un bloc image très actif spatialement a une représentation conséquente dans le domaine transformé. Ses coefficients ont une amplitude et une distribution spectrale importantes. Deux paramètres sont introduits, l'activité spatiale de bloc (bsa) et l'activité spatiale d'image SA (It) . Ils représentent respectivement l'activité spatiale intra-bloc et intra-image. Le contenu analysé est celui de chaque bloc de pixels fnm(x,y)de position à n,m dans l'image (I). L'indicateur (bsan,rn) pour le bloc n,m tient compte des composantes alternatives du bloc de coefficients transformés Fnm(u,v). On a SA(It) est l'expression du contenu global de l'image (ou de la zone d'anaiyse). 1 est calculé par une movenne des activités sDatials (bsa r.,rn des blocs concernes. où (H x W) est le nombre de blocs de pixels de l'image ou de la fenêtre d'analyse. L'indicateur d'activité spatiale global SA(It) permet de qualifier la richesse de l'image. Ce paramètre est affecté par le codage à réduction de débit et par les erreurs de transmission. La comparaison du même paramètre calculé à l'entrée et à la sortie d'un système vidéo numérique contribue à l'identificatIon des dégradations générées par ce système. Ce paramètre est représentatif des dégradations générées par le système, mais sans l'inconvénient d'une compensation entre les défauts. En effet, d'une part le calcul réalisé par blocs élimine les effets de bloc et d'autre part les pertes d'information et les faux contours se traduisent dans la tranformée par des effets de perte des hautes fréquences, et ces deux phénomènes agissent donc dans le même sens. L'évolution temporelle de SA(It) représente les événements apparus dans la séquence vidéo. Par ailleurs, ce paramètre est affecté par le codage à réduction de débit et par les erreurs de transmisslon. La comparaison du même paramètre calculé à l'entrée et à la sortie d'un système vidéo numérique est effectuée sur une image It indexée à l'instant t. Elle contribue à l'identification des dégradations générées par ce système. La formulation mathématLque utilisée pour la comparaison est Le paramètre définitif pour la mesure de la dégradation liée à la perte d'informations est le paramètre MLR MLRM = MaximumM[log tatio(l )] Cet indicateur est calculé sur une séquence de M images. Cette option permet d'avoir une information réduite concernant une période de temps. Elle permet d'avoir une note de qualité sur : T = M x 1/25 seconde, ou plus généralement T = M/NI, NI désignant le nombre d'images par seconde. La procédure de calcul de la perte d'information se présente par exemple comme suit a) sélection d'une image à l'entrée du système, b) calcul de SA(It)e pour cette image, c) indexatIon de cette image avec sa date et son paramètre associé Sêe, d) transmission de paramètres de référence vers la sortie du système, f) calcul de SA(It)s pour cette image, g) calcul du log ratio(I-) pour chaque image It, h) calcul de la note de qualité LITER pour chaque séquence de M images. On notera que le procédé pourrait également être mis en oeuvre en faisant "remonter" vers l'entrée les valeurs de SA calculées pour chaque image de sortie. Un avantage du procédé est que chaque image à l'entrée et à la sortie est caractérisé par la valeur du paramètre SA, ce qui signifie qu'il n'est besoin de transmettre en tant que paramètre Pe qu'une seule donnée par image pour caractériser les dégradations (autres que l'effet de bloc) par opposition à des procédés de l'art antérieur qui nécessitent de transmettre un grand nombre de données permettant d'effectuer une comparaison par différence pixel à pixel entre les images référence et dégradée. ;Q comparaison effectuée à partir du facteur MLR traduit une mesure Qe la qualité des images finales et par conséquent informe sur les performances du système de codage mis en oeuvre lors de la transmission. Afin de valider la pertinence du facteur MLR, nous avons confronté les résultats obtenus aux notes obtenues par des panels de personnes représentatives lors de tests subjectifs. Le fort coefficient de corrélation obtenu entre l'évaluation objective et subjective nous démontre l'efficacité d'un tel paramètre. Dans les séquences vidéo codées à faible débit, la distorsion par effet de bloc se manifeste par des discontinuités du contenu de l'image aux limites des blocs adjacents. Cet artefact est visible sur les zones de faible ou de grande activité spatiale et aussi sur les zones de forte ou de faible luminosité. Il est dû à un codage inter-blocs indépendant et à une quantification grossière du contenu de chaque bloc. Le coefficient continu F(O,O) de la transformée TCD représente la moyenne des valeurs de luminance d'un bloc. Le codage des coefficients continus de deux blocs, adjacents spatialement, à des valeurs sensiblement différentes est une des causes de cette dégradation. Les autres coefficients transformés, qui représentent les composantes alternatives, traduisent les différentes évolutions contenues dans e bloc. Le codage et la quantification grossière de ces coefficients génère également l'effet de bloc dans les zones des contours. Des erreurs dans le flux numérique influent sur la compensation de mouvement. Ceci se traduit par un effet au niveau des blocs ou des macro-blocs. Selon l'inventlon, un paramètre de mesure de l'effet de bloc met en oeuvre une comparaison de deux paramètres calculés de chaque côté du système vidéo. L'indicateur proposé permet d'évaluer l'importance de l'information visuelle erronée apportée par la grille structurée de blocs 8x8 visible dans les images en sortie de chaîne. Les bordures des blocs 8x8 ou les discontinuités de contours sont considérées comme une information perceptible et additionnelle dans l'image. Sachant que le codage est réalisé de manière indépendante pour les blocs adjacents, - les pertes d'information internes aux blocs génèrent des discontinuités de contenu, voire même l'apparition de nouveaux contours sur les frontières. Un procédé adapté à cette dégradation selon l'invention repose sur le calcul de deux paramètres l'activité spatiale globale SA (voir supra) et l'activité spatiale glo paramètre, est donc introduit, Sa,,,,, dans l'équation suivante. SAd (t) SAratio(t) = 100 ) SAratia n'est pas seulement sensible aux contours ou discontinuités aux bords des blocs mais il l'est également au contenu d'origine des blocs. Il est par conséquent nécessaire, pour mesurer spécifiquement la dégradation introduite par- un système vidéo, d'en extraire l'information apportée par le contenu original de l'image. L'indicateur définitif pour la mesure et la détection de l'effet de bloc introduit par un système vidéo est calculé à partir des valeurs SA atia obtenues à l'entrée et à la sortie du-dit système. Ces valeurs sont considérées respectivement comme valeur du signal référence à l'entrée (SAratioe) et valeur du signal dégradé à la sortie (SAratios). L'indicateur de l'effet de bloc BM est défini par l'équation suivante. BMsA = 100. a SAratio s Le procédé de détection de l'effet de bloc peut être mis en oeuvre par exemple comme suit a) sélection d'une image à l'entrée du système, b) calcul de SA(It)e, SA(It)e et SAratno(It)e pour cette image, c) indexation de cette image avec sa date et son paramètre associé, d) transmission de données de référence vers la sortie du système, e) identification de la même image indexée à la sortie du système, f) calcul de SA(It)s, SAf )s et SAratia(It)5 pour la même image à la sortie du système, g) calcul du BMSA pour chaque image Ie, h) calcul de la moyenne de EMSA pour chaque séquence de M' images (avec par exemple M' = M). Cette note qualifie la séquence pour la dégradation en question. La figure 3 démontre l'importance décroissante de BMSA pour 4, 6 et respectivement 9Mbits/s. Cette remarque est valable spécialement dans la dernière partie de la séquence vidéo. En effet, pendant cette période le contenu spatial et temporel de la séquence choisie évolue brutalement. Ceci provoque une baisse de la qualité des images finales due à une réduction de débit. Le codage par quantification des coefficients TCD provoque l'apparition de l'effet de bloc, et ce d'autant plus que le codage est grossier. Les dégradations dans les séquences d'images numériques codées ont des formes spatiales ou temporelles. En effet, des erreurs de codage ou de transmission peuvent affecter le contenu spatial contours, objets et structures contenus dans le signal. Cependant, elles peuvent également entacher 'évolution temporelle de la séquence vidéo : saccades, gel d'images ou déplacements erronés d'objets. Afin dtident Cler es dégradations temporelles qui altèrent les séquences vidéo, nous avons défini un paramètre supplémentaire appelé "Activité Temporelle" TA. Il traduit les évolutions du contenu des séquences d'image dans le temps. Il est sensible aux différents événements qui peuvent survenir et affecter le signal. L'indicateur spatial SA Le signal mono-dimensionnel SA(I ) d'activité spatiale globale de l'image est proportionnel à la richesse du contenu. Dans le cas mono-dimensionnel, la TCD est une transformation qui fait correspondre a un ensemble de M points ;(x) une série de M fréquences. ;Ile se définit par Sur toute fenêtre de M éc.antillons temporels du signal SA(It) est appliquée une transformée en cosinus discrète TCD. Cette nouvelle représentatIon traduit les fréquences temporelles du signal et par conséquent l'activité de la séquence d'image. La première raie F(O) est la composante continue, les autres raies F(u) avec u#o sont les composantes alternatives. La composante continue du spectre renseigne sur a valeur moyenne de la richesse du contenu sur M images. Les fréquences alternatives informent sur les var;aiions temporelles du contenu. Seule l'amplitude des raies est prise en considération. Le spectre TCDT (voir fig.4c) du signal SA (voir fig.4b) informe sur les évolutions temporelles de la séquence vidéo et par conséquent traduit toute dégradation modifiant ces évolutions. (O) par définition représente la moyenne du signal SA(I-) pendant les M images de sortie considérées dans la séquence. Les amplitudes des composantes alternatives F(U)uwo traduisent la manière avec laquelle change le contenu de la vidéo variations rapides (clips vidéo), modifications lentes ("travelling"), plans fixes, etc... Afin d'exploiter les propriétés de l'analyse dans le domaine transformé, un paramètre TA a été défini. Il est calculé pour un groupe de M Images. TA indique l'amplitude des fluctuations de l'activité SA dans la séquence. Il traduit tout événement temporel de la séquence et par conséquent : les coupures et les gels d'images. TA est comparé à un seuil donné S1 établi empiriquement, en dessous duquel on considère qu'on se trouve en présence d'une coupure ou d'un gel d'images. Le procédé de détection des coupures et des gels d'image peut être mis en oeuvre par exemple comme suit a) sélection d'un groupe d'images, b) calcul de SA(It) pour chaque image, c) calcul de la TCDT et de TA, d) - si TA - si TA S2, alors le groupe d'images est une succession de la dernière image non erronée gelée par le décodeur. - si TA > sl, retour à a. Le phénomène de gel d'images est une dégradation caractéristique qui apparaît lors d'erreurs sur la chaîne complète de distribution du signal de télévision. Le paramètre d'activité temporelle TA est applicable dans ce cas de figure. En effet, il détecte cet événement qui est lié à la stratégie du décodeur. Le gel d'images prend deux formes 1. restitution des dernières données images non corrompues, 2. visualisation d'une image uniforme (souvent noire) , et ceci tant que les données reçues sont inexploitables. La famille de codage MPEG introduit la notion de groupe d'images ("Group of Pictures GOP"). Cette structure de la séquence est liée à une syntaxe définissant les algorithmes prédictifs à réduction de débit qui sont appliqués à chaque image (images codées sans prédiction : Intra - I, images prédictives unidirectionnelles P, images prédictives Bidirectionelles B). La particularité de chaque image influe sur l'aspect de l'impact des différentes erreurs qui l'affectent. L'indicateur défini ici traduit les événements temporels qui surviennent durant la séquence vidéo. Il permet par conséquent d'identifier et de comptabiliser les ruptures temporelles des séquences. Par la même occasion il représente un paramètre de qualification de la qualité du service. L'invention ne se limite pas aux modes de réalisation décrits. En particulier, l'effet de bloc pourrait être évalué de manière intrinsèque sur l'image de sortie. Le procédé est alors caractérisé en ce qu'il comporte les étapes suivantes di) déterminer l'activité spatiale globale SAd2 de l'image de sortie dans une fenêtre d'analyse décalée présentant des blocs de points image ou pixels qui sont décalés d'au moins un pixel dans la direction des lignes de l'image et/ou dans une direction perpendiculaire à la direction des lignes de l'image par rapport auxdits blocs de pixels de ladite fenêtre d'analyse, cette détermination mettant en oeuvre les sous-étapes suivantes: i) appliquer ladite transformée par blocs sur chaque bloc décalé (n,m) dudit ensemble de blocs de pixels de la fenêtre d'analyse décalée pour déterminer ies coefficients transformés Fdn,m (;,j), ii) déterminer à partir des coefficients transformés Fdn,m (i,j) des blocs décalés, l'activité spatiale bsad de chaque bloc dudit ensemble de blocs décalés, iii) déterminer l'activité spatiale globale SAd2 de l'ensemble des blocs décalés constituant la fenêtre d'analyse décalée, ei) comparer l'activité spatiale globale SAd2 de l'image de sortie dans la fenêtre d'analyse décalée et l'activité spatiale globale SA2 de l'image de sortie dans la fenêtre d'analyse. L'invention peut être mise en oeuvre à l'aide de tranformées par blocs linéaires et orthogonales telles que la tranformée en cosinus discrète TCD, la transformée en sinus discrète, la transformée de Fourrier, la tranformée de Haar, la transformée d'Hadamar, la transformée de Slant ou bien encore une transformée en ondelettes (DWT) dont les coefficients ont été réorganisés pour obtenir les blocs transformés. REVEND I CAT IONS 1. Procédé d'évaluation de la dégradation d'une image vidéo codée par blocs de points images ou pixels, cette dégradation étant générée par un système de codage et/ou de stockage et/ou de transmission, produisant une Image de sortie à partir d'une image d'entrée, caractérise en ce qu'il comporte les étapes suivantes a) sélectionner - une image d'entrée et déterminer une activité spatiale (SA) de l'image d'entrée dans une fenêtre d'analyse représentant au moins une partie de l'image qui présente un ensemble de dits blocs de points images ou pixels, cette détermination mettant en oeuvre les sous-étapes suivantes i) déterminer pour chaque bloc (n,m) du dit ensemble de blocs de pixels, les coefficients transformés Fn,m (i,y) par une transformée par blocs selon les dits blocs il) déterminer à partir des coefficients transformés Fn,m (i,y) l'activité spatiale bsa de chaque bloc dudit ensemble de blocs de pixels iii) déterminer à partir de l'activité spatiale de bsa de chaque bloc l'activité spatiale globale (SAl) de l'ensemble des blocs de pixels constituant la fenêtre d'analyse b) sélectionner l'image de sortie correspondant à l'image d'entrée et déterminer ladite activité spatiale globale (SA2) de l'image de sortie dans ladite fenêtre d'analyse, en mettant en oeuvre pour l'image de sortie les sous-étapes a i) à a iii) c) comparer l'activité spatiale globale (SA2) de l'image de sortie dans la fenêtre d'analyse et l'activité spatiale globale (SA.) de l'image d'entrée dans la fenêtre d'analyse. 2. Procédé selon la revendication caractérisé en ce que l'activité spatiale osa d'un bloc (n,m) est déterminée par a formule suivante où k(i,j) est un coefficient constant de la pondération pour la composante i,j utilisée et p une constante. 3. Procédé selon la revendication 2 caractérisé en ce que l'activité spatiale SA de l'ensemble des blocs peut être alors déterminée par la formule suivante (HxW) représentant le nombre de blocs de pixels dans la fenêtre d'analyse. 4. Procédé selon la revendication 3 caractérisé en ce que ladite comparaison (étape c) est effectuée à l'aide d'un paramètre LR défini de la manière suivante LR = gj [fi(SA1/SA2) avec i,j E {1,9} avec fa(x,yj=x-y ou f(x, y)=x/y et g1(x)=100|x| ou a2(x)=i00 |log(|x|)| 5. Procédé selon la revendication 4 caractérisé en ce qu'il met en oeuvre le calcul du paramètre LR sur une séquence de M images d'entrée et de M images de sorties correspondants, et en ce qu'il met en oeuvre le calcul d'un paramètre d'évaluation ou note de qualité MLR défini de la manière suivante MLR = Maximum, (LR) 6. Procédé selon une des revenaications précédentes caractérisé en ce qu'il comporte la détermination d'un indicateur d'activité temporelle TA pour un groupe de M images, déterminé de la manière suivante F,, (u), pour u variant de O à M-l, désignant les M coefficients d'une transformée par blocs appliquée à M échantillons temporels de l'activité spatiale SA. 7. Procédé selon la revendication 6 caractérisé en ce qu'il comporte une étape de comparaison de l'indicateur d'activité temporelle TA à un seuil donné S1, de manière à déterminer la présence d'une coupure ou d'un gel d'image. 8. Procédé selon une des revendications précédentes caractérisé en ce qu'il comporte, en vue de la détermination de manière intrinsèque de l'effet de bloc, les étapes suivantes di) déterminer l'activité spatiale globale SAd2 de l'image de sortie dans une fenêtre d'analyse décalée présentant des blocs de points image ou pixels qui sont décalés d'au moins un pixel dans la direction des lignes de l'image et/ou dans une direction perpendiculaire à la direction des lignes de l'image par rapport auxdits blocs de pixels de ladite fenêtre d'analyse, cette détermination mettant en oeuvre les sous-étapes suivantes i) appliquer ladite transformée par blocs sur chaque bloc décalé (n,m) dudit ensemble de blocs de pixels de la fenêtre d'analyse décalée pour déterminer les coeeficients transformés Fdn,m (i,j), ii) déterminer à partir des coefficients transformés n,m , m (i,j) des blocs décalés, l'activité spatiale bsad de chaque bloc dudit ensemble de blocs décalés, iii) déterminer l'activité spatiale globale SAd2 de l'ensemble des blocs décalés constituant la fenêtre d'analyse décalée, ei) comparer l'activité spatiale globale SAd2 de l'image de sortie dans la fenêtre d'analyse décalée et l'activité spatiale globale SA de l'image de sortie dans la fenêtre d'analyse. 9. Procédé selon une des revendications précédentes caractérisé en ce qu'il comporte, en vue de la détermination de l'effet de bloc, les étapes suivantes d2) déterminer l'activité spatiale globale SAdi de l'image d'entrée dans une fenêtre d'analyse décalée présentant des blocs de points image à pixels qui sont décalés d'au moins un pixel dans la direction des lignes de l'image et/ou dans une direction perpendiculaire à la direction des lignes de l'image par rapport auxdits blocs de pixels de ladite fenêtre d'analyse, cette détermination mettant en oeuvre les sous-étapes suivantes i) appliquer une transformée par blocs sur chaque bloc décalé (n,m) dudit ensemble de blocs de pixels de la fenêtre d'analyse décalée pour déterminer les coeeficients transformés Fdn,m (i,y), ii) déterminer à partir des coefficients transformés Fdn,m (i,y) des blocs décalés, l'activité spatiale bsad de chaque bloc dudit ensemble de blocs décalés, iii) déterminer l'activité spatiale globale SAdl de l'ensemble des blocs décalés constituant la fenêtre d'analyse décalée, e) déterminer 1' activité spatiale SAd2 de l'image de sortie dans ladite fenêtre d'analyse décalée en mettant en oeuvre les sous-étapes d i) à d iii) f) comparer d'une part l'activité spatiale SAd2 de l'image de sortie dans la fenêtre d'analyse décalée et l'activité spatiale SA2 de l'image de sortie dans la fenêtre d'analyse, pour évaluer l'effet de bloc dans l'image de sortie et d'autre part l'activité spatiale SAdl de l'image d'entrée dans la fenêtre d'analyse décalée et l'activité spatiale SAl de l'image d'entrée dans la fenêtre d'analyse pour évaluer l'effet de bloc dans l'image d'entrée. 10. Procédé selon la revendication 9 caractérisé en ce que l'étape f met en oeuvre la détermination d'un indicateur BM d'effet de bloc, selon la formule suivante SAd2 / SA2 BM = --------- x 100 SAdl / SA1 11. Procédé selon une des revendications précédentes caractérisé en ce que ladite transformée par blocs est une transformée linéaire et orthogonale, notamment une transformée en cosinus discrète TCD ou une transformée en sinus discrète ou une transformée de Fourier ou de Haar ou de Hadamar ou de Slant. 12. Procédé selon une des revendications 1 à 10 caractérisé en ce que ladite transformée par blocs est une transformée en ondelettes (DWT) dont les coefficients ont été réorganisés pour obtenir des blocs transformés.