FR 2684258 A1 19930528 FR 9114625 A 19911127 PROCEDE D'ESTIMATION ET DE CODAGE HIERARCHISE DU MOUVEMENT DE SEQUENCES D'IMAGES La présente invention concerne un procédé d'estimation et de codage hiérarchisé du mouvement de séquences d' images utilisé notamment dans les chaînes de transmission d'images en vue de réduire le débit d'informations. On a proposé, dans la demande de brevet français N 89 11328 déposée au nom de THOMSON CONSUMER ELECTRONICS, un procédé et un dispositif d'estimation et de codage hiérarchisé du mouvement de séquences d'images. Ce procédé est caractérisé en ce qu'il consiste, à un premier niveau, à découper l'image courante en macroblocs 2p+1 2p+l pixels et à déterminer un premier champ de vecteurs-mouvement associé à cette partition en macroblocs en utilisant un algorithme d'estimation itératif et récursif initialisé par des vecteurs-mouvement estimés pour l'image précédente, puis, à un deuxième niveau à découper en quadrants chacun des macroblocs et à déterminer, pour les blocs qui en résultent, un second champ de vecteurs-mouvement en utilisant le même algorithme d'estimation mais initialisé par des vecteurs du eme champ de mouvement estimé au niveau précédent, puis, à un i niveau, i = 3 à p, à découper chacun des blocs considérés au niveau i-1 en quadrants et à déterminer, pour les blocs qui en ème résultent, un i champ de vecteurs-mouvement en utilisant le même algorithme d'estimation initialisé par des vecteurs du champ de mouvement estimés au niveau précédent, les blocs de taille minimale étant des blocs de 2l+1.2l+1 pixels et à déterminer un champ de vecteurs-mouvement final à partir des p champs de vecteurs-mouvement résultants en choisissant le niveau de partition le moins élevé pour lequel le vecteur-mouvement associé au bloc correspondant conduit à la minimisation d'un critère traduisant les différences de luminance entre blocs se correspondant dans les images successives par les vecteurs de déplacement estimés. Cette méthode d'estimation permet d'augmenter la convergence de l'algorithme d'estimation de mouvement récursif utilisé et permet d'adapter au mieux l'algorithme d'estimation de mouvement au codage de type "Quadtree" qui sera fait ultérieurement du champ de mouvement résultant. Cette méthode d'estimation peut être utilisée, entre autre, dans le cadre du HD-MAC. La présente invention a pour but d'apporter différents perfectionnements au procédé d'estimation et de codage hiérarchisé du mouvement de séquences d' images décrit dans la demande de brevet français N"89 11328. En conséquence, la présente invention a pour objet un procédé d'estimation et de codage hiérarchisé du mouvement de séquences d'images , consistant - à un premier niveau, à découper l'image courante en macroblocs de 2P.2P pixels et à déterminer un premier champ de vecteurs-mouvement, associé à cette partition en macroblocs, en utilisant un algorithme d'estimation itératif et récursif initialisé par des vecteurs-mouvement estimés pour l'image précédente, - puis, à un deuxième niveau, à découper en quadrants chacun des macroblocs et à déterminer, pour les blocs qui en résultent, un second champ de vecteurs-mouvement en utilisant le même algorithme d'estimation mais initialisé par des vecteurs du champ de mouvement estimé au niveau précédent, - puis à un i me niveau, i = 3 à p, à découper chacun des blocs considérés au niveau i-l en quadrants et à déterminer, pour ème les blocs qui en résultent, un i champ de vecteurs-mouvement en utilisant le même algorithme d'estimation initialisé par des vecteurs du champ de mouvement estimés au niveau précédent, les blocs de taille minimale étant des blocs de 2 2p pixels, - et à déterminer un champ de vecteurs-mouvement final à partir des p champs de vecteurs-mouvement résultants, caractérisé en ce que le champs de vecteurs-mouvement est déterminé en utilisant les différences de luminance entre blocs se correspondant dans les images successives par les vecteurs de déplacement estimés des p niveaux de partition correspondant à chaque bloc de taille minimale et en choisissant comme vecteur mouvement associé au bloc de taille minimale celui pour lequel la différence de luminance est minimale sur ce bloc. Cette nouvelle détermination du champ de vecteurs-mouvement final permet de simplifier le système, notamment au niveau de sa réalisation, et d'assurer une meilleure convergence des calculs. D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description faite ci-après de divers modes de réalisation, cette description étant faite avec référence aux figures ci-annexées dans lesquelles - la figure 1 illustre les informations utilisées pour estimer un champ de mouvement selon l'invention, dans le cas d'un bloc courant - la figure 2 illustre le calcul des gradients conformément à la présente invention - la figure 3 illustre les informations utilisées pour estimer un champ de mouvement dans le cas de blocs se trouvant près du bord de l'image ou de la trame - la figure 4 illustre la décomposition d'un macrobloc en blocs plus petits suivant cinq niveaux hiérarchisés - la figure 5 est un organigramme détaillant les différentes étapes du procédé de la présente invention - la figure 6 détaille les étapes de l'estimation de mouvement à un niveau donné selon un premier mode de réalisation de la présente invention, et - la figure 7 détaille les étapes de l'estimation de mouvement à un niveau donné selon un autre mode de réalisation de la présente invention. On décrira tout d'abord un type d'algorithme d'estimation de mouvement qui peut être utilisé dans le cadre de la présente invention. Cet algorithme est dérivé de l'algorithme spatiotemporel décrit notamment par A.N. Netravali et J.D Robbins dans un article intitulé "Motion compensated television coding parut.1" tiré de "The Bell system technical journal" volume 58, N"3, March 1979. L'algorithme de la présente invention permet une estimation de mouvement dite "symétrique". Une estimation symétrique est utilisée principalement pour les blocs courants se trouvant au milieu d'une trame ou image comme représenté sur la figure i. Dans ce cas, l'estimation de mouvement des blocs de la trame courante à l'itération i est faite en utilisant la trame précédente se produisant à l'instant t-l et la trame suivante se produisant à l'instant t+1. Les gradients sont aussi calculés dans les trames précédentes et suivantes et le vecteur de déplacement à l'itération i est alors calculé à partir du vecteur déplacement du bloc courant estimé à 1' itération précédente i-l par les équations suivantes (qui tiennent compte des différences de luminance entre les pixels du bloc courant et les pixels des blocs correspondants, déplacés de ce vecteur estimé précédemment, d'une part dans la trame précédente et d'autre part dans la trame suivante) dans laquelle - Dxi et Dyi ainsi que Dxi-1 et Dyi-1 sont respectivement les composantes horizontales et verticales des vecteurs-mouvement Di et D pour un bloc de pixels respectivement à l'itération i et à l'itération précédente i-l - DFD (z-Di-1, t-l) correspond à la différence inter-trame déplacé pour le pixel z, c'est-à-dire la différence entre la luminance du pixel z dans la trame courante à l'instant t et la luminance du pixel déplacé correspondant à z dans la trame précédente (à l'instant t-l), c'est-à-dire z déplacé du vecteurmouvement Di-1 calculé à l'itération précédente DFD(z-Di-1, t-1) = I (z, t) - I (z - Di-1, t-l) où I(z, t) représente la valeur de luminance d'un pixel z à l'instant T - DFD(z+Di-1, t + 1) représente la différence inter-trame déplacée pour le pixel z, c'est-à-dire la différence entre la luminance du pixel z à la trame courante à l'instant t et la luminance du pixel déplacé correspondant à z dans la trame suivante (à l'instant t+1), c'est-à-dire à z déplacé du vecteur-mouvement symétrique Dii calculé à l'itération précédente - grad. x (z - D i-i t-1), grad.X(z + D i-i t + 1), grad. Cz - D i-i t - 1), grad.Y(Z + Di-1, t + 1) représentent les composantes respectivement horizontales et verticales du gradient de luminance des pixels déplacés, calculées respectivement dans la trame précédente et dans la trame suivante, - et où sg (...) est la notation pour "signe de" la quantité qui suit. Sur la figure 1, on a représenté la position des blocs déplacés dans la trame précédente et dans la trame suivante, correspondant à un bloc courant se trouvant au milieu de l'image pour lequel le vecteur déplacement estimé à l'itération précédente est D 1 de la trame précédente à la trame courante et donc -Di-1 de la trame suivante à la trame courante, en supposant un déplacement symétrique. Si l'on examine un détail du voisinage de z - Di-1 ou z + Di-1, os voit que, en général, le vecteur déplacement conduit à un pixel déplacé z - Di-1 ou z + Di-1 ayant une position qui ne correspond pas à des valeurs entières de x et y. En conséquence, on utilise les valeurs de luminance des quatres points du voisinage de z-D ou z+D 1 , à savoir les points I1, I2, I3 et I4 comme représenté sur la figure 2, pour le calcul des luminances des pixels déplacés et pour le calcul des gradients. Dans ce cas, l'interpolation de la luminance I est réalisée en effectuant une somme pondérée de I1, I2, I3 et I4. Pour le gradient, on utilise, conformément à la présente invention, des formules plus complexes, prenant en compte la position du pixel déplacé Z-D i-i ou z+Di+l par rapport aux points I1, I2, I3 et I4, comme représenté sur la figure 2. Si les calculs sont réalisés avec une précision au quart de pixels, les pixels déplacés peuvent avoir, entre les pixels existants I1, I2, I3, I4, un certain nombre de positions telles que représentées par des points sur la figure 2. Suivant la position du pixel déplacé par rapport au pixel existant, des formules différentes seront utilisées pour le calcul des gradients en x et en y. Ainsi, si le pixel déplacé se trouve dans le voisinage de I1, la valeur du gradient est donnée Gx = (I2 - I1) et Gy= (I1 - I3). Si le pixel déplacé se trouve dans le voisinage de I2 : Gx = (I2 - I1) et Gy = (I2 - I4). Si le pixel déplacé se trouve dans le voisinage de I3 : Gx = (I4 - I3) et Gy = (I1 - I3). Si le pixel déplacé se trouve dans le voisinage de I4 : Gx = (I4 - I3) et Gy = (I2 - I4). Si le pixel déplacé se trouve dans la partie médiane entre I1 et I3 ou entre I2 et I4 : Gx = ((I2 - I1) + (I4 - I3))/2 et Gy = (I1 - I3) ou (I2 - I4). Si le pixel déplacé se trouve dans la partie médiane entre I1 et I2 ou entre 14 et Gx = (I2 - I1) ou (I4 - I3), Gy = ((I1 - I3) + (I2 - I4))/2 Si le pixel déplacé est équidistant de I1, I2, I3, I4 Gx = ((I2 - I1) + (I4 - I3))/2 Gy = ((I1 - I3) + (I2 - I4))/2 Gx et Gy dans les équations ci-dessus correspondent en fait à grad X(Z-Di-l, t - 1) et grad.Y(Z-Di-1, t - 1) avec 11, 12, 13 et I4 à t - 1, de même pour (Z+D , t+1). Toutefois, dans le cas des blocs des bords de l'image, il n'est pas toujours possible de réaliser une estimation symétrique. Un vecteur calculé peut très bien pointer en dehors de l'image sur l'une des deux trames adjacentes à la trame courante. Ainsi, comme représenté sur la figure 3, le pixel Z2 d'un bloc se trouvant sur le bord gauche de l'image pointe en dehors de la trame suivante t + 1 tandis que le pixel Z1 appartenant à un bloc du bord droit de l'image pointe en dehors de la trame précédente se produisant à l'instant t-1. Dans ce cas, afin de disposer de prédicteurs pour les niveaux inférieurs, on réalise une estimation de mouvement non symétrique, c'est-à dire une estimation de mouvement entre la trame courante et la trame suivante ou entre la trame courante et la trame précédente. Les différences de luminance entre blocs ou DFD et les gradients sont calculés sur une seule trame et mis à zéro ou à la même valeur pour l'autre trame afin de pouvoir utiliser le même type d'équation que les équations (1) ci- dessus et conserver la même architecture générale du système. L'algorithme d'estimation de mouvement récursif tel que décrit ci-dessus est donc utilisé dans le procédé d'estimation et de codage hiérarchisé du mouvement de séquences d'images conforme à la présente invention. Pour adapter au mieux cet algorithme d'estimation de mouvement au codage de type "quadtree" qui sera fait ultérieurement du champ de mouvement résultant, les blocs d'un niveau sont le résultat du partage en quatre d'un bloc du niveau immédiatement supérieur, comme représenté sur la figure 3 qui illustre les différentes tailles de blocs résultant de la partition. De plus, la récursivité est obtenue en utilisant le champ de mouvement obtenu à un niveau donné de la hiérarchie pour calculer le champ de mouvement du niveau immédiatement inférieur. Ainsi, pour les calculs à l'intérieur d'un niveau donné, l'estimation de mouvement pour un bloc est indépendante de celle des autres blocs du voisinage au même niveau. On décrira maintenant, avec référence, plus particulièrement à la figure 5, les principales étapes du procédé d'estimation et de codage hiérarchisé du mouvement de séquences d'images conforme à la présente invention. Ces différentes étapes peuvent être mises en oeuvre avec un estimateur de mouvement BRH (pour blocs récursifs hiérarchiques). Les principales étapes du procédé conforme à la présente invention comportent tout d'abord une étape d'estimation de mouvement qui se répète pour les n niveaux suivie d'une étape de convergence qui donne un champ de vecteurs-mouvement final et se terminent enfin par une étape de sélection, cette étape de sélection étant en fait optionnelle suivant l'application concernée. On décrira l'étape d'estimation de mouvement, à savoir, la première étape de la figure 5, selon différents modes de réalisation de l'invention. La taille maximale des blocs d'éléments d'images reconnue pour le codage étant de 128 x 128 pixels et la taille minimale, dans le mode de réalisation représenté, de 8 x 8 pixels mais pouvant aller jusqu'au pixel par partage successifs de chaque bloc en 4 quadrants, le calcul des vecteurs déplacement est effectué selon l'invention d'abord pour les blocs de 128.128 pixels (niveau 1) puis de 64.64 pixels (niveau 2), puis de 32.32 pixels (niveau 3), puis de 16.16 pixels (niveau 4) et enfin de 8.8 pixels (niveau 5), indépendamment du contenu de l'image. Pour l'initialisation pour l'image courante, les N valeurs initiales pour les blocs les plus grands au premier niveau (128.128) sont nulles pour la première image de la séquence, puis égales au N vecteurs-mouvement les plus occurents ou prédicteurs temporels dans le même macrobloc de l'image précédente. Puis les N valeurs initiales recueillies pour un bloc de taille n.n proviennent du niveau précédent, donc d'un bloc de taille 2n.2n. Ces valeurs sont les déplacements estimés pour les N blocs les plus proches de taille supérieure. En général, N est choisi égal à 4, mais d'autres valeurs peuvent être envisagées sans sortir du cadre de la présente invention. Une étape du procédé d'estimation de mouvement à un niveau n est décrite plus en détail ci-après en référence à la figure 6. Dans ce cas, on a décrit un calcul à posteriori des vecteurs déplacements. Comme représenté sur la figure 6, les séquences d'images sont mémorisées dans trois mémoires ou zones mémoires d'images ou de trames qui contiennent donc l'image suivante M1 à t+l, l'image courante M2 à l'instant t et l'image précédente M3 à l'instant t-l ainsi que les valeurs de luminance des 4 pixels entourant le pixel déplacé Z+ à la sortie de M1, la valeur de luminance I(t) = I(z, t) à la sortie de M2 et les valeurs de luminance des 4 pixels entourant le pixel déplacé Z-D 1 à la sortie de M3. Une première étape de calcul par point donne les valeurs des différences inter-trames déplacées et des gradients à partir des valeurs de luminance issues des mémoires et des vecteurs déplacement initiaux, par rapport à l'image précédente à t-1 et par rapport à l'image suivante à t+1. De manière plus spécifique, on calcule dans les blocs référencés 1 et 2 en utilisant les informations stockées dans M1 et M2 - DFDi-1 (t+1), Gradxi-1 (t+1), Grad yi-1 (t+1) - A = DFDi-1 (t+1).signe (Grad xi-1 (t+1)) - B = DFDi-1 (t+1).signe (Grad yi-1 (t+1)) en utilisant les informations stockées dans M2 et M3 - DFDi-1 (t-1), Grad xi-1 (t-1), Grad yi-1 (t-1) - C = DFDi-1 (t-l).signe (Grad xi-1 (t-1)) - D = DFDi 1 (t-l).signe (Grad yi-1 (t-1)) L'étape suivante consiste alors à calculer les composantes des vecteurs déplacement Di et Di selon les formules itératives x y données ci-dessus après avoir calculé les sommes partielles des différences telles que C-A et B-D sur l'ensemble du bloc correspondant de niveau n. Ainsi, comme représenté sur la figure 6, l'on calcule dans le bloc 3 et l'on effectue la correction des prédicteurs D . Lorsque les N valeurs initiales pour chaque bloc courant (N étant en général choisi égal à 4) correspondant à des blocs du voisinage de ce bloc courant mais d'un niveau différent, c' est-àdire de taille immédiatement supérieure, permettent d'initialiser 4 estimations de mouvement séparées pour ce bloc (une seule a été représentée) et donnent lorsque les formules itératives de l'algorithme ont convergé quatre valeurs de vecteurs-mouvement mis à jour, indépendantes, qui complètent les valeurs des vecteurs de niveaux précédents. Les carrés des différences intertrames DFD intermédiaires utilisés pour les mises à jour sont calculés, stokés en mémoire et additionnés pour l'ensemble du bloc de façon à choisir pour chaque bloc à l'issue des 4 calculs, le "meilleur" déplacement parmi les quatre nouveaux vecteurs de déplacement mis à jour, c'est-à-dire celui pour lequel la somme S des différences inter-trames élevées au carré pour ce bloc est minimale. Ce meilleur vecteur constitue alors l'un des vecteurs du champ de mouvement calculé au niveau n et noté DV128 ou DV64 ou DV32 ou DV16 ou DV8 suivant le niveau d'itération. Selon un autre mode de réalisation représenté à la figure 7, on peut réaliser un calcul à priori. Dans ce cas, on applique directement sur les N valeurs initiales utilisées pour chaque bloc courant un critère de décision en utilisant un calcul de DFD et en choisissant comme vecteur celui qui donne la DFD minimale. Sur ce vecteur choisi, on réalise une correction de manière à obtenir les composantes du vecteur déplacement Di et Di avec les x y formules itératives données ci-dessus. De manière plus précise et comme représenté sur la figure 7, dans les blocs 1 et 2, on effectue le même calcul que dans les blocs 1 et 2 de la figure 6. Puis dans les blocs 3', on effectue le calcul des sommes sur chaque bloc 8x8 du bloc courant puis on somme sur le bloc courant, à savoir i-i de i-i (DFD étant la DFD calculées avec les vecteurs D de même pour Grade 1). Ensuite, on effectue la calcul des sommes sur le bloc courant comme dans le cas du bloc 3 de la figure 3 et le calcul de Di et Di, à savoir la correction des prédicteurs D . x Dans les deux cas décrits ci-dessus, la correction du vecteur déplacement a été réalisée en utilisant un calcul de gradient. Toutefois, il est évident pour l'homme de l'art que la correction d'un vecteur mouvement peut être réalisée en effectuant une correspondance de blocs, à savoir un "block matching" en langue anglaise sur une fenêtre limitée, comme cela est décrit dans la littérature notamment dans les articles suivants - J.R Jain et A.K Jain IEEE Trans. ou Comm. COM-29 N 121 "Displacement Measurement and its application in interframe image coding" ; - Y.Niromiga et Y. Ohtsuk IEEE Vol COM-30 Janvier 1982 "A motion compensated interframe coding scheme for television pictures". A l'issue de cette phase d'estimation de mouvement, qui peut être réalisée différemment comme décrit ci-dessus, un champ de vecteurs-mouvement différents correspondant à chaque niveau, c'est-à-dire à chaque taille de blocs, est disponible pour déterminer le champ de mouvement final adapté à la construction de l'arbre de codage tel qu'un arbre de codage de type "quadtree". Une fois les cinq champs de mouvement déterminés, l'étape de convergence consiste donc à déterminer le champ de mouvement final. Conformément à la présente invention, la détermination du champ de mouvement final se fait en examinant pour chaque bloc de taille minimale, à savoir, pour chaque bloc de taille 8 x 8 pixels dans le mode de réalisation représenté, les différents vecteursmouvement DVI qui ont été établis à chaque niveau de partition. Ainsi, pour chaque bloc de taille minimale, on a un vecteur DV128 donnant une DFD 128 établi lors de l'initialisation du processus, à savoir au niveau l puis un vecteur DV64 donnant une DFD 64 établi au niveau 2, un vecteur DV32 donnant une DFD 32 établi au niveau 3, un vecteur DV16 donnant une DFD 16 établi au niveau 4, et un vecteur DV8 donnant une DFD 8 établi au niveau 5. Ainsi, pour ces cinq vecteurs, on compare leur DFD et on choisit le vecteur DVI donnant la DFDI la plus faible comme vecteur affecté à ce bloc de taille minimale dans le champ de vecteur-mouvement final. Cette même sélection est effectuée pour tous les blocs 8x8 d'un macrobloc, comme représenté dans la partie convergence de la figure 5. Comme représenté sur la figure 5, le procédé conforme à la présente invention peut se terminer par une étape de sélection. Cette étape est facultative et en fait a pour but d'homogénéiser les champs de mouvement résultants et d'accroître l'efficacité d'un codage réducteur de débit. Cette étape consiste à réaliser un histogramme des vecteurs présents dans un macrobloc d' images puis à ne garder que les vecteurs majoritaires. Ainsi, on ne garde que les vecteurs-mouvement correspondant à un nombre de pics de l'histogramme supérieur à un seuil donné, ce seuil correspondant à un nombre de vecteurs-mouvement jugé suffisant. Pour les autres vecteurs, soit les vecteurs les moins fréquents, ils sont remplacés par les vecteurs de valeur les plus proches parmi les vecteurs sélectionnés de manière à effectuer une réassignation. Dans le cadre d'une application HD-MAC, cette étape de sélection n'est pas réalisée pour les images paires, mais une réassignation est effectuée sur chaque trame paire des images paires parmi les 9 blocs voisins du bloc considéré sur l'image précédente (impaire). Dans le mode de réalisation représenté à la figure 5, on a réalisé un histogramme sur 13 vecteurs au maximum, avec une réassignation. D'autre part, il est possible d'utiliser le procédé d'estimation ci-dessus en descendant jusqu'à une taille de blocs équivalente au pixel. Dans ce cas, on utilise un algorithme d'estimation non symétrique et on ne réalise pas d'histogramme. REVENDICATIONS 1. Procédé d'estimation et de codage hiérarchisé du mouvement de séquences d'images, consistant - à un premier niveau, à découper l'image courante en macroblocs de 2P.2P pixels et à déterminer un premier champ de vecteurs-mouvement associé à cette partition en macroblocs, en utilisant un algorithme d'estimation itératif et récursif initialisé par des vecteurs-mouvement estimés pour l'image précédente, - puis, à un deuxième niveau, à découper en quadrants chacun des macroblocs et à déterminer, pour les blocs qui en résultent, un second champ de vecteurs-mouvement en utilisant le même algorithme d'estimation mais initialisé par des vecteurs du champ de mouvement estimés au niveau précédent, ème - puis à un i niveau, i = 3 à p, à découper chacun des blocs considérés au niveau i-l en quadrants et à déterminer, pour ème les blocs qui en résultent, un i champ de vecteurs-mouvement en utilisant le même algorithme d'estimation initialisé par des vecteurs du champ de mouvement estimés au niveau précédent, les blocs de taille minimale étant des blocs de 21.21 pixels, - et à déterminer un champ de vecteurs-mouvement final à partir des p champs de vecteurs-mouvement résultants, caractérisé en ce que le champs de vecteurs-mouvement est déterminé en utilisant les différences de luminance entre blocs se correspondant dans les images successives par les vecteurs de déplacement estimés des p niveaux de partition correspondant à chaque bloc de taille minimale,et en choisissant comme vecteurmouvement associé au bloc de taille minimale celui pour lequel la différence de luminance est minimale sur ce bloc. 2. Procédé selon la revendication 1, caractérisé en ce que l'algorithme d'estimation de mouvement utilisé est symétrique ou non symétrique et utilise pour la détermination d'un vecteurmouvement à un rang d'itération i, N vecteurs-mouvement déterminés à l'itération i-l parmi lesquels on choisit le vecteur qui conduit à la minimisation du critère traduisant les différences de luminance entre blocs se correspondant dans les. images successives par les vecteurs de déplacement estimés correpondants, le vecteur choisi étant corrigé en fonction du déplacement des pixels du bloc entre la trame courante et la trame précédente et/ou la trame suivante. 3. Procédé selon la revendication 1, caractérisé en ce que l'algorithme d'estimation de mouvement utilisé est symétrique ou non-symétrique et utilise pour la détermination d'un vecteurmouvement à un rang d'itération i, N vecteurs-mouvement déterminés à l'itération i-l corrigés en fonction du déplacement des pixels du bloc entre la trame courante et la trame précédente et/ou la trame suivante, le vecteur-mouvement choisi étant celui qui conduit à la minimisation du critère traduisant les différences de luminance entre blocs se correspondant dans les images successives par les vecteurs de déplacement estimés correspondants. 4. Procédé selon les revendications 2 ou 3, caractérisé en ce que la correction du vecteur-mouvement est réalisée en diminuant ledit vecteur d'une quantité fonction des différences de luminance intertrames déplacées, c'est-à-dire des sommes de différences de luminance entre les pixels du bloc courant et les pixels des blocs correspondant au bloc courant déplacés du vecteur-mouvement, respectivement dans la trame précédente et/ou dans la trame suivante de la séquence. 5. Procédé selon la revendication 4, caractérisé en ce que la correction du vecteur-mouvement est fonction des composantes horizontale et verticale des gradients de luminance calculés dans les trames précédentes et/ou suivantes. 6. Procédé selon la revendication 5, caractérisé en ce que le gradient est calculé en utilisant les quatre pixels ayant une position correspondant à des valeurs entières de x et y les plus proches du pixel déplacé et en appliquant les formules ci-après - si le pixel déplacé est dans le voisinage de I1 Gx = (I2 - I1) Gy = (I1 - I3) - si le pixel déplacé est dans le voisinage de I2 Gx = (I2 - I1) Gy = (I2 - I4) - si le pixel déplacé est dans le voisinage de I3 Gx = (I4 - I3) Gy = (I1 - I3) - si le pixel déplacé est dans le voisinage de I4 Gx = (I4 - I3) Gy = (12 - I4) - si le pixel déplacé est dans la partie médiane entre I1 et I3 ou I2 et I4 Gx = ((I2 - I1) + (I4 - I3))/2 Gy = (I1 - I3) ou (I2 - I4) - si le pixel déplacé est dans la partie médiane entre I1 et I2 ou I4 et I3 Gx = (I2 - I1) ou (I4 - I3) Gy = ((I1 - I3) + (I2 - I4))/2 - si le pixel déplacé est équidistant de I1, I2, I3, I4 Gx = ((I2 - I1) + (I4 - I3))/2 Gy = ((I1 - I3) + (I2 - I4))/2 7. Procédé selon les revendications 2 ou 3, caractérisé en ce que la correction du vecteur-mouvement est réalisée en effectuant une correspondance de blocs, à savoir un "block matching" en langue anglaise, sur une fenêtre limitée. 8. Procédé selon l'une quelconque des revendications 2 à 6, caractérisé en ce que, lorsque l'algorithme d'estimation de mouvement utilisé est non-symétrique, les différences de luminance entre blocs ainsi que les corrections sont calculées sur la trame précédente ou la trame suivante et mise à zéro ou à la même valeur pour l'autre trame. 9. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que, à la première étape, pour chaque macroblocs, l'algorithme d'estimation de mouvement est mis en oeuvre N fois, initialisé par les N vecteurs-mouvement les plus occurents dans le même macrobloc de l'image précédente. 10. Procédé selon la revendicatin 7, caractérisé en ce que, pour la première image traitée, les N vecteurs-mouvement d'initialisation sont à zéro. 11. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que de plus, on homogénéise le champ de vecteurs-mouvement final en éliminant du champ de vecteursmouvement final les vecteurs les moins fréquents et en les remplaçant par les vecteurs les plus proches, de fréquence plus grande.