FR2884341A1 - Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x - Google Patents
Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x Download PDFInfo
- Publication number
- FR2884341A1 FR2884341A1 FR0503549A FR0503549A FR2884341A1 FR 2884341 A1 FR2884341 A1 FR 2884341A1 FR 0503549 A FR0503549 A FR 0503549A FR 0503549 A FR0503549 A FR 0503549A FR 2884341 A1 FR2884341 A1 FR 2884341A1
- Authority
- FR
- France
- Prior art keywords
- image
- acquired
- time
- layers
- displacements
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000002059 diagnostic imaging Methods 0.000 title description 2
- 230000033001 locomotion Effects 0.000 claims abstract description 12
- 238000006073 displacement reaction Methods 0.000 claims description 53
- 230000002123 temporal effect Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 8
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000002594 fluoroscopy Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000014616 translation Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 241000764238 Isis Species 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/223—Analysis of motion using block-matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10121—Fluoroscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
La présente invention concerne un 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) remarquable 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.Un autre objet de l'invention concerne un système d'amélioration d'une image digitale mettant en oeuvre le procédé suivant l'invention.
Description
2884341 1
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
Claims (10)
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.
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0503549A FR2884341B1 (fr) | 2005-04-11 | 2005-04-11 | Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x |
| US11/401,043 US8121192B2 (en) | 2005-04-11 | 2006-04-10 | Method and system for treating noise in a generated image |
| DE102006017097A DE102006017097A1 (de) | 2005-04-11 | 2006-04-10 | Verfahren und System zur Behandlung des Rauschens in einem erzeugten Bild |
| JP2006108845A JP5072253B2 (ja) | 2005-04-11 | 2006-04-11 | 作成した画像のノイズを処理するための方法及びシステム |
| FR0603195A FR2884342B1 (fr) | 2005-04-11 | 2006-04-11 | Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0503549A FR2884341B1 (fr) | 2005-04-11 | 2005-04-11 | Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2884341A1 true FR2884341A1 (fr) | 2006-10-13 |
| FR2884341B1 FR2884341B1 (fr) | 2008-05-02 |
Family
ID=35355300
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0503549A Expired - Fee Related FR2884341B1 (fr) | 2005-04-11 | 2005-04-11 | Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US8121192B2 (fr) |
| JP (1) | JP5072253B2 (fr) |
| DE (1) | DE102006017097A1 (fr) |
| FR (1) | FR2884341B1 (fr) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11006095B2 (en) | 2015-07-15 | 2021-05-11 | Fyusion, Inc. | Drone based capture of a multi-view interactive digital media |
| US10242474B2 (en) * | 2015-07-15 | 2019-03-26 | Fyusion, Inc. | Artificially rendering images using viewpoint interpolation and extrapolation |
| US10147211B2 (en) | 2015-07-15 | 2018-12-04 | Fyusion, Inc. | Artificially rendering images using viewpoint interpolation and extrapolation |
| US10222932B2 (en) | 2015-07-15 | 2019-03-05 | Fyusion, Inc. | Virtual reality environment based manipulation of multilayered multi-view interactive digital media representations |
| US11095869B2 (en) | 2015-09-22 | 2021-08-17 | Fyusion, Inc. | System and method for generating combined embedded multi-view interactive digital media representations |
| US11783864B2 (en) | 2015-09-22 | 2023-10-10 | Fyusion, Inc. | Integration of audio into a multi-view interactive digital media representation |
| US11202017B2 (en) | 2016-10-06 | 2021-12-14 | Fyusion, Inc. | Live style transfer on a mobile device |
| US10437879B2 (en) | 2017-01-18 | 2019-10-08 | Fyusion, Inc. | Visual search using multi-view interactive digital media representations |
| US10313651B2 (en) | 2017-05-22 | 2019-06-04 | Fyusion, Inc. | Snapshots at predefined intervals or angles |
| US11069147B2 (en) | 2017-06-26 | 2021-07-20 | Fyusion, Inc. | Modification of multi-view interactive digital media representation |
| US10592747B2 (en) | 2018-04-26 | 2020-03-17 | Fyusion, Inc. | Method and apparatus for 3-D auto tagging |
| EP3726459B1 (fr) * | 2019-04-17 | 2022-02-16 | Leica Instruments (Singapore) Pte. Ltd. | Circuit, procédé et programme de réglage de rapport signal-bruit |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1030268A1 (fr) * | 1999-02-18 | 2000-08-23 | Ge Medical Systems Sa | Méthode pour la réduction du bruit de rayons X |
-
2005
- 2005-04-11 FR FR0503549A patent/FR2884341B1/fr not_active Expired - Fee Related
-
2006
- 2006-04-10 US US11/401,043 patent/US8121192B2/en active Active
- 2006-04-10 DE DE102006017097A patent/DE102006017097A1/de not_active Withdrawn
- 2006-04-11 JP JP2006108845A patent/JP5072253B2/ja not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1030268A1 (fr) * | 1999-02-18 | 2000-08-23 | Ge Medical Systems Sa | Méthode pour la réduction du bruit de rayons X |
Non-Patent Citations (4)
| Title |
|---|
| AACH T ET AL: "Bayesian motion estimation for temporally recursive noise reduction in X-ray fluoroscopy", PHILIPS JOURNAL OF RESEARCH, ELSEVIER, AMSTERDAM, NL, vol. 51, no. 2, 1998, pages 231 - 251, XP004126961, ISSN: 0165-5817 * |
| CLOSE R A ET AL: "LAYER DECOMPOSITION OF CORONARY ANGIOGRAMS", PROCEEDINGS OF THE SPIE, SPIE, BELLINGHAM, VA, US, vol. 3979, 2000, pages 1230 - 1234, XP008027880, ISSN: 0277-786X * |
| STUKE I ET AL: "Estimation of multiple motions using block matching and Markov random fields", PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING SPIE-INT. SOC. OPT. ENG USA, vol. 5308, no. 1, 2004, pages 486 - 496, XP002355913, ISSN: 0277-786X * |
| V. AUVRAY, J. LIENARD, P. BOUTHEMY: "Multiresolution parametric estimation of transparent motions and denoising of fluoroscopic images", INT. CONF. ON MEDICAL IMAGE COMPUTING AND COMPUTER ASSISTED INTERVENTION (MICCAI'05), October 2005 (2005-10-01), XP008056201 * |
Also Published As
| Publication number | Publication date |
|---|---|
| JP5072253B2 (ja) | 2012-11-14 |
| JP2006289094A (ja) | 2006-10-26 |
| DE102006017097A1 (de) | 2006-11-02 |
| US20070058880A1 (en) | 2007-03-15 |
| US8121192B2 (en) | 2012-02-21 |
| FR2884341B1 (fr) | 2008-05-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2884341A1 (fr) | Procede et systeme d'amelioration d'une image digitale generee a partir d'un detecteur de rayons x | |
| EP0994442A2 (fr) | Une méthode de poursuite du mouvement d'un organe et de suppression de parasites pour les systèmes d'imagerie par tomographie assistée par ordinateur | |
| EP0936575B1 (fr) | Procédé de traitement d'images pour l'estimation de mouvement dans une séquence d'images | |
| EP0113605B1 (fr) | Procédé et appareil pour former une image radiologique d'un corps | |
| EP0871143B1 (fr) | Système de traitement d'une séquence d'images bruitées, et appareil d'examen médical incluant un tel système | |
| FR2736182A1 (fr) | Procede de traitement d'images pour la reduction du bruit dans une image d'une sequence d'images numeriques et dispositif mettant en oeuvre ce procede | |
| EP0635806A1 (fr) | Procédé de traitement d'images numériques pour déterminer localement le centre et la largeur d'objets en forme de rubans contrastant sur un fond | |
| FR2993768A1 (fr) | Systeme et procede d'imagerie echographique | |
| FR2736181A1 (fr) | Procede de traitement d'images pour la reduction du bruit dans une image d'une sequence d'images numeriques et dispositif mettant en oeuvre ce procede | |
| FR2602602A1 (fr) | Dispositif et procede de correction d'image vue a vue pour le mouvement d'un objet | |
| FR2798804A1 (fr) | Procede et appareil pour reconstruire des donnees d'image acquises par un systeme d'imagerie a rayons x par tomosynthese | |
| FR2978273A1 (fr) | Procede de reduction du bruit dans une sequence d'images fluoroscopiques par filtrage temporel et spatial | |
| FR2924254A1 (fr) | Procede de traitement d'images en radioscopie interventionnelle | |
| FR2848007A1 (fr) | Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal | |
| JPH11501195A (ja) | 画像内のノイズ減少方法 | |
| FR2884948A1 (fr) | Procede et dispositif de reduction du bruit dans une sequence d'images fluoroscopiques | |
| EP0110787A2 (fr) | Procédé et appareil pour former des images en coupe d'un corps | |
| FR2884340A1 (fr) | Procede et dispositif de traitement d'images en angiographie soustraite | |
| EP0101336B1 (fr) | Appareil de génération d'un signal vidéo traité, à contraste amélioré pour radiologie, et procédé pour traiter un tel signal | |
| EP0926626A1 (fr) | Procédé de traitement d'image pour réduire le bruit dans une image d'une séquence de trois images temporelles et appareil d'imagerie médicale mettant en oeuvre un tel procédé | |
| FR2846830A1 (fr) | Procede de filtrage spatio-temporel du bruit fluoroscopique | |
| FR2881941A1 (fr) | Procede de determination des parametres geometriques d'un dispositif d'imagerie par rayon x | |
| FR2848806A1 (fr) | Procede de calibrage d'un appareil d'imagerie radiologique necessitant un nombre limite d'acquisitions | |
| FR2799029A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet en mouvement, en particulier une image tridimensionnelle de vaisseaux du coeur humain | |
| US10068319B2 (en) | Method for noise reduction in an image sequence |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| ST | Notification of lapse |
Effective date: 20141231 |