DE3420576C2 - Anordnung zum Reprojizieren von Bildern aus mehreren eindimensionalen Projektionen in der Computer-Tomographie - Google Patents
Anordnung zum Reprojizieren von Bildern aus mehreren eindimensionalen Projektionen in der Computer-TomographieInfo
- Publication number
- DE3420576C2 DE3420576C2 DE19843420576 DE3420576A DE3420576C2 DE 3420576 C2 DE3420576 C2 DE 3420576C2 DE 19843420576 DE19843420576 DE 19843420576 DE 3420576 A DE3420576 A DE 3420576A DE 3420576 C2 DE3420576 C2 DE 3420576C2
- Authority
- DE
- Germany
- Prior art keywords
- dimensional
- data
- projections
- arrangement
- image
- 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.)
- Expired - Lifetime
Links
- 238000003325 tomography Methods 0.000 title claims description 4
- 230000001131 transforming Effects 0.000 claims description 22
- 230000000875 corresponding Effects 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims 2
- 238000007781 pre-processing Methods 0.000 claims 1
- 238000000034 methods Methods 0.000 description 14
- 238000002591 computed tomography Methods 0.000 description 4
- 239000002184 metals Substances 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 238000000844 transformation Methods 0.000 description 3
- 210000000988 Bone and Bones Anatomy 0.000 description 2
- 239000004793 Polystyrene Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000010586 diagrams Methods 0.000 description 2
- 239000011159 matrix materials Substances 0.000 description 2
- 230000003252 repetitive Effects 0.000 description 2
- 281999990635 Foundations companies 0.000 description 1
- 238000010521 absorption reactions Methods 0.000 description 1
- 238000000205 computational biomodeling Methods 0.000 description 1
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Chemical compound data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgZG9taW5hbnQtYmFzZWxpbmU9ImNlbnRyYWwiIHRleHQtYW5jaG9yPSJzdGFydCIgeD0nMTI1Ljk4MScgeT0nMTU2JyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7ZmlsbDojM0I0MTQzJyA+PHRzcGFuPkZlPC90c3Bhbj48L3RleHQ+Cjwvc3ZnPgo= data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgZG9taW5hbnQtYmFzZWxpbmU9ImNlbnRyYWwiIHRleHQtYW5jaG9yPSJzdGFydCIgeD0nMTguNzk5NicgeT0nNDcuNzk1NScgc3R5bGU9J2ZvbnQtc2l6ZTozOHB4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO2ZpbGw6IzNCNDE0MycgPjx0c3Bhbj5GZTwvdHNwYW4+PC90ZXh0Pgo8L3N2Zz4K [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006011 modification reactions Methods 0.000 description 1
- 230000000051 modifying Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000002250 progressing Effects 0.000 description 1
- 230000003595 spectral Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING; 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; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Description
Die Erfindung betrifft eine Anordnung nach dem Oberbegriff des Anspruches 1 sowie ein Verfahren nach dem Oberbegriff des Anspruches 5.
Die Entwicklung von rechnergesteuerten Röntgen-Tomographiegeräten (CT) hat Abtastvorrichtungen mit kürzeren Datenerfassungs- und Bildrekonstruktionszeiten sowie mit verbesserter Dichte und räumlicher Auflösung ergeben. Die Verbesserungen sind hauptsächlich dadurch erreicht worden, daß komplexere Datenerfassungsanordnungen und eine schnellere Bildkonstruktions-Hardware verwendet wurden. Die Bildqualität ist ferner dadurch verbessert worden, daß Annahmen neu bewertet wurden, die bei der Konstruktion früherer Generationen von CT-Abtastgeräten und beim Einführen von Korrekturen und/oder Verbesserungen in diesen Annahmen innerhalb des Bildrekonstruktionsalgorithmus gemacht wurden. Diese Annahmen wurden anfangs zugrunde gelegt, um die Daten, die durch ein Abtastgerät gewonnen wurden, mit theoretischen Rekonstruktionsalgorithmen kompatibel zu machen.
Ein Beispiel für diese Annahmen betrifft das Spektrum der Röntgenquelle und die Energieabhängigkeit der Dämpfungskoeffizienten von unterschiedlichen Elementen des zu untersuchenden Objektes. Eine wichtige Annahme, die bisher gemacht wurde, um Bilder zu erzeugen, ist die, daß die Quelle monochromatisch ist oder daß die Energieabhängigkeit von Dämpfungskoeffizienten für alle Elemente identisch ist. Es ist bekannt, daß keine dieser beiden Bedingungen erfüllt ist, und es entstehen somit in den gewonnenen Bildern die als polychromatische Artefakte bekannten Erscheinungen. Die Artefakte können als mit "cupping" bezeichnete Vertiefungen und als negative Schlieren zwischen scharfen Objekten, die hohe Dämpfungskoeffizienten haben, identifiziert werden.
Bei einer bekannten Anordnung, z. B. nach US-PS 42 17 641, wird eine sich wiederholende Nach-Rekonstruktionsmethode verwendet, um den Pegel von polychromatischen Artefakten zu, verringern. Andere bekannte Druckschriften beschreiben polychromatische Artefakt-Korrekturtechniken, z. B. US-PS 42 22 104 und 42 23 384 wie auch ein Aufsatz "A Framework for Spectral Artifact Corrections in X-Ray Computed Tomography" von J. Peter Stonstrom et al, erschienen in IEEE Transactions on Biomedical Engineering, Band BME-28, Nr. 2, Februar 1981.
Grundlage dieser bekannten Nach-Rekonstruktions-Korrekturverfahren ist, daß Objekte aus zwei annähernd homogenen Komponenten in bezug auf die Energieabhängigkeit ihrer Dämpfungskoeffizienten bestehen. In biologischen Anwendungsfällen sind die beiden Komponenten Knochen und Weichgewebe. Es wird ein Anfangsbild rekonstruktiert, das polychromatische Korrekturen erster Ordnung für das Majoritätselement, üblicherweise Weichgewebe, enthält. Das Anfangsbild wird Bildelement um Bildelement in Segmente geteilt, um angenäherte Bilder der beiden Komponenten zu erzeugen. Anschließend werden die Pfadlängen durch die beiden Bilder unter Verwendung einer Reprojektionstechnik berechnet. Es werden Fehlerprojektionen aus den Reprojektionen gebildet und den Projektionsdaten hinzuaddiert, die zur Herstellung des Anfangsbildes verwendet wurden. Es wird ein Bild zweiter Ordnung aus den neuen Prolektionsdaten rekonstruiert. Wenn der Pegel der polychromatischen Korrektur ausreichend hoch ist, ist der Algorithmus vollständig. Ist dies nicht der Fall, wird das vorstehend erläuterte Verfahren wiederholt.
Die Anwendung der Reprojektion ist nicht auf polychromatische Korrekturalgorithmen beschränkt. Der Aufsatz "An Algorithm for the Reduction of Metal Clip Artifacts in CT Reconstructions" von G. H. Glover und N. J. Pelc, in Medical Physics, Band 8, Nr. 6, November 1981 beschreibt ein Verfahren zum Entfernen der durch Metallclips verursachten Artefakte unter Verwendung einer Reprojektion als Teil ihres Algorithmus. Der Aufsatz "A Simple Computational Method for Reducing Streak Artifacts in CT Images" von G. Henrich, in Computed Tomography, Band 4, 1981 beschreibt einen Algorithmus, der verwendet werden kann, um Schlieren zu entfernen, z. B. solche, die durch. Teilvolumenartefakte verursacht werden.
Die polychromatischen Metallclip- und Schlieren-Artefakt-Korrekturalgorithmen, die zum Stand der Technik beschrieben worden sind, sind nicht praktisch umgesetzt worden, weil der Reprojektionsschritt außerordentlich zeitaufwendig war. Die bekannten Reprojektionsmethoden waren zu langsam, weil sie auf dem eingeprägten Reprojektionsschritt beruhten, der in den Rekonstruktionsalgorithmen auf der Grundlage von algebraischen Techniken enthalten ist. Die geringe Geschwindigkeit der bekannten Reprojektionsanordnungen und ein Versuch einer Lösung sind in einem Aufsatz "Algorithmus for Fast Back- und Re-Projection in Computed Tomography" von T. M. Peters, erschienen in IEEE Transadions on Nuclear Science, Band NS-28, Nr. 4, August 1981, S. 3641-3647 erläutert. Dieser Aufsatz beschreibt ein Verfahren, das einen modifizierten Rückprojektor zur Erzielung von Reprojektionen verwendet. Das Problem bei dieser Anordnung ist, daß die Modifikationen die Hardware eines Rückprojektors radikal ändern und damit die Anordnung nicht einfach für kommerzielle Anwendungsfälle einsetzbar ist. Die Anordnung macht eine Vorrichtung zum Reversieren des normalen Datenflusses durch den Rückprojektor erforderlich, was zu Reprojektionen am normalen Eingang der Vorrichtung führt. Zusätzlich haben die resultierenden Reprojektionen eine schlechte Qualität und benötigen komplexe Korrekturen, um sie mit einem Artefakt-Korrekturalgorithmus verwenden zu können.
Aus der US 42 05 375 ergibt sich ein Projektions- und Rückprojektionsverfahren, nicht aber ein Reprojektionsverfahren. Der Scanner eines computergesteuerten Tomogra phen dreht um weniger als 180° und liefert trotzdem genügend Bilddaten, die 360° abdecken. Die Absorptionsdaten werden in einer Vielzahl von Ansichten, die einen Winkel θ < 180° abdecken, gesammelt. Die gesammelten Daten werden durch Betrachtung von Winkeln θ1, θ2, . . ., θn sortiert. Die Daten für die nicht durchlaufenen 180° werden aus den sortierten Daten nach einer Fourier-Transformation der sortierten Daten berechnet. Die vollständigen, d. h. die gemessenen und die berechneten Daten werden im Rekonstruktionsvorgang verwendet, um das Bild zu erhalten. Dabei wird eine zweidimensionale Fourier-Transformation vorgeschlagen, um Rekonstruktionsda ten aus den Probenwerten zu gewinnen, die aus dem Funktionsdatenspeicher erhalten werden. Diese Werte werden aus, der Fourier-Transformationseinheit entnommen, die Koordinaten der Werte von Winkel- in Rechteckform umgewandelt und eine zweidimensionale inverse Fourier-Transformationsvorrichtung zur Erzielung der Bilddaten verwendet. Hieraus ergibt sich nicht die Lehre, zweidimensionale Four ier-Transformationen eines Bildes zu verwenden, um Reprojektionen zu gewinnen.
Aus der Druckschrift "Radiology" Band 117, Dezember 1975, Seiten 561-572 (Brooks) ist die Verwendung von Fourier-Transformationen für Rückprojektionen bereits seit 1976 bekannt. Trotzdem sind bisher keine Tatsachen bekannt geworden, die die Verwendung von Fourier-Transformationen zur Erzielung von Reprojektionen für Korrekturzwecke als bekannt oder naheliegend ergeben. Die Anwendung von Reprojektionen für Korrekturzwecke war am Anmeldetag vorliegender Erfindung ebenfalls bekannt. Die Anwendung derartiger Reprojektionen ist jedoch wegen der relativ langen Zeitdauer, die für die Erzeugung und Anwendung von Reprojektionen erforderlich war, nie verwirklicht oder in Betracht gezogen worden, da derartige bekannte Reprojektionssysteme für praktische Anwendungsfälle zu langsam waren. Der Reprojektionsvorgang verwendet "Pseudo-Projektionen", die aus Pixel-Werten konstruiert und dann korrigiert werden, um die zu korrigierenden Artefakte zu beseitigen. Anschließend werden die korrigierten Pseudo-Projektionen gefiltert, rückprojiziert und dem Bild hinzugefügt, so daß das korrigierte Bild zur Anzeige gebracht wird.
Aufgabe der Erfindung ist, die für das Erzeugen und Anwenden von Reprojektionen benötigte Zeit zum Korrigieren von Artefakten; die bei rechnergesteuerten Rönt gen-Tomographiegeräten auftreten, entscheidend zu verkürzen.
Gemäß der Erfindung wird dies mit einer Anordnung mit den Merkmalen des Kennzeichens des Anspruches 1 bzw. mit einem Verfahren mit den Merkmalen des Kennzeichens des Anspruches 5 erreicht. Weitere Ausgestaltungen der Erfindung sind Gegenstand der Unteransprüche.
Mit der erfindungsgemäßen Anordnung bzw. dem erfindungsgemäßen Verfahren wird erreicht, daß die Zeit zum Erzeugen und Anwenden der Reprojektionen unter Verwendung einer Transformation entscheidend verkürzt wird. Die Verwendung der Fourier-Transformationen zur Umwandlung der Bilddaten in "Pseudo-Projektionen" oder Reprojektionen erfolgt im Falle der Erfindung in der Weise, daß die Reprojektionen in eine Korrektureinheit eingeführt werden, wo sie mit den ursprünglichen Projektionen kombiniert, gefiltert und dann zur Erzielung der korrigierten Bilddaten rückprojiziert werden. Die Verwendung eines Fourier-Al gorithmus in einem Rückprojektionsvorgang bei der Erzeugung von Bildern hat keinen Einfluß auf die Verwendung von Fourier-Transformationen nach der Erfindung, um Reprojektionen zur Verwendung bei einem Rückprojektor für Korrekturzwecke zu erzielen. Im Falle der Erfindung wird nicht der inverse Wert des Rückprojektionsvorganges zur Erzielung von Reprojektionen verwendet, sondern Fourier-Transformationen.
Das Problem der Reprojektion zum Korrigieren der Rückprojektion bestand bisher in der verhältnismäßig langen Zeitdauer, die für das Einführen von Reprojektionen erforderlich ist. Durch den Einsatz von Fourier-Transformationen zur Verringerung der Anzahl von Multiplikationen, die für die Reprojektion erzeugt werden müssen, läßt sich die für das Erzeugen und Anwenden der Rückprojektionen erforderliche Zeit verringern. Die Fourier-Transformationen machen in vorteilhafter Weise Gebrauch von sich wiederholenden Sinus- und Cosinus-Funktionen. Fourier-Transformationen stützen sich darauf, daß Signale in eine Vielzahl von Sinus- und Cosinus-Funktionen zerlegt werden können. Die sich wiederholende Art der Sinus- und Cosinus-Funktionen reduziert die Anzahl von Vervielfachungen, die zur Bewertung des Fourier-Integrals erforderlich sind, erheblich. Funktionen von Raum oder Zeit können durch eine Summierung von Sinus- und Cosinus-Wellen (Harmonischen unterschiedlicher Frequenzen) dargestellt werden. Die Amplituden einer jeden Harmonischen werden als Fourier-Koeffizient bezeichnet. Jede zweidimensionale verteilungsdichte Funktion kann als Summe von Sinus- und Cosinus-Wellen ausgedrückt werden, die in verschiedenen Richtungen innerhalb einer Ebene fortschreiten. Die Amplituden der Wellen werden durch die Fourier-Koeffizienten bezeichnet. Grundlage der Fourier-Re projektionen ist, daß Fourier-Koeffizienten der Bilder auf die Fourier-Koeffizienten der Reprojektionen bezogen sind. Die Amplitude der in einem Winkel fortschreitenden Wellen ist gleich dem Fourier-Koeffizienten der Reprojektion bei diesem Winkel.
Mit vorliegender Erfindung werden diese charakteristischen Eigenschaften der Fou rier-Transformation zur Beschleunigung des Korrekturvorganges verwendet und damit die Reprojektionen zu einem brauchbaren Verfahren für das Korrigieren von Bildartefakten gemacht.
Nachstehend wird die Erfindung in Verbindung mit der Zeichnung anhand eines Ausführungsbeispieles erläutet. Es zeigt
Fig. 1 ein Blockschaltbild einer Ausführungsform einer Reprojektionsanordnung nach der Erfindung,
Fig. 2 ein erweiteres Blockschaltbild der Reprojektionsanordnung nach Fig. 1.
Die CT-Abtastanordnung 11 nach Fig. 1 weist ein Portal 12 auf. Das Objekt wird der Strahlung im Portal 12 ausgesetzt und die Strahlung wird, nachdem sie das Objekt durchlaufen hat, erfaßt. Die erfaßten Signale werden in einer Schaltanordnung FEE behandelt. Die elektrischen Signale werden durch einen Prozessor 14 vorverarbeitet. Die Ausgangssignale aus 14 werden Projektionen genannt. Der Ausgang aus 14 wird in die Korrekturvorrichtung 15 eingeführt, die eine polychromatische Korrektur erster Ordnung an den Projektionen ausführt, um das Anfangsbild zu erzeugen.
Der Ausgang der Korrekturvorrichtung 15 wird in den Filter 16 eingeführt. Die gefilterten Projektionen werden durch den Rückprojektor 17 rückprojiziert. Der Ausgang des Rückprojektors wird in Form eines digitalisierten Bildes in die Matrix 18, eingeführt. Die Matrix aus digitalisierten Daten wird verwendet, um auf dem Anzeigegerät 19 Bilder zu erzeugen.
Polychromatische Fehlerkorrekturen werden in der Rückkopplungsschleife 21 vorgenommen, die sich vom Ausgang des Rückprojektors 17 zur polychromatischen Fehlerkorrekturvorrichtung 15 erstreckt. Im polychromatischen Korrekturbetrieb wird der Ausgang der Rückkopplungsschleife 21 mit dem Ausgang aus 14 kombiniert und dann in den Filter 16 eingeführt.
Die Rückkopplungsschleife 21 nimmt den Bildprozessor 22 und den Reprojektor 23 bzw. die Rückkopplungsvorrichtung auf. Der Bildprozessor 22 ist bei einer bevorzugten Ausführungsform der Erfindung so ausgelegt, daß er zwischen Bildelementen von Knochen und Weichgewebe unterscheidet. Der Ausgang aus dem Bildprozessor 22 ist mit dem Reprojektor 23 verbunden.
Eine erweiterte Darstellung des Reprojektors 23 ergibt sich aus Fig. 2. Der Ausgang des Bildprozessors 22 wird in die Vorrichtung 31 eingeführt, die eine zweidimensionale Fourier-Transformation berechnet. Die Werte der Fourier-Transformation, die am Ausgang der Vorrichtung 31 verfügbar sind, werden in der Vorrichtung 32 interpoliert, um radiale Linien der Transformation des verarbeiteten Bildes zu erhalten. Reprojektionen des verarbeiteten Bildes werden durch die Vorrichtung 33 erhalten, die die eindimensionale inverse Fourier-Transformation der radialen Linien der zweidimensionalen Transformation des verarbeiteten Bildes ergibt.
Um das Verständnis der Anordnung nach der Erfindung zu verbes sern, worden nachstehend die mathematischen Grundlagen der Reprojektion erläutert.
Es wird die Funktion f (x, y) betrachtet, die eine Rekonstruk
tion eines Querschnittes eines Gegenstandes und des mit (θ, t)
charakterisierten Pfades wie folgt darstellt:
t = x.cos(θ) + y.sin(θ) (1),
wobei
|t| < ∞ |θ| < π/2.0.
Eine Probe, p (θ, t) der Reprojektion der Objektfunktion längs
des mit (θ, t) charaktererisierten Pfades ergibt sich zu:
wobei δ(z) eine normale δ-Funktion ist, die wie folgt be
schrieben wird:
Die Integration in (2) erfolgt über Streifen mit der Breite Null. Ein Streifen mit Null-Breite kann in (2) eingesetzt werden, indem die δ-Funktion durch den normalisierten Quer schnitt des Streifens ersetzt wird. Die Normalisierung gewähr leistet, daß das Integral des Querschnittes der gewünschten Öffnungsfunktion Eins ist.
Die inverse Fourier-Scheibentheorie wird nachstehend erläutert, um zu zeigen, daß eine inverse Fourier-Transformation einer radialen Linie der zweidimensionalen Fourier-Transformation der Objektfunktion eine Reprojektion der gleichen Objektfunktion ist.
Es sei F (u, v) die zweidimensionale Fourier-Transformation der
Objektfunktion f (x, y). Die Transformation ergibt sich als:
Nunmehr wird F (u, v) in einem Polarkoordinatensystem betrach
tet. Dabei seien ω und ϕ die variablen Größen, die das Polarko
ordinatensystem charakterisieren. Sie sind auf "u" und "v"
bezogen wie folgt:
u = ω.cos(ϕ) (5)
v = ω.sin(ϕ) (6).
F (ϕ, ω) ergibt sich durch Einsetzen von (5) und (6) in (4):
wobei angenommen wurde, daß der Ausdruck F (ϕ, ω) eine Polarko
ordinatendarstellung anstelle des rechtwinkeligen Koordinaten
systems mit der Verwendung von F (u, v) einschließt, und daß "w"
gegeben ist durch:
w = x.cos(ϕ) + y.sin(ϕ) (8).
Für einen festen Wert von ϕ stellt F (ϕ, ω) eine radiale Linie
der zweidimensionalen Fourier-Transformation der Objektfunktion
dar. Es sei die eindimensionale inverse Fourier-Transformation,
g (ϕ, z), einer radialen Linie gegeben durch einen festen Wert
von ϕ. Die Funktion g (ϕ, z) ist gegeben durch:
Gleichung (9) kann angesetzt werden, wenn man (7) in (9)
einsetzt. Dies ergibt:
Nun wird die Reihenfolge der Integration in (10) geändert und
erhalten:
Es läßt sich auf einfache Weise zeigen, daß das innere Integral
(11) sich reduziert auf:
δ(z - x.cos(ϕ) - y.sin(ϕ)) (12),
wobei (12) durch Ersetzen der Definition von "ω" nach (8) erhalten wurden.
Setzt man (12) in (11) ein, erhält man:
Nunmehr wird die Standarddefinition einer Reprojektion nach
(2)-(13) mit θ und "t", ersetzt durch ϕ und "z", verglichen. Daraus
ergibt sich, daß g (ϕ, z) die Reprojektion bei ϕ ist. Diese
Tatsache führt zu folgender Methode, um Reprojektionen eines
Gegenstandes festzustellen:
- 1. Bestimmung der zweidimensionalen Fourier-Transformation einer Objektfunktion.
- 2. Verwendung der eindimensionalen inversen Fourier-Transforma tion einer radialen Linie der Fourier-Transformation der Objektfunktion zur Erzielung einer Reprojektion bei dem die radiale Linie kennzeichnenden Winkel.
Bei einer tatsächlichen Ausführung des oben erläuterten Verfah rens müssen diskrete Fourier-Transformationen (DFT) anstelle der eingeschlossenen kontinuierlichen Fourier-Transformationen verwendet werden. Es wird nunmehr gezeigt, wie das obige Verfahren erweitert werden kann, so daß es in einer Hardware verwirklicht werden kann.
Es wird die diskrete Folge x (i) für i = 0, 1, . . ., N-1 betrach
tet. Die diskrete Fourier-Transformation X (k) für k = 0, 1, . . .,
N-1, ist gegeben durch:
Die inverse diskrete Fourier-Transformation, IDFT, ist gegeben
durch:
Es sei angenommen, daß x (i) Proben in der Zeit von x(t)
darstellt. Der Abstand zwischen den Proben ist durch d_t
gegeben. Es läßt sich unter bestimmten Bedingungen zeigen, daß
X(k) Proben der kontinuierlichen Fourier-Transformation, X(f),
von x(t) darstellt. Der Abstand zwischen den Proben in der
Frequenzdomäne, d_f, kann wie folgt gezeigt werden:
d_f = 1/(N.d_t) (16).
Es sei angenommen, daß x(i) in folgender Weise auf x(t) bezogen
ist:
x(i) = x(i+d_t), i = 0, 1, . . ., N/2
x(i) = x ([i-N].d_t), i = N/2 + 1, . . ., N-1 (17).
Dann ist X(k) auf X(f) (innerhalb eines Skalenfaktors) wie folgt
bezogen:
X(k) = X(k.d_f), k = 0, 1, . . ., N/2
X(k) = X([k-N].d_f), k = N/2 + 1, . . ., N-1 (18).
Die mit (17) und (18) gegebenen Beziehungen sind eine Folge der Periodizitätsanforderungen einer DFT. Aus diesen Gleichungen ergibt sich, daß die negativen Zeit- oder Frequenzdomänenteile den Positionen der positiven Domänen folgen. Die Reihenfolge kann durch Modulieren der Folgen um eine abwechselnde +1/-1 Folge in beiden Domänen reversiert werden.
Die Verwendung der eindimensionalen DFT wird nunmehr so erwei tert, daß sie verwendet werden kann, um die zweidimensionale DFT eines Bildes zu finden.
Es sei angenommen, daß das rekonstruierte Bild in einem Kreis
mit dem Radius R0 enthalten sein kann. Auch sei angenommen, daß
eine NPIC × NPIC Rekonstruktion dieser Objektfunktion gemacht
wird. Die Rekonstruktion, f(i, m), kann auf die originale
Objektfunktion wie folgt bezogen werden:
f(i, m) = f(x, y) (19),
für
x = -R0 + i.DGRID (20a)
y = -R0 + m.DGRID (20b)
DGRID = 2.0.R0/NPIC (21),
und "i" und "m" liegen im Bereich (0, NPIC).
Die zweidimensionale diskrete Fourier-Transformation von f(i,m),
F(k,l) ist definiert durch:
Da der Exponentialausdruck in (22) trennbar ist, kann er wie
folgt ausgedrückt werden:
wobei G(k,m) gegeben ist durch:
Bei einer Prüfung stellt man fest, daß sowohl (23) als auch (24) standardförmige eindimensionale DFTen darstellen. Die mit den Gleichungen (23) und (24) eingeschlossene Methode besteht darin, zuerst die DFTen aller Reihen des Bildes und dann die DFTen der Spalten der DFTen der Reihen zu finden.
Wenn die DFTen anstatt kontinuierlicher Fourier-Transformationen verwendet werden, besteht eine Schwierigkeit darin, daß die resultierenden diskreten Transformationen auf einem rechteckför migen Gitter zur Verfügung stehen. Die gewünschten radialen Linien stellen Proben dar, die in einer Polarkoordinatendarstel lung der Transformation erhalten werden. Es ist somit eine Interpolation notwendig, um die rechtwinkelige Daratellung in eine Polarkoordinatendarstellung umzusetzen. Dabei kann unter Verwendung der Gleichungen (19)-(21) die Gleichung (18) auf zwei Dimensionen erweitert werden, so daß F(k,l) auf F(u,v) bezogen werden kann, um die Konstanten abzuleiten, die durch den Interpolationsvorgang erforderlich sind.
Verwendet man eine zweidimensionale Interpolation, können die in Gleichung (23) gegebenen Proben verwendet werden, um die Proben der Fourier-Transformation längs einer radialen Linie zu finden. Diese Funktion kann dann in eine IDFT geleitet werden, um die Reprojektion der Objektfunktion bei dem Winkel der radialen Linie zu finden.
Es ist bekannt, daß eine DFT unter Verwendung optimaler Techniken als die direkte Anwendung von (14) oder (15) verwirk licht werden kann. Diese Methoden sind als Schnell-Fourier-Trans formationen (FFT) bekannt. Es läßt sich zeigen, daß die FFT in wesentlich kürzerer Zeit durchgeführt werden können als die entsprechenden DFT für Vektoren mit einer großen Anzahl von droben. Die Zeitersparnis, die durch Verwendung einer FFT erreicht wird, reduziert auch entscheidend die Zeit, die erforderlich ist, um Reprojektionen zu erhalten, so daß dieses System für klinisch brauchbare Artefaktkorrekturen verwendet werden kann.
Claims (9)
einer Vorrichtung (14) zum Vorverarbeiten der erfaßten Strahlung, um Daten, zu gewinnen, die eindimensionalen Projektionen der erfaßten Strahlung entsprechen, einer Vorrichtung (16) zum Filtern der Daten der eindimensionalen Projektionen, einer Vorrichtung (17) zum Rückprojizieren der gefilterten Daten, um zweidimensionale digitale Bilddaten zu gewinnen, und einer Rückkopplungsvorrichtung (23) zum Rückkoppeln der zweidimensionalen Bilddaten, zu den Daten, die der Vorrichtung (16) zum Filtern zugeführt werden, dadurch gekennzeichnet, daß die Rückkopplungsvorrichtung (23) aufweist:
- (a) eine Vorrichtung (31) zum Durchführen der zweidimensionalen Fou rier-Transformation, an den zweidimensionalen Bilddaten,
- (b) eine Vorrichtung (32) zur Berechnung der Werte der Fourier-Transformation der Bilddaten längs radialer Linien durch Interpolieren, und
- (c) eine Vorrichtung (33) zur Bildung eindimensionaler inverser Fou rier-Transformationen der radialen Linien, um den Projektionen entsprechende Reprojektionsdaten zu erhalten.
die durch die zweidimensionale Fourier-Transformation gewonnenen Bilddaten interpoliert werden, um Werte längs radialer Linien zu berechnen und
die eindimensionale inverse Fourier-Transformation der radialen Linien gebildet wird, um den Projektionen entsprechende Reprojektionsdaten zu erhalten.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US06/501,832 US4616318A (en) | 1983-06-07 | 1983-06-07 | System for reprojecting images using transform techniques |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| DE3420576A1 DE3420576A1 (de) | 1985-02-07 |
| DE3420576C2 true DE3420576C2 (de) | 1999-04-29 |
Family
ID=23995189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| DE19843420576 Expired - Lifetime DE3420576C2 (de) | 1983-06-07 | 1984-06-01 | Anordnung zum Reprojizieren von Bildern aus mehreren eindimensionalen Projektionen in der Computer-Tomographie |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US4616318A (de) |
| JP (1) | JP2720894B2 (de) |
| DE (1) | DE3420576C2 (de) |
| FR (1) | FR2547439B1 (de) |
| NL (1) | NL192189C (de) |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH02121636A (en) * | 1988-10-31 | 1990-05-09 | Yokogawa Medical Syst Ltd | Image processing method and apparatus for reducing streak |
| US5008822A (en) * | 1988-11-25 | 1991-04-16 | Picker International, Inc. | Combined high speed backprojection and forward projection processor for CT systems |
| US5079697A (en) * | 1989-05-01 | 1992-01-07 | The General Hospital Corporation | Distortion reduction in projection imaging by manipulation of fourier transform of projection sample |
| US5414803A (en) * | 1991-01-11 | 1995-05-09 | Hewlett-Packard Company | Method utilizing frequency domain representations for generating two-dimensional views of three-dimensional objects |
| US5404293A (en) * | 1991-06-11 | 1995-04-04 | The University Of Utah | Cone beam reconstruction using helical data collection paths |
| US5170439A (en) * | 1991-06-11 | 1992-12-08 | Picker International, Inc. | Cone beam reconstruction using combined circle and line orbits |
| US5701360A (en) * | 1993-05-07 | 1997-12-23 | Siemens Aktiengesellschaft | Fourier reconstruction of computer tomography images which represent a selectable region of the examination subject |
| DE4315279C2 (de) * | 1993-05-07 | 1996-02-22 | Siemens Ag | Verfahren zum Erzeugen eines Schichtbildes mit Hilfe eines Computertomographen |
| US5438602A (en) * | 1993-12-23 | 1995-08-01 | General Electric Company | Correction of CT attenuation data using fan beam reprojections |
| AU3341295A (en) * | 1994-08-29 | 1996-03-22 | Eskofot A/S | A method of estimation |
| IL121773D0 (en) * | 1997-09-15 | 1998-02-22 | Elscint Ltd | Method for improving CT images having high attenuation objects |
| US6021215A (en) * | 1997-10-10 | 2000-02-01 | Lucent Technologies, Inc. | Dynamic data visualization |
| US6351548B1 (en) | 1999-06-23 | 2002-02-26 | The Board Of Trustees Of The University Of Illinois | Fast hierarchical reprojection algorithm for tomography |
| US6282257B1 (en) | 1999-06-23 | 2001-08-28 | The Board Of Trustees Of The University Of Illinois | Fast hierarchical backprojection method for imaging |
| US6332035B1 (en) | 1999-06-23 | 2001-12-18 | The Board Of Trustees Of The University Of Illinois | Fast hierarchical reprojection algorithms for 3D radon transforms |
| US6263096B1 (en) | 1999-06-23 | 2001-07-17 | The Board Of Trustees Of The University Of Illinois | Multilevel domain decomposition method for fast reprojection of images |
| US6307911B1 (en) | 1999-06-23 | 2001-10-23 | The Board Of Trustees Of The University Of Illinois | Fast hierarchical backprojection for 3D Radon transform |
| US6420194B1 (en) | 1999-10-12 | 2002-07-16 | Lucent Technologies Inc. | Method for extracting process determinant conditions from a plurality of process signals |
| US6741730B2 (en) * | 2001-08-10 | 2004-05-25 | Visiongate, Inc. | Method and apparatus for three-dimensional imaging in the fourier domain |
| US7209535B2 (en) * | 2003-06-20 | 2007-04-24 | Wisconsin Alumni Research Foundation | Fourier space tomographic image reconstruction method |
| JP5237812B2 (ja) * | 2005-09-06 | 2013-07-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 複数エネルギー窓をもつ計算機断層撮影におけるデータ処理および解析 |
| DE102005051620A1 (de) * | 2005-10-27 | 2007-05-03 | Siemens Ag | Verfahren zur Rekonstruktion einer tomographischen Darstellung eines Objektes |
| EP2049019A4 (de) * | 2006-08-03 | 2011-11-30 | Univ California | Iterative verfahren zur dosisreduzierung und bildvergrösserung in der tomographie |
| US8442353B2 (en) * | 2006-08-03 | 2013-05-14 | The Regents Of The University Of California | Incorporation of mathematical constraints in methods for dose reduction and image enhancement in tomography |
| JP2011510757A (ja) * | 2008-01-30 | 2011-04-07 | ザ リージェンツ オブ ザ ユニヴァーシティー オブ カリフォルニアＴｈｅ Ｒｅｇｅｎｔｓ Ｏｆ Ｔｈｅ Ｕｎｉｖｅｒｓｉｔｙ Ｏｆ Ｃａｌｉｆｏｒｎｉａ | 対象の周囲を動的な束縛として使用することによる、トモグラフィにおける放射線量低減および画質向上 |
| US8280135B2 (en) * | 2009-01-20 | 2012-10-02 | Mayo Foundation For Medical Education And Research | System and method for highly attenuating material artifact reduction in x-ray computed tomography |
| EP2424438A4 (de) * | 2009-04-30 | 2013-04-03 | Univ California | System und verfahren zur schnellen implementierung einer tomographie mit gleichneigung |
| KR102096410B1 (ko) * | 2014-05-02 | 2020-04-03 | 삼성전자주식회사 | 의료 영상 장치 및 그 제어 방법 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4205375A (en) * | 1977-01-31 | 1980-05-27 | Tokyo Shibaura Electric Co., Ltd. | Method and apparatus for performing computed tomography |
| US4222104A (en) * | 1978-11-02 | 1980-09-09 | E M I Limited | Radiography |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4433380A (en) * | 1975-11-25 | 1984-02-21 | Philips Medical Systems, Inc. | Tomographic scanner |
| GB1578442A (en) * | 1976-07-01 | 1980-11-05 | Emi Ltd | Radiography |
| JPS564257B2 (de) * | 1977-06-14 | 1981-01-29 | ||
| US4217641A (en) * | 1978-04-28 | 1980-08-12 | U.S. Philips Corporation | Correction for polychromatic X-ray distortion in CT images |
| JPS6058504B2 (de) * | 1980-11-17 | 1985-12-20 | Tokyo Shibaura Electric Co | |
| JPH0230066B2 (de) * | 1983-02-15 | 1990-07-04 | Tokyo Shibaura Electric Co |
-
1983
- 1983-06-07 US US06/501,832 patent/US4616318A/en not_active Expired - Lifetime
-
1984
- 1984-05-25 NL NL8401684A patent/NL192189C/nl not_active IP Right Cessation
- 1984-06-01 FR FR8408669A patent/FR2547439B1/fr not_active Expired - Fee Related
- 1984-06-01 DE DE19843420576 patent/DE3420576C2/de not_active Expired - Lifetime
- 1984-06-05 JP JP11635484A patent/JP2720894B2/ja not_active Expired - Lifetime
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4205375A (en) * | 1977-01-31 | 1980-05-27 | Tokyo Shibaura Electric Co., Ltd. | Method and apparatus for performing computed tomography |
| US4222104A (en) * | 1978-11-02 | 1980-09-09 | E M I Limited | Radiography |
Non-Patent Citations (1)
| Title |
|---|
| US-Z.: Radiology, Bd.117, Dez.1975, S. 561-572 * |
Also Published As
| Publication number | Publication date |
|---|---|
| NL192189B (nl) | 1996-11-01 |
| NL8401684A (nl) | 1985-01-02 |
| FR2547439A1 (fr) | 1984-12-14 |
| US4616318A (en) | 1986-10-07 |
| DE3420576A1 (de) | 1985-02-07 |
| JPS6041170A (en) | 1985-03-04 |
| NL192189C (nl) | 1997-03-04 |
| JP2720894B2 (ja) | 1998-03-04 |
| FR2547439B1 (fr) | 1991-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Basu et al. | O (n/sup 2/log/sub 2/n) filtered backprojection reconstruction algorithm for tomography | |
| DE69726479T2 (de) | Selektive Verfeinerung von Gitternetzen | |
| DE19851597B4 (de) | Verfahren und Gerät zur Erzeugung verbesserter Bilder von Kranzarterien | |
| Ogniewicz | Discrete voronoi skeletons | |
| DE19916354B4 (de) | Verfahren und Vorrichtung zur Analyse von Bildstrukturen | |
| DE60031050T2 (de) | Segmentation des Rands der Prostata in zwei- und dreidimensionalen Ultraschallbildern | |
| DE19916612B4 (de) | Verfahren und System zur Verbesserung diskreter Pixelbilder | |
| DE69628723T2 (de) | Gradientengestütztes Verfahren zur Erzeugung von Werten für unbekannte Pixel in einem Digitalbild | |
| Baltsavias | Multiphoto geometrically constrained matching | |
| DE602005004694T2 (de) | Verfahren und Vorrichtung für lokal adaptive Bildverarbeitungsfilter | |
| DE69733238T2 (de) | Verfahren und Gerät zur digitalen Substraktionsangiographie | |
| DE69733989T2 (de) | Verfahren zur Extrahierung von Daten und Verfahren zur Verdeckerung von Daten | |
| DE102009003387B4 (de) | Verfahren und System zur Bildrekonstruktion | |
| DE19854917B4 (de) | Verfahren zur Bildrekonstruktion für ein CT-Gerät | |
| DE10357206B4 (de) | Verfahren und Bildbearbeitungssystem zur Segmentierung von Schnittbilddaten | |
| DE2916486C2 (de) | ||
| DE69631225T2 (de) | Bildrekonstruktion von spiralabgetasteten Teilkegelstrahldaten | |
| DE69724906T2 (de) | Numerisches Bildverarbeitungsverfahren zur automatischen Extraktion bandförmiger Objekte | |
| DE69432672T2 (de) | Verfahren zum ausrichten und zum korrigieren des fehlers des kegelstrahls in einem tomographischen abtastsystem mit berechnetem divergentem strahl | |
| DE102004007637B4 (de) | Verfahren zum Erzeugen eines Bildes mit erhöhter Auflösung unter Verwendung einer Mehrzahl von Bildern mit niedriger Auflösung | |
| DE4016245C2 (de) | Verfahren zum Abtasten und Sammeln von Daten | |
| DE60030456T2 (de) | Verfahren und system zur verbesserung von digitalbildern | |
| EP0004258A1 (de) | Computer-Tomograph zur Durchstrahlung eines Körpers | |
| DE2804158C2 (de) | ||
| EP0989520A2 (de) | Computertomographie-Verfahren mit kegelförmigem Strahlenbündel |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 8110 | Request for examination paragraph 44 | ||
| 8125 | Change of the main classification |
Ipc: G06F 15/66 |
|
| D2 | Grant after examination | ||
| 8364 | No opposition during term of opposition | ||
| 8327 | Change in the person/name/address of the patent owner |
Owner name: PICKER INTERNATIONAL, INC., CLEVELAND, OHIO, US |