DE4446819A1 - CT System mit wendelförmiger Abtastung durch Doppelfächerbündel - Google Patents
CT System mit wendelförmiger Abtastung durch DoppelfächerbündelInfo
- Publication number
- DE4446819A1 DE4446819A1 DE4446819A DE4446819A DE4446819A1 DE 4446819 A1 DE4446819 A1 DE 4446819A1 DE 4446819 A DE4446819 A DE 4446819A DE 4446819 A DE4446819 A DE 4446819A DE 4446819 A1 DE4446819 A1 DE 4446819A1
- Authority
- DE
- Germany
- Prior art keywords
- data
- array
- fan
- image
- arrays
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims description 16
- 238000002591 computed tomography Methods 0.000 title abstract description 12
- 230000008569 process Effects 0.000 title description 2
- 238000003491 array Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 238000013519 translation Methods 0.000 description 4
- 238000013170 computed tomography imaging Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000003313 weakening effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Pulmonology (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
Die Erfindung bezieht sich auf eine Computer-Tomographie
(CT)-Bildgebungseinrichtung und insbesondere auf ein Ver
fahren und eine entsprechende Einrichtung zum Rekonstruk
tieren von Bildern aus Daten, die während einer wendelför
migen Abtastung unter Verwendung von Twin- oder Doppelfä
cherbündeln gewonnen sind.
In einem gegenwärtigen Computer-Tomographie-System proji
ziert eine Röntgenstrahlquelle ein fächerförmiges Bündel,
das kollimiert ist, um in eine X-Y-Ebene von einem kartesi
schen Koordinatensystem zu liegen, die die "Bildebene" ge
nannt wird. Das Röntgenbündel tritt durch das abzubildende
Objekt, wie beispielsweise einen medizinischen Patienten,
hindurch und trifft auf eine lineare Anordnung (Array) von
Strahlungsdetektoren. Die Intensität der durchgelassenen
Strahlung ist abhängig von der Schwächung des Röntgen
strahlbündels durch das Objekt, und jeder Detektor erzeugt
ein getrenntes elektrisches Signal, das eine Messung der
Bündelschwächung ist. Die Schwächungsmessungen von allen
Detektoren werden getrennt gewonnen, um das Transmissions-
bzw. Durchlässigkeitsprofil zu erzeugen.
Die Quelle und die lineare Detektoranordnung in einem kon
ventionellen CT-System werden mit einem Gestell in der
Bildebene und um das Objekt herum gedreht, so daß sich der
Winkel, unter dem das Röntgenbündel das Objekt schneidet,
konstant ändert. Eine Gruppe von Röntgenschwächungsmessun
gen aus der Detektoranordnung bei einem gegebenen
Gestellwinkel wird ein "Experiment bzw. View" genannt, und
eine "Abtastung bzw. Scan" des Objektes weist einen Satz
von Views auf, die während einer Umdrehung der Röntgen
quelle und des Detektors bei verschiedenen Gestellwinkeln
gemacht werden. Bei einer axialen Abtastung werden die Da
ten verarbeitet, um ein Bild zu konstruieren, das einer
zweidimensionalen Scheibe durch das Objekt entspricht. Das
vorherrschende Verfahren zum Rekonstruieren eines Bildes
aus einem Satz von Daten wird in der Technik als die gefil
terte Rückprojektions-Technik bezeichnet. Dieses Verfahren
wandelt die Schwächungsmessungen von einer Abtastung in
ganze Zahlen um, die "CT-Zahlen" oder "Hounsfield-Einhei
ten" genannt werden, die zur Steuerung der Helligkeit von
einem entsprechenden Pixel auf einer Kathodenstrahlröhren-
Anzeige verwendet werden.
Um die gesamte Abtastzeit zu verkürzen, wenn zahlreiche
Scheiben gewonnen werden, wird eine sogenannte
"wendelförmige" Abtastung durchgeführt, wo der Patient be
wegt wird, während das Gestell umläuft, um Daten für die
vorgeschriebene Anzahl von Scheiben zu gewinnen. Die wen
delförmige Abtastung führt jedoch gewisse Fehler in bezug
auf die Daten in den gewonnen tomografischen Projektions
sätzen ein. Die Mathematik der tomografischen Rekonstruk
tion geht davon aus, daß der tomografische Projektionssatz
entlang einer konstanten z-Achsen-Scheibenebene gewonnen
wird. Die wendelförmige Abtastbahn weicht deutlich von die
sem Zustand ab, und diese Abweichung, wenn sie nicht kom
pensiert wird, hat Bildartefakte in dem rekonstruierten
Scheibenbild zur Folge. Die Schwere der Bildartefakte hängt
im allgemeinen von der "Wendel-Versetzung" in den Projekti
onsdaten ab, die als die Differenz zwischen den Tischorten
der abgetasteten Daten und dem z-Achsenwert der gewünschten
Scheibenebene gemessen wird. Aus der wendelförmigen
Abtastung resultierende Fehler werden kollektiv als
"Schräg"-Fehler bezeichnet.
Es sind mehrere Verfahren verwendet worden, um Schräg- bzw.
Skew-Fehler bei der wendelförmigen Abtastung zu verklei
nern. Eine erste Lösung ist in der US-PS 5 046 003 mit dem
Titel "Method for Reducing Skew Image Artifacts in Helical
Projection Imaging" beschrieben, wo eine ungleichförmige
Tischbewegung verwendet wird, um die wendelförmig gewon
nenen Projektionen nahe der Scheibenebene zu konzentrieren,
während die Beschleunigungskräfte auf den Patienten be
grenzt werden.
Gemäß der US-PS 5 270 923 mit der Bezeichnung "Computerized
Tomographic Image Recontruction Method for Helical Scan
ning" werden Schräg-Artefakte dadurch verkleinert, daß zwi
schen zwei Halbabtastungen der Daten interpoliert wird, die
jeweils nur 180° plus den Fächerbündelwinkel der Gestell
drehung erfordern. Die Halbabtastungen erfordern weniger
Gestelldrehung und somit weniger Tischbewegung, wodurch die
gesamte wendelförmige Versetzung der Projektionsdaten ver
kleinert wird.
Eine dritte Lösung ist in der US-PS 5 233 518 mit der Be
zeichnung "Extrapolative Reconstruction Method for Helical
Scanning" beschrieben, wonach Schräg-Artefakte durch Inter
polieren und Extrapolieren zwischen zwei Teilprojektionsda
ten von nur 180° der Gestelldrehung verkleinert werden. Die
zwei Teilprojektionssätze erfordern weniger Gestelldrehung
als die obige Halbabtastlösung und dadurch wird die gesamte
wendelförmige Versetzung der Projektionsdaten weiter ver
kleinert. Schräg-Artefakte werden weiter verkleinert durch
Gewichtung der gewonnenen Daten als eine Funktion des Ab
standes von der Scheibenebene, wie es in der US-PS 5 170 346
mit dem Titel "Method for Reducing Patient Translation
Artifacts in Tomographic Imaging" beschrieben ist.
Diese bekannten Verfahren erzeugen gute Bilder, wenn die
Translationsgeschwindigkeit des Patienten mittelmäßig ist,
d. h. nicht schneller als 10 mm pro Gestelldrehung. Weitere
Verkürzungen in der Abtastzeit können nicht erzielt werden
durch Vergrößern der Translationsgeschwindigkeit des Pati
enten, weil signifikate Bildverschlechterungen auftreten.
Es ist Aufgabe der Erfindung, ein Verfahren und eine ent
sprechende Einrichtung der eingangs genannten Art so auszu
gestalten, daß die Geschwindigkeit des Patiententisches
während einer wendelförmigen Abtastung erhöht werden kann,
ohne daß Bewegungsartefakte vergrößert werden.
Erfindungsgemäß wird ein CT-System geschaffen, bei dem zwei
Fächerbündel, die entlang der Achse der Gestelldrehung ver
schoben sind, Projektionsdaten während einer wendelförmigen
Abtastung gewinnen. Dieses System kann approximiert werden
als zwei Detektorreihen, die gleichzeitig Projektionsmes
sungen an unterschiedlichen axialen Orten sammeln. Wenn ein
derartiges System in dem wendelförmigen Modus arbeitet, ge
neriert es verwobene Doppelwendeln im Gegensatz zu einer
einzelnen Wendel von einer üblichen wendelförmigen Fächer
bündelabtastung. Die verwobenen Doppelwendeln, die durch
die zwei Fächerbündel aufgezeichnet werden, erzielen Pro
jektionsdaten, aus denen Bilder in jeder vorgeschriebenen
Scheibe mit verminderter Bildverschlechterung aufgrund der
Patienten-Translation rekonstruiert werden können. Insbe
sondere werden Projektionsraum-Datenarrays aus Projektions
daten gewählt, die durch jedes Fächerbündel gewonnen worden
sind, wobei die Daten in jeder Array gewichtet werden, um
die Translationsbewegung des Patienten zu korrigieren und
Datenredundanzeinflüsse auszugleichen.
Die Erfindung wird nun mit weiteren Merkmalen und Vorteilen
anhand der folgenden Beschreibung und Zeichnung von Ausfüh
rungsbeispielen näher erläutert.
Fig. 1 ist eine Darstellung von einem CT-Bildgebungssy
stem, in dem die Erfindung verwendet werden kann;
Fig. 2 ist ein schematisches Blockdiagramm von dem CT-
Bildgebungssystem;
Fig. 3 ist eine schematische Darstellung von einem doppel
ten Röntgenstrahl-Fächerbündel im Schnitt entlang der z-
Achse;
Fig. 4 ist eine schematische Darstellung von einem Rönt
genstrahl-Fächerbündel, wie es in der Bildebene gesehen
wird;
Fig. 5 ist eine Abbildung von einem vollständigen Satz von
Projektionsdaten, die durch ein einzelnes Fächerbündel ge
wonnen sind;
Fig. 6 ist eine Abbildung von Projektionsdaten, die durch
das Doppel-Röntgenstrahlen-Fächerbündel gemäß Fig. 2 ge
wonnen sind, und zeigt, wie sie kombiniert werden, um Pro
jektionsdaten für eine einzelne Scheibe zu bilden;
Fig. 7 ist ein Blockdiagramm von einem Bild-Rekonstruktor,
der Teil des CT-Bildgebungssystem gemäß Fig. 2 bildet; und
Fig. 8 ist eine Abbildung von Projektionsdaten, die gewon
nen wurden, wenn die wendelförmige Tischgeschwindigkeit re
lativ hoch ist.
Wie in Fig. 3 gezeigt ist, werden zwei Detektorreihen in
einem Doppel- oder Twinfächerbündelsystem verwendet. Das
Röntgen-Fächerbündel ist im Effekt in zwei Fächerbündel ge
spalten, die entlang der z-Drehachse verschoben sind. Wenn
die zwei Fächerbündel als das vordere und das hintere
Bündel bezeichnet werden, ist D der Abstand zwischen der
Mitte von diesen zwei Fächern, wenn er an dem Iso-Zentrum
gemessen wird.
In Fig. 4 ist der Gestellwinkel mit θ bezeichnet. Die De
tektorwinkel und die Viewwinkel von sowohl dem vorderen als
auch dem hinteren Fächerbündel können mit Φ± bzw. θ± be
zeichnet werden, wobei der Viewwinkel von jedem Fächerbün
del durch die entsprechende örtliche Koordinate bezeichnet
ist, die Null ist, wenn das Fächerbündel sich an der zu re
konstruierenden Bildposition befindet.
Wenn der Tischvorschub pro Gestellumdrehung mit Z angenom
men wird, dann ist die axiale Tischposition z und die zu
rekonstruierende Scheibenposition zs. Die Positionen der
vorderen und hinteren Bündel sind z+ bzw. z-. Sie können
ausgedrückt werden als:
Das wendelformige Z± ist eine Funktion von θ± ausgedrückt als
z±(θ±). Es ist bekannt, daß beim Fehlen der Patientenbewe
gung ein Satz von Fächerbündelprojektionen, die über π Ra
dian plus den Winkel des Fächerbündels der Gestelldrehung
gesammelt sind, ausreichend für eine Rekonstruktion ist.
Bezugnehmend auf Fig. 5 kann ein derartiger Projektions
satz in dem Projektionsraum (θ, Φ) beschrieben werden, wo
bei Φ in dem Bereich (-Φm, Φm) und 6 in dem Bereich (-θm,
θm) liegen, wobei Φm die Hälfte des Fächerwinkels und θm=
π/2 + Φm ist. In der Theorie sind die Projektionsdaten in
dem gepunkteten Bereich in Fig. 5 für eine Rekonstruktion
von einem Bild ausreichend.
Eine wendelförmige Twin- oder Doppelabtastung ergibt zwei
Sätze von Projektionsdaten, die als P+(θ+, Φ+) und P-(θ-,
Φ-) bezeichnet sind. Um die Scheibenprofilerweiterung bei
der Rekonstruktion zu minimieren, werden nur diejenigen
Projektionsdaten von einem Datensatz, der bei zs zentriert
und mit θ+ oder θ- in dem Bereich von -θm bis θm gesammelt
ist, verwendet. Diese zwei Datensätze sind in Fig. 6 ge
zeigt, wo die Zentren dieser zwei Datensätze um einen Be
trag Δ=2pπ in der Gestellwinkelrichtung aufgrund der Ver
schiebung der zwei Detektorreihen verschoben sind.
Die Flächen der redundanten Daten für die wendelförmigen
Twinfächerbündelprojektionen sind in Fig. 6 in ähnlicher
Weise schraffiert und mit 1+ bis 5+ bzw. 1- bis 5- bezeich
net. Die Grenzen zwischen diesen Bereichen in den zwei Da
tensätzen sind wie folgt:
Grenze | |
Gleichung | |
1./2. | |
θ = θm - 2Φm - 2Φ | |
2./3. | θ = - Δ + Φm |
3./4. | θ = - Δ + θm - 2Φm - 2Φ |
4./5. | θ = - θm + 2Φm - 2Φ |
1./2. | θ = - θm + 2Φm - 2Φ |
2./3. | θ = Δ - θm |
3./4. | θ = Δ - θm + 2Φm - 2Φ |
4./5. | θ = θm - 2Φm - 2Φ |
Die redundanten Daten können richtig kombiniert werden, um
einen neuen Datensatz an der Scheibenposition zs hervorzu
rufen, aus dem ein Bild rekonstruiert wird. Das Zentrum des
hervorgerufenen Datensatzes könnte irgendwo liegen, ohne
das Resultat zu beeinflussen. Aus Zweckmäßigkeitsgründen
ist das Zentrum so gewählt, daß es mit dem Zentrum des vor
deren Bündels ausgerichtet ist, wie es in Fig. 6 gezeigt
ist. Dieser neue Datensatz ist durch Bereiche 1s bis 5s. Es
sei darauf hingewiesen, daß die Bereiche 2s bis 5s einen
vollständigen Datensatz bilden, aus dem ein Bild an der ge
wünschten Scheibenposition rekonstruiert werden kann.
Der Datenkombinierungsprozeß kann allgemein ausgedrückt
werden:
Ps(θ, Φ, Zs) = w+(θ+, Φ+)P+(θ+, Φ+)
+ w-(θ-, Φ-)P-(θ-, Φ-) (2)
Die Gewichtungen sind gegeben durch:
Für Fig. 6 gilt: Bereiche 2+ und 2- werden kombiniert, um
den Bereich 2s zu erhalten; die Bereiche 3+, 5- und 1- wer
den kombiniert, um den Bereich 3s zu erhalten; die Bereiche
4+ und 4- werden kombiniert, um den Bereich 4s zu erhalten;
und die Bereiche 3-, 5+ und 1+ werden kombiniert, um den
Bereich 5s zu erhalten.
Die Bereiche 3s und 5s werden erhalten, indem drei redun
dante Messungen kombiniert werden. Dies geschieht in der
Weise, daß jeweils zwei redundante Messungen von dem einen
Datensatz mit der redundanten Messung von dem anderen Da
tensatz kombiniert werden, wobei Gleichung (2) verwendet
wird, um zwei kombinierte Projektionen P(θ₁, Φ₁, zs) und
P(θ₂, Φ₂, zs) zu bilden. Weiterhin werden diese zwei kombi
nierten Projektionen wie folgt verknüpft:
P₂(θ+, Φ+, Zs) = a(θ₁, Φ₁) P(θ₁, Φ₁, zs)
+ a(θ₂, Φ₂) P(θ₂, Φ₂, zs) (4)
wobei (θ₁, Φ₁) und (θ₂, Φ₂) die folgende Beziehung haben:
θ₂ = θ₁ + π + 2Φ₁ und Φ₂ = - Φ₁ (5)
Somit ist a(θ, Φ) tatsächlich die Halbtastgewichtungen,
das genügt:
a(θ, Φ) + a(θ + π + 2Φ, - Φ) = 1 (6)
Mit der oben beschriebenen Kombinierungs- bzw. Verknüp
fungsstrategie gibt es zwei Arten von Kombinationen. Die
eine Art enthält die Kombinierung von 5- und 3+, 4- und 4+
und 3- und 5+. Der Rest der Kombinierung bildet die andere
Art. Die letzte wird die erste Kombinierungsart genannt,
während die erste die zweite Kombinierungsart genannt wird.
Für die zweite Kombinierungsart ergibt sich:
θ- = θ+ - Δ und Φ- = Φ+ (7a)
Somit erhält man dann aus Gleichungen (3) und (1)
W2±(θ, Φ) = 1 ± θ/Δ (8a)
Für die erste Kombinierungsart erhält man:
θ- = θ+ - Δ + π + 2Φ+ und Φ- = -Φ+ (7b)
In ähnlicher Weise ergibt sich dann aus den Gleichungen (3)
und (1):
Die resultierenden, kombinierten Gewichtungen, die auf je
den Bereich von P+ (θ, Φ) und P- (θ, Φ) ausgeübt werden,
sind:
Bereich 1+: a+(θ, Φ) w1+(θ, Φ)
Bereich 1-: a-(θ, Φ) w1-(θ, Φ)
Bereich 2+: w1+(θ, Φ)
Bereich 2-: w1-(θ, Φ)
Bereich 3+: a+(θ - Δ, Φ) w2+(θ, Φ) + (1 - a+(θ - Δ, Φ)) w1+(θ, Φ)
Bereich 3-: a-(θ + Δ, Φ) w2-(θ, Φ) + (1 - a-(θ + Δ, Φ)) w1-(θ, Φ)
Bereich 4+: w2+(θ, Φ)
Bereich 4-: w2-(θ, Φ)
Bereich 5+: a+(θ, Φ) w2+(u, F)
Bereich 5-: a-(θ, Φ) w2-(θ, Φ) (9)
Bereich 1-: a-(θ, Φ) w1-(θ, Φ)
Bereich 2+: w1+(θ, Φ)
Bereich 2-: w1-(θ, Φ)
Bereich 3+: a+(θ - Δ, Φ) w2+(θ, Φ) + (1 - a+(θ - Δ, Φ)) w1+(θ, Φ)
Bereich 3-: a-(θ + Δ, Φ) w2-(θ, Φ) + (1 - a-(θ + Δ, Φ)) w1-(θ, Φ)
Bereich 4+: w2+(θ, Φ)
Bereich 4-: w2-(θ, Φ)
Bereich 5+: a+(θ, Φ) w2+(u, F)
Bereich 5-: a-(θ, Φ) w2-(θ, Φ) (9)
wobei:
a+(θ, Φ) = a-(-θ, -Φ) = a(θ, Φ) (10)
Es gibt viele Wege, um die Redundanz-Gewichtung a(θ, Φ) zu
wählen. Als ein Beispiel kann die folgende Gewichtung ver
wendet werden.
a(θ, Φ) = a[x(θ, Φ)] = 3x²(θ, Φ) - 2x²(θ, Φ) (11a)
wobei
In der vorstehenden Ableitung ist angenommen, daß θm = π/2
+ Φm. Da in der Praxis θm größer als dieser Wert sein
kann, kann die Ableitung von ähnlichen Gleichungen für eine
derartige Situation ebenfalls abgeleitet werden.
Die Projektionsdaten in jeder Array P+ und P- sind somit
gewichtet, wie es in Gleichung (9) angegeben ist, und dann
kombiniert bzw. verknüpft, um das gewünschte Bild zs zu
bilden. Diese Verknüpfung kann auf zwei Wegen geschehen.
Die entsprechenden Elemente in jeder richtig gewichteten
Array können addiert werden, um die kombinierte Projekti
onsarray zu bilden, die dann gefiltert und rückprojiziert
wird, um das Bild zu bilden. Alternativ kann jede richtig
gewichtete Array getrennt gefiltert und rückprojiziert wer
den, um getrennte Bilder zu bilden. Die entsprechenden Ele
mente oder Pixels in jeder Bildraumarray können dann ad
diert werden, um das endgültige Bild bei zs zu bilden.
Die vorstehend beschriebene Erfindung kann direkt auf den
Fall angewendet werden, wo die wendelförmige Tischgeschwin
digkeit eine gewisse Geschwindigkeit überschreitet, so daß
(p < 1/2 - Φm/π). Mit der folgenden Transformation ist sie
auch auf den Fall anwendbar, wo die wendelförmige Tischge
schwindigkeit langsamer ist. In diesem Fall wird, wie es in
Fig. 8 gezeigt ist, P+ (θ, Φ) entlang der Richtung von θ
um 2π verschoben. Wenn die periodische Eigenschaft des
Viewwinkels gegeben ist, sind die ursprünglichen Vorderbün
deldaten und ihre verschobene Version äquivalent. Deshalb
kann das ursprüngliche hintere Bündel mit der verschobenen
Version des vorderen Bündels kombiniert werden. Um dies zu
tun, kann das ursprüngliche hintere Bündel als das neue
vordere Bündel und die verschobene Version des vorderen
Bündels als das neue hintere Bündel neu definiert werden.
Die neue Versetzung ist nun definiert als Δ = 2π - 2pπ und
man hat 0 < Δ π 2Φm für 1/2 + Φm/π < π < 1.
Für eine detaillierte Beschreibung der Erfindung wird nun
zunächst auf die Fig. 1 und 2 Bezug genommen, in denen
ein Computer-Tomographie (CT)-Bildgebungssystem 10 gezeigt
ist, das ein Gestell 12 aufweist, das einen CT-Scanner der
"dritten Generation" darstellt. Das Gestell 12 hat eine
Röntgenstrahlquelle 13, das ein Bündel von Röntgenstrahlen
14 in Richtung auf eine Detektoranordnung 16 auf der gegen
überliegenden Seite des Gestells projiziert. Die Detektoran
ordnung (-Array) 16 ist aus zwei Reihen von Detektorele
menten 18 gebildet, die zusammen die projizierten Röntgen
strahlen abtasten, die durch einen medizinischen Patienten
15 hindurchtreten. Jedes Detektorelement 18 erzeugt ein
elektrisches Signal, das die Intensität von einem auftref
fenden Röntgenstrahlbündel und somit die Schwächung des
Bündels darstellt, wenn dieses durch den Patienten hin
durchtritt. Während einer Abtastung, um Röntgenprojektions
daten zu gewinnen, werden das Gestell 12 und die darauf an
gebrachten Komponenten um einen Drehmittelpunkt 19 gedreht.
Die Drehung des Gestells und der Betrieb der Röntgenstrahl
quelle 13 werden durch eine Regeleinrichtung 20 des CT-Sy
stems geregelt. Der Regelmechanismus 20 enthält eine Rönt
gen-Steuerung 22, die Versorgungs- und Zeitsteuersignale an
die Röntgenstrahlquelle 13 und eine Gestellmotor-Steuerung
23 liefert, die die Drehgeschwindigkeit und die Position
des Gestells 12 steuert. Ein Datengewinnungssystem (DAS) 24
in dem Regelmechanismus 20 tastet analoge Daten von den De
tektorelementen 18 ab und wandelt die Daten in digitale Da
ten für eine nachfolgende Bearbeitung um. Ein Bild-
Rekonstruktor 24 empfängt abgetastete und digitalisierte
Röntgenstrahldaten von dem DAS 24 und führt eine Hochge
schwindigkeits-Bildrekonstruktion nach dem erfindungsgemä
ßen Verfahren aus. Das rekonstruierte Bild wird als eine
Eingangsgröße einem Computer 26 zugeführt, der das Bild in
einer Massenspeichervorrichtung 29 speichert.
Der Computer 26 empfängt auch Befehle und Abtastparameter
von einem Operator über eine Konsole 30, die ein Tastenfeld
(Keyboard) aufweist. Eine zugeordnete Kathodenstrahlröhren
anzeige 32 gestattet dem Operator, das rekonstruierte Bild
und andere Daten von dem Computer 26 zu beobachten. Die von
dem Operator zugeführten Befehle und Parameter werden von
dem Computer 26 verwendet, um Steuersignale und Information
an das DAS 24, die Röntgensteuerung 22 und die Gestellmo
tor-Steuerung 23 zu liefern. Zusätzlich arbeitet der Compu
ter 26 als eine Tischmotor-Steuerung, die einen motorisier
ten Tisch 36 steuert, um den Patienten 15 in dem Gestell 12
zu positionieren.
Die vorliegende Erfindung ist in dem Bild-Rekonstruktor 25
implementiert, der in Fig. 7 genauer gezeigt ist. Jedes
Experiment bzw. View von Daten aus dem DAS 24 für das erste
Fächerbündel wird bei 51 empfangen, wo es vorverarbeitet
wird, um es hinsichtlich verschiedener bekannter Fehler zu
korrigieren, wie beispielsweise Bündelhärtung, Versetzungen
(Offsets) und Änderungen in dem Detektor und der Kanalver
stärkung. Weiterhin wird der negative Logarithmus genommen,
um Projektionsdaten zu liefern, die in einer Projektionsda
tenarray 52 gespeichert werden. Die gleiche Vorverarbeitung
wird auf die Abtastdaten des zweiten Fächerbündels bei 53
ausgeübt und sie werden in einer Projektionsdatenarray 54
gespeichert. Es sind die in den Arrays 52 und 54 gespei
cherten Projektionsdaten, die gemäß der Erfindung
kombiniert bzw. verknüpft werden, um ein Scheibenbild zu
erzeugen.
Weiterhin wird gemäß Fig. 7 die rechte Seite der Projekti
onsdaten in der Array 52 für das Bündel 1 ausgelesen, und
die entsprechende Gewichtungsfunktion, die oben in Glei
chung (9) angegeben ist, wird bei 55 ausgeübt. Die gewich
teten Projektionsdaten werden in die entsprechende Stelle
in einer Array 56 geschrieben, und diese gewichteten Pro
jektionsdaten werden in einer üblichen Weise bei 57 gefil
tert und rückprojiziert, um eine Bilddatenarray 58 für das
Bündel 1 zu erzeugen. In ähnlicher Weise wird der entspre
chende Datensatz des Bündels 2 aus der Array 54 ausgelesen,
und die entsprechende Gewichtungsfunktion, die oben in
Gleichung (9) angegeben ist, wird bei 59 ausgeübt. Die ent
stehende gewichtete Projektionsdatenarray 60 wird bei 61
gefiltert und rückprojiziert, um eine Bilddatenarray 62 für
das zweite Bündel zu erzeugen.
Ein Scheibenbild an der Stelle zs wird bei 63 erzeugt, in
dem die zwei Bildarrays 58 und 62 kombiniert bzw. verknüpft
werden. Dies geschieht durch eine Summierstelle 65, die die
Größe von jedem Pixel in Bündel 1 der Array 58 und die
Größe des entsprechenden Pixels in dem Bündel 2 der Array
62 addiert. Die entstehende Scheibenbildarray 63 kann ge
speichert werden für eine spätere Verwendung oder Anzeige
für den Operator.
Es sind jedoch noch weitere Ausführungsbeispiele möglich.
Beispielsweise kann die Filterung und Rückprojektion
zunächst für das eine Bündel und dann für das andere Bündel
erfolgen, anstatt daß dies parallel geschieht, wie es hier
beschrieben wurde. Es ist auch möglich, die gewichteten
Projektionsdaten in den Arrays 56 und 60 zu kombinieren
bzw. zu verknüpfen, bevor die Filterung und Rückprojektion
erfolgt. Die Erfindung ist auch auf CT-Scanner der vierten
Generation anwendbar.
Claims (4)
1. Verfahren zum Erzeugen eines tomografischen Bildes von
einem Objekt aus Projektionsdaten, die in einer wendelför
migen Abtastung gewonnen werden, wobei die Daten von zwei
Röntgenfächerbündeln abgeleitet werden, die entlang einer
z-Achse angeordnet sind und die jeweils eine Reihe von Fä
cherbündelprojektionen bei mehreren Gestellwinkeln θ um die
z-Achse erzeugen und jede Fächerbündelprojektion einen Win
kel Φ einschließt, enthaltend die Schritte:
- a) Gewinnen eines Satzes von Projektionen über Gestell winkeln θ in dem Bereich über 180° + Φ, wenn das abgebil dete Objekt entlang der z-Achse verschoben wird, wobei der Satz von Projektionen eine erste Datenarray bildet, die die Schwächung von Röntgenenergie in dem ersten Fächerbündel durch das Objekt angibt, und eine zweite Datenarray bildet, die die Schwächung der Röntgenenergie in dem zweiten Fä cherbündel durch das Objekt angibt,
- b) Multiplizieren der Daten in der ersten Datenarray mit einem ersten Satz von Gewichtungswerten (55),
- c) Multiplizieren der Daten in der zweiten Datenarray mit einem zweiten Satz von Gewichtungswerten (59),
- d) Kombinieren bzw. Verknüpfen der entsprechenden gewich teten Daten in den ersten und zweiten Datenarrays (65),
- e) Erzeugen des tomografische Bildes (63) aus den ver knüpfen ersten und zweiten Datenarrays.
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß
die gewichteten Daten in den ersten und zweiten Datenarrays
vor der Verknüpfung (65) separat gefiltert und rückproji
ziert werden (56, 60).
3. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß
die ersten und zweiten Datenarrays jeweils in mehrere Be
reiche dividiert werden (Fig. 6) und die gewichtenden
Werte, die in den Schritten b) und c) verwendet werden,
teilweise dadurch bestimmt werden, in welchem der mehreren
Bereiche die multiplizierten Daten liegen.
4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, daß
die ersten und zweiten Datenarrays jeweils in fünf Bereiche
geteilt werden.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/175,703 US5469487A (en) | 1993-12-30 | 1993-12-30 | CT system with twin fan beam helical scan |
Publications (1)
Publication Number | Publication Date |
---|---|
DE4446819A1 true DE4446819A1 (de) | 1995-07-06 |
Family
ID=22641294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE4446819A Withdrawn DE4446819A1 (de) | 1993-12-30 | 1994-12-27 | CT System mit wendelförmiger Abtastung durch Doppelfächerbündel |
Country Status (4)
Country | Link |
---|---|
US (1) | US5469487A (de) |
JP (1) | JPH07250832A (de) |
DE (1) | DE4446819A1 (de) |
IL (1) | IL111982A (de) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5966422A (en) * | 1992-07-20 | 1999-10-12 | Picker Medical Systems, Ltd. | Multiple source CT scanner |
US5606585A (en) * | 1995-12-21 | 1997-02-25 | General Electric Company | Methods and apparatus for multislice helical image reconstruction in a computer tomography system |
US6275562B1 (en) | 1998-04-28 | 2001-08-14 | General Electric Company | Apparatus and methods for performing scalable multislice computed tomography scan |
US6198791B1 (en) | 1998-08-25 | 2001-03-06 | General Electric Company | Scalable multislice imaging system |
US6141398A (en) | 1998-08-25 | 2000-10-31 | General Electric Company | Protocol driven image reconstruction, display, and processing in a multislice imaging system |
WO2000038576A1 (en) | 1998-12-30 | 2000-07-06 | General Electric Company | Image thickness selection for multislice imaging system |
US6421412B1 (en) * | 1998-12-31 | 2002-07-16 | General Electric Company | Dual cardiac CT scanner |
US6466645B1 (en) * | 2000-11-22 | 2002-10-15 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for tube-spit correction |
JP3870105B2 (ja) * | 2002-02-22 | 2007-01-17 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 逆投影方法およびx線ct装置 |
US6816565B1 (en) * | 2003-06-27 | 2004-11-09 | Ge Medical Systems Global Technology Company, Llc | Matched view weighting methods and apparatus using multiple tilted reconstruction planes |
DE102004030550A1 (de) * | 2004-06-24 | 2006-01-19 | Siemens Ag | Bildgebendes Tomographiegerät mit zumindest zwei unter Systemwinkel angeordneten Aufnahmesystemen und Verfahren für ein derartiges Tomographiegerät zur Bestimmung der Systemwinkel der Aufnahmesysteme |
US7929664B2 (en) | 2007-02-13 | 2011-04-19 | Sentinel Scanning Corporation | CT scanning and contraband detection |
CN100462051C (zh) * | 2007-05-24 | 2009-02-18 | 上海交通大学 | 多源螺旋ct并行重建系统 |
WO2010141101A1 (en) | 2009-06-05 | 2010-12-09 | Sentinel Scanning Corporation | Transportation container inspection system and method |
JP6486100B2 (ja) * | 2014-12-22 | 2019-03-20 | キヤノン株式会社 | 画像処理装置、画像処理方法、およびプログラム |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4725963A (en) * | 1985-05-09 | 1988-02-16 | Scientific Measurement Systems I, Ltd. | Method and apparatus for dimensional analysis and flaw detection of continuously produced tubular objects |
US5170346A (en) * | 1989-11-22 | 1992-12-08 | General Electric Company | Method for reducing patient translation artifacts in tomographic imaging |
US5216601A (en) * | 1989-11-22 | 1993-06-01 | General Electric Company | Method for fan beam helical scanning using rebinning |
US5265142A (en) * | 1992-05-08 | 1993-11-23 | General Electric Company | Image reconstruction technique for a computer tomography system |
-
1993
- 1993-12-30 US US08/175,703 patent/US5469487A/en not_active Expired - Lifetime
-
1994
- 1994-12-14 IL IL111982A patent/IL111982A/xx not_active IP Right Cessation
- 1994-12-26 JP JP6321103A patent/JPH07250832A/ja active Pending
- 1994-12-27 DE DE4446819A patent/DE4446819A1/de not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
US5469487A (en) | 1995-11-21 |
IL111982A (en) | 1997-11-20 |
JPH07250832A (ja) | 1995-10-03 |
IL111982A0 (en) | 1995-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69030555T2 (de) | Rechnergesteuerte Tomographie mit schraubenförmiger Abtastung | |
DE4016245C2 (de) | Verfahren zum Abtasten und Sammeln von Daten | |
DE69028999T2 (de) | Rechnergesteuertes Tomographiegerät | |
DE69033923T2 (de) | Spiralabtastrechnertomographie | |
DE4139150C1 (en) | Computer tomograph with part ring formed X=ray source and detector - has double ring system without complementary interpolation | |
DE102007049469B4 (de) | Röntgentomographievorrichtung und Artefaktreduktionsverfahren | |
DE69902326T2 (de) | Bildkorrektur für ein geneigtes Gehäuse in einem Mehrschicht-Computertomograph | |
DE69126243T2 (de) | Computertomograph mit Scaneinrichtung | |
DE69218839T2 (de) | Tomographische Bilderzeugung mittels Schrägstrahlen | |
DE102005009263A1 (de) | Verfahren und System zur Bildgebung unter Verwendung mehrerer versetzter Röntgenstrahlen-Emissionspunkte | |
DE102004054405A1 (de) | Verfahren und Vorrichtung zum Vermindern von Artefakten in Computertomographiebildgebungssystemen | |
DE10211581A1 (de) | Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen Interpolationsalgorithmus | |
DE10251448A1 (de) | Verfahren für die Computertomographie eines periodisch sich bewegenden Untersuchungsobjektes, sowie ein CT-Gerät zur Durchführung dieses Verfahrens | |
DE4446819A1 (de) | CT System mit wendelförmiger Abtastung durch Doppelfächerbündel | |
DE102007039573A1 (de) | Verfahren zur analytischen Rekonstruktion für eine Mehrfachquellen-Inversgeometrie-CT | |
DE102005005839A1 (de) | Verfahren und Vorrichtung zur Verminderung von Artefakten bei Kegelstrahl-CT-Bildrekonstruktionen | |
DE10244180B4 (de) | Verfahren zur Bilderstellung in der Computertomographie eines periodisch bewegten Untersuchungsobjektes und CT-Gerät zur Durchführung des Verfahrens | |
DE10159927B4 (de) | Verfahren zur Bildrekonstruktion für die Computertomographie | |
DE3342075C2 (de) | Einrichtung zum Umordnen von Daten für tomographische Schichtbildern | |
DE19513052A1 (de) | Kalibrierung der Verstärkung eines Detektorkanals mittels Wobbelns des Brennpunktes | |
DE10245116A1 (de) | Verfahren zur Erzeugung eines Bildes mittels eines tomographiefähigen Röntgengeräts mit mehrzeiligem Röntgendetektorarray | |
DE60019537T2 (de) | Halbabtastungsalgorithmus für eine schnelle mehrreihige Detektoranordnung | |
DE69821082T2 (de) | Abbildungsvorrichtung und Verfahren für Computertomographie | |
DE102007021023A1 (de) | Verfahren zur Bilderstellung für die Spiral-CT mit veränderlichem Pitch und CT-Gerät zur Durchführung des Verfahrens | |
DE10127269A1 (de) | Verfahren für die Computertomographie sowie Computertomographie (CT)-Gerät |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8110 | Request for examination paragraph 44 | ||
8139 | Disposal/non-payment of the annual fee |