FR2927719A1 - Procede de traitement d'images obtenues par tomosynthese et dispositif associe - Google Patents
Procede de traitement d'images obtenues par tomosynthese et dispositif associe Download PDFInfo
- Publication number
- FR2927719A1 FR2927719A1 FR0800896A FR0800896A FR2927719A1 FR 2927719 A1 FR2927719 A1 FR 2927719A1 FR 0800896 A FR0800896 A FR 0800896A FR 0800896 A FR0800896 A FR 0800896A FR 2927719 A1 FR2927719 A1 FR 2927719A1
- Authority
- FR
- France
- Prior art keywords
- thickness
- tomosynthesis
- contrast
- image
- limiting
- 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 31
- 238000012545 processing Methods 0.000 title abstract description 13
- 230000007547 defect Effects 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims abstract description 3
- 238000002059 diagnostic imaging Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 4
- 238000003672 processing method Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 2
- 238000000605 extraction Methods 0.000 claims 2
- 230000003902 lesion Effects 0.000 description 19
- 210000000056 organ Anatomy 0.000 description 19
- 238000003745 diagnosis Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 239000002872 contrast media Substances 0.000 description 3
- 238000011282 treatment Methods 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
L'invention concerne un procédé de traitement d'images obtenues par tomosynthèse comprenant : une acquisition (S0) d'une pluralité d'images de projection 2D d'une région d'intérêt d'un patient ; une reconstruction (S1) d'une image numérique 3D à partir des images de projection 2D acquises ; caractérisé en ce qu'il comprend : une détection (S2) d'un objet dans l'image reconstruite ; une estimation (S3) d'une épaisseur limite (T) caractéristique d'un défaut de contraste des voxels pour un diamètre (D) de l'objet ; une estimation (S4) de l'épaisseur (tp) de l'objet ; une comparaison (S5) de l'épaisseur (tp) de l'objet à l'épaisseur limite (T) ; et en ce que si l'épaisseur (tp) de l'objet est inférieure à l'épaisseur limite (T), le procédé comprend en outre une application (S6) au moins aux voxels de l'objet dans l'image reconstruite d'un facteur correctif multiplicatif (Kp) égal au ratio entre l'épaisseur limite (T) et l'épaisseur (tp) de l'objet.
Description
2927719 PROCEDE DE TRAITEMENT D'IMAGES OBTENUES PAR TOMOSYNTHESE ET DISPOSITIF ASSOCIE
DOMAINE TECHNIQUE GENERAL L'invention concerne un procédé d'imagerie médicale de type tomosynthèse. En particulier, elle concerne un procédé de traitement d'images obtenues par tomosynthèse permettant d'améliorer la qualité des images nécessaires au praticien radiologue pour élaborer son diagnostic. i0 ETAT DE LA TECHNIQUE
La tomosynthèse est une modalité d'imagerie médicale dont les principes essentiels sont rappelés ci-dessous et illustrés sur la figure 2. is On a représenté sur la figure 2 schématiquement l'acquisition d'images 2D de l'organe et la reconstruction d'une image 3D de cet organe par tomosynthèse.. Des rayons X, R, issus d'une source S sont émis selon différentes angulations (1, ..., i, ..., n) vers l'organe O. Après avoir traversé l'organe, ils 20 sont détectés par le détecteur Det formant un ensemble d'images de projection Al, ..., Ai,..., An. Il est à noter qu'il y a autant d'images 2D acquises que d'angulations considérées. En tomosynthèse simple énergie le fond des images 2D est supposé localement uniforme autour de la lésion d'intérêt. 25 On note que dans le cas de la tomosynthèse double énergie deux images sont acquises avec deux spectres d'énergies différents pour chaque angulation considérée et que l'on a injecté au préalable dans le patient un produit de contraste. L'acquisition est mise en oeuvre par un détecteur Det situé en regard 30 d'une source de rayons X, par exemple une caméra numérique ou un détecteur solide à base de Silicium amorphe ou de Selenium amorphe. 2 2927719 Une application de la tomosynthèse est la détection et la caractérisation d'une lésion dans un organe, par exemple une lésion cancéreuse. Dans un mode de réalisation de tomosynthèse les images de 5 projection 2D sont utilisées pour reconstruire une image 3D, le praticien radiologue interprète cette image en fonction des différences de contraste observées. Le volume reconstruit contient une pluralité de voxels quantifiant l'atténuation des rayons X par la matière rencontrée. io Dans le cas d'un organe présentant une lésion celle-ci présente une atténuation différente de l'organe. Le résultat sur l'image 3D est une différence de contraste lui permettant de détecter la lésion. Toutefois, un problème intrinsèque à la tomosynthèse est l'ouverture 15 angulaire limitée de l'acquisition des images de projection qui induit en particulier une résolution limitée selon la dimension perpendiculaire au plan du détecteur. On considère ici que le plan du détecteur est selon les axes x, y d'un repère orthonormé et que la dimension perpendiculaire est selon l'axe z de ce même repère. 20 La mauvaise résolution dans la dimension z perpendiculaire au plan du détecteur est due à l'ouverture angulaire couverte par le tube à rayons X limitée en comparaison avec la tomographie. En effet, pour obtenir une résolution isotrope (identique selon toutes les dimensions) il faudrait acquérir continument un ensemble d'images de 25 projection tout autour de l'organe. On comprend toutefois que ceci est en pratique impossible pour l'examen de tomosynthèse. On connaît des procédés permettant d'améliorer la résolution des images acquises pour discriminer des éléments proches. Mais outre ce problème, en deçà d'une certaine épaisseur de lésion 30 (prise sur l'axe z), la valeur des voxels constituant cette lésion dans un plan image passant par la lésion dans le volume reconstruit n'est pas suffisante pour une identification améliorée du type de lésion par le radiologue (lésion bénigne ou maligne).
Ceci a pour effet d'introduire une incertitude sur la caractérisation de la lésion ce qui peut évidemment compromettre le diagnostic.
En d'autres termes l'image reconstruite n'est fiable que pour des lésions présentant une épaisseur minimale relativement à leur taille selon un plan qui est parallèle au détecteur et qui passe par l'objet. PRESENTATION DE L'INVENTION L'invention permet d'améliorer les images obtenues par tomosynthèse et notamment celles obtenues par émission de rayons X avec injection de produit de contraste.
Et plus particulièrement, l'invention permet de s'affranchir du défaut
15 de contraste due à l'ouverture angulaire limitée caractérisant le parcours de la source de rayons X autour de la zone d'intérêt lors de l'acquisition d'images de projection 2D en tomosynthèse.
En particulier, elle évite que la valeur des voxels de l'image 3D ne soit dégradée lorsque l'épaisseur de l'objet à caractériser est inférieure à une
20 épaisseur limite pour une taille d'objet donnée mesurée dans un plan qui est parallèle au détecteur et qui passe par l'objet.
L'invention permet en effet de s'affranchir de cette limite pour obtenir des images 3D contenant des voxels quantifiant la radio-opacité de la zone d'intérêt de manière précise.
25 Ainsi selon un premier aspect, l'invention concerne un procédé de traitement d'images obtenues par tomosynthèse comprenant une acquisition d'une pluralité d'images de projection 2D d'une région d'intérêt d'un patient; une reconstruction d'une image numérique 3D à partir des images de projection 2D acquises.
30 Le procédé de l'invention est caractérisé en ce qu'il comprend : une détection d'un objet dans l'image reconstruite ; une estimation d'une 4 2927719 épaisseur limite caractéristique d'un défaut de contraste des voxels pour un diamètre de l'objet ; une estimation de l'épaisseur de l'objet ; une comparaison de l'épaisseur de l'objet à l'épaisseur limite ; et en ce que si l'épaisseur de l'objet est inférieure à l'épaisseur limite, le procédé comprend 5 en outre une application au moins aux voxels de l'objet dans l'image reconstruite d'un facteur correctif multiplicatif égal au ratio entre l'épaisseur limite et l'épaisseur de l'objet. On obtient donc avec le procédé de l'invention une image 3D de l'organe contenant l'objet à caractériser dont les voxels sont proportionnels à w la grandeur physique caractéristique de la radio-opacité de l'objet. Selon un deuxième aspect, l'invention concerne un système d'imagerie médicale comprenant des moyens apte à mettre en oeuvre un procédé selon le premier aspect de l'invention. Et selon un troisième aspect, l'invention concerne un programme 15 d'ordinateur caractérisé en ce qu'il comprend des instructions machine pour la mise en oeuvre d'un procédé selon le premier aspect de l'invention.
PRESENTATION DES FIGURES
20 D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit laquelle est purement illustrative et non limitative et doit être lue en regard des dessins annexés sur lesquels outre la figure 1 déjà discutée : - la figure 2 illustre schématiquement un système d'imagerie 25 médicale ; - la figure 3 illustre schématiquement des étapes du procédé selon l'invention ; - la figure 4 illustre le contraste volumique de l'objet dans le plan image reconstruit passant par ledit objet en fonction de l'épaisseur du 30 modèle d'objet ; 5 2927719 -la figure 5 illustre la taille du modèle de lésion dans le plan du détecteur en fonction de l'épaisseur limite - la figure 6 illustre schématiquement des sous-étapes du procédé de l'invention ; 5 - la figure 7 illustre le principe d'estimation de l'épaisseur limite.
DESCRIPTION DETAILLEE DE L'INVENTION
Système d'imaqerie médicale io Sur la figure 2, on a illustré, schématiquement, un système d'imagerie médicale 1 permettant l'acquisition d'images 2D de projection pour la reconstruction d'une image 3D d'un organe lésé. Un tel système peut être un appareil mammographique pour la détection et la caractérisation de lésions dans le cas du dépistage, du 15 diagnostic et du traitement du cancer du sein. Le système d'imagerie médical 1 comprend un système d'acquisition d'images 3, un système de traitement d'images 5 et un système d'affichage 4. Le système d'acquisition 3 permet d'acquérir une pluralité d'images de 20 projection 2D d'une région d'intérêt û d'un organe û d'un patient. Le système de traitement 5 est par exemple un ordinateur. Le système de traitement 5 est couplé à des moyens mémoires 6 qui peuvent être intégrés ou séparés du système de traitement 5. Ces moyens peuvent être formés par un disque dur ou tout autre moyen de stockage amovible (un cd-rom, une disquette 25 etc.). Le système d'acquisition d'images 3 est par exemple un système d'acquisition par rayons X ce dernier comprenant tout moyen connu permettant l'émission de rayons X sur l'objet 2 et l'acquisition d'images résultantes. 30 Le système de traitement permet d'une part la reconstruction de l'image 3D ainsi que le traitement d'images tel que nous allons le décrire. 6 2927719 Description générale du procédé de traitement d'image Dans ce qui suit on considère que l'on cherche à analyser l'image d'un organe d'un patient pour éventuellement détecter une lésion. La figure 3 illustre schématiquement les étapes du procédé de 5 traitement d'images. Le procédé de traitement vise à traiter une image 3D représentative de la prise de contraste de l'organe lésée obtenue au moyen d'un procédé de reconstruction de type tomographique à partir d'une pluralité d'images de projection 2D acquises suite à une émission de rayons X sur un organe lésé (émission décrite précédemment). io L'image 3D comprend une pluralité de voxels présentant une valeur de contraste quantifiant la radio-opacité de l'organe lésé. Dans le cas de la tomosynthèse simple énergie, le contraste de l'objet est proportionnel au coefficient d'atténuation linéaire de l'objet. Et dans le cas de la tomosynthèse double énergie, après 15 recombinaison pour chaque angle de projection des images basse et haute énergies pour produire une image de prise de contraste (par exemple au moyen d'une soustraction logarithmique pondérée connue de l'homme de l'art) et reconstruction de l'image 3D à partir de ces images de projection 2D recombinées, le contraste de l'objet est proportionnel à la concentration C du 20 produit de contraste. Le produit de contraste injecté est par exemple de l'iode ou plus généralement un produit iodé. Le procédé comprend une étape SO au cours de laquelle on acquiert une pluralité d'images de projection 2D de l'organe d'un patient. Le procédé comprend également une étape S1 au cours de laquelle 25 on reconstruit une image 3D de l'organe. Le procédé comprend une étape de détection S2 d'un objet dans l'image 3D reconstruite. On entend par détection de l'objet le fait que dans une zone donnée de l'image 3D reconstruite des voxels agglomérés présentent tous une 30 valeur de contraste sensiblement différente de celle des autres voxels de la zone donnée. Les voxels agglomérés constituent l'objet détecté. 7 2927719 Dans une forme de réalisation l'image 3D est constituée d'une pluralité d'images 2D dites plans de coupe, parallèles au plan du détecteur et empilées selon l'axe z perpendiculaire au plan du détecteur, ladite zone donnée est donc elle-même constituée d'une pluralité de zones données 2D 5 présente dans les plans de coupe. L'objet est alors détecté en observant sur une série de plans de coupe, des pixels agglomérés présentant une valeur de contraste sensiblement différente de celle des autres pixels de la zone 2D donnée. L'accumulation selon z des pixels agglomérés sur la série de plans de coupe lo constitue l'objet détecté. Dans la suite de la description, le plan P fera référence à un plan sensiblement parallèle au détecteur et passant par l'objet détecté. Le principe du procédé repose sur une estimation S3 d'une épaisseur limite T à partir de laquelle les valeurs de contraste des voxels ne sont plus 15 proportionnelles à la radio-opacité de l'objet. Le procédé évite que la valeur des voxels de l'image 3D recombinée ne soit dégradée lorsque l'épaisseur de l'objet à caractériser est inférieure à une épaisseur limite pour une taille d'objet donnée mesurée dans un plan qui est parallèle au détecteur et qui passe par l'objet. 20 Sur la figure 4 on a représenté, pour un objet de section carrée de dimensions 3x3 mm dans le plan P, les valeurs de contraste de l'objet dans un plan parallèle au détecteur et passant par l'objet (par exemple le plan P) en fonction de l'épaisseur t de l'objet (exprimée en dixième de millimètres) en considérant un coefficient d'atténuation linéaire constant û 25 caractéristique de la radio-opacité de l'objet. Plus précisément on utilise un objet de forme parallélépipédique de section carrée dans le plan P et on fait varier son épaisseur t (sa dimension selon l'axe z).
Pour obtenir une courbe telle qu'illustrée sur la figure 4, pour chaque 30 objet d'épaisseur t donnée et de forme donnée, on procède à une coupe dans le volume reconstruit passant par l'objet (sensiblement perpendiculaire à l'axe z), par exemple le plan P.
Pour des dimensions d'objet données dans le plan P, on peut tracer les valeurs de contraste en fonction de l'épaisseur de l'objet t.
s On distingue sur la figure 4 deux zones : une zone non linéaire I et une zone linéaire Il. Ces deux zones sont de part et d'autres d'une valeur t=T.
En particulier en tomosynthèse simple énergie :
^ dans la zone linéaire Il, le contraste de l'objet est proportionnel 10 au coefficient d'atténuation linéaire p de l'organe ;
^ dans la zone non linéaire I il a été montré que le contraste est
essentiellement proportionnel à ( M/T. On comprend évidemment que dans la zone non linéaire I il y a une
15 quantification du contraste imparfaite, on peut dire qu'il y alors un défaut de contraste. On comprend en outre que pour rétablir la linéarité il suffit de multiplier les valeurs de contraste dans cette zone par un facteur K=T/t.
Il est à noter que dans le cas de la tomosynthèse à double énergie le coefficient d'atténuation linéaire est à remplacer par la concentration du
20 produit de contraste dans l'organe.
On montre alors de la même manière que
^ dans la zone linéaire Il, le contraste de l'objet est proportionnel à la concentration C du produit de contraste ;
^ dans la zone non linéaire I il a été montré que le contraste est 25 essentiellement proportionnel à (C.t)/T. En conséquence dans le procédé pour un objet de taille D dans le plan P on compare S4 l'épaisseur tp de l'objet à l'épaisseur limite T qui dépend de D et si l'épaisseur tp est inférieure à l'épaisseur limite T on
30 multiplie S6 les valeurs du contraste de l'objet dans l'image 3D reconstruite 9 2927719 par le coefficient Kp=T/tp calculé S4 (ceci étant valable pour la tomosynthèse à simple ou double énergie).
La figure 5 illustre une relation linéaire entre la taille de l'objet dans le s plan P du détecteur et l'épaisseur limite T. Une telle courbe a été obtenue pour une ouverture angulaire comprise entre -20° et +20°. Pour une ouverture angulaire donnée de l'appareil de tomosynthèse, il y a donc relation de proportionnalité entre la taille de l'objet dans le plan P et l'épaisseur limite T. Le coefficient de proportionnalité est stocké avec io l'ouverture angulaire à laquelle il s'applique dans les moyens de traitement du système d'imagerie médicale.
Ainsi, pour une ouverture angulaire donnée, si l'on est capable d'estimer la taille D de l'objet dans le plan P, on peut en déduire l'épaisseur 15 limite T qui va nous permettre de corriger les valeurs de contraste. Il faut avoir construit la courbe donnant le contraste en fonction de l'épaisseur de l'objet (voir figure 4) au moins pour deux valeurs de D afin de pouvoir en déduire le coefficient de proportionnalité entre T et D.
20 Ainsi pour rétablir la proportionnalité entre contraste des voxels de la lésion et coefficient d'atténuation linéaire (ou concentration du produit de contraste) il faut estimer d'une part S3 l'épaisseur limite T et d'autre part S4 l'épaisseur tp de l'objet à caractériser.
25 Estimation de l'épaisseur limite T et de l'épaisseur tp de l'objet à caractériser Pour estimer S3 l'épaisseur limite T on initialise S30 le modèle d'objet (voir figure 7) pour une taille D donnée dans le plan P. Pour ce faire, on peut prendre une hypothèse arbitraire quant à la 30 forme de l'objet : ellipsoïdale, ou parallélépipédique, par exemple.
Le paramètre important étant surtout la dimension selon l'axe z perpendiculaire au plan du détecteur : l'épaisseur t du modèle d'objet. On a représenté sur la figure 7 un modèle théorique d'un volume 50 contenant un modèle d'objet 51 d'épaisseur t.
Sur ce modèle théorique on simule l'acquisition d'une pluralité d'images de projection 2D suivie d'une reconstruction d'images 3D dans les mêmes conditions que les conditions réelles d'un examen radiologique. On obtient alors une reconstruction imparfaite du modèle d'objet et (voir figure 4) on sait que les voxels de l'image reconstruite peuvent ne pas i0 présenter un contraste proportionnel à la radio-opacité de l'objet (à cause de l'ouverture angulaire limitée de la trajectoire de la source rayons X autour de l'objet que l'on image). Pour obtenir l'épaisseur limite T on va répéter S34 les opérations en faisant à chaque fois varier l'épaisseur t de l'objet, et à chaque fois extraire 15 un plan 53 de coupe sensiblement perpendiculaire à l'axe z et passant par l'objet, par exemple le plan P. De ces plans de coupe, on peut alors extraire le contraste de l'objet et tracer une courbe similaire à celle de la figure 4 permettant d'en déduire l'épaisseur limite T pour une taille D de l'objet et l'ouverture angulaire 20 données.
On note que l'on obtient un réseau de courbes pour différentes valeurs de la taille D de l'objet que l'on peut stocker dans les moyens de traitement du système d'imagerie médicale. Une fois l'épaisseur limite T estimée on estime l'épaisseur tp de l'objet à caractériser.
On peut à cet effet utiliser un dispositif permettant d'obtenir les 30 caractéristiques dimensionnelles de l'objet, par exemple en utilisant un système d'imagerie secondaire tel que l'imagerie par résonance magnétique, 25 11 2927719 échographie etc. Si l'on utilise par exemple un appareil d'échographie, en réalisant l'acquisition d'une image de coupe de l'objet dans un plan perpendiculaire au plan x,y du détecteur de l'appareil de tomosynthèse, il devient alors possible de mesurer l'épaisseur tp de l'objet dans ce plan de 5 coupe ultrasonore.
Et l'épaisseur tp de l'objet à caractériser estimée permet ensuite la détermination du facteur correctif Kp.
io On peut également utiliser les caractéristiques dimensionnelles pour donner au modèle d'objet une forme générale qui soit proche de l'objet à caractériser. On peut par exemple considérer le modèle d'objet comme sphérique. Dans ce cas, la mesure de la taille D de l'objet dans le plan P du détecteur 15 nous donne accès directement à une estimation de l'épaisseur tp, avec tp=D. Alternativement, si l'on suppose que l'objet est un ellipsoïde aplati selon l'axe z dont l'épaisseur est liée à la taille D dans le plan P selon un rapport a, la mesure de la taille D de l'objet dans le plan P nous donne accès directement à une estimation de l'épaisseur tp de l'objet, avec tp=a.D. 20 Ainsi d'une part grâce à l'estimation de l'épaisseur limite T et de l'épaisseur tp de l'objet on peut calculer le facteur correctif Kp et corriger le contraste de l'objet dans l'image reconstruite.
25 Le contraste de la lésion dans le cas où une lésion est à caractériser est correctement quantifié et les informations disponibles pour que le praticien fasse son diagnostic sont améliorées.
II a été décrit un exemple dans le cas de l'imagerie médicale mais on 30 peut utiliser un tel procédé de traitement lorsque l'on souhaite caractériser un matériau et notamment détecter ses défauts. 25 12
Claims (8)
1. Procédé de traitement d'images obtenues par tomosynthèse comprenant une acquisition (SO) d'une pluralité d'images de projection 2D d'une région d'intérêt d'un patient; une reconstruction (Si) d'une image numérique 3D à partir des images de projection 2D acquises ; caractérisé en ce qu'il comprend io une détection (S2) d'un objet dans l'image reconstruite ; -une estimation (S3) d'une épaisseur limite (T) caractéristique d'un défaut de contraste des voxels pour un diamètre (D) de l'objet ; une estimation (S4) de l'épaisseur (tp) de l'objet ; - une comparaison (S5) de l'épaisseur (tp) de l'objet à l'épaisseur 15 limite (T) ; et en ce que si l'épaisseur (tv) de l'objet est inférieure à l'épaisseur limite (T), le procédé comprend en outre une application (S6) au moins aux voxels de l'objet dans l'image reconstruite d'un facteur correctif multiplicatif (Ko) égal au ratio 20 entre l'épaisseur limite (T) et l'épaisseur (tp) de l'objet.
2. Procédé selon la revendication 1, caractérisé en ce que le facteur correctif est fonction de l'épaisseur limite (T) à partir de laquelle les voxels ne sont plus à variation linéaire par rapport à la radio-opacité de l'objet détecté.
3. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'estimation (S3) de l'épaisseur limite (T) comprend pour une dimension donnée de l'objet une détermination (S30) d'un volume théorique d'une zone d'intérêt contenant un modèle de l'objet d'épaisseur (t) et de forme donnés 13 2927719 - un calcul (S31) d'un modèle de volume reconstruit issu du volume théorique par reconstruction par tomosynthèse à partir d'une pluralité d'images 2D issues du volume théorique calculé ; une extraction (S32) d'un plan de coupe (Pc) interceptant l'objet s d'épaisseur donnée (t) ; - une répétition (S33) des deux étapes (S30, S31) de calcul et de l'étape d'extraction (S32) pour une pluralité de valeurs d'épaisseur donnée (t) ; - une détermination (S34) de l'épaisseur limite (T) à partir du contraste l0 relevé pour chaque épaisseur considérée.
4. Procédé selon la revendication précédente, caractérisé en ce que le modèle de l'objet est choisi parmi le groupe suivant : sphérique, cubique ou plus généralement ellipsoïdale ou parallélépipédique.
5. Procédé selon l'une des revendications précédentes, caractérisé en ce que les images de projection sont obtenues par tomosynthèse simple énergie. 20
6. Procédé selon l'une des revendications précédentes, caractérisé en ce que les images sont obtenues par tomosynthèse double énergie avec injection d'un produit de contraste..
7. Système d'imagerie médicale comprenant des moyens apte à mettre 25 en oeuvre un procédé selon l'une des revendications là 6.
8. Programme d'ordinateur caractérisé en ce qu'il comprend des instructions machine pour la mise en oeuvre d'un procédé selon l'une des revendications 1 à 6. 15 30
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0800896A FR2927719B1 (fr) | 2008-02-19 | 2008-02-19 | Procede de traitement d'images obtenues par tomosynthese et dispositif associe |
US12/371,968 US8798348B2 (en) | 2008-02-19 | 2009-02-17 | Method for processing images obtained by tomosynthesis and associated device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0800896A FR2927719B1 (fr) | 2008-02-19 | 2008-02-19 | Procede de traitement d'images obtenues par tomosynthese et dispositif associe |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2927719A1 true FR2927719A1 (fr) | 2009-08-21 |
FR2927719B1 FR2927719B1 (fr) | 2010-03-26 |
Family
ID=39722670
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0800896A Active FR2927719B1 (fr) | 2008-02-19 | 2008-02-19 | Procede de traitement d'images obtenues par tomosynthese et dispositif associe |
Country Status (2)
Country | Link |
---|---|
US (1) | US8798348B2 (fr) |
FR (1) | FR2927719B1 (fr) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015217141A1 (de) * | 2015-09-08 | 2017-03-09 | Siemens Healthcare Gmbh | Erzeugen von kontrastverstärkten Bilddaten von zu untersuchendem Brustgewebe |
CN108135558B (zh) | 2015-10-05 | 2021-09-03 | 皇家飞利浦有限公司 | 用于表征身体部分的特征的装置 |
US9943280B2 (en) | 2016-03-07 | 2018-04-17 | General Electric Company | Breast tomosynthesis with flexible compression paddle |
US10893842B2 (en) | 2018-02-08 | 2021-01-19 | Covidien Lp | System and method for pose estimation of an imaging device and for determining the location of a medical device with respect to a target |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070286470A1 (en) * | 2006-06-07 | 2007-12-13 | Sylvain Bernard | Method for the processing of radiological images in tomosynthesis for a detection of radiological signs |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6041952A (ja) | 1983-08-19 | 1985-03-05 | 株式会社東芝 | X線フアントム |
US4792900A (en) * | 1986-11-26 | 1988-12-20 | Picker International, Inc. | Adaptive filter for dual energy radiographic imaging |
US5236363A (en) | 1991-08-09 | 1993-08-17 | General Electric Company | Phantom for simulating an x-ray exam patient |
EP0885616A1 (fr) * | 1997-06-20 | 1998-12-23 | Schering Aktiengesellschaft | Utilisation d'agents de contraste à administration intraveneuse, et appareil pour la mammographie |
US5852647A (en) * | 1997-09-24 | 1998-12-22 | Schick Technologies | Method and apparatus for measuring bone density |
FR2803069B1 (fr) * | 1999-12-28 | 2002-12-13 | Ge Medical Syst Sa | Procede et systeme de compensation de l'epaisseur d'un organe |
FR2814666A1 (fr) * | 2000-07-07 | 2002-04-05 | Ge Med Sys Global Tech Co Llc | Procede et appareil d'examen d'un sein par injection d'un produit de contraste |
JP2002218322A (ja) | 2001-01-16 | 2002-08-02 | Fuji Photo Film Co Ltd | 放射線画像処理方法および装置 |
DE10143131B4 (de) | 2001-09-03 | 2006-03-09 | Siemens Ag | Verfahren zur Ermittlung von Dichte- und Ordnungszahlverteilungen bei radiographischen Untersuchungsverfahren |
US6707878B2 (en) * | 2002-04-15 | 2004-03-16 | General Electric Company | Generalized filtered back-projection reconstruction in digital tomosynthesis |
AU2003252103A1 (en) | 2002-07-23 | 2004-02-09 | Ge Medical Systems Global Technology Company, Llc | Methods and systems for detecting components of plaque |
JP2004208752A (ja) * | 2002-12-27 | 2004-07-29 | Konica Minolta Holdings Inc | 乳房画像撮影装置 |
US6980624B2 (en) * | 2003-11-26 | 2005-12-27 | Ge Medical Systems Global Technology Company, Llc | Non-uniform view weighting tomosynthesis method and apparatus |
EP1743279A4 (fr) | 2004-04-08 | 2010-09-01 | Yeda Res & Dev | Methode de detection, de diagnostic et de pronostic du cancer du poumon faisant appel a trois points temporels |
US7352885B2 (en) * | 2004-09-30 | 2008-04-01 | General Electric Company | Method and system for multi-energy tomosynthesis |
-
2008
- 2008-02-19 FR FR0800896A patent/FR2927719B1/fr active Active
-
2009
- 2009-02-17 US US12/371,968 patent/US8798348B2/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070286470A1 (en) * | 2006-06-07 | 2007-12-13 | Sylvain Bernard | Method for the processing of radiological images in tomosynthesis for a detection of radiological signs |
Non-Patent Citations (3)
Title |
---|
DOBBINS J T ET AL: "Digital x-ray tomosynthesis: current state of the art and clinical potential", PHYSICS IN MEDICINE AND BIOLOGY, TAYLOR AND FRANCIS LTD. LONDON, GB, vol. 48, no. 19, 7 October 2003 (2003-10-07), pages R65 - R106, XP002277335, ISSN: 0031-9155 * |
PUONG SYILVIE ET AL.: "Dual-energy contrast enhanced digital breast tomosynthesis: concept, method and evaluation on phantoms", MEDICAL IMAGING 2007: PHYSICS OF MEDICAL IMAGING, EDITIED BY JIANG HSIEH AND MICHAEL J. FLYNN - PROCEEDINGS OF SPIE, vol. 6510, 2007, XP002497880 * |
WHEELER F.W. ET AL.: "Micro-Calcification Detection in Digital Tomosynthesis Mammography", MEDICAL IMAGING 2006: IMAGE PROCESSING, EDITED BY JOSEPH M. REINHARDT, JOSIEN P. W. PLUIM, PROCEEDINGS OF SPIE VOL. 6144, 13 February 2006 (2006-02-13), pages 614420-1 - 614420-12, XP002497881 * |
Also Published As
Publication number | Publication date |
---|---|
US20090208085A1 (en) | 2009-08-20 |
US8798348B2 (en) | 2014-08-05 |
FR2927719B1 (fr) | 2010-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0925556B1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
EP0389333B1 (fr) | Procédé d'acquisition de données radiologiques relatives à un corps irradié et de reconstruction de structures correspondant à ce corps | |
FR2963976A1 (fr) | Procede de traitement d'images pour la determination de zones suspectes dans une matrice tissulaire, et son utilisation pour la navigation 3d a travers la matrice tissulaire | |
FR2948481A1 (fr) | Procede d'imagerie pour la realisation d'une modelisation en triple energie, et dispositif pour la mise en oeuvre d'un tel procede | |
FR2953964A1 (fr) | Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections | |
FR2830741A1 (fr) | Procede et appareil pour etalonner un systeme d'imagerie | |
FR2954556A1 (fr) | Procede de traitement d'acquisitions de tomosynthese pour obtenir une representation du contenu d'un organe | |
EP3724646B1 (fr) | Procédé, dispositif et programme de traitement d'images de diffraction d'un matériau cristallin | |
EP1852717A1 (fr) | Procédé d'estimation de rayonnement diffusé dans un détecteur bidimensionnel | |
FR2873835A1 (fr) | Procede et dispositif d'imagerie aux rayons x avec produit de contraste permettant une visualisation amelioree | |
FR2797978A1 (fr) | Procede de recalage automatique d'images | |
FR2896607A1 (fr) | Procede et dispositif de controle de la qualite d'images de tomosynthese. | |
FR2927719A1 (fr) | Procede de traitement d'images obtenues par tomosynthese et dispositif associe | |
HUE035491T2 (en) | Transmission Image Reconstruction and Imaging Using Poisson Detector Data | |
FR2792749A1 (fr) | Procede de localisation et de representation tridimensionnelle d'elements d'interet d'un organe | |
EP3161461B1 (fr) | Procédé et système d'analyse d'un objet par diffractométrie utilisant un spectre en diffusion et un spectre en transmission | |
WO2016134903A1 (fr) | Procédé de segmentation de l'image d'un objet reconstruit par reconstruction tridimensionnelle | |
FR2993447A1 (fr) | Procede et systeme de traitement d'images pour l'affichage 3d d'un organe d'un patient | |
EP3851837B1 (fr) | Procédé d'établissement d'une fonction de réponse spectrale d'un système de mesure par rayon x ou gamma | |
EP3385757B1 (fr) | Procédé de traitement d'un spectre d'un rayonnement mesuré par un détecteur | |
EP3987312A1 (fr) | Camera compton et procede d'imagerie 3d campton | |
EP3153888A1 (fr) | Procédé de correction d'un spectre | |
EP0965946B1 (fr) | Procédé d'amelioration du rapport signal/bruit de l'image radiographique d'un objet en mouvement | |
FR2851359A1 (fr) | Procede et dispositif de calibration et de correction de niveaux de gris | |
EP2345007A1 (fr) | Procédé d'estimation de la concentration d'un traceur dans un ensemble de structures tissulaires, support d'information et dispositif correspondants |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 9 |
|
PLFP | Fee payment |
Year of fee payment: 10 |
|
PLFP | Fee payment |
Year of fee payment: 11 |
|
PLFP | Fee payment |
Year of fee payment: 14 |
|
PLFP | Fee payment |
Year of fee payment: 15 |
|
PLFP | Fee payment |
Year of fee payment: 16 |
|
PLFP | Fee payment |
Year of fee payment: 17 |