FR2632427A1 - Procede pour ameliorer sensiblement le temps de reconstruction d'une image dans une tomographie informatisee a angle limite - Google Patents
Procede pour ameliorer sensiblement le temps de reconstruction d'une image dans une tomographie informatisee a angle limite Download PDFInfo
- Publication number
- FR2632427A1 FR2632427A1 FR8906840A FR8906840A FR2632427A1 FR 2632427 A1 FR2632427 A1 FR 2632427A1 FR 8906840 A FR8906840 A FR 8906840A FR 8906840 A FR8906840 A FR 8906840A FR 2632427 A1 FR2632427 A1 FR 2632427A1
- Authority
- FR
- France
- Prior art keywords
- image
- projections
- measured
- missing
- reconstructed
- 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
Links
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/424—Iterative
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Image Processing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Generation (AREA)
Abstract
Pour obtenir une réduction importante du temps d'ordinateur dans la tomographie informatisée à angle limité, l'image reconstruite est décomposée en deux images partielles 22, 24 reconstituées à partir d'une donnée de balayage disponible et d'une donnée de balayage manquante. Comme la première est inchangée pendant le processus de reconstruction itérative des images, elle est reconstituée à partir de projections mesurées 21 une fois seulement au commencement des itérations et utilisée de façon répétitive. L'image composite formée par l'addition des deux images partielles est corrigée par l'information a priori sur l'objet. L'image partielle perfectionnée due aux projections manquantes est calculée et de nouveau additionnée à l'image partielle due à des projections mesurées et corrigée. Ces étapes sont répétées de manière itérative jusqu'à ce qu'un test de convergence donne une image finale. Un dispositif TI radiographique à angle limité utilise la version à espace de projection de l'algorithme de transformée itérative. Application à la reconstruction des images.
Description
La présente invention concerne l'imagerie à angle limité et plus
spécialement la réduction du temps de calcul pour reconstruire une image en utilisant des algorithmes itératifs. Dans certains cas de la tomographie informatisée
(TI) radiographique industrielle, on ne dispose pas de don-
nées de projection complètes relatives à l'objet. Les don-
nées de projection dans une certaine plage angulaire peuvent être manquantes pour des raisons telles qu'une atténuation trop importante ou une inaccessibilité physique. De tels cas sont décrits dans le brevet des Etats-Unis d'Amérique 4 506 327. Un procédé itératif pour reconstruire l'objet dans de telles circonstances est indiqué dans ce brevet ainsi que
dans d'autres demandes de brevet en attente de la demande-
resse et dans des articles publiés par la demanderesse et par V. PerezMendez; le mode opératoire est indiqué en figures 1 et 2. L'algorithme d'itération introduit la donnée de balayage manquante en utilisant une information a priori
sur l'objet. L'image reconstruite est formée dans un mouve-
ment de va et vient entre l'espace-objet et l'espace de pro-
jection, étant corrigée dans l' espace-objet par - 2 -
- 2 - 2632427
l'information a priori relative à cet objet, et dans
l'espace de projection par la donnée de projection mesurée.
Les résultats indiqués dans lesdits articles montrent que l'algorithme d'itération apporte de fait une amélioration importante à l'imagerie à angle limité. En général, il faut environ 5 à 10 itérations pour provoquer la convergence de l'image reconstruite. Dans chaque itération, la majeure partie du temps d'ordinateur est consacrée à la rétroprojection filtrée et aux opérations de projection. Comme l'opération de projection ne doit être effectuée que dans les vues manquantes, le temps d'ordinateur pour l'opération de projection sera beaucoup plus petit pour l'opération de rétroprojection filtrée si le nombre des vues manquantes est une petite fraction du nombre total des vues. Cette condition est généralement remplie
dans la tomographie informatisée radiographique indus-
trielle. Ainsi, le temps total T de l'ordinateur pour reconstruire l'image en utilisant l'algorithme itératif de
la figure 2 est en gros égal au nombre des itérations multi-
plié par n fois le temps d'ordinateur t pour une opération de rétroprojection filtrée: T=nt. Comme la rétroprojection filtrée est l'algorithme pour la reconstruction des images dans une TI radiographique à angle complet, cette équation montre que le temps d'ordinateur pour la reconstruction à
angle limité est environ 5 à 10 fois celui de la reconstruc-
tion à angle complet.
La présente invention a pour objet la mise au
point d'un procédé qui permet d'obtenir une réduction impor-
tante des calculs nécessaires à la reconstruction des images
d'un objet dans une TI à angle limité.
La présente invention a pour autre objet de réduire d'un facteur élevé le temps d'ordinateur pour la TI et la TI radiographique à angle limité lorsque le nombre de vues manquantes est faible par rapport au nombre total des
vues.
-3 - 2632427
La présente invention admet que le mode opératoire
de reconstruction est tel que l'image intermédiaire recons-
truite en utilisant l'algorithme itératif de l'art antérieur
peut être considérée comme composée de deux images par-
tielles, l'une reconstruite à partir de la donnée de projec-
tion mesurée et l'autre à partir des projections manquantes.
Comme la première image reste inchangée pendant les itéra-
tions, elle ne peut être reconstruite qu'une fois au commen-
cement des itérations et utilisée à maintes reprises pendant
celles-ci.
Un procédé de formation d'une image d'un objet par le mode opératoire perfectionné de reconstruction à angle limité comporte les étapes suivantes. L'objet est soumis à
un balayage avec un agent d'imagerie dans une plage angu-
laire limitée et des projections mesurées sont générées à des angles de vision disponibles. Une image partielle de la donnée mesurée est reconstruite à partir des projections
mesurées; cela peut être effectué par rétroprojection fil-
trée, laquelle est une opération linéaire. Une image compo-
site de l'objet est produite en additionnant l'image par-
tielle de la donnée mesurée et l'image partielle de la don-
née calculée qui est reconstruite à partir des projections
manquantes; ces dernières étant initialement mises à zéro.
L'image composite est corrigée en utilisant une information a priori sur l'objet en: 1) mettant à zéro tous les pixels d'image à l'extérieur de la limite connue de l'objet, 2) remettant à une limite supérieure préchoisie tous les pixels ayant une densité dépassant la limite supérieure, et
3) remettant à zéro tout pixel ayant une densité négative.
Les projections manquantes aux angles de vision manquants sont calculées à partir de l'image composite corrigée, intermédiaire. L'étape suivante consiste à reconstruire une
image partielle améliorée de la donnée calculée due aux pro-
jections manquantes venant d'être calculées; cette opéra-
tion peut être une opération de rétroprojection filtrée. Les quatre dernières étapes sont répétées de manière itérative, en commençant par l'addition des deux images partielles, de manière à améliorer progressivement l'image composite en utilisant la même image partielle de la donnée mesurée et l'image partielle successivement améliorée des données cal- culées jusqu'à l'obtention d'une qualité suffisante pour l'image. Avec la diminution du rapport entre le nombre de vues manquantes et le nombre total de vues, et en supposant que le temps d'ordinateur pour la rétroprojection filtrée est approximativement le même que pour la projection, le
gain de temps en matière d'ordinateur augmente.
Une autre caractéristique de l'invention, en par-
ticulier de la TI radiographique à angle limité est que l'information a priori sur l'objet peut comporter la limite extérieure de l'objet dont on se rapproche par l'enveloppe
convexe de l'objet construite à partir des projections mesu-
rées.
La suite de la description se réfère aux figures
annexées qui représentent respectivement: Figure 1, une représentation schématique de la
version à espace de projection d'un algorithme de transfor-
mation itérative de l'art antérieur pour reconstruire un objet à partir d'une donnée à angle limité; Figure 2, un organigramme du mode opératoire de reconstruction des images à angle limité de l'art antérieur de la figure 1; Figure 3, une représentation schématique du mode opératoire itératif modifié de la présente invention qui décompose l'image reconstruite en deux images partielles reconstituées à partir des données de balayage disponibles et manquantes, respectivement; Figure 4, un organigramme du mode opératoire plus efficace pour la reconstruction des images à angle limité de la figure 3; Figure 5, un schéma d'un dispositif de tomographie
_ 5 _ 2632427
informatisée radiographique industrielle; et Figure 6, la reconstruction de l'enveloppe convexe
d'un objet à partir de la donnée de projection mesurée.
Lorsqu'un objet est balayé avec des rayons X dans une plage angulaire limitée, les signaux des détecteurs des rayons X représentent les projections de l'objet dans la plage angulaire limitée. En figure 1, on a représenté les projections mesurées 10 à cinq angles de vision disponibles, et les projections manquantes 11 à trois angles de vision
manquants sont remplies par un algorithme d'itération utili-
sant la limite extérieure de l'objet 12, qui est connue, et une autre information a priori sur l'objet. Les projections manquantes sont mises à l'origine à zéro et l'objet est
reconstruit par une opération de rétroprojection filtrée RF.
Après correction de la densité de l'objet, qui correspond au
niveau du gris de l'image reconstruite, les projections man-
quantes sont calculées à partir de-l'objet reconstruit par
une opération de projection P. Le mode opératoire de recons-
truction d'images de l'art antérieur est indiqué plus com-
plètement dans l'organigramme de la figure 2. En liaison avec les étapes 13-15, les projections mesurées dans une plage angulaire limitée sont acquises (Etape 13). Un ensemble complet de projections de l'objet est constitué de trois projections mesurées plus les projections manquantes à
des angles de vision inaccessibles qui sont mises initiale-
ment à zéro (Etape 14). La densité de l'objet est recons-
truite par rétroprojection filtrée (Etape 15).-L'estimation
initiale de la densité de l'objet, étapes 16 et 17, est cor-
rigée par l'information a priori sur l'objet, c'est-à-dire
par l'étendue et l'emplacement de l'objet, la limite supé-
rieure connue de la densité de l'objet et le fait qu'il n'y-
a pas de densité négative (Etape 17). L'image de l'objet est corrigée, pixel par pixel, en remettant à zéro les pixels
situés à l'extérieur de l'étendue connue dé l'objet, remet-
tant à la limite supérieure les pixels ayant une densité
- 6 - 2632427
dépassant cette limite supérieure et en remettant à zéro les pixels ayant une densité négative (Etape 16). Après l'exécution d'un test de convergence de la donnée, étape 18, les projections manquantes de la densité intermédiaire de l'objet dans les vues manquantes sont calculées (Etape 19). Une seconde itération commence et la premiere estimation des projections manquantes est maintenant fournie
ainsi que les projections mesurées. Une opération de rétro-
projection filtrée est effectuée sur l'ensemble des projec-
tions de manière à reconstruire l'objet. La densité de l'objet est corrigée par l'information a priori, la seconde estimation des projections manquantes est calculée, etc.. En général, l'image reconstruite de l'objet converge après environ 5 à 10 itérations et une densité finale de l'objet
reconstruit ou image reconstruite, étape 20, est sortie.
La transformée itérative modifiée de la figure 3 permet de réduire sensiblement le temps de reconstruction
des images dans une TI à angle limité. Comme la rétroprojec-
tion filtrée est une opération linéaire, les images intermé-
diaires reconstruites en 15 en utilisant les itérations de la figure 2 peuvent être considérées comme constituées de deux images partielles: l'une reconstruite à partir de la
donnée de projection mesurée et l'autre à partir des projec-
tions manquantes. Comme la première reste inchangée pendant les itérations, il faut seulement la reconstruire une fois au commencement des itérations et l'utiliser maintes et maintes fois pendant les itérations. Cette modification conduit à des gains importants en temps d'ordinateur si le nombre des vues contenant des données est élevé par rapport
au nombre des vues manquantes, au détriment du stockage sup-
plémentaire pour l'image partielle de la donnée de projec-
tion mesurée.
Après acquisition des projections mesurées 21, une image partielle 22 de la donnée mesurée, due aux projections
mesurées, est reconstruite par une opération de rétroprojec-
_ 7 _ 2632427
tion filtrée. C'est la seule fois o cela est exécuté; ce mode opératoire ne s'effectue pas dans la boucle d'itération. Les projections manquantes 23 sont initialement
mises à zéro et ensuite calculées à partir de l'image compo-
site corrigée. Une image partielle 24 de donnée calculée due
aux projections manquantes, est reconstruite par rétropro-
jection filtrée. Les deux images partielles 22 et 24 sont additionnées pour produire une image composite -intermédiaire qui est corrigée de la même manière que précédemment, en utilisant une information a priori sur l'objet. Il s'agit
d'une opération non linéaire'et des corrections sont appor-
tées à l'image composite au lieu des images partielles. Les projections manquantes aux angles de vision manquants sont
calculées en ordinateur à partir de l'image composite corri-
gée intermédiaire en utilisant une opération de projection.
Une image partielle perfectionnée 24 de la donnée calculée est calculée une seconde fois par rétroprojection filtrée et
additionnée à l'image partielle 22 de la donnée mesurée sto-
ckée. L'image composite intermédiaire 25 est progressivement améliorée au fur et à mesure de l'itération jusqu'à ce que
l'on obtienne une qualité suffisante pour l'image.
- La figure 4 est un organigramme de la version spa-
tiale de projection de la transformée itérative modifiée et perfectionnée. Les étapes 26 et 27 sont celles o l'objet est soumis à un balayage avec des rayons X à des angles de vision accessibles et o des projections mesurées-dans la plage angulaire limitée sont acquises. L'image partielle de
la donnée mesurée due aux projections mesurées est recons-
truite par rétroprojection filtrée. L'image composite, étape
28, est formée en combinant l'image partielle due aux pro-
jections mesurées et une image partielle due aux projections manquantes. L'image composite, étapes 29 et 30, est corrigée par l'information a priori (30) indiquée sur l'objet en: 1) remettant à zéro les pixels se trouvant -à l'extérieur de l'étendue connue et de la limite extérieure de l'objet, 2)
- 8 -2632427
remettant à la limite supérieure les pixels ayant une den-
sité dépassant la limite supérieure, et 3) remettant à zéro
les pixels ayant une densité négative.
Si un test de convergence n'est pas satisfait à l'étape 31, les étapes 32 et 33 suivantes consistent à cal-
culer les projections manquantes de l'image composite inter-
médiaire et la densité de l'objet dans les vues manquantes (Etape 32), et à reconstruire une autre image partielle de donnée calculée en utilisant les projections manquantes venant d'être calculées et une opération de rétroprojection
filtrée (Etape 33). L'image composite (28) est progressive-
ment améliorée alors que les étapes 29 et 31-33 sont répé-
tées de manière itérative en utilisant la même image par-
tielle dela donnée mesurée et l'image partielle successive-
ment améliorée de la donnée calculée due aux projections manquantes. Lors de la satisfaction du test de convergence en 31, une image composite reconstruite finale est sortie en 34.
La figure 5 représente un dispositif TI radiogra-
phique à angle limité et une réalisation pratique de la pré-
sente invention. Une source 31 de rayons X balaye l'objet 36
dans une plage angulaire limitée; les photons radiogra-
phiques qui pénètrent l'objet sont détectés par un réseau de détecteurs 37 de rayons X. Le balayage peut être effectué soit dans une géométrie à faisceaux parallèles, soit dans
une géométrie à faisceaux en éventail, cette dernière géomé-
trie étant illustrée. Les signaux détectés, qui représentent les projections de l'objet dans la plage angulaire limitée,
sont introduits dans un ordinateur 38. Les données radiogra-
phiques mesurées à angle limité sont entrées dans les algo-
rithmes 39 de reconstruction des images à angle limité. A l'autre entrée 40 de l'ordinateur l'information a priori sur l'objet est présentée, c'est-à-dire la limite extérieure de
l'objet, la limite supérieure et la limite inférieure (géné-
ralement zéro) de sa densité. Les limites supérieure et
9 - 2632427
inférieure de la densité de l'objet peuvent être estimées à partir de sa composition chimique. On peut obtenir la.limite extérieure par sondage, soit mécanique soit optique, ou on peut s'en rapprocher avec l'enveloppe convexe de l'objet construit à partir de la donnée radiographique, comme cela est représenté en 41 en utilisant une technique décrite dans la demande de brevet de la demanderesse n 032 804 déposée le 1er avril 1987. Une autre référence est le document dit TIS Technical Report Numéro 88CRD006, février 1988, de la société dite General Electric: "The Use of Convex Hulls in Limited-Angle Computerized Tomography".(Utilisation des enveloppes convexes dans la tomographie informatisée à angle limité).
La figure 6 représente les principes de construc-
tion de l'enveloppe convexe à partir des projections mesu-
rées. Si la rétroprojection d'une projection mesurée 42 est limitée par les lignes 43 en tirets, on peut montrer que le support de l'objet est complètement contenu dans la partie non zéro de la rétroprojection. Lorsqu'une zone D est construite par l'intersection de toutes les rétroprojections (seules 3 sont représentées), l'objet 36 est entièrement contenu dans la zone D. La zone D est un polygone contenant l'objet, et avec l'augmentation du nombre des angles de vision, la zone D se rapproche de la limite de l'objet si celle-ci -est convexe. Si la limite n'est pas convexe, la
zone D se rapproche dé l'enveloppe convexe de la limite.
En utilisant l' information a priori fournie et la
donnée de balayage à angle limité, l'objet et son image com-
posite sont reconstruits au moyen de l'algorithme de recons-
tructi6n à angle limité de la figure 4. Lors de la satisfac-
tion du test de convergence, l'image composite reconstruite
de l'objet est affichée en 44, figure 5.
La présente invention réduit d'un facteur élevé le
temps d'ordinateur pour la TI à angle limité et la TI radio-
graphique lorsque le nombre des vues manquantes est faible
- 10 - 2632427
par rapport au nombre total de vues. Si l'on suppose que le temps d'ordinateur pour la rétroprojection filtrée est à peu près le même que pour la projection, le temps total d'ordinateur (en unités arbitraires) pour la reconstruction des images en utilisant l'algorithme de transformation modi- fié de la figure 4 est indiqué dans le tableau suivant en fonction du nombre n des itérations et du rapport r entre le nombre de vues manquantes et le nombre total de vues. Le temps d'ordinateur correspondant en utilisant l'algorithme de l'art antérieur de la figure 2 est indiqué à titre de comparaison. On peut voir qu'avec le nouvel algorithme la TI à angle limité devient beaucoup plus souhaitable en matière
de temps d'ordinateur.
TABLEAU 1
n/r 0,05 0,10 0,20
1,5(11) 2,0(11) 3,0(11)
2,0(21) 3,0(21) 5,0(21)
3,0(41) 5,0(41) 9,0(41)
Le rapport technique TIS de la société General
Electric ayant pour numéro 87CRD212, novembre 1987, "Compa-
rison of Three Methods in Limited-Angle Image Reconstruc-
tions" (comparaison de trois méthodes dans la reconstruction des images à angle limité), examine trois algorithmes de reconstrution des images, à savoir le TRA (techniques de reconstruction algébrique), l'entropie maximum et
l'algorithme de transformée itérative de la présente inven-
tion. Il conclut que l'algorithme de transformée itérative offre la combinaison la meilleure de la qualité des images, de la stabilité du bruit et du temps de calcul pour les
reconstructions à angle limité.
La présente invention peut s'appliquer à d'autres
modalités d'imagerie et systèmes TI. Le microscope électro-
nique a un angle de couverture limité et utilise les mêmes - il - 2632427 méthodes de reconstruction. Un système TI à angle limité aux
ultra-sons est utilisé dans la caractérisation des défauts.
L'annexe A sur micro-fiche contient un listing du programme de reconstruction des images à angle limité mais ne comporte pas les sousprogrammes TI standards tels que la projection, la convolution et la rétroprojection qu'on peut
trouver dans le document "Donner Algorithms for Reconstruc-
tion Tomography, R.H. Huesman et al, 1977, Lawrence Berkeley
Laboratory, Université de Californie, Berkeley".
- 12 - 2632427
Claims (10)
1. Procédé pour former une image d'un objet par un mode opératoire perfectionné de reconstruction d'images à angle limité, caractérisé en ce qu'il comprend les étapes consistant à: a) procéder au balayage de l'objet (36) avec un agent d'imagerie dans une plage angulaire limitée et à l'acquisition de projections mesurées à des angles de vision disponibles; b) reconstruire une image partielle de la donnée mesurée à partir des projections mesurées; c) produire une image composite de l'objet en additionnant l'image partielle de la donnée mesurée et une image partielle d'une donnée calculée qui est reconstruite à partir des projections manquantes initialement mises à zéro; d) corriger l'image composite en utilisant une information a priori sur l'objet;
e) calculer à partir d'une image composite corri-
gée intermédiaire les projections manquantes à des angles de vision manquants; f) reconstruire une image partielle perfectionnée
de la donnée calculée due aux projections manquantes calcu-
lées dans l'étape précédente; et g) répéter de manière itérative les étapes c à f et améliorer progressivement l'image composite en utilisant
la même image partielle de la donnée mesurée et l'image par-
tielle successivement améliorée de la donnée calculée jusqu'à ce qu'on atteigne une qualité suffisante pour l'image.
2. Procédé selon la revendication 1, caractérisé en ce que l'agent d'imagerie est choisi dans le groupe
constitué des rayons X, d'ultra-sons et des faisceaux élec-
troniques.
- 13 - 2632427
3. Procédé selon la revendication 1, caractérisé en ce que la correction avec une information a priori dans l'étape d comprend 1) la mise à zéro des pixels de l'image composite se trouvant à l'extérieur d'une limite connue de l'objet, 2) la remise à une limite supérieure pré-choisie des pixels ayant une densité dépassant la limite supérieure,
et 3) la remise a zéro des pixels ayant une densité néga-
tive.
4. Procédé selon la revendication 3, caractérisé en ce que la reconstruction dans les étapes b et f est une
opération de rétroprojection filtrée.-
5. Procédé selon la revendication 4, caractérisé
en ce que le calcul dans l'étape e est une opération de pro-
jection.
6. Procédé selon la revendication 3, caractérisé en ce qu'il comprend en outre l'étape consistant à afficher une image composite finale produite après une multitude d'itérations.
7. Procédé de formation d'une image d'un objet dans un dispositif de tomographie informatisée à rayons X en utilisant un mode opératoire perfectionné de reconstruction d'images à angle limité, caractérisé en ce qu'il comprend
les étapes consistant à: -
a) procéder au balayage de l'objet avec des rayons X dans une plage angulaire limitée et à l'acquisition de projections mesurées aux angles de vision disponibles; b) reconstruire une image partielle d'une donnée mesurée à partir des projections mesurées par une opération de rétroprojection filtrée; c) produire une image composite de l'objet en additionnant l'image partielle de la donnée mesurée et une image partielle d'une donnée calculée qui est reconstruite à
partir des projections manquantes aux angles de vision man-
quants initialement remis.A zéro; d) corriger la densité de l'objet dans l'image
- 14 - 2632427
composition avec une information a priori sur l'objet en 1) mettant à zéro les pixels d'image à l'extérieur d'une limite
connue de l'objet, 2) remettant à une limite supérieure pré-
choisie les pixels ayant une densité dépassant cette limite supérieure, et 3) remettant à zéro les pixels ayant une den- sité négative; e) calculer les projections manquantes à partir
d'une image composite corrigée intermédiaire par une opé-
ration de projection; f) reconstruire une image partielle améliorée d'une donnée calculée par une opération de rétroprojection filtrée à partir des projections manquantes calculées dans l'étape précédente; et g) répéter de manière itérative les étapes c à f et améliorer progressivement l'image composite en utilisant
la même image partielle de la donnée mesurée et l'image par-
tielle successivement améliorée de la donnée calculée.
8. Procédé selon la revendication 7, caractérisé en ce que la correction de la densité de l'objet dans l'étape d(1) comprend l'étape de détermination de la limite de l'objet en se rapprochant d'une enveloppe convexe de
l'objet construite à partir des projections mesurées.
9. Procédé selon la revendication 7, caractérisé en ce qu'il comprend en outre l'étape du test de l'image composite corrigée intermédiaire pour sa convergence et
l'étape de sortie d'une image reconstruite finale.
10. Procédé selon la revendication 9, caractérisé en ce qu'il comprend en outre l'étape d'affichage de l'image finale.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/201,784 US5053958A (en) | 1988-06-06 | 1988-06-06 | Method to reduce image reconstruction time in limited-angle ct systems including using initial reconstruction valves for measured projection data during each iteration |
Publications (1)
Publication Number | Publication Date |
---|---|
FR2632427A1 true FR2632427A1 (fr) | 1989-12-08 |
Family
ID=22747279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR8906840A Withdrawn FR2632427A1 (fr) | 1988-06-06 | 1989-05-25 | Procede pour ameliorer sensiblement le temps de reconstruction d'une image dans une tomographie informatisee a angle limite |
Country Status (6)
Country | Link |
---|---|
US (1) | US5053958A (fr) |
JP (1) | JPH0251788A (fr) |
DE (1) | DE3918354A1 (fr) |
FR (1) | FR2632427A1 (fr) |
GB (1) | GB2221130B (fr) |
IT (1) | IT1229455B (fr) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2648304B1 (fr) * | 1989-06-12 | 1991-08-30 | Commissariat Energie Atomique | Procede de determination d'un espace a partir d'un espace discret connu pour la reconstruction d'images bi ou tridimensionnelles, dispositif de mise en oeuvre et application du procede |
US5233518A (en) * | 1989-11-13 | 1993-08-03 | General Electric Company | Extrapolative reconstruction method for helical scanning |
US5270926A (en) * | 1990-12-21 | 1993-12-14 | General Electric Company | Method and apparatus for reconstructing a three-dimensional computerized tomography (CT) image of an object from incomplete cone beam projection data |
US5396528A (en) * | 1991-06-28 | 1995-03-07 | General Electric Company | Tomographic image reconstruction using cross-plane rays |
DE4139150C1 (en) * | 1991-11-28 | 1993-06-24 | Siemens Ag, 8000 Muenchen, De | Computer tomograph with part ring formed X=ray source and detector - has double ring system without complementary interpolation |
US5414623A (en) * | 1992-05-08 | 1995-05-09 | Iowa State University Research Foundation | Optoelectronic system for implementation of iterative computer tomography algorithms |
US5611026A (en) * | 1992-12-21 | 1997-03-11 | General Electric Company | Combining a priori data with partial scan data to project three dimensional imaging of arbitrary objects with computerized tomography |
US5319693A (en) * | 1992-12-30 | 1994-06-07 | General Electric Company | Three dimensional computerized tomography scanning configuration for imaging large objects with smaller area detectors |
US5446776A (en) * | 1993-08-02 | 1995-08-29 | General Electric Company | Tomography with generation of radon data on polar grid points |
US5461651A (en) * | 1993-11-17 | 1995-10-24 | General Electric Company | Reconstruction of images in cone beam scanning with rectangular detector elements |
FR2721402B1 (fr) * | 1994-06-21 | 1996-08-02 | Commissariat Energie Atomique | Procédé d'obtention par radiographie tangentielle d'une image d'un objet en rotation autour d'un axe. |
US5805098A (en) * | 1996-11-01 | 1998-09-08 | The United States Of America As Represented By The Secretary Of The Army | Method and system for forming image by backprojection |
EP1356433A1 (fr) * | 2001-01-17 | 2003-10-29 | Watkins, Eserick Tanchwel | Procede de reconstruction d'image tomographique |
JP2002290868A (ja) * | 2001-03-27 | 2002-10-04 | Texas Instr Japan Ltd | 周波数変換回路、復調回路及びテレビ受信装置 |
JP2005538344A (ja) * | 2002-07-08 | 2005-12-15 | シデック テクノロジーズ アーベー | 画像化装置および方法 |
US6873677B2 (en) * | 2003-04-18 | 2005-03-29 | Leon Kaufman | Method and device for improving time resolution of an imaging device |
US7309985B2 (en) | 2003-04-23 | 2007-12-18 | Koninklijke Philips Electronics N. V. | Method of reconstructing an MR image |
US20050078861A1 (en) * | 2003-10-10 | 2005-04-14 | Usikov Daniel A. | Tomographic system and method for iteratively processing two-dimensional image data for reconstructing three-dimensional image data |
US7042974B2 (en) * | 2004-05-20 | 2006-05-09 | Eastman Kodak Company | Method for reconstructing tomographic images |
US7856079B2 (en) * | 2004-11-15 | 2010-12-21 | Koninklijke Philips Electronics N.V. | Reconstruction method for computer tomography and computer tomograph |
DE102005053022A1 (de) * | 2005-11-07 | 2007-05-16 | Siemens Ag | Verfahren und Vorrichtung zur räumlichen Darstellung eines Untersuchungsbereichs eines Untersuchungsobjekts |
CN101331516B (zh) * | 2005-12-15 | 2012-01-11 | 皇家飞利浦电子股份有限公司 | 用于多次迭代算法的高级收敛 |
DE102006025759A1 (de) * | 2006-05-31 | 2007-12-06 | Siemens Ag | Verfahren zur Bildrekonstruktion eines Objektes mit Projektionen und Vorrichtung zur Durchführung des Verfahrens |
US7680240B2 (en) * | 2007-03-30 | 2010-03-16 | General Electric Company | Iterative reconstruction of tomographic image data method and system |
US7920670B2 (en) * | 2007-03-30 | 2011-04-05 | General Electric Company | Keyhole computed tomography |
US8135186B2 (en) * | 2008-01-25 | 2012-03-13 | Purdue Research Foundation | Method and system for image reconstruction |
US7885371B2 (en) * | 2008-08-28 | 2011-02-08 | General Electric Company | Method and system for image reconstruction |
US8111893B2 (en) * | 2009-06-09 | 2012-02-07 | Wisconsin Alumni Research Foundation | Method for dynamic prior image constrained image reconstruction |
DE102009047867B4 (de) * | 2009-09-30 | 2016-10-06 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Korrektur von trunkierten Projektionsdaten |
US7916828B1 (en) * | 2010-01-06 | 2011-03-29 | General Electric Company | Method for image construction |
US8238518B2 (en) * | 2010-06-23 | 2012-08-07 | The Institute Of Cancer Research | Radiotherapy system |
US8189735B2 (en) * | 2010-07-22 | 2012-05-29 | General Electric Company | System and method for reconstruction of X-ray images |
DE102011083391A1 (de) * | 2011-09-26 | 2013-03-28 | Siemens Aktiengesellschaft | Verfahren und CT-System zur Erzeugung eines CT-Bilddatensatzes |
CN104182932B (zh) * | 2013-05-27 | 2017-04-12 | 株式会社日立制作所 | Ct 装置、ct 图像系统及ct 图像生成方法 |
CN103745488A (zh) * | 2013-12-31 | 2014-04-23 | 沈阳东软医疗系统有限公司 | 一种计算机断层成像中生成投影数据的方法和装置 |
JP6617368B2 (ja) * | 2015-06-08 | 2019-12-11 | 国立研究開発法人理化学研究所 | 電子源の作製方法 |
CN108010096A (zh) * | 2017-11-30 | 2018-05-08 | 上海联影医疗科技有限公司 | Cbct图像重建方法、装置和cbct设备 |
CN111080740B (zh) * | 2019-12-27 | 2023-06-16 | 上海联影医疗科技股份有限公司 | 一种图像校正方法、装置、设备及介质 |
CN112070856B (zh) * | 2020-09-16 | 2022-08-26 | 重庆师范大学 | 基于非下采样轮廓波变换的有限角c型臂ct图像重建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2613487A1 (fr) * | 1987-04-01 | 1988-10-07 | Gen Electric | Procede pour obtenir une information sur les limites d'un objet dans la tomographie informatisee a angle limite |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4506327A (en) * | 1981-11-23 | 1985-03-19 | General Electric Company | Limited-angle imaging using multiple energy scanning |
US4866614A (en) * | 1985-12-26 | 1989-09-12 | General Electric Company | Ultrasound characterization of 3-dimensional flaws |
DE3716988A1 (de) * | 1986-06-23 | 1988-01-14 | Gen Electric | Verfahren zum rekonstruieren von objekten aus abtastungen mit begrenztem winkel bei der computertomographie |
US4922421A (en) * | 1986-12-15 | 1990-05-01 | General Electric Company | Flaw characterization by multiple angle inspection |
-
1988
- 1988-06-06 US US07/201,784 patent/US5053958A/en not_active Expired - Fee Related
-
1989
- 1989-05-25 FR FR8906840A patent/FR2632427A1/fr not_active Withdrawn
- 1989-06-02 JP JP1139411A patent/JPH0251788A/ja active Pending
- 1989-06-05 GB GB8912867A patent/GB2221130B/en not_active Expired - Fee Related
- 1989-06-06 DE DE3918354A patent/DE3918354A1/de not_active Withdrawn
- 1989-06-06 IT IT8920799A patent/IT1229455B/it active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2613487A1 (fr) * | 1987-04-01 | 1988-10-07 | Gen Electric | Procede pour obtenir une information sur les limites d'un objet dans la tomographie informatisee a angle limite |
Non-Patent Citations (1)
Title |
---|
APPLIED OPTICS vol. 20, no. 3, Février 1981, NEW YORK.USA TAKUSO SATO ET AL: 'tomographic image reconstruction from limited projections using iterative revisions in image and transform spaces' * |
Also Published As
Publication number | Publication date |
---|---|
GB8912867D0 (en) | 1989-07-26 |
DE3918354A1 (de) | 1990-01-11 |
GB2221130B (en) | 1992-07-08 |
IT8920799A0 (it) | 1989-06-06 |
GB2221130A (en) | 1990-01-24 |
IT1229455B (it) | 1991-08-08 |
US5053958A (en) | 1991-10-01 |
JPH0251788A (ja) | 1990-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2632427A1 (fr) | Procede pour ameliorer sensiblement le temps de reconstruction d'une image dans une tomographie informatisee a angle limite | |
JP3373720B2 (ja) | X線断層撮影装置 | |
US7023951B2 (en) | Method and apparatus for reduction of artifacts in computed tomography images | |
US8625870B2 (en) | Method and system for supplementing detail image in successive multi-scale reconstruction | |
US4888693A (en) | Method to obtain object boundary information in limited-angle computerized tomography | |
Kamphuis et al. | Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm | |
JPH07109621B2 (ja) | 不完全な円錐状ビーム投射データから物体の三次元映像を再構成する方法および装置 | |
US11727609B2 (en) | Limited-angle CT reconstruction method based on anisotropic total variation | |
US20030118151A1 (en) | Image reconstruction using multiple X-ray projections | |
US20080165920A1 (en) | Method and apparatus for reduction of metal artifacts in ct imaging | |
FR2872659A1 (fr) | Procede et appareil pour une reconstruction directe en imagerie par tomosynthese | |
FR2736455A1 (fr) | 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 | |
FR2848007A1 (fr) | Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal | |
US4626991A (en) | System for reprojecting images acquired by backprojection | |
JP6505513B2 (ja) | X線コンピュータ断層撮像装置及び医用画像処理装置 | |
CN111127579A (zh) | 一种基于神经网络去锥形束伪影的ct重建方法 | |
US6850587B1 (en) | Reprojection-based three-dimensional image reconstruction | |
US7088800B2 (en) | X-ray CT system and beam-hardening post-processing method | |
US20100054561A1 (en) | System and method for image reconstruction | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
US10417795B2 (en) | Iterative reconstruction with system optics modeling using filters | |
KR100923094B1 (ko) | 트렁케이션 아티팩트를 보정하는 방법 | |
Van Gompel et al. | Elliptical extrapolation of truncated 2D CT projections using Helgason-Ludwig consistency conditions | |
JP4387758B2 (ja) | Spect装置及びspect画像再構成方法 | |
JP2016039902A (ja) | X線コンピュータ断層撮影装置及び医用画像処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |