EP1485873A2 - Procede de tomosynthese a reconstruction multiplan. - Google Patents

Procede de tomosynthese a reconstruction multiplan.

Info

Publication number
EP1485873A2
EP1485873A2 EP03718844A EP03718844A EP1485873A2 EP 1485873 A2 EP1485873 A2 EP 1485873A2 EP 03718844 A EP03718844 A EP 03718844A EP 03718844 A EP03718844 A EP 03718844A EP 1485873 A2 EP1485873 A2 EP 1485873A2
Authority
EP
European Patent Office
Prior art keywords
source
planes
plane
regularization
trajectory
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.)
Withdrawn
Application number
EP03718844A
Other languages
German (de)
English (en)
Inventor
Régis Guillemaud
Pierre Bleuet
Isabelle Magnin
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Centre National de la Recherche Scientifique CNRS
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Centre National de la Recherche Scientifique CNRS
Commissariat a lEnergie Atomique CEA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Centre National de la Recherche Scientifique CNRS, Commissariat a lEnergie Atomique CEA filed Critical Centre National de la Recherche Scientifique CNRS
Publication of EP1485873A2 publication Critical patent/EP1485873A2/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T12/00—Tomographic reconstruction from projections
    • G06T12/20—Inverse problem, i.e. transformations from projection space into object space
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/436—Limited angle

Definitions

  • the invention relates to a method of tomosynthesis by lighting an object using an X-ray source.
  • the invention finds an application in the field of medical imaging, the field of non-destructive testing of objects and, more generally, in any field implementing the reconstruction of objects scrolling past a source of X-rays.
  • Tomosynthesis is a technique which, from a small number of 2D projections (2D for “two-dimensional”) distributed over a limited angular range and acquired on a digital sensor, to reconstruct the 3D volume (3D for “to three dimensions ”) of a studied object.
  • FIG. 1 A system for implementing the tomosynthesis method according to the known art is represented, symbolically, in FIG. 1.
  • a source 1 of X-rays travels along a linear path 2 in front of an object 3.
  • the X-rays 5 which pass through the object 3 are detected by a mobile detector 4 which moves in a direction opposite 6 to the direction of the source 1.
  • a computer system (not shown in the figure) collects information from the mobile detector to acquire, process and reconstruct the 3D image of the object.
  • the 3D reconstruction is relatively fast but provides relatively fuzzy reconstructions due to the low number of projections and the limited angle character of the 2D projections.
  • This patent describes taking into account an object model which makes it possible to strengthen the blood vessels in their transverse direction and to smooth them in their longitudinal direction.
  • the processing carried out is not based on any assumption concerning the acquisition geometry or the type of sampling of the reconstruction volume but on assumptions specific to the object.
  • the article "Fan-Beam reconstruction from a straight Une of source points" (B. D. Smith, IEEE TMI, vol. 12, num. L, 1993) discloses, to reconstruct the 3D volume, to exploit the linear trajectory of the source in a particular geometry consisting of a series of planes P organized in a fan as shown in FIGS. 2A and 2B.
  • a detector plane 4 comprises several parallel lines 7 of detector pixels.
  • each plane P passes through the path 2 of the source and a line 7 of detector pixels.
  • Each plane P is independent of the other planes.
  • Such geometry is only achievable in the. case of a linear trajectory of the source.
  • the reconstruction is of an analytical type. The possibility of reconstructing the 3D volume is studied theoretically with an infinite linear trajectory. An exact reconstruction formula is given under these conditions. By passing to a trajectory of the finished linear source, the quality of the reconstructions turns out to be very insufficient and a treatment subsequent to the reconstruction is then necessary.
  • the method according to the invention does not have these drawbacks.
  • the invention relates to a method of tomosynthesis by lighting an object using an X-ray source having a linear trajectory, the method comprising a step of decomposing the volume of the object into n independent planes forming a fan.
  • the method includes a step of anisotropic regularization on each plane.
  • the tomosynthesis process comprises an additional step of regularization between planes.
  • the tomosynthesis method according to the invention has the advantage of being adapted to the linear path of the source. It is an algebraic processing method particularly well suited to 3D reconstruction with little data. Brief description of the figures
  • FIG. 1 represents a block diagram of a system for implementing the tomosynthesis process
  • FIGS. 2A and 2B represent 2D planes organized in a fan for regularization in which the X-ray source moves along a linear trajectory
  • FIG. 3 represents a tomosynthesis process algorithm according to the invention.
  • the method according to the invention comprises, successively, a step E1 of data acquisition, a step E2 of decomposing the reconstruction volume into independent planes and a step E3 of reconstruction and regularization.
  • Step E1 is a step known in itself during which an X-ray source and a plane detector move linearly, parallel to each other, in opposite directions, as described above (cf. FIG. 1 ).
  • the planar detector consists of a set of lines of detector pixels. With reference to FIG. 1, the space being referenced by the direct trihedron (x, y, z), the trajectory of the source takes place in the direction x and the plane detector is parallel to the plane (x, y).
  • Step E2 performs, from the data acquired during step E1, a decomposition of the volume of the object into n planes P forming a fan.
  • the n planes P are independent of each other.
  • a plane P passes through the path of the source and through a line of pixels of the detector (cf. FIG. 2A).
  • the method according to the invention implements step E3 of reconstruction of the 3D volume.
  • the reconstruction step is associated with an anisotropic regularization step adapted to the geometry of data acquisition, that is to say to the linear trajectory of the source.
  • the anisotropic regularization according to the invention advantageously makes it possible to smooth out certain structures and to enhance others.
  • the term "enhance" should be understood as the action of emphasizing or supporting the contrast of a structure.
  • noise Three types of artifacts are present during reconstruction: noise, limited angle artifacts and source camera shake. These phenomena are distributed in the three directions of space as follows:
  • the anisotropic regularization algorithm according to the invention advantageously makes it possible to treat each of the directions x, y, z differently.
  • the anisotropic regularization on a plane P is carried out independently of the regularization on the other planes.
  • This anisotropic regularization consists of a smoothing type processing in the x direction and an enhancement type processing in the z direction.
  • the regularization on each plane P can be followed by a smoothing type regularization between planes in the direction y.
  • the regularization according to the invention is thus adapted to the anisotropy of the reconstruction artifacts linked to the particular acquisition geometry that is the linear and finite trajectory of the source. Indeed, it is according to the directions z and x that the reconstructed volume presents the least good resolution.
  • the anisotropic regularization implemented in these directions advantageously makes it possible to very significantly improve the quality of the reconstructed images.
  • it is according to the direction y that the acquired data present the best resolution. It follows that a smoothing is sufficient, in the direction y, for improve the quality of reconstructed images.
  • the algebraic regularization algorithm implements the minimization of a function J (f), where the variable f represents the object to be reconstructed.
  • - p is a projection, i.e. the natural logarithm of the ratio between the information acquired by the detector in the absence of the object and the information acquired by the detector in the presence of the object
  • - ⁇ x , ⁇ y , ⁇ z are functions, commonly called potential functions, which determine, according to the respective directions x, y, z, to what extent it is decided or not that an outline of the object to be reconstructed has been detected
  • - ⁇ x , ⁇ y, ⁇ z are weighting factors which estimate, according to the respective directions x, y, z, the difference between the projections and data evaluated a priori (data relating to the object to be reconstructed considered as formed of homogeneous zones separated by clear edges)
  • ⁇ x , ⁇ y, ⁇ z are minimum grayscale heights of the image in the respective directions x, y, z, from which a contour is accepted ( ⁇ x, ⁇ y and ⁇ z are therefore chosen greater
  • 2 is a reconstruction term calculated from the raw acquired data.
  • the other three terms of the equation are calculated reconstruction terms which can be treated independently of each other, by modifying the quantities ⁇ x , ⁇ y , ⁇ z , ⁇ ; ⁇ y, ⁇ z , ⁇ x , ⁇ y , ⁇ z -
  • the acquisition geometry that is to say the finite linear trajectory of the source, which determines the type of regularization and not the object as is the case according to
  • the parameter ⁇ z is chosen greater than ⁇ x and ⁇ y, so as to regularize more strongly in the z direction where the artefacts are strong, because of the linear trajectory of the source.
  • the parameter ⁇ z is chosen to be less than ⁇ x and pour y to take more account of discontinuities in the z direction than in the x and y directions.
  • the parameters ⁇ y and ⁇ y are chosen to be weak because the resolution along the direction y is - a priori - a good quality resolution.
  • the function ⁇ z is preferentially chosen to be concave in order to enhance the reconstructed volume more strongly along the direction z, while the functions ⁇ x and ⁇ y are preferably chosen to be convex.
  • two perpendicular linear paths are used for the source.
  • the reconstruction then takes place on the basis of a square-based pyramid corresponding to the intersection of two sets of perpendicular fan-shaped planes.
  • the method can also comprise a treatment and a correction of the radiation scattered by the object in order to further increase the resolution of the reconstructed image.
  • the lines 7 of detector pixels are substantially parallel to the path of the source.
  • a resampling of the projections is carried out before the reconstruction step.
  • the resampling step can be carried out by interpolation.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Image Processing (AREA)

Abstract

L'invention concerne un procédé de tomosynthíse par éclairage d'un objet à l'aide d'une source de rayonnement X (1) ayant une trajectoire linéaire (2), le procédé comprenant une étape de décomposition en n plans (P) formant éventail entre la trajectoire linéaire (2) et un plan détecteur (4) parallíle à la trajectoire linéaire. Le procédé comprend une étape de régularisation anisotrope sur au moins un plan (P). L'invention s'applique à l'imagerie médicale, au domaine du contrôle non destructif d'objets et, plus généralement, à tout domaine mettant en oeuvre la reconstruction d'objets en défilement devant une source de rayons X.

Description

PROCEDE DE TO OSYNTHESE A RECONSTRUCTION MULTIPLAN
Domaine technique et art antérieur
L'invention concerne un procédé de tomosynthêse par éclairage d'un objet à l'aide d'une source de rayonnement X.
L'invention trouve une application dans le domaine de l'imagerie médicale, le domaine du contrôle non destructif d'objets et, plus généralement, dans tout domaine mettant en œuvre la reconstruction d'objets en défilement devant une source de rayons X.
La tomosynthêse est une technique qui permet, à partir d'un faible nombre de projections 2D (2D pour « à deux dimensions ») réparties sur un domaine angulaire restreint et acquises sur un capteur numérique, de reconstruire le volume 3D (3D pour « à trois dimensions ») d'un objet étudié.
Le nombre de projections 2D et le débattement angulaire étant limités, il est nécessaire d'utiliser des méthodes de reconstruction associées à une opération de régularisation pour obtenir une bonne résolution des reconstructions.
Un système de mise en œuvre de procédé de tomosynthêse selon l'art connu est représenté, de façon symbolique, en figure 1. Une source 1 de rayons X se déplace selon une trajectoire linéaire 2 devant un objet 3. Les rayons X 5 qui traversent l'objet 3 sont détectés par un détecteur mobile 4 qui se déplace dans une direction opposée 6 à la direction de la source 1. Un système informatique (non représenté sur la figure) recueille l'information issue du détecteur mobile pour acquérir, traiter et reconstruire l'image 3D de 1 ' objet .
Plusieurs projections 2D de l'objet à étudier sont acquises pour différents angles de vue. A partir de ces projections, l'algorithme utilisé en tomosynthêse procède soit de manière analytique par rétroprojection de chacune des projections 2D ( cf . « Tomosynthesis : A 3D radiographie imaging technique », D. G. Grant, IEEE TMBE, vol . 19, pp 20 -28, 1972) , soit de manière algébrique éventuellement avec régularisation {cf . « Deterministic Edge Preserving Regularization in Computed Imaging », P. Charbonnier, IEEE TIP, vol . 16, Feb . 1997) .
Dans le cas où l'algorithme procède de manière analytique, la reconstruction 3D est relativement rapide mais fournit des reconstructions relativement floues dues au faible nombre de projections et au caractère à angle limité des projections 2D.
L'utilisation d'un algorithme procédant de manière algébrique avec régularisation permet d'obtenir des reconstructions de meilleure qualité ayant une résolution plus élevée. Les algorithmes qui procèdent de manière algébrique avec régularisation mettent en oeuvre des méthodes basées sur un processus itératif entre projection et reconstruction dans lesquelles il est possible d'introduire de l'information connue a priori sur l'objet. L'hypothèse sous-jacente à ces méthodes est celle d'un objet à étudier formé de zones homogènes séparées par des bords francs . Ainsi le bruit peut-il être lissé, tout en conservant les discontinuités présentes dans l'objet. La régularisation peut être anisotrope et adaptée à l'objet lui-même comme cela est divulgué dans le brevet américain US N°5 930 384 intitulé « Process for the reconstruction of a 3D image with contrast and resolution improvements and application of said process to the production of an atténuation cartography of an object » . Ce brevet décrit la prise en compte d'un modèle d'objet qui permet de renforcer les vaisseaux sanguins dans leur direction transversale et de les lisser suivant leur direction longitudinale. Le traitement effectué n'est basé sur aucune hypothèse concernant la géométrie d'acquisition ou le type d'échantillonnage du volume de reconstruction mais sur des hypothèses propres à l'objet. L'article « Fan-Beam reconstruction from a straight Une of source points » (B . D. Smith, IEEE TMI, vol . 12 , num. l , 1993) divulgue, pour reconstruire le volume 3D, d'exploiter la trajectoire linéaire de la source dans une géométrie particulière constituée d'une série de plans P organisés en éventail comme cela est représenté sur les figures 2A et 2B. Un plan détecteur 4 comprend plusieurs lignes parallèles 7 de pixels détecteurs . Par organisation en éventail des plans P, il faut entendre que chaque plan P passe par la trajectoire 2 de la source et une ligne 7 de pixels détecteurs. Chaque plan P est indépendant des autres plans. Une telle géométrie n'est réalisable que dans le. cas d'une trajectoire linéaire de la source. Dans cet article, la reconstruction est de type analytique. La possibilité de reconstruire le volume 3D est étudiée de façon théorique avec une trajectoire linéaire infinie. Une formule exacte de reconstruction est donnée dans ces conditions. En passant à une trajectoire de la source linéaire finie, la qualité des reconstructions s'avère très insuffisante et un traitement postérieur à la reconstruction est alors nécessaire.
De façon générale, selon l'art connu, en tomosynthêse linéaire, on se trouve en présence d'un problème de déformation verticale des objets suivant l'axe vertical z qui se traduit également par un flou dans la direction de translation de la source. Par ailleurs chacune des directions (x,y,z) est soumise à un bruit d'acquisition.
Le procédé selon l'invention ne présente pas ces inconvénients.
Exposé de l'invention
En effet, l'invention concerne un procédé de tomosynthêse par éclairage d'un objet à l'aide d'une source de rayonnement X ayant une trajectoire linéaire, le procédé comprenant une étape de décomposition du volume de l'objet en n plans indépendants formant éventail. Le procédé comprend une étape de régularisation anisotrope sur chaque plan.
Selon un perfectionnement de l'invention, le procédé de tomosynthêse comprend une étape supplémentaire de régularisation entre plans.
Le procédé de tomosynthêse selon l'invention présente l'avantage d'être adapté à la trajectoire linéaire de la source. C'est un procédé de traitement algébrique particulièrement bien adapté à une reconstruction 3D avec peu de données. Brève description des figures
D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention fait en référence aux figures jointes parmi lesquelles ; la figure 1 représente un schéma de principe de système de mise en œuvre de procédé de tomosynthêse ; les figures 2A et 2B représentent des plans 2D organisés en éventail pour une régularisation dans laquelle la source de rayons X se déplace selon une trajectoire linéaire ; la figure 3 représente un algorithme de procédé de tomosynthêse selon l'invention.
Description détaillée de modes de mise en œuyre de 1' invention
Comme représenté en figure 3, le procédé selon l'invention comprend, successivement, une étape El d'acquisition de données, une étape E2 de décomposition du volume de reconstruction en plans indépendants et une étape E3 de reconstruction et de régularisation.
L'étape El est une étape connue en soi durant laquelle une source de rayons X et un détecteur plan se déplacent linéairement, parallèlement l'un à l'autre, dans des directions opposées, comme cela a été décrit précédemment (cf. figure 1). Le détecteur plan est constitué d'un ensemble de lignes de pixels détecteurs. En référence à la figure 1, l'espace étant référencé par le trièdre direct (x, y, z) , la trajectoire de la source s'effectue selon la direction x et le plan détecteur est parallèle au plan (x, y) .
L'étape E2 effectue, à partir des données acquises lors de l'étape El, une décomposition du volume de l'objet en n plans P formant éventail. Les n plans P sont indépendants les uns des autres. Un plan P passe par la trajectoire de la source et par une ligne de pixels du détecteur (cf. figure 2A) .
Une fois effectuée la décomposition en plans, le procédé selon l'invention met en œuvre l'étape E3 de reconstruction du volume 3D. L'étape de reconstruction est associée à une étape de régularisation anisotrope adaptée à la géométrie d'acquisition des données, c'est-à-dire à la trajectoire linéaire de la source. La régularisation anisotrope selon l'invention permet avantageusement de lisser certaines structures et d'en rehausser d'autres. Le terme « rehausser » doit s'entendre comme l'action de souligner ou d'appuyer le contraste d'une structure.
Trois types d'artefacts sont présents lors de la reconstruction : le bruit, les artefacts d'angle limité et le flou de bougé de la source. Ces phénomènes sont répartis dans les trois directions de l'espace comme suit :
1. Suivant la direction x, il y a du bruit et un flou dû au déplacement de la source.
2. Suivant la direction y, il y a uniquement du bruit (la résolution est alors seulement celle du détecteur) .
3. Suivant la direction z, il y a du bruit et des artefacts d'angle limité à cause de la trajectoire finie de la source. L'algorithme de régularisation anisotrope selon l'invention permet avantageusement de traiter différemment chacune des directions x, y, z.
La reconstruction à l'aide des plans P organisés en éventail permet de réduire considérablement le temps de calcul puisque les traitements se font dans deux dimensions et qu'il est possible d'effectuer le traitement des différents plans en parallèle. La régularisation anisotrope sur un plan P est effectuée indépendamment de la régularisation sur les autres plans. Cette régularisation anisotrope consiste en un traitement de type lissage selon la direction x et un traitement de type rehaussement selon la direction z.
Selon un perfectionnement de l'invention, la régularisation sur chaque plan P peut être suivie d'une régularisation de type lissage entre plans selon la direction y. La régularisation selon l'invention est ainsi adaptée à l' anisotropie des artefacts de reconstruction liés à la géométrie d'acquisition particulière qu'est la trajectoire linéaire et finie de la source. En effet, c'est selon les directions z et x que le volume reconstruit présente la moins bonne résolution. La régularisation anisotrope mise en œuvre selon ces directions permet avantageusement d'améliorer très sensiblement la qualité des images reconstruites. Par contre, c'est selon la direction y que les données acquises présentent la meilleure résolution. Il s'en suit qu'un lissage suffit, selon la direction y, pour améliorer la qualité des images reconstruites.
L'algorithme algébrique de régularisation selon l'invention met en œuvre la minimisation d'une fonction J(f) , où la variable f représente l'objet à reconstruire .
A titre d'exemple non limitatif, la fonction J(f) peut s'écrire :
J (f ) = |p- f + λ2 x (équation ( 1) ) où
- fi,j,k représente l'objet reconstruit aux points de coordonnées i, j, k, - R est une matrice de projection appliqué à l'objet,
- p est une projection, c'est-à-dire le logarithme Népérien du rapport entre 1 * information acquise par le détecteur en l'absence de l'objet et l'information acquise par le détecteur en présence de l'objet, - φx, φy, φz sont des fonctions, communément appelée fonctions potentiel, qui déterminent, selon les directions respectives x, y, z, dans quelle mesure il est décidé ou non qu'un contour de l'objet à reconstruire a été détecté, - λx, λy, λz sont des facteurs de pondération qui estiment, selon les directions respectives x, y, z, la différence entre les projections et des données évaluées a priori (données relatives à l'objet à reconstruire considéré comme formé de zones homogènes séparées par des bords francs) , δx, δy, δz sont des hauteurs minimales de niveau de gris de l'image selon les directions respectives x, y, z, à partir desquelles un contour est accepté (δx, δy et δz sont donc choisis supérieurs à un niveau de bruit) .
Dans l'équation (1), le terme | |p - Rf | | 2 est un terme de reconstruction calculé à partir des données acquises brutes. Les trois autres termes de l'équation sont des termes de reconstruction calculés qui peuvent être traités indépendamment les uns des autres, en modifiant les grandeurs δx, δy, δz, λχ; λy, λz, φx, φy, φz-
Avantageusement, selon l'invention, c'est la géométrie d'acquisition, c'est-à-dire la trajectoire linéaire finie de la source, qui détermine le type de régularisation et non l'objet comme c'est le cas selon
1 ' art connu.
L'adaptation de la régularisation anisotrope à la trajectoire linéaire de la source va maintenant être décrite avec plus de précisions en référence aux grandeurs δx, δy, δz, λx, λy, λz, φx, φy/ φz mentionnées ci-dessus .
Le paramètre λz est choisi supérieur à λx et λy, de manière à régulariser plus fortement suivant la direction z où les artefacts sont forts, à cause de la trajectoire linéaire de la source. Le paramètre δz est choisi inférieur à δx et δy pour prendre en compte plus fortement les discontinuités dans la direction z que dans les directions x et y. Les paramètres λy et δy sont choisis faibles car la résolution suivant la direction y est - a priori - une résolution de bonne qualité. Enfin, la fonction φz est choisie préférentiellement concave pour rehausser le volume reconstruit plus fortement selon la direction z, alors que les fonctions φx et φy sont préférentiellement choisies convexes. Selon un deuxième mode de réalisation de l'invention, deux trajectoires linéaires perpendiculaires sont utilisées pour la source. La reconstruction s'opère alors sur la base d'une pyramide à base carrée correspondant à l'intersection de deux ensembles de plans en éventail perpendiculaires. L'avantage d'une utilisation de deux trajectoires perpendiculaires est de donner plus d'information sur l'objet et, partant, de mieux reconstruire celui-ci.
Quel que soit le mode de réalisation de l'invention, par ailleurs, le procédé peut également comprendre un traitement et une correction du rayonnement diffusé par l'objet afin d'accroître encore la résolution de l'image reconstruite.
Selon le mode de réalisation préférentiel de l'invention décrit ci-dessus, les lignes 7 de pixels détecteurs sont sensiblement parallèles à la trajectoire de la source. Dans le cas où le parallélisme entre la trajectoire de la source et les lignes de pixels détecteurs n'est pas correctement réalisé, par exemple pour des raisons mécaniques, il est procédé à une étape de rééchantillonnage des projections avant l'étape de reconstruction. L'étape de rééchantillonnage peut être réalisé par interpolation.

Claims

REVENDICATIONS
1. Procédé de tomosynthêse par éclairage d'un objet à l'aide d'une source de rayonnement X (1) ayant une trajectoire linéaire (2) , le procédé comprenant une étape de décomposition en n plans (P) formant éventail entre la trajectoire linéaire (2) et un plan détecteur (4) parallèle à la trajectoire linéaire, caractérisé en ce qu'il comprend une étape de régularisation anisotrope sur au moins un plan (P) .
2. Procédé selon la revendication 1, caractérisé en ce que l'étape de régularisation anisotrope sur un plan (P) consiste en un lissage selon la direction (x) de la trajectoire linéaire de la source (2) et en un réhaussement selon une direction perpendiculaire (z) à la trajectoire linéaire de la source (2) et au plan détecteur (4) .
3. Procédé selon la revendication 2, caractérisé en ce qu'il comprend une étape supplémentaire de lissage entre plans (P) selon une direction (y) perpendiculaire à la trajectoire linéaire de la source (2) et parallèle au plan détecteur (4) .
4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la régularisation anisotrope comprend la minimisation d'une fonction J(f) telle que : j(f) -ωr- + *.x où
- fι,j,k représente l'objet reconstruit aux points de coordonnées i, j, k,
- R est une matrice de projection appliqué à l'objet, - p est une projection,
- φx/ φy φz sont des fonctions qui déterminent, selon des directions respectives x, y, z, dans quelle mesure il est décidé ou non qu'un contour de l'objet à reconstruire a été détecté, la direction x étant la direction de la trajectoire linéaire de la source et les directions y et z étant des directions perpendiculaires entre elles qui définissent un plan perpendiculaire au plan détecteur,
- λx, λy, λz sont des facteurs de pondération qui estiment, selon les directions respectives x, y, z, la différence entre les projections et des données évaluées a priori, δx, δy, δz sont des hauteurs minimales selon les directions respectives x, y, z, à partir desquelles un contour est accepté.
5. Procédé de tomosynthêse par éclairage d'un objet à l'aide d'une source de rayonnement X (1), caractérisé en ce qu'il comprend un procédé selon l'une quelconque des revendications 1 à 4 dans lequel la trajectoire de la source est effectuée selon une première direction et un procédé selon l'une quelconque des revendications 1 à 4 dans lequel la trajectoire de la source est effectuée selon une deuxième direction perpendiculaire à la première direction, la régularisation anisotrope s 'opérant sur la base d'une pyramide de plans P correspondant à l'intersection d'un premier ensemble de plans P formant éventail entre la trajectoire linéaire de la source selon la première direction et le plan détecteur (4) et un second ensemble de plans P formant éventail entre la trajectoire linéaire de la source selon la deuxième direction perpendiculaire à la première direction et le plan détecteur (4) .
6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comprend un traitement et une correction d'un rayonnement diffusé par l'objet.
EP03718844A 2002-02-08 2003-02-06 Procede de tomosynthese a reconstruction multiplan. Withdrawn EP1485873A2 (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR0201558A FR2835949B1 (fr) 2002-02-08 2002-02-08 Procede de tomosynthese a reconstruction multiplan
FR0201558 2002-02-08
PCT/FR2003/000374 WO2003067525A2 (fr) 2002-02-08 2003-02-06 Procede de tomosynthese a reconstruction multiplan

Publications (1)

Publication Number Publication Date
EP1485873A2 true EP1485873A2 (fr) 2004-12-15

Family

ID=27620021

Family Applications (1)

Application Number Title Priority Date Filing Date
EP03718844A Withdrawn EP1485873A2 (fr) 2002-02-08 2003-02-06 Procede de tomosynthese a reconstruction multiplan.

Country Status (4)

Country Link
US (1) US8670601B2 (fr)
EP (1) EP1485873A2 (fr)
FR (1) FR2835949B1 (fr)
WO (1) WO2003067525A2 (fr)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2849250B1 (fr) * 2002-12-23 2005-05-13 Commissariat Energie Atomique Procede de reconstruction d'une image radiographique par combinaison de vignettes se recouvrant
WO2005055803A2 (fr) * 2003-12-03 2005-06-23 The General Hospital Corporation D/B/A Massachusetts General Hospital Systeme et procede pour la reconstitution de faisceau conique multisegment en imagerie par tomosynthese
US8538099B2 (en) * 2005-03-23 2013-09-17 General Electric Company Method and system for controlling image reconstruction
US8290225B2 (en) * 2005-12-14 2012-10-16 Koninklijke Philips Electronics N.V. Method and device for relating medical 3D data image viewing planes to each other
DE102011076929A1 (de) * 2011-06-03 2012-12-06 Siemens Ag Verfahren und Vorrichtung zur Darstellung von Volumendaten für eine Untersuchung von Dichteeigenschaften
FR3047339B1 (fr) * 2016-02-01 2018-04-06 Safran Procede de controle non-destructif par redressement
DE102019204765B3 (de) 2019-04-03 2020-06-18 Siemens Healthcare Gmbh Verfahren zur Ermittlung eines dreidimensionalen Tomosynthesedatensatzes, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4054402B2 (ja) * 1997-04-25 2008-02-27 株式会社東芝 Ｘ線断層撮影装置
US5214686A (en) * 1991-12-13 1993-05-25 Wake Forest University Three-dimensional panoramic dental radiography method and apparatus which avoids the subject's spine
DE19509007C2 (de) * 1995-03-13 2001-07-05 Siemens Ag C-Bogen-Röntgendiagnostikgerät zum Erstellen von Schichtaufnahmen
FR2736455B1 (fr) * 1995-07-03 1997-08-08 Commissariat Energie Atomique Procede de reconstruction d'une image 3d avec amelioration du contraste et de la resolution et application de ce procede a la realisation d'une cartographie d'attenuation d'un objet
JP3373720B2 (ja) * 1996-03-25 2003-02-04 株式会社日立メディコ Ｘ線断層撮影装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO03067525A2 *

Also Published As

Publication number Publication date
FR2835949A1 (fr) 2003-08-15
WO2003067525A2 (fr) 2003-08-14
WO2003067525A3 (fr) 2004-03-25
US20050078862A1 (en) 2005-04-14
US8670601B2 (en) 2014-03-11
FR2835949B1 (fr) 2004-07-09

Similar Documents

Publication Publication Date Title
EP0379399B1 (fr) Procédé de calcul et d'exploitation de l'image en projection conique, par exemple au sens des rayons x, d'un objet tridimensionnel echantillonné, et procédé de reconstruction tridimensionnelle d'un objet étudié utilisant ce procédé de calcul
Kang et al. Deep convolutional framelet denosing for low-dose CT via wavelet residual network
FR2800189A1 (fr) Procede de reconstruction multi-resolutions d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique
EP0925556A2 (fr) Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique
EP0777893A1 (fr) Procede de reconstruction d'une image 3d avec amelioration du contraste et de la resolution
WO1998053426A1 (fr) Dispositif electronique de traitement d'images pour la detection de variations dimensionnelles
FR2798804A1 (fr) Procede et appareil pour reconstruire des donnees d'image acquises par un systeme d'imagerie a rayons x par tomosynthese
EP0840252A1 (fr) Procédé de traitement d'image numérique pour l'extraction automatique d'objets en forme de rubans
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
FR2544939A1 (fr) Systeme et procede de reprojection d'images acquises par retroprojection en tomodensitographie
FR2848007A1 (fr) Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal
FR2885717A1 (fr) Procede d'estimation de correction d'une image et dispositif d'angiographie le mettant en oeuvre
Yang et al. Self-supervised cryo-electron tomography volumetric image restoration from single noisy volume with sparsity constraint
EP1485873A2 (fr) Procede de tomosynthese a reconstruction multiplan.
FR2848806A1 (fr) Procede de calibrage d'un appareil d'imagerie radiologique necessitant un nombre limite d'acquisitions
Kong et al. D-star: Diffusion-based sparse tomographic angular recovery for isotropic-resolution photoacoustic imaging
US8184887B2 (en) System and method for image reconstruction
WO2005001775A1 (fr) Procede de reconstruction d'une image tomographique par une methode analytique comprenant une modelisation amelioree du mouvement de l'objet
Arun et al. Despeckling of OCT images using DT-CWT based fusion technique
EP1296286A1 (fr) Procédé permettant de suivre la déformation sur trois dimensions d'un organe déformable
Sunnegårdh Iterative filtered backprojection methods for helical cone-beam CT
EP3411851A1 (fr) Procédé de contrôle non-destructif par redressement
CN119090758A (zh) 一种基于自监督学习的冷冻电子断层扫描图像去噪方法
FR2823345A1 (fr) Procede d'amelioration de la qualite d'une image radiographique tridimensionnelle d'un objet et dispositif radiographique correspondant
FR3032543A1 (fr) Procede de determination d'un axe de rotation d'un objet en tomographie et procede de caracterisation par tomographie

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20040719

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PT SE SI SK TR

RIN1 Information on inventor provided before grant (corrected)

Inventor name: GUILLEMAUD, REGIS

Inventor name: BLEUET, PIERRE

Inventor name: MAGNIN, ISABELLE

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: COMMISSARIAT A L'ENERGIE ATOMIQUE

Owner name: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE (CNRS

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: COMMISSARIAT A L'ENERGIE ATOMIQUE ET AUX ENERGIES

Owner name: CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE (CNRS

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

INTG Intention to grant announced

Effective date: 20160901

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20170112