FR 3073068 A1 20190503 FR 1760144 A 20171027 L'invention concerne un procédé de restauration d'images d'une séquence d'images, et un dispositif mettant en oeuvre ledit procédé. Certaines applications utilisant des images ou des séquences d'images nécessitent une très haute qualité d'images. C'est le cas par exemple d'applications de surveillance où une faible qualité d'images peut entraîner de mauvaises interprétations d'une image et, par exemple, des déclenchements de fausses alarmes. La restauration d'images consiste à améliorer une qualité d'une image en lui appliquant diverses techniques de traitement d'images telles que des techniques de réduction ou de suppression de bruit, d'amélioration de contours, d'amélioration de contraste, ... Une image d'une séquence d'images (i.e. d'une séquence vidéo), est un cas particulier d'image puisque cette image est généralement proche visuellement d'images temporellement voisines dans la séquence d'images. Un dispositif d'acquisition de séquences d'images, tel qu'une caméra vidéo, possède en effet généralement une fréquence d'acquisition d'images suffisamment élevée pour que des corrélations temporelles subsistent entre des images successives d'une séquence d'images. Utilisées à bon escient, ces corrélations temporelles peuvent permettre à une image de bénéficier d'améliorations apportées à une ou plusieurs images voisines grâce à des méthodes de restauration d'images. Toutefois, identifier des corrélations temporelles existant entre des images n'est pas toujours une tâche facile. C'est notamment le cas lorsqu'une séquence d'images représente au moins un objet en mouvement sur un fond lui-même en mouvement, le fond pouvant alors être considéré comme un objet. En effet, afin de pouvoir bénéficier efficacement des corrélations temporelles existant entre deux images, il est souhaitable de prendre en compte les mouvements de l'objet et les mouvements du fond. Des méthodes connues de restauration d'images comportant des objets en mouvement comprennent une phase d'analyse de mouvement. Cette phase d'analyse de mouvement permet de mettre en correspondance des pixels corrélés entre deux images. Une amélioration apportée à un premier pixel d'une première image peut alors bénéficier à un second pixel d'une seconde image corrélé avec le premier pixel. La phase d'analyse de mouvement utilise généralement des techniques de flot optique comprenant une estimation d'un champ dense de mouvement dans lequel chaque pixel 2 (ou au moins un sous-ensemble de pixels) d'une image à restaurer est associé, par un vecteur de mouvement, à un pixel d'une image de référence par rapport à laquelle le mouvement est déterminé. Ces techniques de restauration d'images sont considérées comme insuffisamment fiables pour permettre une restauration d'images de qualité. De 5 plus, une estimation d'un champ dense de mouvement a un coût calculatoire important. Un tel coût calculatoire n'est pas compatible avec certains systèmes embarqués ou portables susceptibles de mettre en oeuvre des méthodes de restauration d'images, telles que des caméras vidéo, des jumelles numériques, des lunettes à réalité augmentée, ... On préfère donc en général remplacer, lors de la phase d'analyse de mouvement, les 10 techniques d'estimation d'un champ dense de mouvement par des techniques d'analyse de mouvement ayant un coût calculatoire moins élevé, telle que les techniques d'estimation de mouvement global. Les techniques d'estimation de mouvement global sont particulièrement efficaces lorsqu'une hypothèse de stationnarité des scènes représentées par la séquence d'images est vérifiée. 15 Par ailleurs, suite à l'analyse de mouvement, une même méthode de restauration d'images est généralement appliquée à chaque pixel de l'image à restaurer. Or, tous les pixels d'une image n'ont pas un intérêt égal. Il est courant en effet qu'une image représente un objet et un fond. Obtenir un rendu visuel détaillé de l'objet est utile, par exemple, lorsque l'objet doit être identifié précisément, alors que, dans le meilleur des 20 cas, le fond n'est utile que pour situer l'objet dans un contexte. Il est connu des méthodes de restauration d'images dans lesquelles la phase d'analyse de mouvement est suivie d'une segmentation de l'image. La segmentation de l'image permet de diviser l'image en sous-parties homogènes suivant un critère prédéfini. Connaissant les sous-parties constituant l'image, il est possible d'appliquer un traitement d'images plus ou moins 25 efficace (et donc plus ou moins complexe) en fonction de l'intérêt de chaque sous- partie. Lorsque la segmentation suit une estimation d'un champ dense de mouvement, cette segmentation peut s'appuyer sur le champ dense obtenu pour diviser l'image en sous-parties homogènes au sens du mouvement. Cette méthode de segmentation peut toutefois donner des résultats approximatifs, notamment lorsque le champ dense de 30 mouvement est bruité. D'autres procédés de segmentation existent comme par exemple les procédés de segmentation basés contours actifs. Toutefois, ces procédés sont généralement destinés à des images fixes et, lorsqu'ils sont appliqués à des images d'une séquence d'images, ils utilisent peu (ou pas) les corrélations temporelles existant entre les images. 3 Il est souhaitable de pallier ces inconvénients de l'état de la technique. Il est notamment souhaitable de proposer un procédé et un dispositif de restauration d'images adapté aux images comprenant des objets en mouvement. Il est souhaitable de plus que ce procédé et ce dispositif soient adaptés aux images où le fond 5 de l'image est lui-même en mouvement. Il est enfin souhaitable que ledit procédé ait un coût calculatoire faible et puisse être mis en oeuvre par un système ayant de faibles capacités calculatoires tel qu'un système embarqué ou portable. Selon un premier aspect de l'invention, l'invention concerne un procédé de restauration d'images d'une séquence d'images comprenant, lorsqu'il est appliqué à une 10 première image de la séquence d'images : estimer une information représentative d'un mouvement global d'un arrière-plan de la première image par rapport à une deuxième image ; compenser le mouvement global de l'arrière-plan dans la deuxième image en utilisant ladite information représentative du mouvement global de l'arrière-plan afin 15 d'obtenir une version recalée de la deuxième image, dite deuxième image recalée ; obtenir un contour d'un objet de la première image en appliquant un procédé de segmentation, ledit procédé de segmentation étant itératif et comprenant, lors d'une itération, une modification d'un contour de l'objet dans la première image obtenue lors d'une itération précédente dudit procédé de segmentation, dit contour précédent, de 20 sorte à obtenir un contour de l'objet dans la première image, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent, un contour final de l'objet étant obtenu lorsqu'une condition prédéfinie d'arrêt dudit procédé de segmentation est remplie, le coût d'un contour de l'objet dans la première image étant une somme entre une première valeur représentative d'une énergie interne audit contour 25 et une seconde valeur représentative d'une énergie externe audit contour, l'énergie externe audit contour étant fonction d'au moins une énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image recalée et d'une énergie, dite énergie de contour, correspondant à une somme de valeurs de modules de gradient calculées pour des pixels d'un deuxième ensemble de pixels 30 appartenant au contour courant de l'objet ; une valeur représentative de l'énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image étant calculée sous forme d'une somme de différences entre des valeurs représentatives de pixels d'un premier ensemble de pixels de la première image appartenant au contour courant et des valeurs représentatives de pixels situés aux 4 mêmes positions spatiales que les pixels du premier ensemble de pixels dans la seconde image ; utiliser le contour de l'objet ainsi obtenu pour estimer une information représentative d'un mouvement global de l'objet; et, 5 appliquer à la première image une méthode de restauration d'images utilisant les informations représentatives du mouvement global de l'arrière-plan et du mouvement global de l'objet estimées. Utiliser une énergie externe intégrant une énergie dépendante d'un mouvement global de l'objet permet de prendre en compte le mouvement de l'objet dans la 10 segmentation du contour de l'objet. On tire ainsi profit des corrélations existant entre deux images pour améliorer la segmentation. Dans un mode de réalisation, pour calculer la valeur représentative de l'énergie interne au contour courant, une dérivée première locale et une dérivée seconde locale du contour sont calculées pour des pixels d'un troisième ensemble de pixels de l'image 15 courante appartenant au contour courant de l'objet, ladite valeur représentative de l'énergie interne étant fonction desdites dérivées calculées. Dans un mode de réalisation, les premier, deuxième et troisième ensembles de pixels sont identiques, et chaque ensemble comprend au moins une sous-partie des pixels de l'image courante appartenant au contour courant de l'objet. 20 Dans un mode de réalisation, lors d'une première itération dudit procédé, un contour initial de l'objet dans l'image courante est obtenu à partir d'un contour final obtenu lors d'une application du procédé de segmentation à l'image de référence ou d'un contour spécifié par un opérateur dans l'image de référence. Dans un mode de réalisation, lors de chaque estimation d'une information 25 représentative d'un mouvement global, une information représentative de la forme et de la position de l'objet est obtenue, ladite information représentative de la forme et de la position de l'objet étant utilisée pour masquer des pixels ne devant pas être pris en compte dans ladite estimation. Dans un mode de réalisation, suite à l'estimation de ladite information 30 représentative du mouvement global de l'objet, dite première information, un filtrage est appliqué à ladite première information afin de garantir des variations régulières du mouvement de l'objet entre deux images successives de la séquence d'images, ledit filtrage comprenant les étapes suivantes : déterminer une première matrice permettant d'estimer un mouvement de l'objet dans un référentiel centré sur un barycentre de 5 l'objet dans la première image et une deuxième matrice permettant d'estimer un mouvement de l'objet dans un référentiel centré sur un barycentre de l'objet dans la deuxième image recalée ; utiliser la première et la deuxième matrices pour calculer une information représentative du mouvement de l'objet, dite seconde information, à partir 5 de ladite première information; utiliser la seconde information pour obtenir une troisième matrice représentative de composantes de translation du mouvement de l'objet ; utiliser la seconde information et la troisième matrice pour obtenir une quatrième matrice représentative de composantes du mouvement de l'objet autres que les composantes de translation ; obtenir une version filtrée de la troisième matrice, dite 10 troisième matrice filtrée, en calculant une somme pondérée entre la troisième matrice et une troisième matrice filtrée obtenue lors de la mise en oeuvre dudit procédé sur la deuxième image ; obtenir une version filtrée de la quatrième matrice, dite quatrième matrice filtrée, en calculant une somme pondérée entre la quatrième matrice et une quatrième matrice filtrée obtenue lors de la mise en oeuvre du procédé sur la deuxième 15 image ; et, obtenir une information représentative d'un mouvement global filtré de l'objet en utilisant la première et la deuxième matrices, la troisième matrice filtrée et la quatrième matrice filtrée. Dans un mode de réalisation, la seconde information est calculée de la manière suivante : 20 dHicobjet = 1/171. dHkobjet. où Vk est la première matrice, Vk_i est la deuxième matrice, dHk°bietla première information et delet la seconde information. Dans un mode de réalisation, la troisième matrice est calculée de la manière suivante : 25 Hk = ApproxT(Vjl. dHkObjet où Hf, est la troisième matrice et ApproxT(X) est une approximation en translation de la matrice homographique X. Dans un mode de réalisation, la quatrième matrice est calculée de la manière 30 suivante : _ Ht-i dwuobjet , "k k - v k - ik - v k-1 OÙ Hikl est la quatrième matrice. Dans un mode de réalisation, la troisième matrice filtrée est calculée de la manière suivante : 6 utFilt ut Filt "k a- "k-1 + (1 - a)Hkt tFilt t Filt où Hk est la troisième matrice filtrée, Hk_i est la troisième matrice filtrée obtenue lors de la mise en oeuvre dudit procédé sur la deuxième image et a est une constante prédéfinie comprise entre « 0 » et « 1 ». 5 Dans un mode de réalisation, la quatrième matrice filtrée est calculée de la manière suivante : Filt Hkh = I + (1- le)Hlk2 où Hik2Filt est la quatrième matrice filtrée, I est une matrice identité et /3 est une constante prédéfinie comprise entre « 0 » et « 1 ». 10 Dans un mode de réalisation, l'information représentative d'un mouvement global filtré de l'objet est calculée de la manière suivante : ,4"objetFut = vk. Hi t,Fut.Hik2Fut. "'"k Selon un deuxième aspect de l'invention, l'invention concerne un dispositif de restauration d'images d'une séquence d'images comprenant : 15 des moyens d'estimation pour estimer une information représentative d'un mouvement global d'un arrière-plan d'une première image de la séquence d'images par rapport à une deuxième image de la séquence d'images ; des moyens de compensation de mouvement pour compenser le mouvement global de l'arrière-plan dans la deuxième image en utilisant ladite information 20 représentative du mouvement global de l'arrière-plan afin d'obtenir une version recalée de la deuxième image, dite deuxième image recalée ; des moyens d'obtention permettant d'obtenir un contour d'un objet dans la première image recalée, comprenant des moyens de modification d'un contour de l'objet dans la première image obtenue lors d'une itération précédente, dit contour 25 précédent, permettant d'obtenir un contour de l'objet dans la première image, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent, un contour final de l'objet étant obtenu lorsqu'une condition prédéfinie d'arrêt est remplie, le coût d'un contour de l'objet dans la première image étant une somme entre une première valeur représentative d'une énergie interne audit contour et 30 une seconde valeur représentative d'une énergie externe audit contour, l'énergie externe audit contour étant fonction d'au moins une énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image recalée et d'une énergie, dite énergie de contour, correspondant à une somme de valeurs de modules de gradient calculées pour des pixels d'un deuxième ensemble de pixels appartenant au contour 7 courant de l'objet ; une valeur représentative de l'énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image étant calculée sous forme d'une somme de différences entre des valeurs représentatives de pixels d'un premier ensemble de pixels de la première image appartenant au contour courant et des valeurs 5 représentatives de pixels situés aux mêmes positions spatiales que les pixels du premier ensemble de pixels dans la seconde image ; des moyens d'estimation de mouvement utilisant le contour de l'objet dans la première image pour estimer une information représentative d'un mouvement global de l'objet entre la première image et la deuxième image ; et, 10 des moyens de restauration d'images utilisant les informations représentatives du mouvement global de l'arrière-plan et du mouvement global de l'objet estiméesSelon un troisième aspect de l'invention, l'invention concerne un programme d'ordinateur, comprenant des instructions pour mettre en oeuvre, par un dispositif, le procédé selon le premier aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif 15 Selon un quatrième aspect de l'invention, l'invention concerne des moyens de stockage, stockant un programme d'ordinateur comprenant des instructions pour mettre en oeuvre, par un dispositif, le procédé selon le premier aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif. Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, 20 apparaîtront plus clairement à la lecture de la description suivante d'un exemple de réalisation, ladite description étant faite en relation avec les dessins joints, parmi lesquels : - la Fig. 1 illustre schématiquement un exemple de contexte dans lequel l'invention peut être mise en oeuvre ; 25 - la Fig. 2 illustre schématiquement un exemple d'image contenant un objet et un fond ; - la Fig. 3 illustre schématiquement un exemple d'architecture matérielle d'un module de traitement compris dans un système d'acquisition d'images ; - la Fig. 4 illustre schématiquement un procédé de restauration d'images selon 30 l'invention ; - la Fig. 5 illustre schématiquement un procédé de segmentation d'images selon l'invention ; - la Fig. 6 illustre schématiquement un procédé de calcul d'un coût d'un contour d'un objet compris dans le procédé de segmentation selon l'invention ; 8 - la Fig. 7 illustre schématiquement un procédé d'estimation de mouvement global ; et, - la Fig. 8 illustre schématiquement un procédé de filtrage de mouvement. L'invention est décrite par la suite dans un contexte où un système de visualisation 5 comprend un dispositif d'acquisition d'images, un module de traitement et un dispositif de visualisation d'images. L'invention peut toutefois être mise en oeuvre dans un contexte où le dispositif d'acquisition d'images, le module de traitement et le dispositif d'affichage sont séparés et distants géographiquement. Dans ce cas, le dispositif d'acquisition d'images, le module de traitement et le dispositif de visualisation 10 d'images comprennent des moyens de communication pour communiquer entre eux. Par ailleurs, le procédé selon l'invention s'appuie sur un procédé de segmentation basé contours actifs basé contours. Nous montrons par la suite que d'autres types de procédés de segmentation basés contours actifs peuvent être utilisés tels que par exemple, les procédés de segmentation basés contours actifs basés région, les procédés 15 de segmentation basés contours actifs implicites basés ensembles de niveaux (« level sets » en terminologie anglo-saxonne),... De plus, les images utilisées dans le cadre de l'invention sont essentiellement des images monochromes où chaque pixel d'une image ne possède qu'une composante. L'invention peut toutefois être appliquée à des images multi-composantes dans 20 lesquelles chaque pixel d'une image possède une pluralité de composantes. La Fig. 1 illustre schématiquement un exemple de contexte dans lequel l'invention peut être mise en oeuvre. Dans la Fig. 1, une scène 1 est observée par un système de visualisation 5 suivant un champ optique 3. La scène 1 comprend un objet 4. Le système de visualisation 5 25 comprend un dispositif d'acquisition d'images 51, un module de traitement 52 apte à mettre en oeuvre un procédé de restauration d'images et une méthode de segmentation selon l'invention et un dispositif de visualisation d'images 53. Le dispositif d'acquisition d'images 51 comprend un ensemble optique et un capteur d'images tel que par exemple un capteur CCD (« Charge-Coupled Device » en terminologie anglo- 30 saxonne, dispositif à transfert de charge) ou un capteur CMOS (« Complementarity metal-oxide-semiconductor », en terminologie anglo-saxonne, semi-conducteur métal-oxyde complémentaire). Le capteur d'images fournit une séquence d'images 2 représentatives de la scène 1 selon une fréquence temporelle prédéfinie F (par exemple de F = 25 à 50 images par secondes) au module de traitement 52. Après traitement par 9 le module de traitement 52, le module de traitement 52 fournit des images améliorées au dispositif de visualisation d'images 53. Le dispositif de visualisation d'images 53 est par exemple un écran. Dans un mode de réalisation, les images fournies par le dispositif d'acquisition 5 d'images 51 sont des images monochromes. Dans un mode de réalisation, les images fournies par le dispositif de visualisation d'images 51 sont des images multi-composantes. La Fig. 2 illustre schématiquement un exemple d'image contenant un objet et un fond. 10 L'exemple décrit en relation avec la Fig. 2 représente une image 2 fournie par le système de visualisation 5. On retrouve dans l'image 2, l'objet 4 (ici un véhicule) se déplaçant sur un fond (ici un paysage). On peut imaginer que le système de visualisation 5 est mobile de sorte que le fond possède un mouvement. La Fig. 3 illustre schématiquement un exemple d'architecture matérielle du 15 module de traitement 52 compris dans le système de visualisation 5. Selon l'exemple d'architecture matérielle représenté à la Fig. 3, le module de traitement 52 comprend alors, reliés par un bus de communication 520 : un processeur ou CPU (« Central Processing Unit » en anglais) 521 ; une mémoire vive RAM (« Random Access Memory » en anglais) 522 ; une mémoire morte ROM (« Read Only 20 Memory » en anglais) 523 ; une unité de stockage telle qu'un disque dur ou un lecteur de support de stockage, tel qu'un lecteur de cartes SD (« Secure Digital » en anglais) 524 ; au moins une interface de communication 525 permettant au module de traitement 52 de communiquer avec le dispositif d'acquisition d'images 51 et/ou le dispositif de visualisation d'images 53. 25 Dans un mode de réalisation dans lequel le dispositif d'acquisition d'images 51, le module de traitement 52 et dispositif de visualisation 53 sont séparés et distants, le dispositif d'acquisition d'images 51 et le dispositif de visualisation 53 comprennent aussi une interface de communication pouvant communiquer avec l'interface de communication 525 par l'intermédiaire d'un réseau tel qu'un réseau sans fils. 30 Le processeur 521 est capable d'exécuter des instructions chargées dans la RAM 522 à partir de la ROM 523, d'une mémoire externe (non représentée), d'un support de stockage (tel qu'une carte SD), ou d'un réseau de communication. Lorsque le module de traitement 52 est mis sous tension, le processeur 521 est capable de lire de la RAM 522 des instructions et de les exécuter. Ces instructions forment un programme 10 d'ordinateur causant la mise en oeuvre, par le processeur 521, de tout ou partie du procédé décrit ci-après en relation avec les Figs. 4 à 8. Le procédé décrit ci-après en relation avec les Figs. 4 à 8 peut être implémenté sous forme logicielle par exécution d'un ensemble d'instructions par une machine 5 programmable, par exemple un DSP (« Digital Signal Processor » en anglais),un microcontrôleur ou un GPU (processeur graphique, « Graphics Processing Unit » en terminologie anglo-saxonne), ou être implémenté sous forme matérielle par une machine ou un composant dédié, par exemple un FPGA (« Field-Programmable Gate Array » en anglais) ou un ASIC (« Application-Specific Integrated Circuit » en anglais). 10 La Fig. 4 illustre schématiquement un procédé de restauration d'images selon l'invention. Le procédé décrit en relation avec la Fig. 4 est un procédé itératif, mis en oeuvre par le module de traitement 52 sur chaque image d'une séquence d'images fournie par le dispositif d'acquisition d'images 51, excepté la première image de la séquence 15 d'images. Par la suite, nous appelons image courante, une image en cours de traitement par le module de traitement 52 et nous notons l'image courante Ik, où k représente un indice d'une image. L'indice k indique que l'image Ik apparaît dans la séquence d'images à un instant Tc, + k. i. Tc, correspond au début de la séquence d'images (et donc à la première image de la séquence) et T = F. Soit Ik_1, une image, dite image 20 précédente, précédant immédiatement l'image courante Ik dans la séquence d'images fournie par le dispositif d'acquisition d'images 51. Dans une étape 41, le module de traitement 52 estime une information représentative d'un mouvement global d'un arrière-plan de l'image courante Ik (ou mouvement du fond) par rapport à l'image précédente Ik_1. L'image précédente Ik_i 25 est alors une image de référence pour l'image courante Ik pour l'estimation de l'information représentative du mouvement du fond. Cette étape est mise en oeuvre par une méthode d'estimation de mouvement global. Une estimation de mouvement global fait l'hypothèse qu'un ensemble de pixels d'une même image est animé d'un même mouvement. Ce mouvement peut être simple, tel qu'un mouvement de translation ou de 30 rotation, ou complexe représenté, par exemple, par une transformation affine ou une homographie. Une homographie est une transformation de coordonnées projective à huit paramètres. Dans un mode de réalisation, le module de traitement considère que le mouvement du fond entre deux images successives de la séquence d'images est représenté par une homographie. Soit (x,y) des coordonnées d'un pixel 11 Pfc."`i appartenant au fond dans l'image courante Ik et (x ',y ') des coordonnées du même pixel pci-ei appartenant au fond dans l'image précédente Ik_1. L'estimation de mouvement global mise en oeuvre lors de l'étape 41 consiste à déterminer les huit paramètres d'une homographie permettant de transformer les coordonnées (x ',y ') de 5 chaque pixel pfc.07-1,1 dans l'image précédente Ik_1 en coordonnées (x,y) d'un pixel eond dans l'image courante Ik. En déterminant les huit paramètres de l'homographie, on détermine une information représentative d'un mouvement global du fond entre l'image précédente Ik_i et l'image courante Ik. Nous notons dlirc°nd l'homographie représentative du mouvement global du fond entre l'image précédente Ik_i et l'image 10 Ik. Dans une étape 42, le module de traitement 52 compense le mouvement du fond dans l'image précédente Ik_i afin d'obtenir une image précédente recalée Iirceci. Pour ce faire, le module de traitement 52 applique l'homographie derd trouvée lors de l'étape 41 à l'ensemble des pixels de l'image précédente Ik_i. 15 Dans une étape 43, le module de traitement 52 obtient un contour C de l'objet 4 en appliquant un procédé de segmentation à l'image courante Ik. Nous décrivons par la suite en relation avec les Figs. 5 et 6 un procédé de segmentation selon l'invention. Dans une étape 44, le module de traitement 52 estime une information représentative d'un mouvement global de l'objet 4, entre l'image courante et une image 20 précédente. Comme lors de l'étape 41, le module de traitement 52 considère que le mouvement global de l'objet 4 est représenté par une homographie, notée dHk°bi Dans un mode de réalisation, l'homographie dHk°hiet représentative du mouvement global de l'objet 4 est obtenue en utilisant l'image courante Ik et l'image précédente Ik_i et en prenant en compte le mouvement du fond derd mesuré lors de l'étape 41. Soit 25 (x,y) des coordonnées d'un pixel Clet appartenant à l'objet 4 dans l'image courante Ik et (x ',y ') des coordonnées du même pixel Pnet appartenant à l'objet 4 dans l'image précédente Ik_1. L'estimation de mouvement global mise en oeuvre lors de l'étape 44 comprend une détermination des huit paramètres d'une homographie dl I illcjesurée permettant de transformer les coordonnées (x,y) de chaque pixel Clet dans l'image 30 courante Ik en coordonnées (x',y') du pixel Pnet dans l'image précédente Ik_i. On peut écrire : uF deesurée = ijObjet d ond ueilk - "'"Ic et et L'homographie dHk°bi représentative est alors obtenue de la manière suivante : 12 dHkObjet = deesurée .(dendyl Dans un mode de réalisation, l'homographie dHObjetk représentative du mouvement global de l'objet 4 est mesurée entre une image recalée et une image non recalée, ce qui permet de ne pas faire intervenir l'homographie dllfc°12c1 Par exemple, 5 l'homographie dHic°biet est mesurée entre l'image précédente recalée Iirceci et l'image courante Ik. Dans une étape 45, le module de traitement 52 applique une méthode de restauration d'images à l'image courante Ik. Dans un mode de réalisation, la méthode de restauration d'images appliquée utilise les informations représentatives du 10 mouvement global du fond et les informations représentatives du mouvement global de l'objet 4 estimées lors des étapes 41 et 44 pour mettre en correspondance les pixels de l'image courante Ik et les pixels de l'image précédente Soit Pk un pixel de l'image courante Ik et Pk_1 un pixel de l'image précédente mis en correspondance avec le pixel Pk en utilisant les homographies dlirc°nd et dHkObjet Le pixel Pk (respectivement 15 le pixel Pk_i) possède un nombre entier positif non nul N, de composantes CiPk (respectivement C:3k-i), i E [1; Ne]. Dans ce mode de réalisation, la valeur Ci1,3k de chaque composante de chaque pixel Pk de l'image Ik est remplacée par une somme pondérée calculée de la manière suivante : CpkrPk _L w rPk-i i = 30 vNI w. rP - I Pk = N vI j=1"".1 Wk Wk-1 20 où Wk et Wk_i_ sont des valeurs de pondération prédéfinie. Par exemple, les pondérations Wk et Wk_i_ peuvent être telles que Wk Dans un mode de réalisation, la méthode de restauration d'images appliquée utilise une fenêtre d'images comprenant un nombre N1 d'images précédant l'image courante Ik. Les pixels de chacune des images sont mis en correspondance en utilisant 25 les informations représentatives du mouvement du fond et du mouvement de l'objet 4 obtenues pour chacune des images de la fenêtre d'images. Dans ce mode de réalisation, la valeur CiPk de chaque composante de chaque pixel Pk de l'image Ik est remplacée par une somme pondérée calculée de la manière suivante : 13 où chaque W est une valeur de pondération prédéfinie. On note que la méthode de restauration utilisant une pondération impliquant deux images est un cas particulier de la méthode de restauration basée sur une fenêtre d'images où N1 = 2. Dans un mode de réalisation, seuls les pixels de l'image courante 1k appartenant 5 à l'objet 4 sont restaurés. Dans un mode de réalisation, plusieurs méthodes de restauration sont utilisées en fonction de l'objet auquel appartient chaque pixel de l'image 1k. Une première méthode de restauration est appliquée aux pixels appartenant au fond de l'image courante 1k et une seconde méthode de restauration est appliquée aux pixels appartenant à l'objet 4 10 dans l'image courante 1k. La première méthode de restauration est par exemple la méthode de restauration utilisant une fenêtre d'images dans laquelle la fenêtre d'images comprend deux images. La seconde méthode de restauration est par exemple la méthode de restauration utilisant une fenêtre d'images dans laquelle la fenêtre d'images comprend cinq images. 15 La Fig. 7 illustre schématiquement un procédé d'estimation de mouvement global. Le procédé de la Fig. 7 est appliqué par le module de traitement 52 lors de l'étape 41 pour déterminer le mouvement global du fond (i.e. l'homographie derd) et lors de l'étape 44 pour déterminer le mouvement global de l'objet 4 (i.e. l'homographie 20 dHObjet). Lors de l'étape 41, seul le mouvement des pixels correspondant au fond est recherché. Le mouvement des pixels correspondant à l'objet 4 ne doit pas être pris en compte. Dans une étape 410, le module de traitement obtient une position et une forme de 25 l'objet 4 dans l'image 1k. Dans un mode de réalisation, la forme et la position de l'objet 4 dans l'image 1k sont données par une liste ordonnée de pixels, appelés points de contrôle, appartenant à un contour C de l'objet 4. La liste ordonnée de points de contrôle peut comprendre tous les pixels appartenant au contour C de l'objet 4 ou un sous-ensemble de pixels du 30 contour C permettant d'obtenir une bonne approximation du contour C. Un parcours de la liste ordonnée de points de contrôle permet d'obtenir le contour C. Dans un mode de réalisation, lors de l'étape 410, le module de traitement 52 fait une hypothèse de mouvement faible de l'objet 4 entre deux images successives de la séquence d'images. Comme dit plus haut, le procédé décrit en relation avec la Fig. 4 est 14 itératif de sorte que lorsque l'image courante Ik est en cours de traitement par ce procédé, l'image précédente Ik_i a déjà été traitée par ce procédé. La position et la forme de l'objet 4 dans l'image Ik_i sont donc connues. En s'appuyant sur l'hypothèse de mouvement faible, le module de traitement 52 considère que la forme et la position 5 de l'objet 4 dans l'image Ik sont identiques à la forme et la position de l'objet 4 dans l'image Ik_1. Le module de traitement 52 réutilise donc la liste ordonnée de points de contrôle de l'image Ik_i pour définir le contour C dans l'image Ik. La première image de la séquence Ik=0 est un cas particulier puisque cette image n'est précédée d'aucune autre image. Dans un mode de réalisation, la position et la 10 forme de l'objet 4 dans la première image Io de la séquence d'images sont données par un opérateur. Pour ce faire, l'opérateur peut détourer l'objet 4 à l'aide d'un dispositif de pointage, tel qu'une souris, sur le dispositif d'affichage 53, qui dans ce cas est un écran tactile. Le fond de la première image est considéré comme ayant un mouvement nul. Par conséquent le mouvement du fond n'a pas à être compensé pour la première 15 image de la séquence. Le procédé décrit en relation avec la Fig. 4 n'est pas appliqué à la première image Io de la séquence, mais cette image Io est utilisée dans ce procédé pour déterminer le mouvement du fond et le mouvement de l'objet 4 entre cette image /0 et l'image Ik_i qui la suit dans la séquence d'images. Dans une étape 411, la position et la forme de l'objet dans l'image courante Ik 20 étant connues, le module de traitement masque chaque pixel de l'image courante Ik appartenant à l'objet 4, i.e. le module de traitement masque chaque pixel appartenant au contour ou interne au contour de l'objet 4. Dans un mode de réalisation, le module de traitement 52 associe chaque pixel de l'image courante Ik à une première valeur de masque lorsque ledit pixel est masqué et à une seconde valeur de masque lorsque ledit 25 pixel n'est pas masqué. La première valeur de masque est par exemple la valeur « 1 » et la seconde valeur de masque est par exemple la valeur « 0 ». Dans une étape 412, le module de traitement 52 estime le mouvement global du fond entre l'image courante Ik et l'image précédente Ik_i (i.e. le module de traitement 52 estime l'homographie derd). Lors de cette estimation seuls les pixels de l'image 30 Ik non masqués (i.e. les pixels associés à la deuxième valeur de masque) sont pris en compte. De plus, seuls les pixels de l'image précédente Ik-1 non masqués sont pris en compte, en utilisant le masque obtenu lors de l'application de la méthode décrite en relation avec la Fig. 4 à l'image précédente Ik_1. Dans un mode de réalisation, la détermination des huit paramètres de l'homographie derd utilise la méthode 15 d'ajustement projectif (« projective fit » en terminologie anglo-saxonne) ou la méthode des flux projectifs (« projective flow » en terminologie anglo-saxonne) décrite dans l'article « Orbite vidéo du groupe projectif : une approche simple de l'estimation sans caractéristiques des paramètres (« vidéo Orbits of the Projective Group : a simple 5 approach to featureless estimation of paramètres » en terminologie anglo-saxonne) », Steve Mann and Rosaling W. Picard, IEEE Tr. On Image Processing, Vol. 6, N°. 9, September 1997. Lors de l'étape 44, seul le mouvement des pixels correspondant à l'objet 4 est recherché. Le mouvement des pixels correspondant au fond ne doit pas être pris en 10 compte. Lors de l'étape 44, le module de traitement met en oeuvre l'étape 410 d'obtention de la position et de la forme de l'objet 4. La position et la forme de l'objet 4 dans l'image courante Ik sont obtenues grâce au résultat de l'étape 43 Lors de l'étape 411, le module de traitement 52 masque chaque pixel de l'image 15 courante Ik.appartenant au fond, i.e. n'appartenant pas à l'objet 4. Lors de l'étape 412, le module de traitement 52 estime le mouvement de l'objet entre l'image courante Ik et l'image précédente Ik_i (i.e. le module de traitement estime 1' homographie dHknesurée puis en déduit l'homographie dHkObjet, ) Lors de cette estimation seuls les pixels de l'image courante Ik et de l'image précédente Ik_i non 20 masqués sont pris en compte. De nouveau, la détermination des huit paramètres de l'homographie dirknesurée utilise la méthode d'ajustement projectif ou la méthode des flux projectifs. Dans un mode de réalisation, lors de la mise en oeuvre de l'étape 410 lors de l'étape 41, le module de traitement 52 fait une hypothèse de mouvement continu de 25 l'objet 4 dans la séquence d'images. L'hypothèse de mouvement continu implique que le mouvement de l'objet 4 entre l'image courante Ik et l'image précédente Ik_i est le même que le mouvement de l'objet 4 entre l'image précédente Ik_i et une image Ik_2 précédent l'image précédente Ik_1. Le procédé décrit en relation avec la Fig. 4 étant itératif, lors du traitement de l'image courante Ik, le mouvement de l'objet 4 entre 30 l'image précédente Ik_i et l'image précédente Ik_2 est connu. Par ailleurs la position et la forme de l'objet 4 dans l'image précédente Ik_i sont elles aussi connues. La position et la forme de l'objet 4 dans l'image courante Ik peuvent donc être trouvées en utilisant eisurée = * dHkObjiet une homographie dHfc12 dl-eV représentant le mouvement de l'objet 4 entre l'image précédente Ik_i et l'image précédente Ik_2. L'homographie 16 Cl frkneisurée est une combinaison d'une homographie dlircmi'd représentative du mouvement du fond entre l'image précédente Ik_i et l'image précédente Ik_2 et d'une homographie dilk°k1 «let représentative du mouvement de l'objet 4 entre l'image précédente Ik_i et l'image précédente Ik_2. L'homographie dllikneisurée est appliquée à 5 l'objet 4 dans l'image précédente Ik_1. Plus précisément l'homographie die appliquée aux points de contrôle de la liste ordonnée de points de contrôle représentant le contour de l'objet 4 dans l'image précédente Ik_i afin d'obtenir la liste ordonnée de points de contrôle représentant le contour de l'objet 4 dans l'image Ik. L'hypothèse de mouvement continu est applicable lors de la mise en oeuvre de l'étape 410 lors de l'étape 10 44. Dans un mode de réalisation, pour prendre en compte que les hypothèses de mouvement faible et de mouvement continu ne permettent d'obtenir qu'une approximation de la forme et la position de l'objet 4 dans l'image courante Ik, une dilatation est appliquée au contour de l'objet 4 dans l'image courante Ik . La dilatation 15 est obtenue par exemple en utilisant une méthode de morphologie mathématique. Dans d'autres modes de réalisation, d'autres méthodes d'estimation des paramètres d'une homographie connues peuvent être utilisées. Dans d'autres modes de réalisation, le module de traitement considère que le mouvement entre deux images successives de la séquence d'images est représenté par 20 d'autres modèles de mouvement telles qu'une translation, une rotation, une transformation affine ou une transformation bilinéaire. Dans un mode de réalisation, préalablement à chaque estimation de mouvement global (du fond ou de l'objet 4), chaque image impliquée dans l'estimation de mouvement global est interpolée au demi, au quart ou au huitième de pixel. De cette 25 manière, la précision de l'estimation de mouvement global est améliorée. La méthode d'ajustement projectif (respectivement la méthode des flux projectifs) consiste à trouver parmi un ensemble de paramètres d'un modèle de mouvement (ici les 8 paramètres d'une homographie), les paramètres du modèle de mouvement minimisant une métrique représentative d'une erreur entre un mouvement réel d'un objet dans une 30 image et un mouvement de l'objet représenté par le modèle de mouvement. Dans la méthode d'ajustement projectif (respectivement la méthode des flux projectifs), chaque combinaison possible de paramètres du modèle de mouvement est testée. Une telle méthode exhaustive de recherche des paramètres du modèle de mouvement peut avoir un coût calculatoire important. Il est possible de réduire le coût calculatoire de la eisurée est 17 méthode d'ajustement projectif (respectivement la méthode des flux projectifs) en utilisant, par exemple, un algorithme de descente de gradient plutôt qu'une recherche exhaustive. Toutefois, un problème connu des méthodes de descente de gradient est que, lorsque la métrique à minimiser possède plusieurs minimums locaux, la méthode de 5 descente de gradient peut converger vers un minimum local qui n'est pas un minimum global, i.e. qui n'est pas la valeur minimale que peut prendre la métrique. Une méthode permettant d'assurer une convergence rapide vers le minimum global de la métrique consiste à initialiser la méthode de descente de gradient avec une valeur proche du minimum global recherché. Dans un mode de réalisation, la recherche exhaustive des 10 paramètres du modèle de mouvement de la méthode d'ajustement projectif (respectivement la méthode des flux projectifs) est remplacée par une méthode de descente de gradient lors de la mise en oeuvre de l'étape 412 lors des étapes 41 et 44. Lors de la mise en oeuvre de l'étape 41 sur l'image Ik, la méthode de descente de gradient est initialisée à une valeur dlircmi'd représentative du mouvement du fond 15 trouvé pour l'image Ik_1. Plus précisément, les huit paramètres de l'homographie dlircmi'd représentative du mouvement du fond entre l'image précédente Ik-1 et l'image précédente Ik_2 sont utilisés pour initialiser les huit paramètres de l'homographie die.ona représentative du mouvement du fond entre l'image courante Ik et l'image précédente Ik_i dans la méthode de descente de gradient. 20 «le De même, lors de la mise en oeuvre de l'étape 44 sur l'image Ik, la méthode de descente de gradient est initialisée à une valeur(dHkHObjF otnd * decond). Plus précisément, les huit paramètres de l'homographie (dHkObj iik et * ci )sont utilisés pour initialiser les huit paramètres de l'homographie dlliknesurée représentative du mouvement de l'objet 4 mesuré entre l'image courante Ik et l'image précédente Ik_i dans la méthode 25 de descente de gradient. L'homographie el°biet est ensuite déduite de l'homographie dtpknesurée et de l'homographie dend. Dans un mode de réalisation, suite à l'estimation de l'information représentative du mouvement de l'objet 4 entre l'image courante Ik et l'image précédente recalée Ik_i , (i.e. suite à l'estimation des huit paramètres de l'homographie dH°biet ) l'information 30 représentative du mouvement estimé est filtrée afin de garantir des variations régulières du mouvement de l'objet entre deux images successives de la séquence d'images. La méthode décrite en relation avec la Fig. 4 est particulièrement efficace lorsqu'une hypothèse de mouvement inertiel de l'objet 4 est vérifiée. Lorsque les informations de mouvement estimées de l'objet 4 sont trop variables, il est préférable de corriger ces 18 informations afin qu'elles se rapprochent d'informations de mouvement compatibles avec l'hypothèse de mouvement inertiel. La Fig. 8 décrit un procédé de filtrage de mouvement permettant de corriger des informations de mouvement. Le procédé décrit en relation avec la Fig. 4 n'utilise alors plus les informations de mouvement estimées 5 mais les informations de mouvement estimées corrigées. La Fig. 8 illustre schématiquement un procédé de filtrage de mouvement. Le procédé décrit en relation avec la Fig. 8 permet de garantir des variations régulières du mouvement de l'objet 4 entre deux images successives de la séquence d'images. 10 Dans une étape 800, le module de traitement 52 détermine une matrice de passage Vk (respectivement une matrice de passage Vk_i) permettant d'estimer un mouvement de l'objet dans un référentiel centré sur un barycentre de l'objet 4 dans l'image courante Ik (respectivement dans l'image précédente recalée Iircecl). 15 Vk =[1 Vk-1 = 0 0 Xb kb-11 0 1 yk 1 0 1 0 0 Xb 0 1 v k-1 0 1 où (xe, ye) (respectivement (xe_1, ye_1)) sont des coordonnées du barycentre de l'objet 4 dans l'image courante /k (respectivement dans l'image précédente recalée fic'-ec1). 20 Dans une étape 801, le module de traitement 52 calcule une information dH°k1et représentative du mouvement de l'objet 4 dans le référentiel centré sur le barycentre de l'objet 4. Tx - Ty "Objet =uObjet u k .1-vilk - v = - Les coefficients de la matrice dejet notés « . » sont des coefficients qui ne sont 25 pas utilisés par la suite du fait d'une approximation mise en oeuvre lors d'une étape suivante 802. Dans l'étape 802, le module de traitement 52 obtient une matrice Hf, représentative de composantes de translation du mouvement de l'objet 4 entre l'image courante Ik et l'image précédente recalée Iirceci (représenté par l'homographie dHkObjet) 30 dans le référentiel centré sur le barycentre de l'objet 4 de la manière suivante : 19 1 0 Tx Hk 1 = ApproxT(Vjl.dHobjet Vk-1) = 0 1 Ty 0 0 1 où Tx et Ty sont des paramètres d'une translation et ApproxT(X) est une approximation en translation de la matrice homographique X Dans une étape 803, le module de traitement 52 obtient une matrice flik2 5 représentative de composantes du mouvement de l'objet 4 entre l'image courante Ik et l'image précédente recalée /in autres que les composantes de translation de la manière suivante : He = H itc-i .vk_i dHkObjet .vk_i Les composantes du mouvement de l'objet 4 autres que les composantes de 10 translation peuvent être par exemple des composantes de rotation, de zoom, etc. Dans une étape 804, le module de traitement 52 filtre les composantes de translation du mouvement de l'objet 4 entre l'image courante Ik et l'image précédente recalée Iirceci comme suit : Hit,Fitt = a. Hitc_iFitt + (1 - a)Hit, 15 où Hitc_i_Filt est une matrice représentative des composantes de translation filtrées du mouvement de l'objet 4 entre l'image précédente Ik-1 et l'image précédente recalée Iircec2 et a est une constante prédéfinie comprise entre « 0 » et « 1 ». Dans un mode de réalisation, a = 0.8. Dans une étape 805, le module de traitement 52 filtre les composantes du 20 mouvement de l'objet 4 entre l'image courante Ik et l'image précédente recalée Iircec2 autres que les composantes de translation de la manière suivante : HFilt ik2 = p./ + (1 - le)flik2 où HeFilt est une matrice représentative des composantes filtrées du mouvement de l'objet 4 autres que les composantes de translation, / est une matrice identité de taille 25 3 x 3 et le est une constante prédéfinie comprise entre « 0 » et « 1 ». Dans un mode de réalisation le = 0.5. Dans une étape 806, le module de traitement 52 détermine une information Filt AuObjet représentative du mouvement global filtré de l'objet 4 de la manière suivante : 30 dH° bjetFilt utFilt uhFilt =11k.iik .iik .vk Dans le mode de réalisation dans lequel le mouvement estimé est filtré afin de garantir des variations faibles du mouvement de l'objet entre deux images successives de la séquence d'images, le mouvement global filtré de l'objet 4 représenté par 20 Filt l'homographie filtrée clHObjet k remplace le mouvement global de l'objet 4 représenté par l'homographie dilic°biet lors de la restauration de l'image courante Ik. La Fig. 5 illustre schématiquement un procédé de segmentation d'images selon l'invention. 5 Le procédé décrit en relation avec la Fig. 5 est mis en oeuvre lors de l'étape 43. Dans un mode de réalisation, le procédé de segmentation d'images mis en oeuvre lors de l'étape 43 est un procédé de segmentation basé contours actifs. Un principe des procédés de segmentation basés contours actifs est de définir un contour initial au voisinage d'un objet puis de modifier itérativement ce contour pour qu'il épouse au 10 mieux la forme de l'objet. A chaque itération, le contour de l'objet obtenu lors d'une itération précédente, dit contour précédent, est modifié de sorte à obtenir un contour de l'objet, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent. Dans les procédés de segmentation basés contours actifs, le coût d'un contour est fonction d'une énergie interne et d'une énergie externe du contour. 15 Nous donnons par la suite des exemples de méthodes de calcul de valeurs représentatives d'une énergie interne et d'une énergie externe d'un contour. Un contour final de l'objet est obtenu lorsqu'une condition prédéfinie d'arrêt du procédé de segmentation basé contours actifs est remplie. Une condition d'arrêt peut par exemple être un nombre maximum d'itérations ou une obtention d'un écart entre deux coûts de 20 contours obtenus dans deux itérations successives inférieures à un seuil prédéfini. On note que plus le contour initial est proche du contour réel de l'objet, plus le procédé de segmentation basé contours actifs converge rapidement vers un contour proche du contour réel de l'objet. Un choix judicieux d'une position et d'une forme du contour initial permet donc d'améliorer des performances du procédé de segmentation basé 25 contours actifs. Dans un mode de réalisation le procédé de segmentation basé contours actifs est basé contour. Dans une étape 431, le module de traitement 52 obtient un contour C initial de l'objet 4 dans l'image courante Ik. Dans un mode de réalisation, lors de l'étape 431, le module de traitement fait 30 l'hypothèse de mouvement faible de l'objet 4 entre l'image courante Ik et l'image précédente Ik_1. Dans ce cas, comme lors de la mise en oeuvre de l'étape 411 lors de l'étape 44, le module de traitement 52 réutilise la liste ordonnée de points de contrôle déterminée lors de la mise en oeuvre du procédé décrit en relation avec la Fig. 4 sur 21 l'image précédente Ik_i pour obtenir le contour C initial de l'objet dans l'image courante/k. Dans une étape 432, le module de traitement 52 calcule un coût du contour courant C en appliquant un procédé que nous décrivons par la suite en relation avec la Fig. 6. 5 Lors de la première itération du procédé de segmentation basé contours actifs décrit en relation avec la Fig. 5, le contour courant C est le contour C initial. Dans une étape 433, le module de traitement 52 vérifie si une condition d'arrêt du procédé de segmentation basé contours actifs est remplie. Dans un mode de réalisation, ledit procédé itératif s'arrête lorsqu'un nombre d'itérations du procédé de segmentation 10 basé contours actifs atteint un nombre d'itérations maximal. Lorsque la condition d'arrêt est remplie, le procédé de segmentation basé contours actifs prend fin lors de l'étape 434 et le module de traitement 52 met en oeuvre l'étape 44 déjà expliquée. Lorsque la condition d'arrêt n'est pas remplie, le module de traitement 52 met en 15 oeuvre une étape 435. Lors de l'étape 435, le module de traitement 52 met en oeuvre une procédure de raffinement du contour C de l'objet 4 obtenu lors de l'itération précédente du procédé de segmentation basé contour actif. Lors de l'étape 435, le module de traitement 52 modifie le contour C de l'objet 4 obtenu lors de l'itération précédente du procédé de segmentation basé contour actif, dit contour précédent, de sorte à obtenir un 20 contour C de l'objet, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent. La modification du contour C utilise, par exemple, une méthode décrite dans l'article « snakes : Active Contour Models » Michael Kass, Andrew Witkin, Demetri Terzopoulos, International Journal of Computer Vision, 321331 (1988), 1987 Kluwer Academic Publishers, Boston. 25 L'étape 435 est suivie de l'étape 432. Dans un mode de réalisation, lors de l'étape 431, le module de traitement fait l'hypothèse de mouvement continu de l'objet 4 entre l'image courante Ik et l'image précédente Ik_1. Dans ce cas, le module de traitement 52 déplace les points de contrôle de la liste ordonnée de points de contrôle pour obtenir le contour C initial de l'objet 4 30 dans l'image courante Ik. Ces points de contrôle, déterminée lors de la mise en oeuvre du procédé décrit en relation avec la Fig. 4 sur l'image précédente * ciH: te/k_11;osnodnt déplacés du mouvement de l'objet 4 représenté par l'homographie dHkObj Dans un mode de réalisation, les points de contrôle de la liste ordonnée de points de contrôle sont déplacés du mouvement filtré de l'objet 4 représenté par l'homographie 22 d uFond- 'ObjetFilt pour obtenir le contour C initial de l'objet 4 dans l'image courante Ik. La Fig. 6 illustre schématiquement un procédé de calcul d'un coût d'un contour d'un objet compris dans le procédé de segmentation selon l'invention. 5 Le procédé décrit en relation avec la Fig. 6 est mis en oeuvre lors de l'étape 432. Dans une étape 4321, le module de traitement 52 calcule une énergie interne Eint du contour C de la manière suivante : Eint = ç a. I Cici I + 2 I 2 i=1 i=1 OÙ a et b sont des constantes prédéfinies égales par exemple à la valeur « 0.01 », 10 N est un nombre de points de contrôle dans la liste de points de contrôle représentant la courbe C, PCi est le i-ème point de contrôle de la liste de points de contrôle représentant la courbe C dans l'image courante Ik, ciçci est une dérivée première locale de la courbe C dans l'image courante Ik calculée au niveau du point de contrôle PCi, et Cp"c,i est une dérivée seconde locale de la courbe C dans l'image courante Ik calculée au niveau du 15 point de contrôle PCi. Dans une étape 4322, le module de traitement 52 calcule une énergie externe Eext du contour C de la manière suivante : Eext = -(Wedge - Eedge Wmvt-Emvt) où Wedge et Wmvt sont des constantes prédéfinies par exemple égales à la valeur 20 « 1 ». Eedge est une énergie, dite énergie de contour, calculée sur une image de module de gradient I kgrad obtenue à partir de l'image courante Ik : Eedge = ic lgrad (PCi) i=1 où 1/grad (PCi) est une valeur de module de gradient de l'image /kg correspondant à la position du point de contrôle PC,. 25 On note que différentes méthodes de calcul d'une image de module de gradient sont applicables ici. Pour obtenir l'image de de module gradient I kgrad on peut par exemple appliquer à chaque pixel de l'image Ik : - une combinaison linéaire de pixels adjacents dudit pixel, chaque pixel adjacent étant pondéré par un poids, la somme desdits poids étant égale à 30 zéro, puis calculer l'amplitude (i.e. le module) de cette combinaison linéaire ; - un filtre de Sobel ; rad - un filtre de Canny ; - Dans un mode de réalisation, l'image Igraa n'est pas calculée, et les valeurs de module de gradient utilisés dans le calcul de l'énergie de contour Eedge sont calculées 5 uniquement aux positions des N points de contrôle PCi. Em't est une énergie dépendante du mouvement de l'objet 4 entre l'image courante /k et l'image recalée précédente N Em't = I'k(PCi)-'ke1(PCi)It=1 où Ik (PCi)est une valeur d'un pixel de l'image courante Ik correspondant au point 10 de contrôle PC, et Iirt_eci(PC,) est une valeur d'un pixel de l'image ITt_eci situé à la même position que le pixel de l'image Ik correspondant au point de contrôle PCi. Dans une étape 4323, le module de traitement 52 calcule le coût J du contour courant C de la manière suivante : J = Eext Eint 15 On constate donc que le mouvement de l'objet 4 est pris en compte dans le procédé de segmentation selon l'invention ce qui permet d'obtenir une meilleure segmentation de l'objet. La minimisation du coût Jpermet de maximiser Em't et Eedge sur les points de contrôle du contour, afin de favoriser les zones à forts gradients spatiaux et/ou temporels. 20 Le principe de l'invention reste le même dans le cas d'une utilisation d'un autre type de procédés de segmentation basés contours actifs que les procédés de segmentation basés contours actifs basé contour. Chaque procédé de segmentation basé contours actifs comprend une estimation d'une énergie externe Eext. Toutefois, les procédés de segmentation basés contours actifs étant adaptés à des images fixes, ils ne 25 prennent pas en compte les mouvements dans une séquence d'images lors de la segmentation. L'invention permet de prendre en compte ces mouvements en intégrant une énergie représentative du mouvement dans l'estimation de l'énergie externe Eext. Ce principe s'applique aux énergies externes Eext calculées dans le cadre des procédés de segmentation basés contours actifs basés régions et des procédés de segmentation 30 basés contours actifs basés ensembles de niveaux. Jusque-là nous avons considéré des images ne comprenant qu'un objet. L'invention s'applique lorsque les images de la séquence d'images comprennent une pluralité d'objets. Lors des étapes 41 et 42, chaque objet est masqué lors de l'estimation 23 24 et la compensation du mouvement du fond dans une image. Les étapes 43, 44 et 45 sont mises en oeuvre indépendamment sur chaque objet. Par ailleurs, jusque-là, nous avons considéré que l'objet 4 était rigide et que par conséquent la forme apparente de l'objet était approximativement constante. Dans un 5 cas réel, suivant les mouvements de l'objet et/ou de la caméra, l'objet peut être vu selon des angles de vue différents ce qui peut provoquer des déformations dans la forme apparente de l'objet. Dans un mode de réalisation, lorsqu'une variation de la forme de l'objet sur une pluralité d'images successives de la séquence d'images dépasse un seuil prédéfini, le module de traitement 52 considère que l'objet apparaissant dans les images 10 a changé. Dans ce cas, lorsque le module de traitement 52 détecte un changement d'objet, il considère qu'une nouvelle séquence d'images démarre et invite l'opérateur à détourer à nouveau l'objet. Dans un autre mode de réalisation, le module de traitement 52 applique le procédé de segmentation décrit en relation avec la Fig. 5 quelles que soient les variations de la forme de l'objet, sans demander d'intervention d'un 15 opérateur. Dans un mode de réalisation, lorsque les images fournies par le dispositif d'acquisition d'images sont des images multi-composantes, le module de traitement 52 applique le procédé de restauration décrit en relation avec la Fig. 4 à chaque composante indépendamment. Chaque composante peut ensuite être affichée indépendamment ou 20 en combinaison avec une ou plusieurs autres composantes sur le dispositif de visualisation 53. Dans un mode de réalisation, lorsque les images fournies par le dispositif d'acquisition d'images sont des images multi-composantes, le module de traitement 52 applique le procédé de restauration décrit en relation avec la Fig. 4 à au moins une des 25 composantes, ou à au moins une composante calculée à partir des composantes disponibles dans les images. Seule l'étape de restauration 45 est appliquée à chaque composante indépendamment en utilisant les informations représentatives du mouvement du fond et du mouvement de l'objet 4 obtenues lors des étapes précédentes. Chaque composante peut ensuite être affichée indépendamment ou en combinaison 30 avec une ou plusieurs autres composantes sur le dispositif de visualisation 53. Par exemple, lorsque les images multi-composantes comprennent une composante de luminance et deux composantes de chrominance, le procédé de restauration décrit en relation avec la Fig. 4 est appliqué uniquement à la composante de luminance, l'étape REVENDICATIONS1) Procédé de restauration d'images d'une séquence d'images, caractérisé en ce que le procédé comprend lorsqu'il est appliqué à une première image de la séquence d'images : estimer (41) une information représentative d'un mouvement global d'un arrière-plan de la première image par rapport à une deuxième image ; compenser (42) le mouvement global de l'arrière-plan dans la deuxième image en utilisant ladite information représentative du mouvement global de l'arrière-plan afin d'obtenir une version recalée de la deuxième image, dite deuxième image recalée ; obtenir (43) un contour d'un objet de la première image en appliquant un procédé de segmentation, ledit procédé de segmentation étant itératif et comprenant, lors d'une itération, une modification (435) d'un contour de l'objet dans la première image obtenue lors d'une itération précédente dudit procédé de segmentation, dit contour précédent, de sorte à obtenir un contour de l'objet dans la première image, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent, un contour final de l'objet étant obtenu lorsqu'une condition prédéfinie d'arrêt dudit procédé de segmentation est remplie (433), le coût d'un contour de l'objet dans la première image étant une somme (4323) entre une première valeur représentative d'une énergie interne audit contour et une seconde valeur représentative d'une énergie externe audit contour, l'énergie externe audit contour étant fonction d'au moins une énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image recalée et d'une énergie, dite énergie de contour, correspondant à une somme de valeurs de modules de gradient calculées pour des pixels d'un deuxième ensemble de pixels appartenant au contour courant de l'objet ; une valeur représentative de l'énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image étant calculée sous forme d'une somme de différences entre des valeurs représentatives de pixels d'un premier ensemble de pixels de la première image appartenant au contour courant et des valeurs représentatives de pixels situés aux mêmes positions spatiales que les pixels du premier ensemble de pixels dans la seconde image ; utiliser le contour de l'objet ainsi obtenu pour estimer (44) une information représentative d'un mouvement global de l'objet ; et, 26 appliquer (45) à la première image une méthode de restauration d'images utilisant les informations représentatives du mouvement global de l'arrière-plan et du mouvement global de l'objet estimées. 5 2) Procédé selon la revendication 1, caractérisé en ce que pour calculer la valeur représentative de l'énergie interne au contour courant, une dérivée première locale et une dérivée seconde locale du contour sont calculées pour des pixels d'un troisième ensemble de pixels de l'image courante appartenant au contour courant de l'objet, ladite valeur représentative de l'énergie interne étant fonction desdites dérivées calculées. 10 3) Procédé selon la revendication 2, caractérisé en ce que les premier, deuxième et troisième ensembles de pixels sont identiques, et chaque ensemble comprend au moins une sous-partie des pixels de l'image courante appartenant au contour courant de 1 ' obj et. 15 4) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lors d'une première itération dudit procédé, un contour initial de l'objet dans l'image courante est obtenu (431) à partir d'un contour final obtenu lors d'une application du procédé de segmentation à l'image de référence ou d'un contour spécifié 20 par un opérateur dans l'image de référence. 5) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lors de chaque estimation d'une information représentative d'un mouvement global, une information représentative de la forme et de la position de l'objet est 25 obtenue, ladite information représentative de la forme et de la position de l'objet étant utilisée pour masquer des pixels ne devant pas être pris en compte dans ladite estimation. 6) Procédé selon la revendication 5, caractérisé en ce que suite à l'estimation de 30 ladite information représentative du mouvement global de l'objet, dite première information, un filtrage est appliqué à ladite première information afin de garantir des variations régulières du mouvement de l'objet entre deux images successives de la séquence d'images, ledit filtrage comprenant les étapes suivantes : 27 déterminer (800) une première matrice permettant d'estimer un mouvement de l'objet dans un référentiel centré sur un barycentre de l'objet dans la première image recalée et une deuxième matrice permettant d'estimer un mouvement de l'objet dans un référentiel centré sur un barycentre de l'objet dans la deuxième image recalée ; 5 utiliser la première et la deuxième matrices pour calculer (801) une information représentative du mouvement de l'objet dans ledit référentiel, dite seconde information, à partir de ladite première information ; utiliser la seconde information pour obtenir (802) une troisième matrice représentative de composantes de translation du mouvement de l'objet ; 10 utiliser la seconde information et la troisième matrice pour obtenir (803) une quatrième matrice représentative de composantes du mouvement de l'objet autres que les composantes de translation ; obtenir (804) une version filtrée de la troisième matrice, dite troisième matrice filtrée, en calculant une somme pondérée entre la troisième matrice et une troisième 15 matrice filtrée obtenue lors de la mise en oeuvre dudit procédé sur la deuxième image ; obtenir (805) une version filtrée de la quatrième matrice, dite quatrième matrice filtrée, en calculant une somme pondérée entre la quatrième matrice et une quatrième matrice filtrée obtenue lors de la mise en oeuvre du procédé sur la deuxième image ; et, obtenir une information représentative d'un mouvement global filtré de l'objet en 20 utilisant la première et la deuxième matrices, la troisième matrice filtrée et la quatrième matrice filtrée. 7) Procédé selon la revendication 6, caractérisé en ce que la seconde information est calculée de la manière suivante : 25 dHicobjet = dHkobjet. où Vk est la première matrice, Vk_i est la deuxième matrice, dHkObjetla première information et dHkObjet la seconde information. 8) Procédé selon la revendication 7, caractérisé en ce que la troisième matrice est 30 calculée de la manière suivante : Hk = ApproxT(Vk-1. dHobjet vk_i) où Hf, est la troisième matrice et ApproxT(X) est une approximation en translation d'une matrice homographique X 28 9) Procédé selon la revendication 8, caractérisé en ce que la quatrième matrice est calculée de la manière suivante : = Ht-1.11-1. u0b jet u "k k v k k - v k-1 où Hi,' est la quatrième matrice. 5 10) Procédé selon la revendication 9, caractérisé en ce que la troisième matrice filtrée est calculée de la manière suivante : Hit,Fut = a. Hitc_1Fut + (1 - a)Hit, où 1-1f,Filt est la troisième matrice filtrée, Hitc_i_Filt est la troisième matrice filtrée 10 obtenue lors de la mise en oeuvre dudit procédé sur la deuxième image et a est une constante prédéfinie comprise entre « 0 » et « 1 ». 11) Procédé selon la revendication 10, caractérisé en ce que la quatrième matrice filtrée est calculée de la manière suivante : 15 Hik,Fut = I + (1 _ e)Hik, où Hik2Filt est la quatrième matrice filtrée, I est une matrice identité et le est une constante prédéfinie comprise entre « 0 » et « 1 ». 12) Procédé selon les revendications 10 et 11, caractérisé en ce que l'information 20 représentative d'un mouvement global filtré de l'objet est calculée de la manière suivante : dH° bjetFilt iltFilt uhFilt = Vk-iik - 11k - vk_i. 13) Dispositif de restauration d'images d'une séquence d'images, caractérisé en 25 ce que le dispositif comprend : des moyens d'estimation (41) pour estimer une information représentative d'un mouvement global d'un arrière-plan d'une première image de la séquence d'images par rapport à une deuxième image de la séquence d'images ; des moyens de compensation de mouvement (42) pour compenser le mouvement 30 global de l'arrière-plan dans la deuxième image en utilisant ladite information représentative du mouvement global de l'arrière-plan afin d'obtenir une version recalée de la deuxième image, dite deuxième image recalée ; des moyens d'obtention (43) permettant d'obtenir un contour d'un objet dans la première image recalée, comprenant des moyens de modification (435) d'un contour de 29 l'objet dans la première image obtenu lors d'une itération précédente, dit contour précédent, permettant d'obtenir un contour de l'objet dans la première image, dit contour courant, tel qu'un coût du contour courant soit inférieur à un coût du contour précédent, un contour final de l'objet étant obtenu lorsqu'une condition prédéfinie 5 d'arrêt est remplie (433), le coût d'un contour de l'objet dans la première image étant une somme (4323) entre une première valeur représentative d'une énergie interne audit contour et une seconde valeur représentative d'une énergie externe audit contour, l'énergie externe audit contour étant fonction d'au moins une énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image recalée et 10 d'une énergie, dite énergie de contour, correspondant à une somme de valeurs de modules de gradient calculées pour des pixels d'un deuxième ensemble de pixels appartenant au contour courant de l'objet ; une valeur représentative de l'énergie dépendante d'un mouvement global de l'objet entre la première image et la seconde image étant calculée sous forme d'une somme de différences entre des valeurs 15 représentatives de pixels d'un premier ensemble de pixels de la première image appartenant au contour courant et des valeurs représentatives de pixels situés aux mêmes positions spatiales que les pixels du premier ensemble de pixels dans la seconde image ; des moyens d'estimation de mouvement (44) utilisant le contour de l'objet dans 20 la première image pour estimer une information représentative d'un mouvement global de l'objet entre la première image et la deuxième images ; et, des moyens de restauration d'images (45) utilisant les informations représentatives du mouvement global de l'arrière-plan et du mouvement global de l'objet estimées. 25 14) Programme d'ordinateur, caractérisé en ce qu'il comprend des instructions pour mettre en oeuvre, par un dispositif (52), le procédé selon l'une quelconque des revendications 1 à 12, lorsque ledit programme est exécuté par un processeur (521) dudit dispositif (52). 30 15) Moyens de stockage, caractérisés en ce qu'ils stockent un programme d'ordinateur comprenant des instructions pour mettre en oeuvre, par un dispositif (52), le procédé selon l'une quelconque des revendications 1 à 12, lorsque ledit programme est exécuté par un processeur (521) dudit dispositif (52).