WO1996017327A1

WO1996017327A1 - Animation par interpolation

Info

Publication number: WO1996017327A1
Application number: PCT/FR1995/001553
Authority: WO
Inventors: Serge Beucher; Fernand Meyer; Jean Serra
Original assignee: Association Pour La Recherche Et Le Developpement Des Methodes Et Processus Industriels (Armines)
Priority date: 1994-11-25
Filing date: 1995-11-23
Publication date: 1996-06-06
Also published as: FR2727543B1; FR2727543A1

Abstract

Le procédé permet de générer N-1 images mosaïques intermédiaires entre deux images successives d'une même scène. Pour cela, on mémorise l'image initiale et l'image finale sous forme de pixels adressables affectés chacun à une des M classes; on compare deux à deux les pixels correspondants des deux images et on identifie ceux des pixels correspondants des deux images qui appartiennent à la même classe, constituant une partie commune, de façon à constituer au plus M ensembles X, Y, ..., chacun constitué de pixels d'une même classe; et on affecte, dans chaque image intermédiaire générée, les pixels de la partie variable chacun à une des M classes, suivant la valeur d'une des distances entre le pixel respectif et des régions prédéterminées de l'image, constituées par la partie commune ou des parties appartenant à une classe prédéterminée dans une ou des images pré-existantes.

Description

ANIMATION PAR INTERPOLATION

La présente invention a pour objet un procédé de génération d'images intermédiaires par interpolation entre deux images disponibles successives d'une même scène. Le terme "image" doit être considéré dans un sens large et comme désignant notamment des représentations d'une scène réelle ou virtuelle, mais aussi toute représentation à au moins deux dimensions sous forme de points élémentaires et pouvant être affectés chacun d'une adresse et d'un paramètre d'affectation à une classe parmi M classes. L'invention est applicable dans de nombreux domaines, et en particulier dans celui de la compression d'images animées, de codage d'images et de restitution d'images manquant dans une séquence.

En matière de compression d'images, on a déjà proposé de constituer des "images mosaïques" constituées par partition du plan en des cellules dont chacune est constituée de pixels homogènes, appartenant à une classe parmi plusieurs classes affectées chacune d'un ou plusieurs paramètres d'identification. La compression s'effectue alors en ne retenant pour transmission qu'une image sur plusieurs. Lors de la réception, il faut reconstituer une séquence ayant le même nombre d'images qu'à l'origine. L'invention vise dans ce cas particulier - qui n'est nullement exclusif - à reconstituer les images intermédiaires omises en tenant compte des modifications intervenues entre les deux images successives transmises. Une autre application est la restauration de films de cinéma, dont certains plans sont détériorés. Une autre application encore est la réalisation de dessins animés, en ne dessinant que certains des plans et en interpolant entre ces plans pour obtenir des images intermédiaires. Ces applications ne sont nullement limitati- ves. De façon générale, l'invention vise à fournir un procédé de génération de N-l images mosaïques intermédiaires (U i « _I._N-_I entre deux images successives U₀ et U_N d'une même scène, fractionnées en cellules appartenant chacune à une classe homogène parmi M classes, N étant une puissance de 2, en assurant une progressivité élevée de transformation pour passer de 1 ' image U₀ à 1 ' image U_N et cela en ne mettant en oeuvre que des opérations relativement simples.

Dans ce but, l'invention propose notamment un procédé suivant lequel :

(a) on mémorise l'image initiale et l'image finale sous forme de pixels adressables affectés chacun à une des M classes,

(b) on compare deux à deux les pixels correspondants des deux images et on identifie ceux des pixels correspondants des deux images qui appartiennent à la même classe, consti¬ tuant une partie commune, de façon à constituer au plus M ensembles X, Y, ..., chacun constitué de pixels d'une même classe ; et (c) on affecte, dans chaque image intermédiaire générée, les pixels de la partie variable chacun à une des M classes, suivant la valeur d'une fonction de la distance entre le pixel respectif et des régions prédéterminées de l'image, constituées par la partie commune ou des parties appartenant à une classe prédéterminée dans une ou des images pré¬ existantes.

Une solution particulière consiste à effectuer, au cours des étapes (c) , des dichotomies successives. La première étape (c) permet de déterminer, pour chaque pixel des parties variables, celui des M ensembles auquel il a la distance géodésique la plus courte et d'affecter le pixel à cet ensemble. Les étapes suivantes utilisent l'image ainsi créée pour recommencer l'opération.

Du fait que cela conduit à un doublement du nombre de points interpolés à chaque étape, le procédé impose, sauf à admettre une partition irrégulière dans le temps, une valeur de N égale à une puissance de 2.

Dans la plupart des cas, il y a non seulement évolution de la scène, mais déplacement global de la portion de la scène représentée, notamment dans le cas d'une scène réelle observée par une caméra effectuant un mouvement panoramique. Dans ce cas, il est avantageux d'évaluer le vecteur mouve¬ ment de la scène dans l'image, par un processus qui peut être l'un de ceux bien connus dans le domaine de la télévi- sion haute définition à diffusion avec compression tempo¬ relle, et d'effectuer la comparaison de l'étape (b) après avoir imposé des déplacements de sens inverse tenant compte de ce vecteur mouvement et permettant d'augmenter les parties communes. D'autres caractéristiques de l'invention apparaîtront à la lecture de la description qui suit de modes particuliers de réalisation, donnés à titre d'exemples non limitatifs. La description se réfère aux figures qui l'accompagnent, dans lesquelles : - les figures 1 et 2 sont des schémas montrant respec¬ tivement une image mosaïque initiale A et une image mosaïque finale B, entre lesquelles doivent être interpolées des images intermédiaires ;

- les figures 3 et 4 montrent, en noir, les parties communes à l'image initiale et à l'image finale, qui sont également communes à toutes les images de la séquence et, hachurées, les parties variables avec l'étiquette qu'elles ont dans l'image A (figure 3) ou dans l'image B (figure 4) ;

- la figure 4 montre les parties communes affectées de l'étiquette qu'elles ont dans toutes les images ;

- la figure 6 est un schéma destiné à montrer des premières distances géodésiques, désignées par d^B _irι, de deux points x et y appartenant à la partie variable, égales, pour x au plus court chemin constitué de points ayant tous la même étiquette que x dans l'image B et reliant x à un point quelconque des parties communes de même étiquette que x ;

- la figure 7, similaire à la figure 6, montre des secondes distances d*_in de deux points x et y dans l'ensemble A ; - la figure 8, similaire aux figures 6 et 7, montre des troisièmes distances géodésiques d^A _out des points x et y, c'est-à-dire pour x la longueur du plus court chemin reliant x à un point d'étiquette différente dans A et constitué uniquement de points ayant la même étiquette que x dans 1 ' image A ;

- la figure 9 montre des lignes de niveau correspondant aux frontières d'une facette dans trois interpolés ;

- les figures 10, 11 et 12 montrent trois interpolés complets entre A et B, correspondant à des distances respectives de 2/10, 4/10 et 6/10 de l'image A, les facettes correspondantes étant indiquées par des hachures de même sens ;

- les figures 13 et 14 sont des schémas montrant les profils de deux fonctions à niveaux de gris entre lesquelles doivent être interpolées des fonctions intermédiaires, le sous-graphe et le sur-graphe étant représentés avec des hachures différentes, correspondant à des étiquettes diffé¬ rentes ; et la figure 15 montre, en noir, la partie variable des sous-graphes et sur-graphes des figures 13 et 14, ainsi que le SKIZ des parties communes, qui constitue le graphe de la fonction interpolée médiane entre les fonctions à niveaux de gris.

On décrira tout d'abord l'invention dans son application à des images dites "mosaïques", dans lesquelles chaque point d'image ou pixel est affecté d'une étiquette indiquant la classe, parmi M classes, à laquelle appartient le point, au lieu d'un paramètre représentant directement une grandeur physique. II existe déjà de nombreuses méthodes de segmentation pemettant de produire, à partir d'une image, une telle image mosaïque. Chaque classe de la mosaïque représente une partie de l'image présentant certains caractères d'homogénéité, en ce qui concerne par exemple le niveau de gris, la couleur ou la texture. Certaines techniques de compression d'images utilisent cette représentation pour coder les images : les contours de la mosaïque sont codés et transmis séparément, puis un modèle de texture est transmis pour chaque classe de la mosaïque. De cette manière le décodeur peut générer d'abord l'image mosaïque puis "habiller" chacune des facettes ou cellules.

Un cas particulièrement simple des images mosaïque est l'image binaire, qui ne comporte que deux ensembles, dont l'un peut être regardé comme le fond et l'autre comme la forme.

La mise en oeuvre de l'invention nécessite le calcul de distances géodésiques dont on rappellera la définition.

Dans un ensemble X de points x. et x₂, la distance géodésique entre x, et x₂ dans X est la longueur du plus court chemin entre Xj et x₂ qui soit contenu dans X. Si x. et x₂ appartiennent à des composantes connexes distinctes de X, il n'existe aucun chemin inclus dans X reliant x. et x₂, on convient que la distance géodésique entre x, et x₂ vaut +•. Si X couvre tout l'espace, la distance géodésique se réduit à la distance simple ou euclidienne.

Dans un ensemble X incluant un point x et un ensemble Y, la distance géodésique dans X entre x et Y est égal à la plus courte distance géodésique dans X entre x et un point appartenant à Y. Elle se confond avec la distance eucli- dienne si une droite constitue le plus court chemin répon¬ dant à cette condition.

Un mode simple de mise en oeuvre de l'invention, lorsqu'on accepte certaines contraintes, fait intervenir la génération d'une image mosaïque dérivée qu'on désignera par le terme de squelette par zone d'influence ou SKIZ. Si on considère une famille d'ensembles

W = (X_j) _lβl>M, chaque ensemble ayant une étiquette différente et pouvant être constitué par des facettes ou cellules regroupant les points affectés à la même classe, le SKIZ de la famille d'ensembles W est une nouvelle image mosaïque formée par une autre famille d'ensembles :

Chaque ensemble Z, est l'ensemble des points du plan qui sont plus proches de l'ensemble Xj que de tout autre ensemble X_kf où k est différent de i. L'ensemble Z_s contient l'ensemble X_t et porte la même étiquette. On écrira :

Z- ⁼ - Z_W(Xi| De façon générale, l'interpolation entre deux images U₀=A et U_N=B revient à chercher une famille d'interpolés :

L'interpolation se fait en deux phases :

(a) on recherche les parties supposées communes à tous les interpolés en partant de l'hypothèse que ce sont les parties communes à A et B. Tout point qui n'appartient pas à une partie commune sera affecté à une classe au cours de la phase (b) ;

(b) on affecte chaque point des parties variables en fonction de la valeur des fonctions distance géodésique entre ce point et des régions prédéfinies de l'image, telles que parties communes ou parties possédant une étiquette prédéterminée ; c'est à ce stade qu'intervient le SKIZ, dans un mode simplifié de réalisation.

CAS OENERAT. DE DEUX TMAOES MOSAÏQUES (Y COMP TS RTUATRI,S.

Au cours de l'étape (a) on compare les images A et B.

Par exemple les figures 1 et 2 montrent des images A et B dans lesquelles les pixels sont répartis entre quatre classes. Dans l'interpolation, on traite les facettes pour lesquelles il y a une correspondance bi-univoque entre les facettes ou cellules A_λ de l'image A et B_A de l'image B (λ désignant une même étiquette dans les deux images) . Les facettes qui apparaissent ou disparaissent sont traitées par des processus classiques. Sur les figures 1 et 2, des lignes en tirets indiquent la segmentation de celle des images qui n'est pas représen¬ tée ; elles font apparaître les parties communes, ayant la même étiquette dans les deux images mosaïques. Elles conserveront la même étiquette dans toutes les images interpolées. On notera par :

- AŒ, l'image mosaïque constituée par :

- les parties communes, représentées en noir sur les figures 3 et 3 et sur la figure 5 par le même rem¬ plissage que sur les figures 1 et 2 ; elles seront désignées par (AfïB.j pour les points d'étiquette λ ; les parties variables, représentées sur les figures 3 et 4 de la même façon que sur les figures 1 et 2. Les parties communes sont conservées dans toutes les images intermédiaires 1,..., N-l.

MÉTHODE SIMPLIFI E = utilisation du SKIZ

La construction de l'interpolé médian i = N/2 à l'aide du SKIZ est simple et rapide. On effectue le SKIZ de la famille des parties fixes munies chacune de son étiquette.

Cette mosaïque milieu M(A,B) constitue l'interpolé médian.

L'algorithme pour la construction du SKIZ d'une famille d'ensembles peut être celui décrit par Luc Vincent dans "Algorithmes morphologiques à base de files d'attente et de lacets. Extension aux graphes", Thèse de l'Ecole des Mines de Paris, 14 mai 1990.

La construction des interpolés intermédiaires peut s'effectuer directement à partir des images A et B ou en séquence en utilisant chaque fois un interpolé intermédiaire déjà disponible.

Dans le second cas, la séquence de N-l images résultant des N étapes de déformation pour passer de U₀ = A à U_N ≈ B implique des itérations successives, d'abord : U_N/2 = Moyenne (U₀,U_N) puis ι^' U_N/ = Moyenne (U₀,U_N/2)

C U_3N/4 = Moyenne (U_N/2,U_N)

et pour chaque itération d'ordre i :

U_N/2 ^l = Moyenne (U₀.U_M/2 ^,_1) Si N est une puissance de 2, on obtient une séquence régulière d'images interpolées.

A titre d'exemple, les figures 4 et 5 montrent les deux premiers de sept interpolés entre les images des figures 1 et 2, les limites des facettes de la figure 1 étant montrées en tirets.

UTILISATION DE PARAMETRES DTSTANOES (.EODESTOTTES PT.TT.q COMPLETS :

La phase (b) peut également être effectuée en utilisant les distances géodésiques et des images "différence" entre les mosaïques U₀=A et U_N=B, par un procédé n'introduisant aucune contrainte sur le nombre des interpolés et s'appuyant sur des distances géodésiques supplémentaires.

L'ensemble "différence" entre A et B est l'ensemble des points d'étiquette différente dans A et B (en noir sur la figure 5) . Cette différence sera désignée par (A/B) et on notera (A/B)_λ les points ayant l'étiquette λ dans la mosaïque A sans l'avoir dans la mosaïque B. Chaque point des parties variables possède donc une étiquette différente dans les mosaïques A et B.

On définit de plus, en chaque point x, deux distances géodésiques pour chaque mosaïque A et B. Une première distance géodésique pour chaque point, appartenant à un ensemble d'étiquette λ, est la distance au point le plus proche qui n'a pas la même étiquette dans la même mosaïque. Etant donné qu'il s'agit d'une distance, dans la même mosaïque A ou B, à une partie qui ne possède pas la même étiquette, on appellera cette distance : d^A ₀-_t pour un point dans la mosaïque A, d^B _out pour un point dans la mosaïque B. On peut considérer que d*_out par exemple est la distance géodésique du point d'étiquette à A_λ dans

La figure 8 montre la distance d^A _out pour deux points x et y appartenant à deux ensembles différents dans la mosaïque A.

Une seconde distance géodésique, pour chaque point appartenant à un ensemble d'étiquette λ, est définie pour chaque point n'appartenant pas aux parties communes At_3. Elle sera désignée par : d^A _in dans la mosaïque A d^B _in dans la mosaïque B

La distance d^A _ιn pour un point qui a l'étiquette λ dans A est la distance géodésique au point le plus proche ayant la même étiquette λ dans B, c'est-à-dire la distance géodésique de ce point à (AfB)_λ. Elle serait nulle pour les parties communes. La figure 8 montre, à titre d'exemple, les distances d^A _in pour deux points x et y. De même, la distance d^B _in pour un point qui a l'étiquette λ dans B est la distance géodésique à (AnB)_i.

Les figures 6 et 7 montrent les distances d^B _in et d^A _in pour les mêmes points x et y que ci-dessus.

On obtient ainsi, pour chacune des mosaïques A et B, deux fonctions distances : d^A _in et d^A _out pour A, d^B _ln et d^B _out pour B, calculées pour chaque point.

On utilise alors un ou plusieurs de critères d'affecta¬ tion fondés chacun sur une comparaison entre deux des valeurs d ci-dessus. Les couples de valeurs dj et d₂ que l'on compare peuvent être parmi les suivants, pour un même point x :

à. (1;

d^A _in d^B _in d^Bo_. d\_u- d^A _in d^A _ou- d^B _{i n} d^B _out

Le premier terme de chaque couple constitue une indica¬ tion sur la distance à la mosaïque A, le second sur la distance à la mosaïque B.

Application d'un -.---I -r-r,-.-.--

Si une des distances ci-dessus est particulièrement représentative à elle seule, on peut utiliser une règle simple.

Dans l'interpolé médian on affecte, à un point x : - l'étiquette qu'il a dans A si d₂(x) > d. (x)

- l'étiquette qu'il a dans B si d₂(x) ≤ d, (x)

Plus généralement, dans un interpolé autre que l'inter¬ polé médian, on est amené à effectuer une pondération. Elle peut se faire en comparant la valeur prise par la fonction : f(d,,d₂) = d₂ / (d, + d₂) à un seuil α compris entre 0 et 1, dont la valeur croît avec le rapport entre le rapport des distances temporelles de l'interpolé aux mosaïques initiale et finale. Si la valeur de f(d-,d₂) est inférieure à α, on affecte à x l'étiquette qu'il a dans la mosaïque A et inversement.

La figure 9 montre, à titre d'exemple, des lignes de niveau (correspondant à une comparaison avec trois seuils différents) pour trois interpolés de l'une seulement des facettes de l'image ; ces lignes de niveau délimiteront la facette dans les interpolés. Les figures 10 à 12 montrent trois images interpolées successives.

Annlica ion d'une pondérât-ion de plusieurs critères Chacun des couples (d.,d₂) ci-dessus permet de construire une séquence d'interpolés différente et toute some pondérée des résultats permet de construire des interpolés encore différents.

Dans de nombreux cas, de bons résultats en matière de progressivité sont obtenus en faisant simplement la moyenne de deux fonctionnelles, par exemple : l /2 [d^A _in , d^A _out ) , (d^B _in , d^B _out ) ]

c'est-à-dire d'une comparaison , sur une somme pondérée faisant intervenir deux fois d_in et d_out.

EXTENS ION ATTX FONCTIONS N.TMÉR TOTIES

Jusqu'ici il n'a été fait référence qu'à des images mosaïques dont les facettes ont implicitement la même étiquette : elles peuvent être regardées comme à deux dimen¬ sions. Des images à plus de deux dimensions spatiales peuvent cependant être traitées de la même façon.

Le procédé suivant l'invention est également utilisable dans le cas d'une image à deux - ou plus généralement n dimensions - ayant également des niveaux de gris, où chaque pixel, de coordonnées x, y dans le cas de deux dimensions, est affecté d'une fonction à niveaux de gris g (x,y) . Une telle fonction peut être considérée, dans une image à deux dimensions géométriques, comme une image mosaïque particu- lière dans un espace à trois dimensions x,y,z, où z est le niveau de gris.

Dans ce cas, on affectera à tout point de l'espace tridimensionnel (x,y,z) une étiquette égale à 1 si g (x,y) ≤ z et 0 dans le cas contraire. A titre d'exemple les figures 13 et 14 montrent deux fonctions, représentant deux images mosaïques, entre lesquelles on veut interpoler. Dans un but de simplifica¬ tion, seule une variation en x, représentant la variation spatiale, est figurée. Les points où g(x)fc z constituent le sur-graphe (au-dessus de la courbe en traits pleins, qui correspond à l'image A) et les points où g(x) ≤ z consti¬ tuent le sous-graphe, montré différemment sur les figures. Les interpolés dans les parties variables peuvent être déterminés de la même façon que dans le cas simplement bidimensionnel.

Une solution simple lorsque l'on souhaite déterminer des interpolés dont le nombre est une puissance de deux, plus un, consiste à adopter la méthode du SKIZ déjà décrite. Le SKIZ conduit d'abord à l'interpolé médian, indiqué sur la figure 15 ; les interpolés intermédiaires peuvent être ensuite déterminés successivement à partir de ceux déjà disponibles.

Dans Le cas où on accepte de n'utiliser que le SKIZ appliqué à des distances euclidiennes, la recherche d'un interpolé médian entre deux fonctions gl et g2 de variables x et y peut s'effectuer de la même façon suivante.

On recherche tout d'abord les frontières haute zo = inf (gl,g2) de la partie basse commune et w₀ = sup (gl,g2) de la partie haute commune. Elles délimitent la zone montrée en noirsur la figure 15. Puis on applique un algorithme de dilatation des parties communes qui réduit progressivement l'intervalle entre les frontières jusqu'à ce qu'elle coïncident .

L'algorithme de dilatation utilise un élément structu- rant dont la forme dépendra de la trame de répartition des pixels. Si cette trame est hexagonale, on utilise avanta¬ geusement une prisme hexagonal de hauteur unité ou une pyramide à base hexagonale centrée sur le point à partir duquel s'effectue la dilatation. Classiquement, on applique les fonctions de morphologie: DILATER (z₀) pour faire monter la frontière z₀ ERODER (w₀) pour faire descendre w₀ Il faut limiter la croissance de la fonction dilatée et la décroissance de la fonction érodée. Pour cela on peut utiliser deux contraintes.

L'une est que la fonction basse dilatée ne doit pas dépasser la valeur de la fonction haute pour les mêmes valeurs x,y. Elle est donc finalement z'_Q : z'o = inf [dilate (z₀) , w₀] L'autre est que la fonction haute érodée ne doit pas être inférieure à la fonction basse initiale au même point : w'o = sup [(erode (w_σ) , z₀] La même contrainte est répétée, à chaque étape. Par exemple, à l'itération d'ordre i : z'i = inf [dilate (z⁹^.), sup [erode(W j..) , z'j-,]] w'i = sup [erode(W i,.) , inf [dilate(z' j_.) , Wj.,]]

COMPENSATION DE MOTTVEMENT

En règle générale, les images initiale A et finale B représentent la même scène à deux instants successifs.

On connaît par ailleurs des procédés d'estimation de mouvement d'images suivantes, utilisés notamment pour déterminer le vecteur mouvement moyen lors du codage d'une image de télévision numérique, permettant de réduire le débit d'informations à transmettre. Dans le cas de la présente invention, il est possible, une fois une image de départ ou d'arrivée segmentée, d'estimer le vecteur mouve¬ ment de chaque facette.

Pour interpoler l'image mosaïque qui se trouve à une distance temporelle α de l'image mosaïque initiale A et

1 - α de l'image temporelle B, on applique à chaque facette pour laquelle le vecteur mouvement est T de A à B, une translation αT dans l'image A et une translation (α-1) (T) a la facette dans la mosaïque B. Plus fréquemment, on déterminera un seul vecteur mouvement T pour toutes les facettes constituant un même ensemble d'étiquettes λ.

Au cours de l'une et l'autre des translations, un pixel x peut avoir plusieurs antécédents d'étiquettes différentes. Dans ce cas on choisit l'étiquette permettant la plus grande intersection avec la translation depuis l'autre image mosaïque terminale.

On obtient ainsi deux images mosaïques, l'une à partir de A et l'autre à partir de B. Chacune de ces réalisations peut avoir des trous, dus à la présence de points qui n'ont pas été couverts à la suite des translations. Ces trous peuvent être remplis en appliquant la fonction SKIZ. Les formes des facettes sont alors légèrement modifiées. De plus, il peut y avoir une modification due à ce que le mouvement réel pour passer de A à B n'est pas une simple translation.

On applique ensuite, aux réalisations OA (α-1) B une interpolation du genre décrit ci-dessus, mais appliquée à A et B. HABILLAGE DES FACETTES PAR DES TEXTURES

Les facettes peuvent présenter des textures particuliè¬ res, identifiables par une information pour chaque ensemble dans les images initiale A et finale B.

Il faut également calculer dans ce cas 1 ' interpolation des textures pour les images interpolés.

Dans le cas d'une interpolation avec translation, une stratégie d'habillage d'une facette d'un interpolé par une texture (par exemple lignes parallèles, quadrillage, etc..) est la suivante. On déplace la texture d'une facette de l'ouvrage A sur la mosaïque par la translation OT pour constituer la mosaïque à distance α de A. La translation est validée si le pixel de destination possède, dans l'image mosaïque à distance α, la même étiquette que dans l'image de dépar . Des zones périphériques de certaines facettes de la mosaïque peuvent alors ne pas se voir affecter une texture par cette opération. On leur attribue alors les niveaux de gris des points les plus proches de même étiquette.

On obtient ainsi une première réalisation d'une image codée à distance α de A et (1- α) de B.

Une seconde image interpolée peut être obtenue à partir de l'image B, cette fois par translation égale à (α-1) T.

L'image codée finale peut alors être obtenue en prenant le barycentre entre les deux réalisations ainsi obtenues.

Claims

REVENDICATIONS

1. Procédé de génération de N-l images mosaïques intermédiaires ( j), ,,^., entre deux images successives U₀ = A et U_N = B d'une même scène, fractionnées en cellules appartenant chacune à une classe homogène parmi M classes, suivant lequel :

(b) on compare deux à deux les pixels correspondants des deux images et on identifie ceux des pixels correspondants des deux images qui appartiennent à la même classe, consti¬ tuant une partie commune, de façon à constituer au plus M ensembles X, Y, ..., chacun constitué de pixels d'une même classe ; et

(c) on affecte, dans chaque image intermédiaire générée, les pixels _.de la partie variable chacun à une des M classes, suivant la valeur d'une des distances entre le pixel respectif et des région prédéterminées de l'image, consti¬ tuées par la partie commune ou des parties appartenant à une classe prédéterminée dans une ou des images pré-existantes.

2. Procédé selon la revendication 1, caractérisé en ce que, au cours des étapes (c) , on effectue des dichotomies successives, la première étape (c) permettant de déterminer, pour chaque pixel des parties variables, celui des M ensembles ou classes auquel il a la distance géodésique la plus courte et d'affecter le pixel à cet ensemble et les étapes suivantes utilisent l'image ainsi créée pour recom- mercer l'opération.

3. Procédé selon la revendication 1, caractérisé en ce qu'on affecte successivement les pixels de la partie vo dans l'image médiane U_(N/2), dans les images intermédiaires U_(N/4) et ^U _{(3 I} ^et ainsi de suite jusqu'à U₂, ... , U^, par affectation de chaque pixel de la partie variable des deux images adjacentes déjà disponibles à celui des M ensembles consti¬ tuant la partie commune de ces deux images avec lequel il a la distance géodésique la plus courte.

4. Procédé selon la revendication 1, caractérisé en ce que, au cours de la phase (c) :

(Cl) on détermine, pour chaque point, dans l'image initiale A et dans l'image finale B, la distance géodésicrue respective d^A _out ou d^B _out de ce point au point le plus proche appartenant à une classe différente, et/ou (C2) on détermine, pour chaque point, dans l'image initiale A et dans l'image finale B, la distance géodésique d^A _ιn ou d^B _ιn au point le plus proche ayant la même étiquette dans l'autre image, puis

(C3) pour affecter chaque point à un ensemble dans une image interpolée on compare les deux termes, pour ce point, d'au moins un des deux couples suivants : d^A _in et d^B _in d^B _out et d^A _out d^A _in et d^A _out d^B _{l n} et d^B _out

5. Procédé selon la revendication 4, caractérisé en ce que, pour les interpolés autres que l'interpolé médian, on effectue une pondération des deux termes avant comparaison.

6. Procédé selon la revendication 4 ou 5, caractérisé en ce que, au cours de l'étape (Cl), on effectue une pondéra¬ tion entre plusieurs comparaisons.

7. Procédé selon la revendication 6, caractérisé en ce que la pondération consiste à faire les moyennes des comparaisons entre : d^A _in et d^A _out d^B _ln et d^B _out

8. Procédé selon l'une quelconque des revendications 1 à 7, suivant lequel on estime le vecteur mouvement T de la scène de l'image U₀ à l'image U_N et, avant détermination d'une image intermédiaire située à une fraction α de l'intervalle de temps entre U₀ et U_N, on applique une translation αT à l'image mosaïque U₀ et (l-α)T à l'image U_N.

9. Procédé selon l'une quelconque des revendications 1 à 7 , caractérisé en ce que l'étiquette est prévue pour différencier des facettes suivant au moins un paramètre radiométrique de l'image, tel que luminance, chrominance et texture.

10. Procédé selon les revendications 2 et 9, caractérisé en ce que les opérations de dichotomie successives sont effectuées en utilisant également la différentiation fournie par l'étiquette.