DE2806219A1 - Vorrichtung fuer die tomographie - Google Patents

Vorrichtung fuer die tomographie

Info

Publication number
DE2806219A1
DE2806219A1 DE19782806219 DE2806219A DE2806219A1 DE 2806219 A1 DE2806219 A1 DE 2806219A1 DE 19782806219 DE19782806219 DE 19782806219 DE 2806219 A DE2806219 A DE 2806219A DE 2806219 A1 DE2806219 A1 DE 2806219A1
Authority
DE
Germany
Prior art keywords
radiation
data
layer section
fourier transform
unit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
DE19782806219
Other languages
English (en)
Other versions
DE2806219C2 (de
Inventor
Tamon Inouye
Susumu Nishida
Kiyoto Saito
Tadatoki Yoshida
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.)
Toshiba Corp
Original Assignee
Tokyo Shibaura Electric Co Ltd
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 Tokyo Shibaura Electric Co Ltd filed Critical Tokyo Shibaura Electric Co Ltd
Publication of DE2806219A1 publication Critical patent/DE2806219A1/de
Application granted granted Critical
Publication of DE2806219C2 publication Critical patent/DE2806219C2/de
Expired legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/421Filtered back projection [FBP]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Algebra (AREA)
  • Pulmonology (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

Henkel, Kern, Feuer & HSnzel Patentanwälte
-4-
MShlstraße 37 TOKYO SHIBAURA ELECTRIC CO., LTD., D-8000München80
Kawasaki-shi, Japan Tel.: 089/982085-87
Telex: 0529802 hnkld Telegramme: ellipsoid
14. Feb. 1978
Vorrichtung für die Tomographie
Die Erfindung betrifft eine Vorrichtung für die Tomographie unter Verwendung von durchdringender Strahlung, mit einer Strahlungsquelle zur Abstrahlung einer durchdringenden Strahlung in fächerförmig aufgespreiztem Zustand in das Innere eines Scheiben- bzw. Schichtabschnitts eines Untersuchungsobjekts, mit einem Detektor mit Detektor- bzw. Meßelementen, um die durch den Schichtabschnitt hindurchgegangene Strahlung im wesentlichen in eine Anzahl von Strahlungsbündeln aufzuteilen, die jeweils unter einem Winkel auseinanderlaufen, und dabei die Intensität der Strahlungsbündel zu messen, mit einer Verschiebungseinheit zum Bewegen der Strahlungsbündel zwecks Abtastung des Schichtabschnitts, mit einer Datenverarbeitungseinheit zur Rekonstruktion eines Bilds des Schichtabschnitts auf der Grundlage der vom Detektor gelieferten Meßdaten und mit einer Anzeigeeinheit zur Wiedergabe des Bilds auf der Grundlage der Ergebnisse der Rekonstruktion.
Bl/eg -5-
809833/1070
Eine Vorrichtung für die Tomographie unter Verwendung durchdringender Strahlung, d.h. Röntgen- oder Gammastrahlung, ist bekannt. Typische Beispiele für eine derartige Vorrichtung sind verschiedene Arten von rechnergestützten Tomographievorrichtungen, bei denen das Bild eines Schichtabschnitts eines Untersuchungsobjekts durch Datenverarbeitung mittels eines elektronischen Rechners wiederaufgebaut und wiedergegeben wird. Bei einer solchen rechnergestützten Tomographievorrichtung sind auf den Übertragungsweg der Strahlung gestützte Projektionsdaten erforderlich, die in der Weise gewonnen werden können, daß das Strahlungsbündel aus allen Richtungen in das Untersuchungsobjekt projiziert wird. Für die Gewinnung des gewünschten Tomogramms müssen daher die Strahlungsquelle und mindestens ein Detektor zur Messung der Intensität der durch den Schichtabschnitt hindurchtretenden Strahlung über mindestens 180° um das Untersuchungsobjekt herumgedreht werden, um die erforderlichen Projektionsdaten zu gewinnen. Zur Durchführung der erwähnten Messung ist bei einer solchen Vorrichtung ein kreisförmiger Hohlraum vorgesehen, um den herum die Strahlungsquelle und der Detektor montiert sind, wobei zusätzlich ein Drehmechanismus vorgesehen sein muß, um den rechnergestützten Tomographieabtaster über mindestens 180 um das Untersuchungsobjekt herumzudrehen. Dabei befindet sich das Untersuchungsobjekt in einer Ruhelage im zentralen Hohlraum. Infolgedessen dürfen die Abmessungen des darzustellenden Schichtabschnitts nur vergleichsweise klein sein, damit sie nicht durch die Größe des Hohlraums und/oder des Untersuchungsobjekts beeinflußt werden. Für die Darstellung eines großen Schichtabschnitts muß bei einer bisherigen Vorrichtung ein sehr großer Drehmechanismus vorgesehen sein. Der Durchmesser des zentralen Hohlraums beträgt bis zu 80 cm bei einer Länge von etwa 25 bis 45 cm, so daß selbst dann, wenn der
-6-
809833/1070
zu untersuchende und wiederzugebende Schichtabschnitt vergleichsweise klein ist bzw. eine kurze Länge auf der Längsachse des Untersuchungsobjekts besitzt, die unmittelbare Darstellung dieses Abschnitts unmöglich sein kann.
Aufgabe der Erfindung ist damit die Schaffung einer Vorrichtung für die Tomographie unter Verwendung durchdringender Strahlung, mit deren Hilfe ohne weiteres ein großer Schichtabschnitt eines Untersuchungsobjekts bzw. ein Teilbereich dieses Schichtabschnitts unter Vermeidung der bei den bisherigen Vorrichtungen aufgetretenen Mangel dargestellt bzw. wiedergegeben werden kann.
Diese Aufgabe wird bei einer Vorrichtung der eingangs definierten Art erfindungsgemäß dadurch gelöst, daß die Verschiebungseinheit so ausgebildet ist, daß sie zumindest eine einmalige praktisch lineare Abtastung jedes Strahlungsbündels in einer Richtung quer zum Schichtabschnitt zuläßt, wobei diese lineare Abtastung durchführbar ist, während eine entsprechende von mehreren, im Bereich der fächerförmigen Ausbreitung der Strahlungsbündel liegenden Richtungen eingehalten wird.
Da bei dieser Vorrichtung die Strahlungsbündel bzw. die Durchdringungsstrahlen, die in fächerförmig aufgespreiztem Zustand von der Strahlungsquelle ausgehen, praktisch linear in Querrrichtung zum Schichtabschnitt verschoben werden, wird eine Linearabtastung dieses Abschnitts durchgeführt, während eine entsprechende Projektionsrichtung von mehreren Projektionsrichtungen eingehalten wird, die innerhalb des Bereichs der fächerförmigen Ausbreitung der Strahlen liegen. Wenn daher
-7-
809833/1070
der Ausbreitungswinkel der Strahlung mit OC/ bezeichnet wird, können die Projektionsdaten, die auf der Basis der vom Detektor gelieferten Meßdaten berechnet werden, praktisch dieselben sein wie die Projektionsdaten, die dann gewonnen werden, wenn der Abtaster mit einer Strahlungsquelle, die eine Vielzahl von prallel zueinander liegenden Einzelstrahlen für die Parallelabtastung projiziert, über einen Winkel von °0 ° um das Untersuchungsobjekt gedreht wird und die Intensitäten der Strahlungsbündel in einer entsprechenden bzw. zweckmäßigen Position im Verlauf der Drehbewegung gemessen werden. Der genannte Wert 06 kann entsprechend der Art und der Größe des Schichtabschnitts, der erforderlichen Schärfe der Bildwiedergabe und dergleichen bestimmt werden, wobei dieser Wert in der Praxis mit kleiner als 180° gewählt wird, vorzugsweise mit höchstens etwa 20 bis 30°.
Da bei der erfindungsgemäßen Vorrichtung ein wiederzugebendes Bild nur durch praktisch lineare Verschiebung der Strahlungsbündel gegenüber dem Untersuchungsobjekt abgenommen werden kann, läßt sich die Konstruktion des Abtasters vereinfachen, und die erforderlichen Projektionsdaten für die Wiedergabe eines Bilds der gewählten Schicht auf der Längsachse des Untersuchungsobjekts sowie eines Teilbereichs dieser Schicht lassen sich in kurzer Zeit gewinnen.
In bevorzugter Ausführungsform der Erfindung weist die Bewegungs- bzw. Antriebseinheit eine Zahnstange und ein in diese eingreifendes Ritzel auf, wobei die Zahnstange längs der Bewegungsrichtung eines Rahmens montiert ist, auf dem die Strahlungsquelle und der Detektor angeordnet sind. Infolge dieser Anordnung lassen sich die Strahlungsbündel und der Detektor mittels einesvsich bekannten, einfachen Mechanismus verschieben.
-8-809833/1070
In spezieller Ausführungsform der Erfindung wird als Strahlungsquelle ein Röntgenstrahier oder ein Röntgen-Schmalbündelgenerator verwendet, bei dem ein bleistiftförmiges bzw. schmales Röntgenstrahlungsbündel in Fächerform geschwenkt wird. Bei dieser Konstruktion kann das gewünschte Bild mittels einer an sich bekannten Strahlungsquelle gewonnen werden, wobei der Röntgen-Schmalbündelgenerator geeignet ist, die Strahlungsdosis zu verringern, welcher das üntersuchungsobjekt bzw. der menschliche Körper ausgesetzt ist.
In weiterer bevorzugter Ausführungsform der Erfindung umfaßt die Datenverarbeitungseinheit für den Aufbau bzw. die Wiederbildung des Bilds des Schichtabschnitts auf der Grundlage der vom Detektor übermittelten Meßdaten eine Projektionsdaten-Rechnereinheit zur Berechnung der Projektionsdaten bezüglich des übertragungswegs der Strahlungsbündel aus den vom Detektor gelieferten Meßdaten und eine Grunddaten-Recheneinheit mit einer ersten Rechen- bzw. Rechnereinheit zum Berechnen der eindimensionalen Fourierschen Transformierten der Projektionsdaten für die jeweiligen Projektionsrichtungen der Strahlungsbündel anhand der genannten Projektionsdaten und mit einer zweiten Rechnereinheit zur Berechnung der eindimensionalen Fourierschen Transformierten bezüglich der über den vollen Umfang des Schichtabschnitts gewählten Richtungen auf der Grundlage der vorgenannten, durch die erste Rechnereinheit berechneten Fourierschen Transformierten, welche die im Ausbreitungsbereich der Strahlungsbündel enthaltenen Strahlungsrichtungen betrifft. Weiterhin ist eine Bildaufbaueinheit mit einer dritten Rechnereinheit zur Berechnung der zweidimensionalen Fourierschen Transformierten der Absorptionskoeffizienten-Verteilung des Schichtabschnitts in Abhängigkeit von der eindimensionalen Fourierschen Transformierten, die von der Grunddaten-Verarbeitungseinheit übermittelt wird, und mit einer invertierenden Fourier-
-9-
809833/1070
sehen Transformationseinheit vorgesehen, welche die Absorptionskoeffizienten-Verteilung des zu untersuchenden Schichtabschnitts unter Zugrundelegung der zweidimensionalen invertierten Fourierschen Transformation für die genannte, durch die dritte Rechnereinheit berechnete zweidimensionale Transformierte berechnet. Unter Verwendung der vorstehend beschriebenen Datenverarbeitungseinheit kann die Absorptionskoeffizienten-Verteilung des Schichtabschnitts auf der Grundlage der Daten berechnet werden, die mit Hilfe des Detektors im Bereich des Winkels OO gewonnen wurden, um die berechneten Ergebnisse in der Anzeigeeinheit wiederzugeben.
Im folgenden sind bevorzugte Ausführungsformen der Erfindung anhand der beigefügten Zeichnung näher erläutert. Es zeigen:
Fig. 1 eine schematische Darstellung einer Röntgenstrahlung verwendenden Tomographievorrichtung gemäß der Erfindung,
Fig. 2 eine schematische Darstellung eines durch Linearverschiebung eines fächerförmigen Strahlungsbündels gemäß Fig. 1 erzielten Zustande,
Fig. 3 eine Darstellung zur Veranschaulichung eines Verfahrens für die Bildung eines Projektionsprofils mittels der zum Schichtabschnitt projizierten Röntgenstrahlung,
Fig. 4 eine graphische Darstellung der Periodizität des reellen Teils in einer anhand des Projektionsprofils berechneten eindimensionalen Fourierschen Transformierten, wobei der reelle Teil als Funktion des Projektionswinkels wiedergegeben ist,
-10-
809833/1070
Fig. 5 eine graphische Darstellung der Periodizität des imaginären Teils in der anhand des Projektionsprofils berechneten eindimensionalen Fourierschen Transformierten, wobei der imaginäre Teil als Funktion des Projektionswinkels dargestellt ist,
Fig. 6 ein Blockschaltbild zur Veranschaulichung des genauen Aufbaus der Datenverarbeitungseinheit,
Fig. 7 einen Schnitt durch einen Röntgen-Schmalbündelgenerator und
Fig. 8 eine schematische Darstellung einer Möglichkeit zur Erzeugung eines im wesentlichen fächerförmigen Röntgenstrahlungsbündels vom Schmalbündelgenerator .
Fig. 1 veranschaulicht schematisch Aufbau und Arbeitsweise einer Ausführungsform der Erfindung unter Verwendung von Röntgenstrahlung als durchdringende Strahlung.
Die Anordnung gemäß Fig. 1 umfaßt eine als Strahlungsquelle dienende Röntgenröhre 10, einen Strahlungsquellen-Kollimator 12 zur Umwandlung der von der Röhre 10 ausgestrahlten Röntgenstrahlung in ein fächerförmiges Röntgenstrahlungsbündel, ein mittels der erfindungsgemäßen Vorrichtung zu untersuchendes Objekt 16, z.B. einen menschlichen Körper, einen der Röntgenröhre 10 gegenüberliegend angeordneten Röntgenstrahlungs-Detektor 18, wobei sich das Untersuchungsobjekt 16 zwischen Detektor 18 und Röhre 10 befindet, einen die Röhre 10 und den Detektor 18 tragenden und längs des Untersuchungsobjekts verschiebbaren Rahmen 20 jeweils an einem von mehreren, den
-11-809833/1070
Detektor 18 bildenden Detektorelementen 24 angebrachte Detektor-Kollimatoren 22, eine das Untersuchungsobjekt 16 unterstützende Auflage bzw. Liege 26, eine Verschiebungseinheit 28 zur Linearverschiebung des Rahmens 20 parallel zu dem zu untersuchenden Schichtabschnitt, und eine Datenabnahmeeinheit 30 zum Sammeln der von den Detektorelementen 24 übermittelten Meßdaten.
Eine Röntgenstrahlungsquellen-Regeleinheit 32 zur Regelung bzw. Steuerung der Arbeitsweise der Röntgenröhre 10 umfaßt einen Hochspannungsgenerator 34 zur Ansteuerung der Röhre 10 und eine Hochspannungsquellen-Steuereinheit 36 zur Steuerung der Röntgenstrahlungserzeugung. Die durch eine Rahmenantriebs-Steuereinheit 38 gesteuerte Verschiebungseinheit 28 vermag den Rahmen 20 quer zur Zeichnungsebene zu verschieben. Bei der dargestellten Ausführungsform ist der Rahmen 20 mittels einer an seiner Unterseite befestigten Zahnstange 40 und eines Ritzels 32 verschiebbar, das mit der Zahnstange 40 kämmt und durch die Steuereinheit 38 angetrieben wird. Bei der Verschiebung des Rahmens 20 sammelt die Datenabnahmeeinheit 30 alle Daten bezüglich der Intensität der das Untersuchungsobjekt durchdringenden Röntgenstrahlung, die von den einzelnen Detektor- bzw. Meßelementen 24 in vorbestimmten Positionen gewonnen werden, um diese Daten sodann zu einer Datenverarbeitungseinheit 44 zu übermitteln. Die Strahlungsquellen-Steuereinheit 32, die Rahmenantriebs-Steuereinheit 38 und die Datenverarbeitungseinheit 44 werden in Abhängigkeit von Zeitsteuer- bzw. Taktsignalen betätigt, die von einer Taktsteuereinheit 46 abgegeben werden. Die Datenverarbeitungseinheit 44 berechnet die Daten bezüglich der Strahlungsintensität, die nach entsprechender Taktsteuerung gewonnen und von der Datenabnahmeeinheit 30 übermittelt werden, um Projektionsdaten bezüglich ausgewählter Projektionsrichtungen innerhalb des Ausbreitungswinkelbereichs des fächerförmigen Strahlungsbündels 14 abzuleiten. Diese Einheit
-12-
809833/1070
44 wird dabei so betrieben, daß sie die Grunddaten, d.h. die eindimensionale Fouriersehe Transformierte auf der Grundlage der genannten Projektionsdaten bezüglich praktisch aller Richtungen berechnet, die für den Aufbau des Bilds der zu untersuchenden Schicht erforderlich sind und weiterhin das Bild aus den genannten Grunddaten aufbaut. Das Bild kann an der Anzeigeeinheit 48 entsprechend den Daten für das wiedergebildete Bild dargestellt werden. Die Projektionsdaten und Grunddaten werden später noch näher erläutert werden.
Fig. 2 veranschaulicht einen Zustand, in welchem sich das von der Röhre 10 emittierte, fächerförmige Röntgenstrahlungsbündel 14 quer zu seiner Fächerebene und parallel zur Zeichnungsebene bewegt. Zur Vereinfachung der Darstellung sind in Fig. 2 nur drei Positionen 14a, 14b und 14c des Fächerbündels dargestellt.
Gemäß Fig. 1 sind zahlreiche Detektor- bzw. Meßelemente 24 auf einer Kreisbogenlinie um die Röhre 10 herum und letzterer gegenüberliegend angeordnet, wobei jedes Element 24 einen vorbestimmten Teil des Fächerbündels 14 in den verschiedenen Projektionsrichtungen empfängt. Ersichtlicherweise besteht somit das Fächerbündel 14 praktisch aus Röntgenstrahlungsbündeln, deren Zahl derjenigen der Detektor- bzw. Meßelemente entspricht, die im Bereich des Winkels OC um die Röhre 10 herum angeordnet sind. In Fig.1 ist beispielsweise ein Röntgenstrahl 15 dargestellt, welcher auf das Meßelement 24a in einer Position fällt, die unter einem Winkel φ zu einer lotrechten, senkrecht zur Auflage 26 verlaufenden Linie 50 liegt. Der Röntgenstrahl 15 ist in den fächerförmigen Röntgenstrahlungsbündeln 14a bis 14c gemäß Fig. 2 mit 15a, 15b bzw. 15c bezeichnet. Aus Fig. 2 geht ohne weiteres hervor, daß bei Bewegung des fächerförmigen Röntgenstrahlungsbündels durch Verschiebung des Rahmens 20 das in Fig. 2 nicht dargestellte, auf der Auflage bzw. Liege ruhende Untersuchungsobjekt 16 durch die
809833/1070 ~13~
Röntgenstrahlen 15, die unter dem Winkel φ zur lotrechten Linie 50 (Fig. 1) liegen, parallel abgetastet wird. Alle anderen Röntgenstrahlen, die von den anderen Meßelementen 24 empfangen werden, tasten zusammen mit dem genannten Röntgenstrahl 15 das Untersuchungsobjekt parallel ab, so daß dann, wenn die Röntgenstrahlen mit den verschiedenen Projektionsrichtungen in den verschiedenen Stellungen nur während einer praktisch linearen Bewegung des Rahmens 20 abgegriffen werden, die Daten bezüglich der Intensitäten der Röntgenstrahlen bzw. -bündel gewonnen werden können, ebenso wie dies der Fall ist, wenn der parallel zur Zeichnungsebene liegende Schichtabschnitt unter Drehung durch eine Anzahl paralleler Röntgenstrahlenbündel unter dem Winkel OU in der Zeichnungsebene abgetastet wird.
Die folgende Beschreibung bezieht sich auf die Grundtheorie sowie auf die Einheiten bzw. Geräte für die Wiederbildung des Schichtabschnitts des Untersuchungsobjekts 16 anhand der Meßdaten, die von der Datenabnahmeeinheit 30 zur Datenverarbeitungseinheit 44 übertragen werden.
Die Meßdaten von einem Detektor- bzw. Meßelement 24 geben die Intensität eines Röntgenstrahlungsbündels an, welche durch den integralen Wert des Röntgenstrahlungs-Absorptionskoeffizienten in jedem Übertragungsweg bestimmt wird. Bei der erfindungsgemäßen Vorrichtung werden die durch den integralen Wert bzw. die integrale Größe dieses Absorptionskoeffizienten bestimmten Projektionsdaten anfänglich auf der Grundlage der Meßdaten berechnet. Anschließend wird die Absorptionskoeffizienten-Verteilung im zu untersuchenden Schichtabschnitt, d.h. das Bild dieses Schichtabschnitts, berechnet bzw. wiederhergestellt. Wenn die Absorptionskoeffizienten-Verteilung nach dem beschriebenen Verfahren berechnet werden kann, läßt sich das Bild auf der Anzeigeeinheit auf der Grundlage dieser Verteilung in an sich bekannter Weise sichtbar wiedergeben. Bei der dargestellten Ausführungsform umfaßt die Anzeigeeinheit einen
809833/107 0
-14-
Schreiber zum Ausdrucken der Daten des aufgebauten bzw. wiedergebildeten Bilds.
Für die Bildrekonstruktion sind bereits verschiedene Verfahren vorgeschlagen worden, beispielsweise das Umkehrverfahren, das Filter-Rückprojektionsverfahren und das Fouriersche Transformations- bzw. Umsetzverfahren. Die genannten Verfahren unterscheiden sich bezüglich der Wiedergabe oder bezüglich der praktischen Datenverarbeitung, während sie andererseits im arithmetischen Sinne nicht voneinander verschieden sind, weshalb im folgenden ein Ausführungsbeispiel unter Anwendung des Fourierschen Umsetz- bzw. TransformationsVerfahrens beschrieben ist.
Nachstehend ist eine Theorie des Fourierschen Transformations-Verfahrens anhand von Fig. 3 erläutert. Fig. 3 veranschaulicht eine zu rekonstruierende Schnittschicht bzw. einen Schichtabschnitt 60 des Untersuchungsobjekts, in einer Lage parallel zur Zeichnungsebene. Auf der Oberfläche des Schichtabschnitts 60 ist ein rechtwinkeliges Koordinatensystem (x, y) mit einem Ursprungsort 0 auf dem Schichtabschnitt festgelegt. Bei 62 ist ein durch den Schichtabschnitt hindurchgehender Röntgenstrahl angedeutet. Weiterhin ist in Fig. 3 ein rechtwinkeliges Koordinatensystem (X, Y) mit dem gemeinsamen Ursprungsort O dargestellt, dessen X-Achse senkrecht zur übertragungsrichtung des Röntgenstrahl 62 verläuft und einen Winkel θ mit der x-Achse bildet. Die Y-Achse erstreckt sich daher in dieselbe Richtung wie der Röntgenstrahl. Da θ mittels der Projektionsrichtung des Strahls 62 variiert werden kann, wird die Größe von θ für den Röntgenstrahl 62 zur Vereinfachung der Beschreibung auf ©o festgelegt- Der lotrechte Abstand von einer Position, in welcher sich der Röntgenstrahl 62 über die X-Achse bewegt, d.h. der Abstand vom Ursprungsort O zum Röntgenstrahl 62, wird
-15-
809833/1070
mit Xq festgelegt. Der Ausdruck f(χ, y) in Fig. 3 bedeutet einen Absorptionskoeffizienten an einer wahlfreien Stelle (x, y) auf der Schichtebene 60.
Der auf dem übertragungsweg des Röntgenstrahls 62, welcher den Schichtabschnitt 60 durchdrungen hat, gezogene Balken g(X0, ©0) gibt die Projektionsdaten an, die anhand der vom nicht dargestellten Meßelement gelieferten Meßdaten bezüglich des Übertragungswegs 64 des Röntgenstrahls 62 berechnet wurden.
In Fig. 3 ist nur eine durch die beiden Werte bzw. Größen von θ~ bzw. X0 bestimmte Projektionsdateneinheit g(XQ, θ0) für einen Röntgenstrahl 62 veranschaulicht. Bei der praktisch linearen Verschiebung des Rahmens 20 gemäß Fig. 1 wird jedoch eine Vielzahl von Projektionsdaten für den übertragungsweg parallel zum Röntgenstrahl 62 gewonnen. Gemäß Fig. 3 kann die Größe der Projektionsdateneinheit anhand der Funktion von θ und X ermittelt werden. Infolgedessen wird diese Einheit auf übliche Weise als g(X, Θ) ausgedrückt. Die jeweiligen Größen müssen als graphische Balken zusammen mit der genannten Größe g(X(w 6Q) in dem eine Querachse X und eine lotrechte Achse g(X, Θ) umfassenden Koordinatensystem angegeben werden, doch sind zur Vereinfachung der Darstellung die einzelnen Balken für g(X, Θ) weggelassen. Vielmehr veranschaulicht Fig. 3 lediglich eine Kurve, die durch Verbindung der oberen Enden der einzelnen Balken erhalten wird, d.h. das Projektionsprofil 66 des Schichtabschnitts 60, im folgenden auch einfach als "Projektion des Schichtabschnittsir bezeichnet.
Die eindimensionale Fourier-Transformierte G( W , Θ) bezüglich der Achse X des erwähnten Projektionsprofils 66 berechnet sich nach folgender Gleichung:
G(ü), Θ) « /_"g(X. 8)e~iüjXdX (1)
-16-809833/1070
In dieser Gleichung bedeuten g(X, Θ) das erwähnte Projektionsprofil 66 und it) eine Raumfrequenz, die üblicherweise für die Ableitung der Fourier-Transformierten eingeführt wird. Es sei angenommen, daß die Absorptionskoeffizienten-Verteilung f(x, y) der Schichtabschnittebene 60 bekannt ist. Eine zweidimensionale Fourier-Transformierte F(^ '7J) für x un<^ Y wird anhand der Daten von f(x, y) berechnet. Die Berechnung von F(^ , Tj ) auf der Grundlage von f(x, y) erfolgt nach folgender Gleichung:
n)
O " fix, y)e-iax+X]y)dxäy (2)
in welcher £ und Tj (jeweils) eine Raumfrequenz bedeuten und neue Veränderliche darstellen, die zur Ableitung der zweidimensionalen Fourier-Transformierten eingeführt werden.
Anschließend wird F(^ ,n ) in F( W , Θ), als Polarkoordinate dargestellt, unter Benutzung der folgenden Gleichung (3) für und 71 von F(^ , η ) umgewandelt:
= cos©
T) = sine (3)
Auf der Grundlage eines Lehrsatzes dahingehend, daß die eindimensionale Fourier-Transformierte der Projektion (Projektionsprofil) g(X, Θ) eines bestimmten Rekonstruktionsbilds f(x, y) gleich dem Zentralschnitt ist, d.h. dem Schnitt mit dem Ursprungspunkt, der durch Schneiden der zweidimensionalen Fourierschen Transformierten dieses Bilds f(x, y) unter einem entsprechenden Winkel erhalten wird, läßt sich die folgende Gleichung (4) ableiten:
-17-
809833/1070
F(w cos9, ω sine) = G(üj, θ) (4)
Der genannten Zentralschnitt ist ein Schnitt durch den Ursprungsort des Frequenzbereichs, d.h. eine Position, die der Gleichung (JU=O genügt.
Wenn daher der nach Gleichung (1) erhaltene Ausdruck G(U/, θ) unter Heranziehung von Gleichung (3) in das Koordinatensystem ( 4» r 77 ) umgewandelt werden kann, ist er gleich dem Ausdruck F (c , Tj ) , der durch Umwandlung von F ( Cü , Θ) in das Koordinatensystem ( £ ,η ) erhalten wird, wie dies aus Gleichung (4) hervorgeht.
Aus den obigen Gleichungen geht hervor, daß G( OJ, Θ), d.h. F( ζ f 77 ) , für die verschiedenen Werte bzw. Größen von θ , die entsprechend in einem Bereich, welcher den gesamten Umfang von 360° des Untersuchungsobjekts 16 umfaßt, und die nach einem noch näher zu beschreibenden Verfahren berechnet werden, auf der Grundlage der verschiedenen Projektionsdaten g(X, Θ) berechnet werden kann, die anhand der von allen Detektor- bzw. Meßelementen 24 und durch linearen Antrieb zur Verschiebung des fächerförmigen Röntgenstrahlungsbündels 14 gemäß Fig. 1 und 2 gewonnenen Intensitätsdaten berechnet werden.
Der Absorptionskoeffizient und folglich die Absorptionskoeffizienten-Verteilung an jedem Punkt (x, y) des Schichtabschnitts wird dadurch abgeleitet, daß die zweiüimensionale invertierte Transformation auf F ( £ , η ) nach folgender Gleichung angewandt wird:
y) -
4ir
809833/1070
-18-
Aus den obigen Gleichungen geht hervor, daß bei dem beschriebenen Bildrekonstruktionsverfahren die Intensitätsdaten für das Röntgenstrahlungsbündel bei den verschiedenen Werten bzw. Größen von ©,über den vollen Umfang von 360 herausgegriffen, unbedingt notwendig sind, während jedoch die gewünschten bzw. angestrebten Intensitätsdaten über 360°, d.h. das Projektionsprofil 66 (Fig. 3), durch lineare Verschiebung der Röntgenröhre 10 und der Meßelemente 24 auf beschriebene Weise nicht gewonnen werden können. In diesem Fall können nur die im Bereich des Ausbreitungs- bzw. Aufspreizwinkels ÖO des Fächerstrahls 14 gemessenen Intensitätsdaten erhalten werden.
Bei der erfindungsgemäßen Vorrichtung werden jedoch die Projektionsdaten g(X, Θ) und die eindimensionale Fourier-Transformierte G(GJ, Θ) (abhängig von Gleichung (1)) innerhalb des genannten Winkelbereichs OCgemessen und wiederum auf der Grundlage der Intensitätsdaten bezüglich θ im obigen Bereich des Winkels OU berechnet, worauf die Gleichung zur Berechnung von G( U)1 Q) für die vollen bzw. alle Richtungen auf der Grundlage von G( £J , Θ) für die Größe von θ innerhalb des Bereichs des Winkels θ bzw. Oc/ bestimmt wird. Die Größe G( Ul , Θ) bezüglich aller Richtungen (360°) wird nach der obigen Gleichung berechnet, um F( ξ , 71 ) entsprechend G(U,', Θ) zu berechnen, worauf schließlich das Bild des Schichtabschnitts auf der Grundlage von Gleichung (5) rekonstruiert und wiedergegeben wird.
Die obige Gleichung zur Berechnung von G( U, Θ) für alle Richtungen bzw. den Gesamtumfang kann unter Heranziehung des Umstands aufgestellt werden, daß G( U, Θ) eine periodische Funktion von θ ist. Diese Tatsache ist nachstehend noch näher erläutert.
-19-
809833/1070
Gleichung (2) kann dann zunächst wie folgt umgeschrieben werden:
, η) = /."/." f(xr y)[σο8(ξχ+ηγ)-1 sin(ζχ+ny)]dxdy
Wie aus Gleichung (6) hervorgeht, umfaßt die Funktion F(^ ,Tf) den reellen Teil ReF (£ , η ) sowie den imaginären Teil ImF(^, η ) .
Dementsprechend gilt folgendes:
ReF(ξ,η) = /_«/_" f(x,y)cos(Cx+ny)dxdy
-ce* —oo
00 „ 00
ImF(C»n) =-/_„/.„ f(x,y)sin(ξχ+ηγ)dxdy (7)
ReF(-C,-n) = /_^/_" f(x,y)cos(-ξχ-ny)dxdy = /_*/„" f(x,y)cos(ξχ+ηγ)dxdy
ξ,η) (8)
Ebenso gilt:
-/_"/_".f(x,y) sin (-ξχ-ηγ)αχαγ fJ^fJZ f (x,y) sin (ξχ+ηγ) dxdy
C,η) (9)
-20-
809833/1070
Andererseits werden die folgenden Gleichungen von Gleichung (3) abgeleitet:
-ξ = -ω cos9 = ω cos(e+ir) -η - -ω sin9 = ω sin(9+Tr)
Infolgedessen können Gleichungen (8) und (9) mit den Größen <aJ und θ wie folgt ausgedrückt werden:
ReF(U cO8(9+tt), ω sin(0+Tr))
= ReF(ü) cos9, ω sin9)
ImF(ω cos(0+u)f ω 3±η(θ+π))
β -ImG(ω COS0, ω sin6)
Der Ausdruck bzw. die Größe G(OJ , Θ) kann auch mit dem reellen Teil ReG(CJ, Θ) und dem imaginären Teil ImG(OJ, Θ) ausgedrückt werden, so daß sich anhand von Gleichung (4) die folgenden Gleichungen (10) und (11) wie folgt ausdrücken lassen:
ReG(ω, θ+ιτ) * ReF(io cos(9+ir), ω sin(0+n)) « ReF(u CO89, ω sin9) - ReG(u),Θ) (10)
Auf ähnliche Weise qilt:
ImG(ü), θ+π) ■ ImF(u) cos(0+n), ω ein(e+n)) - -ImF(ω cose, ω βΐηβ)
Ο) (11)
-21-
809833/1070
Ersichtlicherweise besitzt somit die Größe für die Gleichheit W von G(U , Θ) die Periodizität bezüglich der Größe Θ. Wie aus Gleichung (10) hervorgeht, ist der reelle Teil ReG{ (J , Θ) eine periodische Funktion von Θ, d.h. 'If , und wie ebenfalls aus Gleichung (11) hervorgeht, ist ImG ( t»J , Θ) eine periodische Funktion von Θ, dessen Vorzeichen bei jeder Periode It* umgekehrt und durch eine Änderung entsprechend 2 Tf* auf das ursprüngliche Vorzeichen zurückgeführt werden kann. Fig. 4 veranschaulicht die Beziehung zwischen ReG(Co, Θ) und Θ. Fig. 5 verdeutlicht dagegen die Beziehung zwischen ImG ( U) , Θ) und Θ.
Da die Funktion G(CJ, Θ) die erwähnte Periodizität besitzt, können ReG(UI, Θ) und ImG(CJ, Θ) jeweils in Form einer unendlichen Reihe entwickelt werden:
oo
ReG(ü),9) = Σ [An(u)cos 2ηθ + Bn(w)sin 2ηθ] (12)
n=0
OO
ImG(W,θ) =« Σ [Cn(ü))cos(2n+1)9 + Dn (ω) sin (2n+l) θ] ...(13) n=0
ReG ( CJ, θ) und ImG(CJ, θ) lassen sich mit ausreichender Genauigkeit berechnen, indem die Ausdrücke addiert bzw. summiert werden, bis eine Ausdruckszahl η eine im Hinblick auf die gewünschte Genauigkeit entsprechend gewählte Größe erreicht, weil die Koeffizienten An ((J ) , Bn ((J ) , Cn ((J ) und Dn ( (J ) mit sich vergrößernder Ausdruckszahl η (Gliedzahl) an Null angenähert werden können.
Die folgende Gleichung (14) wird anstelle von Gleichungen (12) und (13) benutzt.
-22-
809833/1070
ReG((is,θ) =f Σ [An(u)cos 2ηθ + Βη(ω)βχη 2ηθ] η=0
ΙπΚ5(ω,θ) « Σ [Cn((o)cos(2n+l)e + Dn(üj)sin(2n+1) θ] ...(14) η=0
Zur Auswertung bzw. Berechnung der erwähnten Koeffizienten wird eine Zahl K definiert, die nach folgender Gleichung gewählt wird:
K > 2(N+1)
Sodann werden die Four ier-Tr ans formierten G ( u/, Θ-) , G(CJfG-)
G(Cd, 9 ) der Projektionsprofile g(CJ, Q1), g((J, ©2)
g(6J, ÖjJ bezüglich der Winkel θ.. , θ» ... θ, auf der Grundlage von Gleichung (1) berechnet. Die Winkel Q1, Θ-, ... Θ, werden
gewählt ' ^ K im Winkelbereich von 0 < θ \ 0^<% Infolgedessen werden die Größen O1, Θ« ... Θ, ihrerseits an der rechten Seite von Gleichung (14) eingesetzt, während die Größen des reellen Teils ReG (6J, Q^), ReG(UJ, Q^i ... ReG(CJ, 9fc) sowie die imaginären Teile ImG((J, θ^, ImG(CJ, θ2) ... ImG(cJ, efc) von G(CJ, θ.), G(Cj, θ2) ... G((J, ek) für die verschiedenen Größen von θ wiederum an der linken Seite dieser Gleichung eingesetzt werden, um lineare Simultangleichungen für die Koeffizienten An((J) bis Dn((J) zu erhalten, so daß die verschiedenen, vorstehend angegebenen Koeffizienten durch Lösung der angegebenen Gleichungen berechnet werden können. Nach dem Einsetzen der Größen für die Koeffizienten in Gleichung (14) kann folglich die Größe bzw. der Ausdruck G (td, Θ) für die gewünschten Werte UJ und θ in Abhängigkeit von Gleichung (14) berechnet werden. Jede so berechnete Größe G(CJ, θ) wird für
-23-
809833/1070
Grunddaten für die Wiederherstellung bzw. Rekonstruktion des Bilds des untersuchten Schichtabschnitts benutzt. Die Bauteile der Datenverarbeitungseinheit 44 gemäß Fig. 1 umfassen eine Projektionsdaten-Rechnereinheit 70, eine Grunddaten-Rechnereinheit 72 und eine Bildrekonstruktionseinheit 74. Die Datenverarbeitungseinheit 44 umfaßt die für den einwandfreien Betrieb der genannten Einheiten 70, 72 und 74 erforderlichen Schaltkreise, etwa Speicher, Torschaltungen und dergleichen, doch sind die zusätzlichen bzw. Hilfseinrichtungen oder Steuereinrichtungen in Fig. 1 mit Rücksicht auf die Darstellung der wesentlichen Grundeinheiten 70, 72 und 74 weggelassen.
Die Projektionsdaten-Rechnereinheit 70 dient zur Aufnahme der Meßdaten von der Datenabnahmeeinheit 30 und zur Berechnung der Projektionsdaten G(X, Θ) entsprechend den Meßdaten. Die Grunddaten-Rechnereinheit 72 empfängt die Projektionsdaten g(X, Θ) zur Klassifizierung dieser Daten entsprechend jeder Projektionsrichtung, d.h. entsprechend jedem Projektionswinkel, und sie berechnet die Fourier-Transformierte G( U), Θ) entsprechend dem Winkel θ auf der Grundlage der Projektionsdaten. Die Grunddaten-Rechnereinheit 72 berechnet die Größe von G(CJ, Θ) für jeden Winkel θ in allen Richtungen nach dem oben beschriebenen Verfahren auf der Grundlage der durch die tatsächliche Messung gewonnenen Daten der Fourier-Transformierten G(UJ, Θ). Die Bildrekonstruktionseinheit berechnet F ( ξ ,77 ) anhand der durch tatsächliche Messung und nach obiger Berechnung erhaltenen Transformierten G(CJ, θ), und sie berechnet sodann die zweidimensionale inverse bzw. invertierte Fourier-Transformierte von F (£, 77 ) zur Wiederherstellung von F(x, y ), d.h. des Bilds des Scheiben- bzw. Schichtabschnitts des Untersuchungsobjekts.
-24-
809833/1070
44 In Fig. 6 ist die DatenverarbeitungseinheitVnäher dargestellt. Gemäß Fig. 6 umfaßt die Grunddaten-Rechnereinheit 72 eine erste Rechnereinheit 72a zur Berechnung einer ersten eindimensionalen Fourierschen Transformierten anhand der durch tatsächliche Messung gewonnenen Projektionsdaten sowie eine zweite Rechnereinheit 72b zur Berechnung einer zweiten eindimensionalen Fourierschen Transformierten, d.h. der oben genannten Grunddaten, in Abhängigkeit von der eindimensionalen Fourier-Transformierten, d.h. einer ersten eindimensionalen Transformierten, auf der Grundlage, daß die erste eindimensionale Fourier-Transformierte eine Periodizität von 360° zur Projektionsrichtung besitzt. Fig. 6 veranschaulicht auch die Bildrekonstruktionseinheit 74 mit einer dritten Rechnereinheit 74a zur Berechnung der zweidimensionalen Fourier-Transformierten der Absorptionskoeffizienten-Verteilung des Schichtabschnitts 60 anhand der ersten und der zweiten eindimensionalen Fourier-Transformierten, welche durch die Grunddaten-Rechnereinheit 72 berechnet worden sind, und einer zweidimensionalen Fourierschen Inversionstransformationseinheit 74b zur Berechnung der Absorptionskoeffizienten-Verteilung des Schichtabschnitts 60 anhand der zweidimensionalen Fourier-Transformierten .
Bei der vorstehend beschriebenen Ausführungsform wird die von der Röntgenstrahlungsquelle 10 emittierte Röntgenstrahlung durch den Strahlungsquellen-Kollimator in ein fächerförmiges Röntgenstrahlungsbündel 14 umgewandelt, das auf das Untersuchungsobjekt bzw. den menschlichen Körper 16 projiziert wird. Ein praktisch fächerförmiges Röntgenstrahlungsbündel kann auch auf andere Weise erhalten werden, beispielsweise mittels eines Röntgen-Schmalbündelgenerators (auch als Röntgenmikrostrahl-Generator bezeichnet), in^dem das Strahlungsbündel fächerförmig verschwenkt wird (auch als Röntgenmikrostrahl- bzw. -bündel bezeichnet).
-25-
809833/1070
Der erwähnte Röntgen-Schmalbündelgenerator 80 ist in Fig. 7 veranschaulicht. Ein von einem Elektronenrohr 82 emittierter Elektronenstrahl EB kann durch das von einer an eine Anode 84 angelegten Hochspannung erzeugte elektrische Feld beschleunigt und dann durch eine Richtspule 86, eine Fokussierspule 88 und andere Elemente in eine vorbestimmte Form gebracht werden, um zu einer Ablenklinse bzw. Ablenkspule 90 übertragen zu werden. Die Ablenkspannung, die unter der Steuerung durch einen nicht dargestellten elektronischen Rechner nach einem vorbestimmten Programm variiert wird, wird an die Ablenkspule 90 angelegt, wobei durch diese Spannung ein Auftreffpunkt G des Elektronenstrahls EB auf einer Fangelektrode 92 bestimmt wird. Die Fangelektrode 92 besteht aus einem für Röntgenstrahlenerzeugung geeigneten Schwermetall. Die Fangelektrode 92 ist schräg zur Achslinie der praktisch miteinander fluchtenden Elemente vom Elektronenrohr 82 zur Ablenkspule 90 ausgerichtet. Beim Auftreffen des Elektronenstrahls EB auf die Fangelektrode wird Röntgenstrahlung vom Auftreffpunkt G emittiert, wobei ein Teil der emittierten Röntgenstrahlung durch ein in einem Gehäuse 94 vorgesehenes Feinloch PH als Röntgen-Schmalbündel bzw. -Bleistiftstrahl PB nach außen abgeleitet und auf das Untersuchungsobjekt 16 gerichtet wird.
Fig. 8 ist eine schematische Darstellung eines Systems zur Erzeugung eines praktisch fächerförmigen Röntgenstrahlungsbündels bei Verwendung des Schmalbündelgenerators 80. In Fig. 8 sind nur die für die Erläuterung dieses Systems erforderlichen Teile sowie das Untersuchungsobjekt 16 dargestellt, wobei die Fangelektrode 92 als Seitenansicht von der rechten Seite des Generators 80 gemäß Fig. 7 her veranschaulicht ist. Die Fangelektrode 92 nimmt eine schräge Lage zum Elektronenrohr 82 ein. Wenn der Elektronenstrahl EB auf den Punkt G gemäß Fig. 7 und 8 geworfen wird, wird ein Teil der Röntgenstrahlung längs der Linie G-Ga geworfen und in Form eines Bleistiftstrahls bzw.
-26-
809833/1070
Schmalbündels XG auf das Untersuchungsobjekt 16 projiziert. Wenn der Auftreffpunkt durch Verschiebung bzw. Richtung des Elektronenstrahls EB nach links durch entsprechende Ansteuerung der Ablenklinse bzw. -spule 90 gemäß Fig. 8 auf einen Punkt F verschoben wird, wird ein Schmalbündel X„
erzielt, das sich in Richtung F-Fa durch das Feinloch PH bewegt. Wird andererseits der Auftreffpunkt auf einen Punkt H verschoben, so bewegt sich das Schmalbündel X„ in der Richtung H-Ha. Wenn daher der Auftreffpunkt wiederholt im Bereich von F bis H verschoben wird, verlagert sich das Schmalbündel fortlaufend im Bereich des Winkels S , wodurch sich dieselbe Wirkung ergibt wie bei der Bestrahlung des Untersuchungsobjekts 16 mit einem fächerförmigen Röntgenstrahlungsbündel. Infolgedessen kann auf die Strahlungsquelle 10 und auf den Kollimator 12 gemäß Fig. 1 verzichtet werden, und wenn der Schmalbündelgenerator 80 so angeordnet ist, daß die Schwenkebene des Schmalbündels mit der Zeichnungsebene übereinstimmt, kann eine Wiedergabe des Scheiben- bzw. Schichtabschnitts des Untersuchungsobjekts 16 auf die in Fig. 1 bis 6 veranschaulichte Weise erzielt werden.
Bei Verwendung des Röntgen-Schmalbündelgenerators kann die Röntgenstrahlungsdosis, welcher der menschliche Körper ausgesetzt ist, ganz erheblich verringert werden.
Bei der vorstehend beschriebenen Ausführungsform wird Röntgenstrahlung als durchdringende Strahlung benutzt. Die Strahlung ist jedoch nicht auf Röntgenstrahlung beschränkt, da gegebenenfalls auch Gammastrahlung anwendbar ist.
Mit den vorstehend beschriebenen Ausführungsformen der Erfindung läßt sich also beispielsweise das Tomogramm eines zu untersuchenden menschlichen Körpers einfach durch Linearverschiebung des Strahlungsbündels ohne Drehung desselben um den
-27-
809833/1070
menschlichen Körper bzw. das Untersuchungsobjekt herum erzielen. Infolgedessen kann durch einfache Linearbewegung des Strahlungsbündels bzw. der Strahlungsbündel auf vorher erläuterte Weise ohne weiteres die Bildrekonstruktion und -wiedergabe eines großen Scheiben- bzw. Schichtabschnitts oder eines langen Schichtabschnitts durchgeführt werden, beispielsweise eines in Axial- oder Längsrichtung des Untersuchungsobjekts verlaufenden Gesamtschichtabschnitts oder eines Teilbereichs desselben. Die Bildrekonstruktion und -wiedergabe eines derart großen bzw. langen Schichtabschnitts ist bei den bisherigen Tomographie-Vorrichtungen nicht möglich. Mit Hilfe der erfindungsgemäßen Vorrichtung kann ein dreidimensionales Bild (size) eines Abbildungspunkts im Patienten auf der Grundlage der Daten hergestellt werden, die von verschiedenen Tomogrammen praktisch senkrecht zur Axialrichtung des Untersuchungsobjekts und verschiedenen anderen, praktisch parallel zur Axialrichtung des Untersuchungsobjekts geführten Tomogrammen gewonnen wurden. Anhand der Größe des gemessenen Abbildungspunkts kann der untersuchende Mediziner ein genaues Bild der betreffenden Stelle gewinnen und entsprechende Maßnahmen zur richtigen Behandlung einleiten. Erfindungsgemäß läßt sich weiterhin das gewünschte Tomogrcmm durch einfache Linearverschiebung des Abtasters herstellen, so daß der Aufbau des Abtasters im Vergleich zur bisherigen Vorrichtung vereinfacht werden kann. Hierdurch werden hohe Zuverlässigkeit und geringe Herstellungskosten für die Vorrichtung gewährleistet. Da außerdem die für die Herstellung des Tomogramms nötigen Daten bei einer ausschließlich linearen Verschiebung bzw. Bewegung des Abtasters gesammelt werden können, wird die Datengewinnungszeit erheblich verkürzt, so daß auch Unscharfen im rekonstruierten Bild verhindert werden können, weil die Bewegung des Patienten bzw. Untersuchungsobjekts während dieser kurzen Datengewinnungszeit im allgemeinen sehr gering ist.
809833/1070

Claims (7)

  1. PATENTANSPRÜCHE
    ,J Vorrichtung für die Tomographie unter Verwendung von durchdringender Strahlung, mit einer Strahlungsquelle zur Abstrahlung einer durchdringenden Strahlung in fächerförmig aufgespreiztem Zustand in das Innere eines Scheiben- bzw. Schichtabschnitts eines Untersuchungsobjekts, mit einem Detektor mit Detektor- bzw. Meßelementen, um die durch den Schichtabschnitt hindurchgegangene Strahlung im wesentlichen in eine Anzahl von Strahlungsbündeln aufzuteilen, die jeweils unter einem Winkel auseinanderlaufen, und dabei die Intensität der Strahlungsbündel zu messen, mit einer Verschiebungseinheit zum Bewegen der Strahlungsbündel zwecks Abtastung des Schichtabschnitts, mit einer Datenverarbeitungseinheit zur Rekonstruktion eines Bilds des Schichtabschnitts auf der Grundlage der vom Detektor gelieferten Meßdaten und mit einer Anzeigeeinheit zur Wiedergabe des Bilds auf der (Grundlage der Ergebnisse der Rekonstruktion, dadurch gekennzeichnet , daß die Verschiebungseinheit (28) so ausgebildet ist, daß sie zumindest eine einmalige praktisch lineare Abtastung jedes Strahlungsbündels (15)
    Bl/eg
    -2-
    809833/1070
    ORIGINAL INSPECTED
    in einer Richtung quer zum Schichtabschnitt (60) zuläßt, wobei diese lineare Abtastung durchführbar ist, während eine entsprechende von mehreren, im Bereich der fäc herförmigen Ausbreitung der Strahlungsbündel liegenden Richtungen eingehalten wird.
  2. 2. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, daß die Verschiebungseinheit (28) eine Zahnstange (40) und ein damit kämmendes Ritzel (42) aufweist, das in Bewegungsrichtung der Verschiebungseinheit an einem Rahmen (20) derselben vorgesehen ist, auf dem die Strahlungsquelle (10) und der Detektor (24) montiert sind.
  3. 3. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet , daß die Strahlungsquelle (10) eine Röntgenstrahlröhre ist.
  4. 4. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet , daß die Strahlungsquelle (10) ein Röntgen-Schmalbündelgenerator (80) zum Verschwenken eines bleistiftförmigen Röntgen-Schmalbündels in der Weise ist, daß das Schmalbündel aufgrund seiner Schwenkbewegung insgesamt eine Fächerform erhält.
  5. 5. Vorrichtung nach Anspruch 1, dadurch gekennzeichnet, daß die Datenverarbeitungseinheit (44) eine Projektionsdaten-Rechnereinheit (70) zur Berechnung der Projektionsdaten bezüglich des übertragungswegs (64) jedes Strahlungsbündels auf der Basis der vom Detektor (24) übermittelten Meßdaten aufweist.
  6. 6. Vorrichtung nach Anspruch 5, dadurch gekennzeichnet , daß die Datenverarbeitungseinheit (44)
    -3-
    809833/1070
    eine Grunddaten-Rechnereinheit (72) aufweist, die ihrerseits eine erste Rechnereinrichtung (72a), um auf der Grundlage der von der Projektionsdaten-Rechnereinheit (70) gelieferten Projektionsdaten eindimensionale Fourier-Transformierte der Projektionsdaten bezüglich der betreffenden Abstrahlrichtung(en) des Strahlungsbündels (15) zu berechnen, sowie eine zweite Rechnereinrichtung (72b) umfaßt, um die eindimensionale Fourier-Transformierte bezüglich der über den Gesamtumfang des Schichtabschnitts gewählten Röntgenstrahlrichtungen auf der Grundlage der Fourier-Transformierten für die Strahlungsrichtungen im Ausbreitungsbereich der Strahlungsbündel zu berechnen, wobei diese Fourier-Transformierte durch die erste Rechnereinrichtung berechnet wird.
  7. 7. Vorrichtung nach Anspruch 6, dadurch gekennzeichnet, daß die Datenverarbeitungseinheit (44) eine Bildrekonstruktions-Einrichtung (74) mit einer dritten Rechnereinrichtung (74a) aufweist, um die zweidimensionale Fourier-Transformierte der Absorptionskoeffizienten-Verteilung des Scheiben- bzw. Schichtabschnitts (60) auf der Grundlage der von der Grunddaten-Rechnereinrichtung (72) gelieferten eindimensionalen Fourier-Transformierten zu berechnen, und weiterhin eine Fouriersche Inversionstransformations-Einrichtung (74b) umfaßt, um die Absorptionskoeffizienten-Verteilung des Schichtabschnitts (60) durch Anwendung der zweidimensionalen inversen bzw. invertierten Fourier-Transformierten auf die von der Rechnereinrichtung berechnete zweidimensionale Fourier-Transformierte zu berechnen.
    -4-
    809833/1070
DE2806219A 1977-02-14 1978-02-14 Computer-Tomograph Expired DE2806219C2 (de)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP52014116A JPS582695B2 (ja) 1977-02-14 1977-02-14 透過性放射線による検査装置

Publications (2)

Publication Number Publication Date
DE2806219A1 true DE2806219A1 (de) 1978-08-17
DE2806219C2 DE2806219C2 (de) 1987-12-23

Family

ID=11852136

Family Applications (1)

Application Number Title Priority Date Filing Date
DE2806219A Expired DE2806219C2 (de) 1977-02-14 1978-02-14 Computer-Tomograph

Country Status (5)

Country Link
US (1) US4282438A (de)
JP (1) JPS582695B2 (de)
DE (1) DE2806219C2 (de)
FR (1) FR2380593B1 (de)
GB (1) GB1598058A (de)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4196352A (en) * 1978-04-28 1980-04-01 General Electric Company Multiple purpose high speed tomographic x-ray scanner
DE3010780A1 (de) * 1980-03-20 1981-09-24 Siemens AG, 1000 Berlin und 8000 München Strahlendiagnostikeinrichtung
JPS56136529A (en) * 1980-03-28 1981-10-24 Tokyo Shibaura Electric Co Apparatus for reconstituting image
US4559597A (en) * 1982-07-07 1985-12-17 Clayton Foundation For Research Three-dimensional time-of-flight positron emission camera system
JPS59111738A (ja) * 1982-12-16 1984-06-28 株式会社東芝 X線断層撮影装置
JPS6140678A (ja) * 1984-08-01 1986-02-26 Hitachi Medical Corp 画像処理装置
JPS61209641A (ja) * 1985-03-15 1986-09-17 株式会社東芝 X線ct装置
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
WO2001069500A1 (en) * 2000-03-10 2001-09-20 Medorder, Inc. Method and system for accessing healthcare information using an anatomic user interface
US6568343B1 (en) * 2001-11-05 2003-05-27 Grant Roy Hill Protective enclosure for watercraft hulls
US8535337B2 (en) * 2010-04-26 2013-09-17 David Chang Pedicle screw insertion system and method
EP2724147A1 (de) * 2011-06-21 2014-04-30 Koninklijke Philips N.V. Bildgebungsvorrichtung
JP2014042564A (ja) * 2012-08-24 2014-03-13 Sony Corp 画像処理装置、画像処理方法、および画像処理システム
US9418415B2 (en) * 2014-02-05 2016-08-16 Shimadzu Corporation Trabecular bone analyzer
US9490099B2 (en) 2014-08-20 2016-11-08 Wisconsin Alumni Research Foundation System and method for multi-source X-ray-based imaging
CN116664714B (zh) * 2023-07-26 2023-10-20 济南汉江光电科技有限公司 一种基于x射线微束传输模型的ct算法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2442009A1 (de) * 1973-08-31 1975-04-03 Emi Ltd Geraet zur untersuchung eines koerpers mittels durchdringender strahlung
DE2520539A1 (de) * 1974-05-08 1975-11-13 Emi Ltd Vorrichtung und verfahren zum betrieb eines radiologischen geraetes
DE2538517A1 (de) * 1974-08-28 1976-03-11 Emi Ltd Radiologisches geraet

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2667585A (en) * 1951-02-15 1954-01-26 Hartford Nat Bank & Trust Co Device for producing screening images of body sections
US3101407A (en) * 1959-04-09 1963-08-20 Jr John Daniel Shipman Fluoroscope system utilizing an image storage tube
US3922552A (en) * 1974-02-15 1975-11-25 Robert S Ledley Diagnostic X-ray systems
US3965353A (en) * 1974-12-06 1976-06-22 Albert Macovski Cross-sectional X-ray emission imaging system
US4031395A (en) * 1975-02-21 1977-06-21 Emi Limited Radiography
US3976885A (en) * 1975-03-18 1976-08-24 Picker Corporation Tomography system having nonconcurrent, compound axial scanning
JPS51127689A (en) * 1975-04-28 1976-11-06 Jeol Ltd X-ray crossing layer image photogrphing device
US4045672A (en) * 1975-09-11 1977-08-30 Nihon Denshi Kabushiki Kaisha Apparatus for tomography comprising a pin hole for forming a microbeam of x-rays
DE2613809B2 (de) * 1976-03-31 1979-01-04 Siemens Ag, 1000 Berlin Und 8000 Muenchen Röntgenschichtgerät zur Herstellung von Transversal-Schichtbildern
US4138721A (en) * 1976-11-11 1979-02-06 Board of Trustees of The Lelane Standard Junior University Limited scan angle fan beam computerized tomography

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2442009A1 (de) * 1973-08-31 1975-04-03 Emi Ltd Geraet zur untersuchung eines koerpers mittels durchdringender strahlung
DE2520539A1 (de) * 1974-05-08 1975-11-13 Emi Ltd Vorrichtung und verfahren zum betrieb eines radiologischen geraetes
DE2538517A1 (de) * 1974-08-28 1976-03-11 Emi Ltd Radiologisches geraet

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Electronics, 14.Okt.76, H.21, S.89-94 *

Also Published As

Publication number Publication date
JPS5399892A (en) 1978-08-31
US4282438A (en) 1981-08-04
FR2380593B1 (fr) 1985-07-26
FR2380593A1 (fr) 1978-09-08
JPS582695B2 (ja) 1983-01-18
DE2806219C2 (de) 1987-12-23
GB1598058A (en) 1981-09-16

Similar Documents

Publication Publication Date Title
DE102006015358B4 (de) Fokus/Detektor-System einer Röntgenapparatur zur Erzeugung von Phasenkontrastaufnahmen, zugehöriges Röntgen-System sowie Speichermedium und Verfahren zur Erzeugung tomographischer Aufnahmen
EP0989520B1 (de) Computertomographie-Verfahren mit kegelförmigem Strahlenbündel
DE69114932T2 (de) Rechnergesteuertes Tomographiesystem.
DE69310295T2 (de) Verbesserter röntgenstrahlenraumcomputertomographscanner
DE69033921T2 (de) Rechnergesteurtes tomographisches Bildrekonstruktionsverfahren für Spiralabtasten
DE60304786T2 (de) Röntgen-Computertomograph
DE112005002171B4 (de) Teilchenstrahl-Bestrahlungsverfahren und dafür verwendete Teilchenstrahl-Bestrahlungsvorrichtung
EP0028431B1 (de) Anordnung zur Ermittlung der Streustrahlungsdichteverteilung in einem ebenen Untersuchungsbereich
DE102005048891B4 (de) Verfahren zur Kalibrierung eines CT-Systems mit mindestens zwei winkelversetzt zueinander angeordneten Fokus/Detektorsystemen und Computertomographie-System
DE2806219A1 (de) Vorrichtung fuer die tomographie
DE102006017290A1 (de) Fokus/Detektor-System einer Röntgenapparatur zur Erzeugung von Phasenkontrastaufnahmen
DE2613809A1 (de) Roentgenschichtgeraet zur herstellung von transversal-schichtbildern
EP0153786A2 (de) Röntgengerät
DE102006046034A1 (de) Röntgen-CT-System zur Erzeugung projektiver und tomographischer Phasenkontrastaufnahmen
DE102006015355A1 (de) Fokus/Detektor-System einer Röntgenapparatur zur Erzeugung von Phasenkontrastaufnahmen
DE3517101C1 (de) Vorrichtung zur digitalen Subtraktions-Angiographie im Energiesubstraktions-Modus
DE2717349A1 (de) Roentgenschichtgeraet zur herstellung von transversalschichtbildern
DE102005031893A1 (de) Verfahren zur Fokusjustage in einem CT
DE102007014829B3 (de) Verfahren zur Streustrahlungskorrektur in bildgebenden Röntgengeräten sowie Röntgenbildgebungssystem
DD140202A5 (de) Verfahren und anordnung zur verminderung von bildfehlern in computer-tomographiebildern
DE2723401A1 (de) Schichtgeraet zur herstellung von transversalschichtbildern
DE2648503A1 (de) Radiographisches geraet
EP0037008A1 (de) Strahlendiagnostikeinrichtung
DE102011080263B4 (de) Vorrichtung zur Röntgenbildgebung für Projektionsradiographie und Computertomographie sowie ein entsprechend ausgebildetes Verfahren zur Röntgenbildgebung
DE2745390C2 (de) Röntgensichtgerät für die Herstellung von Transversalschichtbildern

Legal Events

Date Code Title Description
OAP Request for examination filed
8128 New person/name/address of the agent

Representative=s name: HENKEL, G., DR.PHIL. FEILER, L., DR.RER.NAT. HAENZ

D2 Grant after examination
8363 Opposition against the patent
8327 Change in the person/name/address of the patent owner

Owner name: KABUSHIKI KAISHA TOSHIBA, KAWASAKI, KANAGAWA, JP

8366 Restricted maintained after opposition proceedings
8325 Change of the main classification

Ipc: A61B 6/03

8381 Inventor (new situation)

Free format text: NISHIDA, SUSUMU, YOKOHAMA, JP INOUYE, TAMON, KAWASAKI, JP YOSHIDA, TADATOKI, TOKIO/TOKYO, JP SAITO,KIYOTO, YOKOHAMA, JP

8305 Restricted maintenance of patent after opposition
D4 Patent maintained restricted
8339 Ceased/non-payment of the annual fee