Brevet FR2884341A1 — publié le 20061013

TITRE Procédé et système d'amélioration d'une image digitale générée à partir d'un détecteur de rayons X. DOMAINE DE L'INVENTION La présente invention concerne le domaine de l'imagerie radiologique à rayons X qui permet de visualiser, par exemple, des organes d'un patient. L'invention concerne plus particulièrement le traitement d'images fluoroscopiques permettant d'atténuer le bruit des images fluoroscopiques pour les rendre plus nettes; néanmoins, l'invention pourra également être appliqué au traitement d'images radiographiques par exemple. ARRIERE PLAN DE L'INVENTION Dans le domaine de l'imagerie médicale, il est bien connu d'utiliser des images fluoroscopiques afin d'aider le guidage d'instruments chirurgicaux au cours d'une opération. Ces images fluoroscopiques sont acquises par un dispositif d'imagerie classiquement constitué d'une source de rayons X positionnée en vis à vis d'un récepteur d'image, ladite source de rayons X et le récepteur d'image étant aptes à être entraînés en rotation autour d'au moins un axe, usuellement trois axes, des moyens de contrôle, des moyens d'acquisition, des moyens de visualisation des images et des moyens de commande. Le patient est positionné sur une table susceptible de se déplacer dans les trois translations possibles associées à un espace donné, c'est-à-dire longitudinalement, latéralement et verticalement, de telle sorte que la partie du corps du patient à examiner et/ou à traiter s'étende entre la source de rayons X et le récepteur d'image. Cette mobilité de la table et de la source de rayons X et du récepteur d'image permet à un praticien de faire l'acquisition d'images pour n'importe quelle partie du corps d'un patient couché sur la table. Ainsi, il est usuel d'utiliser des images fluoroscopiques en deux dimensions obtenues par l'irradiation du patient par de faibles 2884341 2 doses de rayons X, un agent de contraste étant de préférence préalablement injecté, lors d'une intervention pour guider l'instrument dans l'organe du patient à traiter. Les informations associées à ces images fluoroscopiques peuvent être avantageusement introduites dans des images reconstruites en trois dimensions pour améliorer le guidage des instruments chirurgicaux. De manière alternative, des images en trois dimensions acquises peuvent être projetées sur les images fluoroscopiques en deux dimensions acquises au cours de l'intervention. Contrairement aux images radiographiques acquises par le dispositif d'imagerie en émettant des rayons X de fortes doses permettant d'obtenir des images de bonne qualité avec un niveau de bruit faible, c'est-à-dire présentant un rapport signal sur bruit élevé, les images fluoroscopiques obtenues avec des doses de rayons X plus faibles comportent un niveau de bruit plus élevé, c'est-à-dire présentant un rapport signal sur bruit faible, et sont donc de moins bonne qualité ce qui est susceptible de gêner le bon déroulement de l'intervention chirurgicale. En effet, le bruit enregistré par le détecteur de rayons X et apparaissant sur les images est d'origine quantique et dépend de la racine carré du nombre de photons détectés par pixel. Lorsque la dose de rayons X est diminuée, le bruit décroît moins rapidement que ladite dose de sorte que le ratio bruit/signal augmente. De plus au bruit quantique vient s'ajouter le mouvement de fond lié notamment à la respiration du patient, aux déplacements des instruments chirurgicaux et aux mouvements de la table sur laquelle est placée le patient. Afin de supprimer le bruit des images acquises, on pourrait appliquer un filtre temporel dans l'hypothèse d'images immobiles; Toutefois, les images acquises étant mobiles en fluoroscopie, l'application d'un simple filtre temporel se traduit par un mouvement flou et/ou une perte de contraste des objets mobiles, outre un pic de bruit. Les procédés de traitement d'images en fluoroscopie consistent à appliquer un filtre dit à mouvement compensé. La plupart des filtres de l'art antérieur font appel à un critère de distinction entre une variation due au bruit et une variation due au mouvement; néanmoins, l'arrêt du filtrage 2884341 3 provoque la réapparition du bruit, ce qui se traduit sur les images par des traînées de bruit derrière les objets mobiles. Afin de remédier à ces inconvénients, on a déjà imaginé des procédés de traitement d'une séquence d'images fluoroscopiques améliorant la qualité des images visualisées. C'est le cas, par exemple, de la demande de brevet français FR 2 790 123 déposée par la société GE Medical Systems qui décrit un procédé dans lequel, pour chaque image courante acquise, on détermine le déplacement de ladite image courante par rapport à l'image précédente acquise dans le plan d'acquisition des images, on élabore une image filtrée précédente décalée, en décalant spatialement l'image filtrée précédente, et on élabore l'image filtrée courante par la moyenne pondérée entre l'image courante acquise et l'image filtrée précédente décalée. Ce type de procédé présente l'inconvénient de ne pas procurer des résultats satisfaisants. En effet, les images fluoroscopiques se présentent comme une superposition de couches d'images transparentes de sorte qu'il est impossible d'identifier clairement les pixels et les objets physiques. Ainsi, même si le mouvement global de la séquence d'image peut être déterminé, on ne peut pas séparer les différentes couches des objets pour les filtrer temporellement indépendamment les unes des autres, ce qui induit une perte de contraste de l'image et une traînée de bruit derrière les objets mobiles. BREVE DESCRIPTION DE L'INVENTION L'un des buts de l'invention est donc de remédier à ces inconvénients en proposant un procédé d'amélioration d'une image digitale d'un objet générée à partir d'un détecteur de rayons X recevant des rayons X émis par une source permettant d'atténuer efficacement le bruit de l'image pour la rendre plus nette. Conformément à l'invention, il est proposé un procédé d'amélioration d'une image digitale d'un objet définie par des pixels, acquise à l'instant t et 2884341 4 générée à partir d'un détecteur de rayons X recevant des rayons X émis par une source remarquable en ce qu'il consiste au moins à déterminer une image dite prédite de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-i, i étant un nombre entier positif supérieur ou égal à 1, et des déplacements de chacune des couches constituant l'image acquise à l'instant t, le nombre de couches étant préalablement fixé et les déplacements de chacune desdites couches étant préalablement déterminés, et à générer une image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. Selon une caractéristique générale de l'invention, pour chaque image acquise à l'instant t, on fixe le nombre N de couches dont est constituée l'image acquise à l'instant t, puis on détermine les déplacements vn de chacune des N couches à partir des images acquises aux instants t-i, puis on génère N images primaires correspondant aux images acquises aux instants t-i, dont les pixels sont déplacés respectivement de la somme des N déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-N, de la somme de N-1 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-1-N, de la somme de N-2 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-2-N et ainsi de suite jusqu'à chaque déplacement vn, pour l'image acquise à l'instant t-1,, puis on génère une image dite prédite de l'objet à l'instant t correspondant à la somme en signes alternés des N images primaires, puis on génère l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. De préférence, on fixe le nombre N de couches, dont est constituée l'image acquise à l'instant t, à N=2, le procédé d'amélioration consiste alors à déterminer les déplacements v1 et v2 de chacune des deux couches à partir des images acquises aux instants t-1 et t-2, à générer deux images primaires, une première image primaire correspondant à l'image acquise à l'instant t-1 dont les pixels sont déplacés du déplacement v1 et une seconde image primaire correspondant à l'image acquise à l'instant t-1 dont les pixels sont déplacés du déplacement v2, à générer une image intermédiaire 2884341 5 correspondant à la somme des deux images primaires, à générer une troisième image primaire correspondant à l'image acquise à l'instant t-2 dont les pixels sont déplacés suivant la somme des déplacements v1 et v2 de chacune des couches, à générer une image dite prédite de l'objet à l'instant t correspondant à l'image intermédiaire à laquelle est soustraite la troisième image secondaire, et finalement à générer l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. Quel que soit le nombre de couches, le procédé suivant l'invention comporte une étape d'application d'un filtre temporel d'une part à l'image prédite et d'autre part à l'image de l'objet acquise à l'instant t préalablement à la génération de l'image visualisée. De préférence, la génération de l'image visualisée est obtenue par l'addition du tiers de l'intensité de chaque pixel de l'image prédite avec les deux tiers de l'intensité de chaque pixel de l'image de l'objet acquise à l'instant t. Afin d'améliorer encore la netteté de l'image visualisée, chaque image acquise à l'instant t est divisée en au moins deux cadres et les étapes de détermination du nombre N de couches, de détermination des déplacements v de chacune des N couches, de génération de N images primaires, de génération de l'image prédite et de génération de l'image visualisée sont appliquées indépendamment pour chaque cadre. Bien évidemment, l'image visualisée est obtenue par l'assemblage de chacun des cadres après traitement. Un autre objet de l'invention concerne un système pour l'amélioration d'une image digitale d'un objet définie par dés pixels, acquise à l'instant t et générée par un détecteur de rayons X recevant des rayons X émis par une source remarquable en ce qu'il comporte des moyens de détermination d'une image dite prédite de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-i, i étant un nombre entier positif supérieur ou égal à 1, et des déplacements de chacune des couches, le nombre de couches dont est constituée l'image acquise à l'instant t étant préalablement fixé et les déplacements de chacune desdites couches étant préalablement 2884341 6 déterminés, et des moyens de génération de l'image visualisée correspondant à la somme de l'image prédite et de l'image de l'objet acquise à l'instant t. Selon une caractéristique essentielle de l'invention, les moyens de détermination de l'image prédite consistent en des moyens de détermination des déplacements vn de chacune des N couches à partir des images acquises aux instants t-i, des moyens de génération de N images primaires correspondant aux images acquises aux instants t-i, dont les pixels sont déplacés respectivement de la somme des N déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-N, de la somme de N-1 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-1-N, de la somme de N-2 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-2-N et ainsi de suite jusqu'à chaque déplacement vn, pour l'image acquise à l'instant t-1, et des moyens de génération de l'image dite prédite de l'objet à l'instant t. BREVE DESCRIPTION DES DESSINS D'autres avantages et caractéristiques ressortiront mieux de la description qui va suivre, de plusieurs variantes d'exécution, données à titre d'exemples non limitatifs, du procédé et du système d'amélioration d'une image digitale d'un objet acquise à partir d'un détecteur de rayons X conforme à l'invention, à partir des dessins annexés sur lesquels: - la figure 1 est une représentation schématique du système permettant de mettre en oeuvre le procédé conforme à l'invention, - la figure 2 est une représentation schématique de la décomposition en deux couches d'une image fluoroscopique acquise à l'instant t et son traitement selon le procédé conforme à l'invention, - la figure 3 est un diagramme schématique du procédé de traitement 30 d'une image décomposée en deux couches suivant l'invention. 2884341 7 DESCRIPTION DETAILLEE DE L'INVENTION En référence à la figure 1, le dispositif d'imagerie pour l'acquisition d'images fluoroscopiques, notamment, est classiquement constitué d'un récepteur d'image numérique 1, une source de rayons X 2 émettant des rayons X sur le récepteur d'image 1, ledit récepteur d'image et la source de rayons X 2 étant respectivement positionnés aux extrémités d'un bras en forme de C ou de U non représenté sur la figure qui pivote autour de trois axes. De plus, l'objet 3 à radiographier, usuellement un patient, est placé entre la source de rayons X 2 et le récepteur d'image 1. Par ailleurs, le dispositif d'imagerie à rayons X comprend classiquement un collimateur 4 ajustable positionné à la sortie de la source de rayons X 2. Le dispositif d'imagerie comprend, par ailleurs, des moyens de traitement 5 de l'image acquise par le récepteur d'image 1, lesdits moyens de traitement 5 comportant au moins un processeur 6, et des moyens de visualisation 7 connectés aux moyens de traitement 5 pour visualiser l'image acquise après traitement de l'image pour atténuer le bruit de cette dernière. Le système d'amélioration de l'image acquise à l'instant t, ladite image se présentant sous la forme d'une image digitale, c'est-à-dire définie par des pixels, et transmise aux moyens de traitement 5 comporte une mémoire dans laquelle sont enregistrées au fur et à mesure les images fluoroscopiques et des algorithmes de traitement de ces images introduits dans le processeur 6 pour permettre leur filtrage. En référence à la figure 2, le système d'amélioration de l'image acquise à l'instant t comporte des moyens 101 de traitements de l'image acquise à l'instant t dite lam (t) constitués de moyens 102 de détermination des déplacements vl et v2 de chacune des deux couches qui constituent l'image et de moyens 103 de détermination d'une image dite prédite lpred de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-1 et t-2, et des déplacements v1 et v2 de chacune des couches. Le système comporte, par ailleurs, des moyens 104 de génération de l'image visualisée correspondant à la somme 2884341 8 pondérée de l'image prédite Ipred (t) et de l'image de l'objet acquise à l'instant t lacq (t). On notera que les moyens 103 de détermination d'une image dite prédite Ipred de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-1 et t-2 et les moyens 102 de détermination des déplacements v1 et v2 des couches n'ont pas besoin de générer des images de chacune des couches. En effet, si l'on considère qu'une image I est constituée de deux couches Il et 12 se déplaçant respectivement suivant des déplacements v1 et v2, on peut écrire la différence D sous la forme suivante: D =I(p,t) I(p+v1,t 1) I(p+v2,t 1)+I(p+v1 +v2,t 2) dans laquelle p correspond à la position du pixel et t correspond au temps présent. En décomposant l'image I suivant deux couches Il et 12, nous obtenons 15 alors la relation suivante: D=11(p,t) I1(p+vl,t 1) I1(p+v2, t 1)+11(p+v, + v2,t 2) +I2(p,t) I2(p+v,,t l) I2(p+v2,t 1)+I2(p+v, +v2,t 2) Soit D = (Il(p,t) Il(p+v1,t 1)) (Il((p+v2),t 1) I1((p+v2)+v1,t 2)) +(I2(p, t) I2(p+ v2,t 1)) (I2(p+v1,t 1) I2((p+ v1)+ v2,t 2)) Les couches Il et 12 se déplaçant respectivement suivant des déplacements vl et v2, chacun des quatre termes entre parenthèses dans l'équation précédente est nul et la différence D est nulle. Par conséquent, conformément à la première équation, nous avons: I(p,t) =I(p+vl,t 1)+I(p+v2,t 1) I(p+v1 +v2,t 2) II apparaît ainsi que l'on peut déterminer lpred (t) en fonction des images acquises aux instants t-1 et t-2 et des déplacements v1 et v2 des deux couches sans pour autant avoir besoin de générer des images desdites couches. Les moyens 102 de détermination des déplacements v1 et v2 des couches consistent en un algorithme tel que décrit dans la publication Estimation of multiple motions using block-matching and Markov random fields , Ingo Stuke, Til Aach, Erhardt Barth et Cicero Mota, ln Visual Communications and Image Processing 2004, S. Panchanathan and B. Vasudev, EDS., IS&TISPIE 16th Annual Symposium Electronic Imaging, San Jose, California, USA, January 18-22, 2004. On observera que, dans cet algorithme, les déplacements v1 et v2 des couches sont considérés comme constant dans le temps ce qui ne correspond pas aux mouvements anatomiques du coeur, du diaphragme, etc... Selon une variante d'exécution du système de traitement d'une image suivant l'invention, les moyens de détermination des déplacements consistent en un algorithme permettant de déterminer les déplacements des couches suivant une fonction affine à 6 paramètres par couche, à savoir deux translations, deux rotations et deux homothéties. Ainsi, le vecteur déplacement à chaque point (x,y) pour la couche i peut s'écrire sous la forme suivante: vx(x,y)=c11+a=1.x+a;2.y et v.ry (x,y)=cy2 +ai,3.x+a. r,4'y Par conséquent, pour la totalité de l'image, il convient de minimiser la 20 fonction suivante à 12 paramètres: J(191,192)= Er(x,y, v1(x,y),v2(x,y)) 2 (x,y Z dans laquelle Zs indique le cadre de l'image et r(x,y,v1(x,y), v2(x,y)) est donné par: r(x, y,v1,v2)= I(x+v1x +v2x, y+vly + v2y,t -1)+ I(x,y,t +1) -I(x+v1x, y+v1y,t)-I(x+v2x y+v2y,t)= 0 La minimisation de la fonction à 12 paramètres est obtenue par application de la méthode de gradients conjugués sur la fonction linéarisée et appliquée sur trois images successives décomposées en pyramide gaussienne. Les incréments linéaires des paramètres sont alors déterminés par la méthode de GaussNewton dans chaque niveau de résolution de la 2884341 10 pyramide, en commençant par le niveau de plus basse résolution et en propageant les paramètres obtenus vers les niveaux de résolution plus élevés. Par ailleurs, les moyens 103 de détermination d'une image dite prédite lpred de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-1 et t-2 comportent, en référence à la figure 3, des moyens 105 de génération de deux images primaires Ipriml et Iprirri2 correspondant à l'image acquise à l'instant t-1 dont les pixels sont déplacés respectivement du déplacement vl et du déplacement v2 de chacune des couches et des moyens 106 de génération d'une image intermédiaire limer correspondant à la somme des deux images primaires Ipriml et Ipnm2. Lesdits moyens 101 de détermination d'une image dite prédite lpred de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-1 et t-2 comportent, par ailleurs, des moyens 107 de génération d'une troisième image primaire Iprim3 correspondant à l'image acquise à l'instant t-2 dont les pixels sont déplacés suivant la somme v1+ v2 des déplacements v1 et v2 de chacune des couches et des moyens 108 de génération de l'image dite prédite Ipred (t) de l'objet à l'instant t correspondant à l'image intermédiaire linter à laquelle est soustraite la troisième image primaire Iprim3. Les moyens 104 de génération de l'image visualisée Lis comporte deux multiplieurs 109 et 109', un premier multiplieur 109 pour pondérer par un facteur, de préférence 1/3, l'image acquise à l'instant t lacq (t) et un second multiplieur 109' pour multiplier par le complément à un du facteur précédent, c'est-à-dire 2/3, l'image prédite Ipred, pour additionner l'intensité ainsi pondérée de chaque pixel de l'image prédite Ipred avec l'intensité pondérée de chaque pixel de l'image de l'objet acquise à l'instant t Iacq. L'image visualisée Isis présente ainsi un bruit atténué améliorant la qualité de l'image. Selon une variante d'exécution du système suivant l'invention, ce dernier comporte des moyens pour diviser chaque image acquise en au moins deux cadres distincts, les images de chaque cadre étant traitées indépendamment l'une de l'autre de la même manière que précédemment. Afin d'obtenir une plus grande atténuation du bruit, les cadres sont de préférence des petits cadres réguliers qui sont assemblés à la fin du traitement de chaque cadre pour reconstituer l'image visualisée Lis. Il va de soi que la décomposition de l'image acquise à l'instant t peut être effectuée suivant plus de deux couches. En effet, conformément aux enseignements de la publication Estimation of multiple motions using block-matching and Markov random fields , Ingo Stuke, Til Aach, Erhardt Barth et Cicero Mota, ln Visual Communications and Image Processing 2004, S. Panchanathan and B. Vasudev, EDS., IS&T/SPIE 16th Annual Symposium Electronic Imaging, San Jose, California, USA, January 18-22, 2004, l'équation I(p,t)=I(p+v t-1)+I(p+v2,t 1) I(p+vl +v2,t 2) peut être généralisée pour N couches sous la forme suivante: I(p,t)=( 1) NI(p,v, +.. .+vN,t N) +( 1)N-'E{pickingN 1 differentvelocitiesxmongthev,} I(v;, +...+ v,N l,t 1 N) +. +( 1)' {pickingN k differentvelocitiesimongthevl}I(v,.l +...+v,Nk, t k) +. + E{Everyv, }I(v; , t -1) Le système suivant l'invention peut ainsi être généralisé pour la décomposition de l'image acquise à l'instant t lacq (t) en N couches. Le système comporte alors des moyens de détermination d'une image dite prédite de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-i, i étant un nombre entier positif supérieur ou égal à N, et des déplacements de chacune des couches, et des moyens de génération de l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. Les moyens de détermination de l'image prédite consiste en des moyens de détermination des déplacements v de chacune des N couches à partir des images acquises aux instants t-i, des moyens de génération de N images primaires correspondant aux images acquises aux instants t-i, dont les pixels sont déplacés des 2884341 12 différentes combinaisons des déplacements v de chacune des N couches, et des moyens de génération d'une image prédite correspondant à la somme en signes alternés des N images primaires, . Enfin, il est bien évident que les exemples que l'on vient de donner ne 5 sont que des illustrations particulières en aucun cas limitatives quant aux domaines d'application de l'invention. 2884341 13 REVENDICATIONS 1 Procédé d'amélioration d'une image digitale d'un objet définie par des pixels, acquise à l'instant t et générée à partir d'un détecteur de rayons X (1) recevant des rayons X émis par une source (2) caractérisé en ce qu'il consiste au moins dans les étapes suivantes de: -détermination d'une image dite prédite de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-i, i étant un nombre entier positif supérieur ou égal à 1, et des déplacements de chacune des couches constituant l'image acquise à l'instant t, le nombre de couches étant préalablement fixé et les déplacements de chacune desdites couches étant préalablement déterminés, - génération de l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t, afin d'atténuer le bruit de l'image de l'objet acquise à l'instant t. 2 Procédé suivant la revendication précédente, caractérisé en ce qu'il consiste au moins dans les étapes suivantes de: - fixation du nombre N de couches dont est constituée l'image acquise à l'instant t, -détermination des déplacements vn de chacune des N couches à partir des images acquises aux instants t-i, - génération de N images primaires correspondant aux images acquises aux instants t-i, dont les pixels sont déplacés respectivement de la somme des N déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-N, de la somme de N-1 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-1-N, de la somme de N-2 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-2-N et ainsi de suite jusqu'à chaque déplacement vn, pour l'image acquise à l'instant t-1, - génération d'une image dite prédite de l'objet à l'instant t correspondant à la somme en signes alternés des N images primaires, 2884341 14 - génération de l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. 3 Procédé suivant l'un quelconque des revendications 1 ou 2, caractérisé en ce qu'il comporte au moins les étapes suivantes de: - fixation du nombre N de couches, dont est constituée l'image acquise à l'instant t, à N=2 - détermination des déplacements v1 et v2 de chacune des deux couches à partir des images acquises aux instants t-1 et t-2, - génération de deux images primaires, une première image primaire correspondant à l'image acquise à l'instant t-1 dont les pixels sont déplacés du déplacement vl et une seconde image primaire correspondant à l'image acquise à l'instant t-1 dont les pixels sont déplacés du déplacement v2, - génération d'une image intermédiaire correspondant à la somme des deux images primaires, génération d'une troisième image primaire correspondant à l'image acquise à l'instant t-2 dont les pixels sont déplacés suivant la somme des déplacements v1 et v2 de chacune des couches, - génération d'une image dite prédite de l'objet à l'instant t correspondant à l'image intermédiaire à laquelle est soustraite la troisième 20 image secondaire, - génération de l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. 4 Procédé suivant l'une quelconque des revendications 1 à 3, caractérisé en ce qu'il comporte une étape d'application d'un filtre temporel à l'image visualisée résultant de la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. Procédé suivant la revendication 4 caractérisé en ce que la génération de l'image visualisée est obtenue par l'addition du tiers de l'intensité de chaque pixel de l'image prédite avec les deux tiers de l'intensité de chaque pixel de l'image de l'objet acquise à l'instant t. 6 Procédé suivant l'une quelconque des revendications 2 à 5, caractérisé en ce que chaque image acquise est divisée en au moins deux cadres et en ce que les étapes, de détermination des déplacements vn de chacune des N couches, de génération de N images primaires, de génération de l'image prédite et de génération de l'image visualisée sont appliquées indépendamment pour chaque cadre. 7 Système pour l'amélioration d'une image digitale d'un objet définie par des pixels, acquise à l'instant t et générée par un détecteur de rayons X (1) recevant des rayons X émis par une source (2), caractérisé en ce qu'il comporte des moyens (103) de détermination d'une image dite prédite de l'objet à l'instant t en fonction des images dudit objet acquises aux instants t-i, i étant un nombre entier positif supérieur ou égal à 1, des moyens (102) de détermination des déplacements de chacune des couches, le nombre de couches dont est constituée l'image acquise à l'instant t étant préalablement fixé et les déplacements de chacune desdites couches étant préalablement déterminés, et des moyens (104) de génération de l'image visualisée correspondant à la somme pondérée de l'image prédite et de l'image de l'objet acquise à l'instant t. 8 Système suivant la revendication 7, caractérisé en ce que les moyens (103) de détermination de l'image prédite consistent en des moyens (105, 106,107) de génération de N images primaires correspondant aux images acquises aux instants t-i, dont les pixels sont déplacés respectivement de la somme des N déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-N, de la somme de N-1 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-1-N, de la somme de N-2 déplacements vn de chacune des N couches, pour l'image acquise à l'instant t-2-N et ainsi de suite jusqu'à chaque déplacement vn, pour l'image acquise à l'instant t-1, et en des moyens (108) de génération de l'image dite prédite de l'objet à l'instant t. 9 Système suivant la revendication 8, caractérisé en ce que les moyens (104) de génération de l'image visualisée comporte des moyens (109, 109') pour additionner le tiers de l'intensité de chaque pixel de l'image prédite avec les deux tiers de l'intensité de chaque pixel de l'image de l'objet acquise à l'instant t. 2884341 16 Système suivant l'une quelconque des revendications 7 à 9, caractérisé en ce qu'il comporte des moyens pour diviser chaque image acquise en au moins deux cadres. 11 Application du procédé suivant l'une quelconque des 5 revendications 1 à 6 au traitement d'images fluoroscopiques.