FR 2673795 A1 19920911 FR 9102605 A 19910305 Procédé et dispositif d'interpolation adapté de signaux sous-bandes pour le codage compatible inter-image de signaux de télévision La présente invention concerne un procédé et un dispositif d'interpolation adapté de signaux sous-bandes pour le codage compatible inter-image de signaux de télévision. Elle s'applique notamment aux systèmes de transmission et de réception de signaux vidéo numériques. Dans les systèmes de transmission d'images de télévision les multiplications des formats de résolution employés pour la télévision classique entrelacée, le vidéotéléphone, la télévision progressive, et l'avènement des systèmes de télévision à haute résolution entrelacée ou progressive rendent nécessaire l'élaboration de systèmes de codage compatible, quelque soit les sens ascendant ou descendant des résolutions. Dans la compatibilité ascendante un récepteur travaillant dans un format vidéo déterminé doit pouvoir recevoir et visualiser un signal transmis suivant un format plus élevé, c'est-à-dire un format selon lequel la résolution du signal est plus grande. C'est le cas par exemple, d'un signal de télévision haute définition progressive connu sous l'abréviation HDTV qui possède un format plus élevé qu'un signal de télévision progressif connu sous l'abréviation EDP et du signal EDP qui possède lui-même un format plus élevé que celui d'un signal de télévision entrelacé classique. Inversement dans le cas d'une compatibilité descendante un récepteur doit pouvoir recevoir et visualiser un signal de résolution moindre que son format de travail. Naturellement les diversités des résolutions entraînent une diversité du dispositif de codage et de décodage à mettre en oeuvre et compliquent considérablement la gestion des communications. Une solution à ces problèmes consiste à effectuer au niveau des codeurs d'émission un découpage de chaque image à transmettre en sous-bandes suivant une rnême structure arbores cente par filtrage et décimation de la bande de fréquence spatiale des signaux à transmettre quelque soit leur famille d'appartenance, de façon à coder indépendamment les signaux à chaque sous-bande avant de les transmettre par multiplexage sur le canal de transmission. En réception les décodeurs effectuent le démultiplexage des signaux codés reçus dans chaque sous-bande en décodant les signaux reçus relatifs à chaque sous-bande en fonction de la résolution des signaux transmis et de la résolution propre des récepteurs.Malheureusement ce procédé perd rapidement de son efficacité lorsque le débit sur le canal de transmission diminue et il apparaît que pour des débits inférieurs à 1,4 bits par pixel la qualité visuelle des images observées n'est plus suffisante et il faut alors introduire un codage inter-image supplémentaire tenant compte de la corrélation temporelle naturelle des séquences d'images. Comme décrit par exemple dans la demande de certificat d'addition nO 90 04816 déposée le 13 avriI 1990 au nom de la Demanderesse, le codage inter-image supplémentaire consiste à effectuer un codage différentiel avec compensation de mouvement au niveau de chaque sous-bande.Cette opération de compensation de mouvement doit être effectuée au niveau des sous-bandes et non au niveau de l'image pleine définition (c'est-à-dire l'image non décomposée) de manière à assurer un codage compatible des images et à éviter tout phénomène de dérive au niveau d'un décodeur compatible. Dans ces conditions, l'estimation de mouvement est effectuée sur les images non décomposées suivant une technique connue sous la désignation "block matching" : l'image est découpée en blocs et un vecteur mouvement par bloc est estimé par un calcul de corrélation. Cette technique permet d'obtenir des vecteurs mouvements entiers, c'est-à-dire dont les composantes sont entières. On peut généraliser aisément pour estimer des vecteurs à la précision demi-pixel en interpolant d'abord les images par exemple avec un filtre bilinéaire.Ces vecteurs mouvements sont ensuite divisés par le taux de sous échantillonnage de la sous-bande considérée pour servir de base à la compensation de mouvement : ce procédé consiste, au lieu de retirer le pixel en vis à vis codé de la sous-bande précédente de la valeur du pixel à coder, à retirer la valeur d'un pixel voisin dont les coordonnées sont calculées en fonction de ces vecteurs mouvements. Ainsi, lorsque les composantes du vecteur mouvement sont entières, le problème se réduit uniquement à un problème d'adressage mémoire (il suffit d'aller chercher un pixel existant). Par contre, lorsque le vecteur mouvement n'est pas entier la valeur du pixel recherché doit être interpolée. La demande de certificat d'addition précitée 90 04516 propose une simple interpolation bilinéaire utilisant les points voisins dans la sous-bande.Néanmoins, cette solution ne suffit pas car elle conduit à une erreur de prédiction beaucoup trop forte dans certaines bandes, et donc réduit les performances globales du système de codage. Le but de l'invention est de pallier les inconvénients précités en proposant une interpolation plus efficace, tout en conservant les propriétés de compatibilité. A cet effet, l'invention a pour objet un procédé d'interpolation adapté de signaux sous-bandes pour le codage et le décodage compatible inter-image de signaux de télévision de résolutions différentes, les signaux étant transmis sur un canal de transmission entre au moins un codeur d d'émission et un décodeur de réception relié à un récepteur, dans lequel le codage consiste à découper en sous-bandes, suivant une même structure arborescente par filtrage et décimation, l'image pleine bande à transmettre et coder indépendamment les signaux sous bandes obtenus et dans lequel le décodage consiste à décoder les signaux reçus de chaque sous-bande en fonction de leur résolution et de celle propre au récepteur, caractérisé en ce qu'il consiste à effectuer au niveau du codeur d'émission un codage inter-image compensé en mouvement dans chaque sous-bande, à partir des vecteurs mouvements estimés sur l'image non décomposée et remis à l'échelle en fonction du taux de sous échantillonnage de la bande considérée, et à effectuer dans chaque sous bande, lorsque le vecteur mouvement (déduit de l'estimation sur l'image non décomposée) ne correspond pas à un nombre entier de pixels, un calcul d'interpolation du pixel correspondant en tenant compte non seulement du signal contenu dans la sous bande elle-même mais également des signaux des autres sous-bandes. L'invention a également pour objet un dispositif pour la mise en oeuvre du procédé précité. Le procédé selon l'invention a pour principal avantage qu'il permet de mieux tenir compte de la corrélation temporelle des signaux sous-bandes tout en assurant la conservation des propriétés de compatibilité du système de codage-décodage. Grâce à l'invention les calculs d'interpolation sont beaucoup plus efficaces et conduisent à des erreurs de prédiction dans les bandes beaucoup plus faibles comparées à celles obtenues par simple interpolation bilinéaire. La transmission est donc beaucoup moins coûteuse ce qui accroît les performances. D'autres caractéristiques et avantages de l'invention apparaîtront ci-après à l'aide de la description qui suit faite en regard des dessins annexés qui représentent - La figure 1 une représentation de résolutions de différents standards de télévision représentés dans un espace spatio-temporel de Fourier à deux dimensions. - La figure 2 un mode d'organisation de cellules de décomposition en sous-bandes mises en oeuvre par l'invention. - La figure 3 une localisation de quatre sous-bandes composées au moyen d'une cellule du type de celle qui est représentée à la figure 2 dans un espace spatio-temporel de Fourier à deux dimensions. - La figure 4 un arbre minimal de décomposition en trois formats d'une image de télévision dans un espace spatio-temporel à deux dimensions de Fourier. - La figure 5 une localisation selon l'invention des parties du signal à traiter de manière indépendante dans l'arbre de décomposition de la figure 4. - La figure 6 une structure d'arbre optimal à 16 bandes déduit de l'arbre de la figure 4 permettant d'obtenir une décorrélation suffisante et une quantification scalaire efficace des signaux. - Les figures 7a et 7b un mode de réalisation d'un codeur sous-bandes mis en oeuvre par l'invention. - La figure 8 un mode de décomposition monodimensionnelle servant de base à la présentation d'un calcul de filtre d'interpolation sous-bande. - La figure 9 un mode de réalisation d'une cellule de décomposition équivalente au mode de décomposition de Ia figure 8. - La figure 10 des schémas d'opérateurs fonctionnels pour montrer les équivalences fonctionnelles existant entre les opérateurs des modes de réalisation des figures 8 et 9. - La figure 11 un mode de réalisation d'un interpolateur sous-bandes selon l'invention. - La figure 12 une variante de réalisation de l'interpolateur de la figure 11. - la figure 13 une décomposition sur un axe fréquentiel d'une image en 4 sous-bandes - Les figures 14a à 14d des graphes représentant les contributions des filtres aux bandes 1, 2, 3, 4. - La figure 15 une modification du mode de mise en oeuvre de la figure 11. - La figure 16 un mode de réalisation d'un interpolateur de points dans une sous-bande convenant à un estimateur de mouvement au demi-pixel. - La figure 17 une description des interdictions d'utilisation d'informations dans certaines bandes pour l'interpolation d'une bande donnée pour respecter les contraintes de compatibilité. Dans la mise en oeuvre de l'invention, le spectre du signal de télévision à coder est décomposé en sous-bandes en respectant la hiérarchie des résolutions imposées par les diffé rents standards de télévision en usage actuellement, de la manière représentée à la figure 1 où les différents niveaux de résolution d'une hiérarchie d'image définie par différents standards sont représentés dans un espace de Fourier à deux dimensions dans lequel sont consignées les fréquences spatiales de l'image mesurées en cycle par largeur et hauteur d'image. Sur la figure 1 seuls les standards HDP, EDP et VT désignant respectivement le format pour la haute définition progressive définie par des images de 1250 lignes transmises à la cadence de 50 images par seconde (1250/50/1:1), EDP désignant le format de définition étendue progressive (625/50/I: 1) définie par des images de 625 lignes transmises à la cadence de 50 images par seconde, et VT désignant le format de l'image à 312 lignes progressives du vidéo téléphone (312/50/1:1) ont été représentés et il n'est considéré dans ce qui suit que les signaux vidéo aux formats progressifs des standards HDP, EDP et VT, l'application du système et du procédé selon l'invention aux signaux entrelacés pouvant être réalisée en utilisant les signaux progressifs équivalents de la façon dont cela est décrit par exemple dans la demande de brevet principal 89 15175 à laquelle est rattachée la demande de certificat d'addition précitée. Dans l'exemple de la figure 1, pour effectuer un codage-décodage compatible entre les trois standards représentés, chaque récepteur doit être capable de décoder et de visualiser, quelque soit son niveau propre de résolution, le signal qu'il reçoit et qui a un format quelconque parmi les 3 formats VT, EDP ou HDP. Ceci ramène le système de codage compatible à un système de codage multirésolution. Pour effectuer un codage compatible inter-image de signaux vidéo, le procédé selon l'invention effectue une décomposition hiérarchique et compatible en sous-bandes des images à transmettre. Les signaux sous-bandes obtenus sont quantifiés linéairement puis codés à l'aide de codes à longueurs variables avec une description par plages de zéros.Le procédé selon l'invention emploie des cellules de décomposition de base, du type de celle qui est représentée à la figure 2 avec le sys tème de reconstruction associé. Une telle cellule est formée par la mise en cascade suivant un arbre prédéfini de cellules élémentaires bidimensionnelles découpant les spectres du signal en quatre bandes, comme décrit à la figure 3. Chaque cellule élémentaire bidimensionnelle est elle-même formée par la mise en cascade de trois cellules élémentaires de décomposition monodimensionnelle selon un arbre à deux étages.Le premier étage, formé par les filtres 1 et 2 de fonctions de transfert H0(Z), et H1(Z), filtre les colonnes de l'image et le second étage, formé par les filtres 3 à 6 de fonctions de transfert H0(Z) et H1(Z), filtre les lignes des deux images sous-bandes obtenues du premier étage formé par les filtres I et 2. La mise en cascade de telles cellules permet de découper le spectre spatial de l'image selon une structure d'arbre à quatre branches délimitant chacune quatre bandes de fréquence telles que représentées à la figure 3.L'image qui est ainsi décomposée peut être reconstruite dans un décodeur de réception au moyen de cellules d'interpolation formées par la mise en cascade de trois cellules d'interpolation monodimensionnelles formées des filtres 7 à 12 selon un arbre également à deux étages. La reconstruction presque parfaite de l'image est obtenue en utilisant des filtres connus sous l'abréviation CQF du terme anglo-saxon "conjugate quadrature filter" ou encore des filtres connus sous l'abréviation QMF de "quadrature mirror filter". En utilisant ce principe de décomposition et de reconstruction, l'arbre minimal qui autorise un codage compatible entre les trois formats VT, EDP et HDP du spectre représenté à la figure 1 est un arbre à sept bandes tel que représenté dans l'espace fréquentiel de Fourier de la figure 4. Dans cette décomposition l'exigence de compatibilité se traduit simplement par un traitement indépendant des trois parties du spectre, de la partie basse fréquence VT (bande 4) du spectre EDP, de sa partie complémentaire CEDP (bandes 5, 6, 7) pour former le spectre EDP et de la partie complémentaire CHDP (bandes 1, 2, 3) du spectre EDP pour former le spectre total HDP que cet arbre engendre de la façon représentée à la figure 5. En particulier cette découpe suggère que toute prise en compte du mouvement de l'image doit être faite de façon indépendante dans chacune de ces trois parties. Cependant en pratique il apparaît qu'une découpe au moyen d'un arbre à sept bandes ne conduit pas à une décorrélation suffisante des signaux pour permettre une quantification scalaire efficace, et que des découpes supplémentaires sont nécessaires. Un bon compromis peut cependant être trouvé au moyen d'un arbre à seize bandes de la façon représentée à la figure 6. Comme décrit précédemment, pour atteindre par exemple les débits de l'ordre de 0,8 bits par seconde il est nécessaire de prendre en compte la corrélation temporelle existant dans les séquences d'image, sans détruire les propriétés de compatibilité du système, c'est-à-dire en respectant les traitements indépendants des parties VT, CEDP et C11DP du spectre de la figure 5, pour assurer un codage compatible, c > est-à-dire éviter tout phénomène de dérive au niveau des décodeurs compatibles. Pour obtenir ce résultat la compensation du mouvement et le calcul de l'erreur de prédiction sont effectués au niveau des sous-bandes ou tout du moins au niveau des parties VT, CEDP et CHDP. Cependant dans ces traitements il faut prendre garde de ne pas faire apparaître un phénomène de dérive sur les images décodées. En effet, en prenant par exemple le cas d'un décodeur EDP, comme celui-ci ne peut utiliser que l'information contenue dans les bandes de basses fréquences, partie VT et partie complémentaire CEDP du signal HDP de résolution supérieure transmis, si au niveau du codeur la compensation de mouvement est effectuée en utilisant toute l'information du signal HDP (parties VT, CEDP, CHDP), le décodeur EDP n'aura plus accès à l'information C11DP et dans ce cas ne pourra plus effectuer un traitement similaire à celui qui a été effectué au niveau du codeur. De ce fait un phénomène de dérive peut apparaître sur les images décodées.Une estimation de mouvement n'affectant pas les propriétés de compatibilité du système peut être faite soit au niveau de l'image, entière, soit au niveau des sous-ban des. Comme le mouvement physique est le même dans chaque image sous-bande, il apparaît préférable d'utiliser un seul ensemble de vecteurs mouvement plutôt que d'estimer le mouvement dans chaque bande, ceci permet d'éviter d'une part une certaine redondance entre les vecteurs mouvement de chaque bande et d'autre part, un codage assez compliqué des vecteurs correspondants pour limiter le débit d'informations requis pour leur transmission. Selon l'invention l'estimation de mouvement a aussi lieu sur des images sources et non comme c'est généralement le cas entre l'image source à coder et l'image précédente décodée.Comme décrit précédemment l'algorithme d'estimation de mouvement utilisé est un algorithme de type "block matching" suivant lequel les images sont découpées en blocs de NxN pixels de façon à estimer un vecteur mouvement par bloc par intercorré lation entre deux images successives, dans ce cas le vecteur mouvement correspond au pic de corrélation obtenu. Le vecteur mouvement ainsi obtenu est transmis au décodeur. Un mode de réalisation d'un codeur correspondant est représenté aux figures 7a et 7b. Le codeur représenté à la figure 7a comprend un banc de filtres de décomposition en n sous-bandes 13 du type de ceux décrits précédemment, chaque sous-bande obtenue du filtre de décomposition 13 étant appliquée à un circuit de codage différentiel formé d'un circuit de quanti fication 14. et d'un circuit de déquantification 15.. Chaque I I pixel quantifié d'une sous-bande est ainsi appliqué à l'entrée d'un dispositif de codage à longueur variable 16. dont la sortie est reliée à une entrée d'un multiplexeur 17. La sortie du multiplexeur 17 est couplée à un canal de transmission 18 par l'intermédiaire d'une mémoire tampon 19.Le circuit de la figure 7a comprend également un dispositif d'estimation de mouvement 20 fournissant des vecteurs mouvement V. à une entrée de chaque circuit de déquantification 15i au travers d'un dispositif de mise à l'échelle 21. Il comprend également des bancs de filtre de reconstruction 22. Un mode de réalisation d'un circuit de codage différentiel formé des circuits de quantification 14. et des circuits de déquantification 15. est représenté à la figure 7b. Sur cette figure le bloc de quantification 14. comprend un dispositif de formation de blocs de NxN pixels 23 couplé à un quantificateur 24 par l'intermédiaire d'un dispositif de commutation 25, ainsi qu'un dispositif sommateur 26. Le dispositif de déquantification 15. comprend un déquantificateur 27 couplé à un bloc de compensation de mouvement 28 par l'intermédiaire d'un circuit de commutation 29 et d'une mémoire 30. Les blocs de commutation 25 et 29 effectuent une commutation entre les modes de fonctionnement inter et intra-image.Dans le mode de fonctionnement "inter-image", le dispositif de codage différentiel retire à chaque pixel à coder x(k,l,n) dans la sous-bande à l'instant n la valeur codée du pixel qui lui correspond interpolé ou non en fonction du vecteur mouvement dans la sous-bande précédente à l'instant n-1 pour coder ensuite cette différence. Suivant le schéma de la figure 7b chaque pixel x(k,l,n) d'adresse k,l dans l'image n est appliqué à l'entrée du codeur sur la deuxième entrée du circuit soustracteur 26.Le circuit soustracteur 26 soustrait la valeur du pixel x(k,l,n) dans l'image à l'instant n de la valeur décodée d'un signal reconstitué x(k+vx,l+vy,n-l) d'instant n-l fourni par le bloc de compensation de mouvement 28 et calculé à partir du vecteur mouvement de composante Vx et Vy dans les directions x et y de l'image. Les composantes Vx et Vy sont fournies par l'estimateur de mouvement 20 après avoir été mis à l'échelle par le dispositif de mise à l'échelle 21 en fonction du taux de sous échantillonnage de la sous-bande considérée. L'image est reconstituée à l'entrée de la mémoire 30 grâce à un circuit additionneur 31 qui additionne le signal reconstitué x au signal d(k,l,n) obtenu à la sortie du circuit de dé quantification 27. L'estimation par "block-matching" qui est ainsi obtenue conduit à estimer un vecteur mouvement avec une précision qui est celle du pixel, c'est-à-dire que les coordonnées du vecteur estimé sont de la forme vx=l pixels, vy = k pixels, I et k étant des entiers. Cependant pour accroître les performances globales du codage il est aussi possible et même préférable d'utiliser une estimation au demi pixel, 1 et k étant dans ce cas des multiples de 0,5. Dans ces conditions la compensation du mouvement est effectuée dans une bande donnée sur des blocs correspondants à l'estimation mais remis à l'échelle par homotétie en fonction du facteur de sous échantillonnage de la bande.A titre d'exemple, en se plaçant dans le cas d'une estimation au demi pixel sur des blocs de dimensions NxN définis au niveau de l'image originale et en désignant par v = (vx, vy) le mouvement estimé sur un bloc, vx et vy étant des multiples de 1/2, le bloc correspondant dans la bande 1 de l'image a une dimension de N: 2xN: 2 pixels et il lui est appliqué un vecteur mouvement (V1=Vlx, Vly) tel que Vîx = Vx/2 et V2y = Vy/2, Vîx et Vly étant alors des multiples de 1/4. Pour le cas de la bande 4 les blocs correspondants ont alors une taille de N/4 par N/4 et les vecteurs appliqués ont comme composante Vx/4 et Vy/4. La compensation de mouvement est alors au quart de pixel dans la bande 1 et un huitième de pixel dans la bande 4.Lorsque le vecteur mouvement obtenu dans une bande est entier, (par exemple pour une valeur estimée du mouvement V= (2,0), le mouvement appliqué dans la bande 1 est alors égal à V(1,0) entier) la compensation du mouvement se réduit à un simple décalage car il suffit d'aller rechercher des points déjà existants dans la bande. Par contre, si le mouvement appliqué dans la bande a l'une de ces composantes non entière, il faut dans ce cas aller chercher un point qui n'existe pas physiquement et qu'il faut trouver par interpolation. L'interpolation peut être réalisée au moyen d'une simple interpolation bilinéaire, au 1/2, 1/4, 1/8 ou 1/16 de pixel, des signaux sous-bandes.Mais si cette solution présente I'avantage de la simplicité puisque pour calculer un point manquant il suffit de prendre en compte les points les plus proches qui ltentourent, elle présente toutefois l'inconvé- nient de ne pas prendre en compte la manière dont les sous-ban des sont obtenues et particulièrement le fait qu'elles souffrent d'aliasing. Cependant cette interpolation donne des résultats relativement convenables dans la bande basse fréquence, par exemple dans la bande 13 de la figure 6 ou dans la bande 4 de la figure 5, mais son application sur les autres bandes est loin d'être efficace car l'erreur de prédiction résultante est alors beaucoup plus importante et nuit aux performances globales du système de codage. Un mode d'interpolation beaucoup plus efficace est décrit ci-après à l'aide des schémas de décomposition en sous-bandes monodimensionnelles représentés aux figures 8 à 12 étant bien entendu que l'extension à des traitements bidimensionnels est immédiate du fait qu'un traitement bidimensionnel peut être réalisé de façon séparable en séparant les lignes et les colonnes dans l'image. Le système de décomposition monodimensionnel qui est représenté à la figure 8 conduit à une décomposition de l'image en quatre sous-bandes numérotées 1, 2, 3 et 4 disposées sur un axe fréquentiel de la façon représentée à la figure 13. Dans la configuration de la figure 8 le signal S1(z) de la sous-bande 1 est obtenu, après trois décimations successives par deux d'un signal pleine bande X(z) appliqué à l'entrée du système de décomposition, à la sortie de trois filtres passe-bas 31, 32 et 33 de fonction de transfert Ho(z) reliés en série. Le signal S2 (z) de la sous-bande 2 est obtenu également après trois décimations successives par deux du signal pleine bande X(z), à la sortie d'un filtre passe haut 34 de fonction de transfert H1(z) couplé à la sortie du filtre passe-bas 32 précédemment utilisé pour former la sous-bande 1. Le signal S3 (z) de la sous-bande 3 est obtenu après deux décimations par deux seulement à la sortie d'un filtre passe haut 35 de fonction de transfert Hl(z) couplé à la sortie du filtre 31 déjà utilisé pour l'obtention des sous-bandes 1 et 2.Enfin la sous-bande 4 est obtenue après une seule décimation par deux à la sortie d'un filtre passe haut 36 de fonction de transfert Hl(z) Une recomposition du signal pleine bande à partir des signaux sous-bandes S1(z), S2(z), S3(z) et S4(z) est obtenue au moyen de filtres de reconstruction respectivement passe bas et passe haut de fonction de transfert Go(z) et G1(z) référencés de 37 à 42. Une première reconstruction a lieu par addition des signaux S1 (Z) et S2 (z), après qu'ils aient subi un suréchantillonnage par deux et un filtrage respectivement par deux filtres 37 et 38 de fonction de transfert Gg(z) et G1(z).Le résultat de l'addition obtenu est à nouveau additionné au signal S3 (z) après qu'ils aient subi un suréchantillonnage par deux et un filtrage respectivement par deux filtres 39 et 40 de fonction de transfert Go(z) et G1(z) pour obtenir une deuxième reconstruction. Enfin la reconstruction complète est réalisée en ajoutant au signal reconstitué par les filtres 39 et 40 le signal S4(z), après suréchantillonnage par deux et filtrage respectivement par des filtres 41 et 42 de fonction de transfert Go(z) et G1(z) du résultat de l'addition du signal reconstitué et du signal S4(z). Les filtres de décomposition HoJ Hi et de reconstruction Go et G1 sont sur la figure 8 des filtres à réponse impulsionnelle finie de type CQF ou QMF qui sont les abréviations anglo-saxonnes de "Conjugate Quadrature Filters" et "Quadrature Mirror Filter" pour satisfaire les propriétés de reconstruction parfaite du signal à la sortie des filtres de reconstruction. En utilisant les équivalences des transformées en z de la figure 10 qui reproduisent la propriété H(ZN)= nH(n)Z Nn) où H(n) est la réponse impul )= n11(n)Z , ou sionnelle d'un filtre H, le schéma de la figure 8 peut être transformé comme le montre la figure 9 par un schéma équivalent où les opérateurs de décimation et ceux d'interpolation sont regroupés. Dans la figure 9, Y1(z), Y2 (z), Y3(z) et Y4(z) désignent les signaux sous-bandes S1(z), S2(z), S3(z) et S4(z) définis précédemment avant décimation et Hl(z), H2(z), H3 (z), 11(z), G1(z), G2(z), G3(z) et G4(z) désignent les fonctions de transfert équivalentes à l'association des filtres 31 à 42 de la figure 8 utilisés dans l'obtention des sous-bandes 1 à 4. Sur la figure 9 Hi(z) et Gi(z) sont définies par les produits suivants H1(z)=H0(z).H0(z2) H0(z4 H2(z)=H0(z).H0(Z2).H1(z4) H3(z)=H0(z).H1(z) H4(z)=H1(z) G1(z)=G0(z) .G0(z22) .G0(z44) G2(z)=G0(z).G0(z ).G1(z ) G3(z)=G0(z).G1(z) G4(z)=G1(z). En admettant qu'un vecteur mouvement V, unidimensionnel et entier, estimé au pixel entier, ait été estimé entre un signal x(n) et un signal xp(n) représentant les images dans un espace à une dimension aux instants courant (n) et précédent p, l'estimation parfaite est obtenue lorsque la relation x(n-v)=x (n) est vérifiée. p Dans ce cas x(n-v) peut être pris comme prédiction de la valeur xp (n) . Les vecteurs mouvements appliqués aux bandes I, 2, 3 et 4 ont respectivement pour valeurs v1=v/8, v2=v/8, v3=v/4, v4=v/2. Le problème revient alors à estimer une prédiction des signaux sp1(n), sp2(n), sp3(n), sp4(n). Dans le cas de la bande 1 par exemple, il est clair que si vl est entier, c'est-à-dire multiple de 8, alors la prédiction du signal sp1(n) est A s (n) = sl(n-vl) pi Si la relation x(n-v) = xp(n) est vérifiée alors Apl(n) = spl(n) strictement. Cependant dans le cas où v n'est plus un multiple de 8 il n'est plus possible de prendre le signal s1(n-v) comme prédiction puisque le point correspondant n'existe plus et que ce point doit alors être interpolé. Etant donné qu'en application de la relation précédente la relation y1 (n-v) = Yp1 (n) est toujours vérifiée du fait de l'invariance par translation, l'interpolation idéale du signal s1 (n-Vl) et par suite la prédiction optimale du signal sp1 (n), peuvent être décrites par la relation s1(n-vl) = sp1(n) = y1(n-V) qui met en évidence le fait que l'interpolation optimale des signaux sous-bande S. ne peut être obtenue qu'en recalculant les signaux Y. Ainsi l'interpolation optimale est celle qui pour une bande donnée, tient compte de l'information contenue dans toutes les autres bandes.Pour cela un dispositif de reconstruction de l'image dans une bande donnée peut être réalisé de la manière représentée à la figure 11 à l'aide de filtres d'interpolation en tenant compte en plus des informations relatives à une bande donnée des informations des autres bandes utilisées. Ainsi sur la figure 11 les filtres d'interpolation d'une bande i formés de filtres Gi (z) et Hí(z) travaillent sur les informations relatives à la bande i elle-même et les filtres G:J (z) Hi(z) tiennent compte de la contribution des autres bandes j sur la bande i. L'interpolation est réalisée par la somme des fonctions de transferts G. (z) Hi(z) des filtres d'interpolation. J i Naturellement ce calcul n'est à effectuer qu'au point de la bande où le vecteur de mouvement v. n'est pas entier. En désignant par Kil (z) = la contribution de la bande i sur elle-mêmeet par Kij(z) = G.(z)xH.(z) la contribution de la bande j sur la bande J i i et en notant par N. le taux de sous échantillonnage d'une bande j, il paraît clair que les filtres d'interpolation décrits précédemment doivent vérifier les conditions Kii(O) = 1 K.. (nN1) = 0 pour n différent de 0 K.. (nN.) = 0 pour tout n et tout j différents de i, du iJ J fait qu'aux conditions limites lorsque le vecteur v. est entier le point interpolé doit conserver la même valeur.En utilisant les notations Kij et K.. le schéma de la figure 11 se ramène au iJ schéma de la figure 12 dans laquelle le décalage V est ramené devant chacun des filtres K.. et K.., la fonction de transfert globale du filtre interpolateur étant égale à la somme de toutes les contributions K. . Dans ces conditions le calcul de la contribution d'une bande j à l'interpolation d'un point d'une bande i peut être fait en prenant l'une des N. phases de la décomposition polyphase à l'ordre N. d'un filtre K.. et le numéro de la phase à prendre en compte est alors déterminé par le décalage V.A titre d'exemple, si K(n) est la réponse impulsionnelle d'un filtre K.., une décomposition polyphase d'ordre N. de ce filtre permet de définir N. phases K. de noj où j est un entier allant de 0 à J J N-l tel que K. = (n) =K Nj); Ainsi pour un ordre de décomposition N j = 4, la décomposition polyphase d'un filtre comportant n+1 échantillons de phase centrée sur l'échantillon d'ordre 0 donnera pour la suite ordonnée des échantillons (-6, -5...0,1...0) les échantillons suivants. Phase 9 j K.(-4), K.(0), K.(4), I'origine étant placée sur J J J l'échantillon 0 Phase 1 K.(-3), Ka (1), Kj(3), I'origine étant placée sur J J J l'échantillon 1 Phase 2 Kj(-6), Kj(-2), K. (2), K. (6), I'origine étant placée J J sur l'échantillon 2 et Phase 3 K.(-5), K.(-î), Kj(3), I'origine étant placée sur J J J l'échantillon 3. Ainsi si V est un multiple de N. la phase 0 du filtre est utilisée et rien n'est modifié, par contre, dans le cas contraire, le numéro de phase n à utiliser est défini par n = Nj-R où R est le reste de la division de V par N.. Comme déjà souligné précédemment l'interpolation d'un point d'une bande donnée devrait pour être parfaite utiliser l'information contenue dans toutes les autres bandes, la contribution de ces bandes étant déterminée par les filtres K... Du fait de la contrainte de compatibilité il n'est pas possible d'utiliser toutes les bandes mais en fait comme le montre les figures 14a, 14b, 14c et 14d la contribution de toutes les bandes n'apparaît pas indispensable. Il apparaît sur les courbes représentées que ce sont essentiellement les sous-bandes adjacentes à une sous-bande considérée qui contribuent le plus à l'interpolation de la sous bande considérée et que les contributions des autres sous-bandes peuvent être considérées comme négligeables. Ainsi, il est possible de se contenter d'utiliser les bandes adjacentes à une bande considérée en ne conservant qu'une interpolation pratiquement optimale tout en satisfaisant les contraintes de compatibilité. Dans les exemples précédents, il a été supposé que les vecteurs mouvements estimés sur un signal avant décomposition étaient entiers mais une généralisation au cas d'une estimation ayant une précision au demi pixel est tout aussi envisageable comme le montre le schéma de la figure 15. En effet le vecteur mouvement V devient un multiple de 1/2 si bien que le décalage V de la figure 11 peut alors conduire sur un point n'existant pas. Il suffit alors de remplacer ce bloc de décalage par le schéma de la figure 15 où le décalage de 2V devient entier et où il suffit d'utiliser le filtre B (z) interpolateur ayant servi à estimer le mouvement au demi pixel. Ainsi le schéma d'interpolation d'un point d'une bande de la figure 12 peut être ramené à celui représenté à la figure 16 dans le cas d'une estimation au demi pixel. Si l'interpolation optimale d'un point d'une sous-bande doit en théorie utiliser l'information contenue dans toutes les bandes au voisinage des points correspondants au point interpolé, ce voisinage étant déterminé par la longueur des filtres d'interpolation Kit, l'application directe de ce procédé à un système de codage pourrait détruire les propriétés de compatibilité puisque par exemple des signaux de bande EDP devraient pouvoir utiliser l'information contenue dans la partie complémentaire CHDP. Dès lors il paraît indispensable de limiter l'interpolation à l'utilisation de certaines bandes adjacentes à la bande considérée seulement. En fait, en considérant l'arbre optimal à 16 bandes de la figure 6 l'utilisation des bandes adjacentes peut être limitée de la façon suivante. Pour les bandes numérotées 13, 14, 15 et 16 l'interpolation n'aura à utiliser que l'information des bandes 13, 14, 15 et 16. Pour les bandes 2, 3, 4 l'interpolation devra être limitée à l'utilisation des bandes 2, 3, 4, 13, 14, 15 et 16. Pour les bandes 1, 5, 6, 7, 8, 9, 10, 11, 12 I'interpolation de ces bandes pourra utiliser les informations contenues dans toutes les autres bandes. Ces conditions de compatibilité sont résumées sur la figure 17 où les flèches barrées décrivent les interdictions d'utilisation de l'information contenue dans les bandes pour l'interpolation d'une bande de donnée. En fait comme cela apparaît sur les courbes des figures 14a à 14c ces limitations n'entraînent pas de pertes de performances puisqu'au bord des bandes immédiatement adjacentes les contributions à l'interpolation de la bande considérée sont négligeables.Les seuls problèmes qui pourraient se poser concernent l'interpolation des bandes situées à la frontière de deux domaines, par exemple, dans le cas de la bande 3 qui ne peut pas utiliser I'information de la bande li qui lui est directement adjacente et contribue donc théoriquement de manière non négligeable. Néanmoins ces contraintes n' amènent pas en pratique de pertes visibles des performances. Par ailleurs la limitation à l'utilisation des bandes directement adjacentes vient limiter considérablement le nombre d'opérations à effectuer et les accès mémoires, et réduit ainsi les coûts des réalisations matérielles. Le tableau récapitulatif suivant donne les numéros des bandes utilisées pour interpoler les différentes bandes de l'arbre de décomposition décrit à la figure 6. Numéro de bande Bandes adjacentes utilisées 1 1, 7, 8, 10, 12 2 2, 4, 14, 16 3 3, 4, 15, 16 4 4, 2, 3 5 5, 6r 7 6 6, 2, 5, 8 7 7, 1, 5, 8 8 8, 1, 4, 7, 6 9 9, 10, 11 10 10, 1, 9, 12 11 i1, 3, 9, 12 12 12, 1, 4, 10, il 13 13, 14, 15 14 14, 16, 13 15 15, 13, 16 16 16, 14,15 Une autre simplification peut encore être effectuée pour limiter le nombre d'opérations et les accès mémoires, en effet les phases des filtres d'interpolation K.. à utiliser ont des longueurs qui sont de l'ordre de la longueur des filtres de décomposition en sous-bandes 110 et H1. Les filtres sous-bandes couramment utilisés ayant des longueurs de l'ordre de 16 ceci représente un nombre d'opérations important à effectuer mais surtout des accès mémoires considérables, d'autant plus que ces accès sont aléatoires puisqu'ils dépendent de vecteurs mouvements En pratique une limitation des longueurs des phases des filtres à 8 coefficients, quatre coefficients de part et d'autre du coefficient central, n'entraîne pas de pertes d'efficacité du point de vue visuel du système de codage. Les accès mémoires se réduisent alors à des masques de 8x8 contre 16x16. REVENDICATIONS 1. Procédé d'interpolation adapté de signaux sous-bandes pour le codage et le décodage compatible inter-image de signaux de télévision de résolutions différentes, les signaux étant transmis sur un canal de transmission (18) entre au moins un codeur d'émission (13. . .22) et un décodeur de réception relié à un récepteur, dans lequel le codage consiste à découper en sous-bandes (i,..6 ; 31...36), suivant une même structure arborescente par filtrage et décimation, l'image pleine bande à transmettre et coder indépendamment les signaux sous-bandes obtenus et dans lequel le décodage consiste à décoder les signaux reçus de chaque sous-bande en fonction de leur résolution et de celle propre au récepteur, caractérisé en ce qutil consiste à effectuer au niveau du codeur d'émission un codage inter-image compensé en mouvement dans chaque sous-bande, à partir des vecteurs mouvements estimés sur l'image non décomposée et remis à l'échelle en fonction du taux de sous échantillonnage de la bande considérée, et à effectuer dans chaque sous bande, lorsque le vecteur mouvement (déduit de l'estimation sur l'image non décomposée) ne correspond pas à un nombre entier de pixels, un calcul d'interpolation du pixel correspondant en tenant compte non seulement du signal contenu dans la sous bande elle-même mais également des signaux des autres sous-bandes. 2. Procédé selon la revendication 1, caractérisé en ce que le calcul d'interpolation consiste à déterminer le pixel interpolé dans une sous-bande donnée au moyen de filtres (Kil) tenant compte de l'information contenue dans la sous-bande eIle-même, et de filtres (Kü) tenant compte de l'information contenue dans les autres sous-bandes. 3. Procédé selon la revendication 2, caractérisé en ce que les filtres d'interpolation (Kil et Ki;) opèrent par ij décomposition polyphase de leur réponse impulsionnelle, la phase des filtres utilisée étant déterminée par le vecteur mouvement transmis par le codeur. 4. Procédé selon l'une quelconque des revendications là 3, caractérisé en ce que le calcul d'interpolation est effectué dans chaque sous-bande à partir de la sous-bande elle-même sur les signaux relatifs à la sous-bande et des sous-bandes adjacentes. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que le calcul d'interpolation est effectué avec une précision déterminée par le taux de sous-échantillonnage de la sous-bande considérée et de la précision d'estimation du vecteur mouvement égale à une fraction de la distance séparant deux pixels consécutifs de l'image 6. Procédé selon la revendication 5 caractérisé en ce que pour assurer la compatibilité du codage entre standards de télévision type VT, EDP et HDP, le calcul d'interpolation dans une sous-bande est effectué à partir de l'information contenue dans cette sous-bande et uniquement à partir d'informations contenues dans certaines sous-bandes adjacentes déterminées en fonction de la position spectrale de la sous-bande correspondant au pixel interpolé. 7. Procédé selon les revendications 4 et 6 caractérisé en ce qu'il consiste lorsque l'estimation du mouvement sur les images non décomposées est effectuée au demi pixel, à multiplier par deux l'ordre d'interpolation des sous-bandes, en remplaçant les filtres d'interpolation au pixel entier K.. (z) par des filtres de fonctions de transfert B (z). K1j (z 2) dans lesquelles B (z) est la fonction de transfert d'un filtre interpolateur d'une estimation au demi pixel, et en utilisant pour phase de chaque filtre B(z)Kij(z) deux fois la longueur du vecteur mouvement V. 8. Procédé selon l'une quelconque des revendications 3 à 7 caractérisé en ce que la phase des filtres est limitée à un nombre de coefficients inférieurs ou égaux à 8. 9. Dispositif pour la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce qu'il comprend pour chaque sous-bande un codeur inter-image commandé à partir d'un dispositif d'estimation de mouvement. 10. Dispositif selon la revendication 9, caractérisé en ce que chaque codeur d'une sous-bande comprend un dispositif quantificateur (14i) couplé à un dispositif de codage différentiel par l'intermédiaire d'un dispositif d'aiguillage (25) pour quantifier chaque pixel fourni par un bloc directement lorsque le mode de codage intra-image a été choisi ou quantifier sa différence par rapport à son homologue corrigé par des vecteurs mouvement dans l'image précédente lorsque le mode de codage inter-image a été choisi. 11. Dispositif selon la revendication 10, caractérisé en ce que chaque dispositif de codage différentiel est couplé à un dispositif d'estimation de mouvement (20) des points de l'image par l'intermédiaire d'un dispositif de compensation de mouvement (96.). 12. Dispositif selon la revendication 11, caractérisé en ce que chaque dispositif de compensation de mouvement (28) est couplé au dispositif d'estimation de mouvement (20) par l'intermédiaire d'un dispositif de mise à l'échelle (21) des vecteurs mouvement fournis par le dispositif d'estimation de mouvement (20), de manière que l'estimation de mouvement soit faite au niveau de l'image non décomposée alors que la compensation de mouvement a lieu sur les sous-bandes. 13. Dispositif selon l'une quelconque des revendications 9 à 12, caractérisé en ce que chaque décodeur d'une sous-bande comporte des filtres à décomposition polyphase pour effectuer l'interpolation du vecteur mouvement dans la sous-bande en tenant compte des signaux des autres sous-bandes, l'ordre de chaque filtre étant déterminé par le vecteur mouvement reçu du codeur.