DE10211581A1 - Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen Interpolationsalgorithmus - Google Patents
Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen InterpolationsalgorithmusInfo
- Publication number
- DE10211581A1 DE10211581A1 DE10211581A DE10211581A DE10211581A1 DE 10211581 A1 DE10211581 A1 DE 10211581A1 DE 10211581 A DE10211581 A DE 10211581A DE 10211581 A DE10211581 A DE 10211581A DE 10211581 A1 DE10211581 A1 DE 10211581A1
- Authority
- DE
- Germany
- Prior art keywords
- detector
- line
- conjugate
- reconstruction plane
- spiral
- 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.)
- Ceased
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 title description 6
- 238000003384 imaging method Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims description 15
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000013170 computed tomography imaging Methods 0.000 abstract description 13
- 230000006870 function Effects 0.000 description 17
- 238000005259 measurement Methods 0.000 description 13
- 230000005855 radiation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 229910052704 radon Inorganic materials 0.000 description 3
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 230000000007 visual effect 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]
-
- 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
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
Es werden Verfahren und eine Vorrichtung zur Verringerung von Bildartefakten bei der Rekonstruktion eines Bilds mit einer Mehrfachschnitt-Computertomographieabbildungsabtasteinrichtung (Mehrfachschnitt-TT-Abbildungsabtasteinrichtung) (10) bereitgestellt. Ein Objekt (22) wird spiralförmig abgetastet, um eine Vielzahl von Projektionsansichten eines Objekts zu erhalten (52), es wird eine Rekonstruktionsebene (70) definiert (54), es wird eine konjugierte Abtastlinie (72) definiert (56), die die Rekonstruktionsebene nicht kreuzt, es wird eine Gewichtungsfunktion auf die konjugierten Abtastwerte angewendet (58), um Bildartefakte zu verringern, und das Bild wird rekonstruiert (60), nachdem die gewichteten Daten gefiltert und rückprojiziert sind. Das Verfahren verringert potentielle Bildartefakte, ohne daß zusätzliche Hardware verwendet oder ersetzt werden muß.
Description
Die Erfindung bezieht sich allgemein auf Verfahren und eine
Vorrichtung zur CT-Abbildung und andere
Strahlungsabbildungssysteme und insbesondere auf eine
Verwendung eines verallgemeinerten spiralenförmigen
Interpolationsalgorithmus.
Bei zumindest einigen "Computertomographie"-
Abbildungssystemkonfigurationen (CT-
Abbildungssystemkonfigurationen) projiziert eine
Röntgenquelle ein fächerförmiges Strahlenbündel, das
kollimiert wird, damit es in einer X-Y-Ebene eines
kartesischen Koordinatensystems liegt, auf die im
allgemeinen als eine "Abbildungsebene" Bezug genommen wird.
Das Röntgenstrahlenbündel geht durch ein abgebildetes
Objekt wie beispielsweise einen Patienten hindurch. Das
Strahlenbündel trifft nach seiner Abschwächung durch das
Objekt auf eine nachstehend als Array bezeichnete
regelmäßige Anordnung von nachstehend als
Strahlungsdetektoren bezeichneten
Strahlungserfassungseinrichtungen auf. Die Intensität der
bei einem Detektorarray empfangenen Strahlung des
abgeschwächten Strahlenbündels hängt von der Abschwächung
des Röntgenstrahlenbündels durch das Objekt ab. Jedes
Detektorelement des Arrays erzeugt ein getrenntes
elektrisches Signal, bei dem es sich um eine Messung der
Strahlenbündelabschwächung an dem Ort des Detektors
handelt. Die Abschwächungsmessungen von allen Detektoren
werden zur Erzeugung eines Übertragungsprofils getrennt
erfaßt.
Bei bekannten CT-Systemen der dritten Generation werden die
Röntgenquelle und das Detektorarray mit einem Portal in der
Abbildungsebene und um das abzubildende Objekt gedreht, so
daß der Winkel, bei dem das Röntgenstrahlenbündel das
Objekt kreuzt, sich ständig ändert. Röntgenquellen umfassen
typischerweise Röntgenröhren, die das Röntgenstrahlenbündel
bei einem Brennfleck abstrahlen. Röntgendetektoren umfassen
typischerweise einen Kollimator zur Kollimation von bei dem
Detektor empfangenen Röntgenstrahlenbündeln, einen an den
Kollimator angrenzenden Szintillator und an den
Szintillator angrenzende Fotodetektoren. Auf eine Gruppe
von Röntgenabschwächungsmessungen, d. h. Projektionsdaten,
von dem Detektorarray bei einem Portalwinkel wird als eine
"Ansicht" Bezug genommen. Eine "Abtastung" des Objekts
umfaßt einen Satz von während einer Umdrehung der
Röntgenquelle und des Detektors bei unterschiedlichen
Portalwinkeln oder Ansichtswinkeln ausgebildeten Ansichten.
Bei einer axialen Abtastung werden die Projektionsdaten zum
Aufbau eines Bilds, das einem durch das Objekt genommenen
zweidimensionalen Schnitt entspricht, verarbeitet. Auf ein
Verfahren zur Rekonstruktion eines Bilds aus einem Satz von
Projektionsdaten wird in dem Fachgebiet als das Verfahren
der gefilterten Rückprojektion Bezug genommen. Dieser
Prozeß wandelt die Abschwächungsmessungen aus einer
Abtastung in als "CT-Zahlen" oder "Hounsfield-Einheiten"
bezeichnete ganze Zahlen, die zur Steuerung der Helligkeit
eines entsprechenden Bildelements bei einer
Kathodenstrahlröhrenanzeige verwendet werden.
Zur Verringerung der für mehrere Schnitte erforderlichen
gesamten Abtastzeit kann eine "spiralenförmige" Abtastung
durchgeführt werden. Zur Durchführung einer
"spiralenförmigen" Abtastung wird der Patient synchron zu
der Drehung des Portals in einer Richtung der z-Achse
bewegt, während die Daten für eine vorgeschriebene Anzahl
von Schnitten erfaßt werden. Ein derartiges System erzeugt
aus einer spiralenförmigen Abtastung des
Fächerstrahlenbündels eine einzelne Spirale. Die durch das
Fächerstrahlenbündel ausgearbeitete Spirale ergibt
Projektionsdaten, aus denen Bilder in jedem
vorgeschriebenen Schnitt rekonstruiert werden können.
Bei einigen CT-Abbildungssystemen ist das Detektorarray
segmentiert, so daß eine Vielzahl von quasi parallelen
Schnitten von Projektionsdaten erfaßt und verarbeitet wird,
um eine Vielzahl von mehreren Schnitten durch ein Volumen
entsprechenden Bildern aufzubauen. Auf derartige CT-
Abbildungssysteme wird als "Mehrfachschnitt"-Systeme Bezug
genommen. Mehrfachschnitt-Systeme stellen im Vergleich zu
einem Einfachschnitt-System mehr Daten zur Rekonstruktion
eines Bilds bereit. Eine Gewichtung von konjugierten
Abtastwerten bei Mehrfachschnitt-Systemen ähnlich der
Gewichtung von Abtastwerten bei Einfachschnitt-Systemen
führt zu einer Störung bei der Bildrekonstruktion, wenn
eine große Anzahl von Schnitten berücksichtigt wird.
Für ein Einfachschnitt-System werden zwei Messungen zur
Interpolation/Extrapolation von aus einer einzelnen Drehung
des Portals gesammelten Daten auf eine
"Rekonstruktionsebene" (POR) verwendet. Bei
unterschiedlichen Quellenpositionen erfaßte Messungen sind
als "konjugierte Messungen" bekannt. Die konjugierten
Messungen werden zur Erzeugung einer konjugierten
Abtastlinie verwendet. Mehrfachschnitt-CT-Abbildungssysteme
erhöhen die Flexibilität bei der Wahl einer Zeile oder
einer Vielzahl von Zeilen, die am nächsten an der POR
befindliche Daten bereitstellen. Diese Flexibilität bei der
Wahl einer Zeile kann an beiden Enden eines gegebenen
Strahls vorhanden sein, d. h. für beide konjugierten
Messungen für den gleichen Strahl durch den Patienten.
Teilungen, für die konjugierte Messungen verfügbar sind,
sind als HQ-Teilungen bekannt.
Die POR-Linie kreuzt die konjugierte Abtastlinie an einem
Schnittpunkt. An dem Schnittpunkt wird ein Gewicht
nichtdeterministisch, da das Gewicht bei der POR-Linie 1
sein sollte und das Gewicht für die konjugierte Abtastlinie
0 sein sollte. Bekannte Algorithmen schalten die
Interpolationspaare dynamisch um, so daß zwei sich am
nächsten an der POR befindende Abtastwerte für die
Interpolation verwendet werden, um die Nicht-Deterministik
zu vermeiden. An dem Umschaltpunkt tritt jedoch eine
Unstetigkeit bei den Gewichten auf. Diese Unstetigkeit
führt nach dem Filtern und der Rückprojektion zu einem
periodischen Riffelmuster von dem mittleren Bereich des
rekonstruierten Bilds zu der äußeren Kante des
Gesichtsfelds. Typischerweise wird als ein
Kompensationsverfahren eine Verlauftechnik zur Überbrückung
von Lücken bei den Gewichten verwendet. Dennoch können
Bildartefakte nicht vermieden werden, da es sich bei der
Filteroperation bei der Bildrekonstruktion im wesentlichen
um einen Ableitungsoperator handelt, der etwaige kleine
Diskrepanzen bei den interpolierten Abtastwerten verstärkt.
Es werden Verfahren und eine Vorrichtung für ein
Mehrfachschnitt-Computertomographieabbildungssystem
(Mehrfachschnitt-CT-Abbildungssystem) zur Verringerung von
Bildartefakten bei der Rekonstruktion eines Bilds
bereitgestellt. Bei einem beispielhaften
Ausführungsbeispiel des Verfahrens wird ein Objekt
spiralenförmig abgetastet, um Projektionsdaten für eine
Vielzahl von Projektionsansichten des Objekts zu erhalten,
es wird eine Rekonstruktionsebene definiert, es wird eine
konjugierte Abtastlinie definiert, die die
Rekonstruktionsebene nicht kreuzt, es wird eine
Gewichtungsfunktion auf die konjugierten Abtastwerte
angewendet, um Bildartefakte zu verringern, und das Bild
wird rekonstruiert, nachdem die gewichteten Daten gefiltert
und rückprojiziert sind.
Bei einer weiteren Ausgestaltung umfaßt ein
Abbildungssystem einen Computer, ein Portal mit einem
Detektorarray, eine Röntgenquelle zur Ausstrahlung eines
Röntgenstrahlenbündels zu dem Detektorarray hin, und das
Abbildungssystem erfaßt Projektionsdaten für eine Vielzahl
von Projektionsansichten des Objekts, definiert eine
Rekonstruktionsebene, definiert eine konjugierte
Abtastlinie, die die Rekonstruktionsebene nicht kreuzt,
wendet eine Gewichtungsfunktion auf die konjugierten
Abtastwerte an, um Bildartefakte zu verringern, und
rekonstruiert das Bild, nachdem die gewichteten Daten
gefiltert und rückprojiziert sind.
Bei einer weiteren Ausgestaltung ist eine
Verarbeitungseinrichtung in dem Abbildungssystem dazu
programmiert, Projektionsdaten für eine Vielzahl von
Projektionsansichten des Objekts zu erfassen, eine
Rekonstruktionsebene zu definieren, eine konjugierte
Abtastlinie derart zu definieren, daß sie die
Rekonstruktionsebene nicht kreuzt, eine Gewichtungsfunktion
auf die konjugierten Abtastwerte anzuwenden, um
Bildartefakte zu verringern, und das Bild zu
rekonstruieren, nachdem die gewichteten Daten gefiltert und
rückprojiziert sind.
Bei noch einer weiteren Ausgestaltung wird ein von einem
Computer lesbarer Träger in dem Abbildungssystem
bereitgestellt, der eine Vielzahl von Datensätzen von
Projektionsdaten umfaßt, die zur Definition einer
Rekonstruktionsebene verwendet werden. Darüber hinaus ist
ein Datensatz von konjugierten Abtastwerten auf dem von
einem Computer lesbaren Träger gespeichert. Ein sich auf
dem von einem Computer lesbaren Träger befindendes Programm
verwendet eine Vielzahl von Regeln zur Definition einer
konjugierten Abtastlinie auf der Grundlage des Datensatzes
von konjugierten Abtastwerten, die die Rekonstruktionsebene
nicht kreuzt. Ferner umfaßt das Programm eine Vielzahl von
Regeln zur Bestimmung einer Gewichtungsfunktion, die zur
Verringerung von Bildartefakten auf die konjugierten
Abtastwerte angewendet wird, und zur Rekonstruktion des
Bilds.
Fig. 1 zeigt eine bildliche Ansicht eines CT-
Abbildungssystems;
Fig. 2 zeigt ein schematisches Blockschaltbild des in der
Fig. 1 veranschaulichten Systems;
Fig. 3 zeigt ein Flußdiagramm, das die durch das CT-System
zur Verwendung eines verallgemeinerten spiralenförmigen
Interpolationsalgorithmus ausgeführten Schritte
veranschaulicht;
Fig. 4 zeigt eine Veranschaulichung der POR und
konjugierter Abtastwerte zu der POR, die von dem
verallgemeinerten spiralenförmigen
Interpolationsalgorithmus verwendet werden;
Fig. 5 zeigt eine Veranschaulichung einer stückweise
linearen POR;
Fig. 6 zeigt eine Veranschaulichung einer allgemein
kurvenförmigen POR; und
Fig. 7 zeigt eine Veranschaulichung konjugierter
Abtastbereiche für eine Spiralenteilung von 7 : 1.
Mit Bezug auf die Fig. 1 und 2 ist stellvertretend für
eine CT-Abtasteinrichtung bzw. einen CT-Abtaster der
"dritten Generation" ein
Computertomographieabbildungssystem (CT-Abbildungssystem)
10 als ein Portal 12 umfassend gezeigt. Das Portal 12 weist
eine Röntgenquelle 14 auf, die ein Strahlenbündel von
Röntgenstrahlen 16 zu einem Detektorarray 18 auf der
gegenüberliegenden Seite des Portals 12 hin projiziert. Das
Detektorarray 18 ist durch Detektorelemente 20 gebildet,
die zusammen die durch ein Objekt 22 wie beispielsweise
einen medizinischen Patienten hindurchgehenden projizierten
Röntgenstrahlen erfassen. Jedes Detektorelement 20 erzeugt
ein elektrisches Signal, das die Intensität eines
auftreffenden Röntgenstrahlenbündels und folglich die
Abschwächung des Strahlenbündels, während es durch das
Objekt oder den Patienten 22 hindurchgeht, darstellt.
Während einer Abtastung zur Erfassung von
Röntgenprojektionsdaten drehen sich das Portal 12 und die
daran angebrachten Komponenten um einen Drehmittelpunkt 24.
Bei einem Ausführungsbeispiel und wie es in der Fig. 2
gezeigt ist, sind Detektorelemente 20 in einer Zeile
angeordnet, so daß während einer Abtastung einem einzelnen
Bildschnitt entsprechende Projektionsdaten erfaßt werden.
Bei einem weiteren Ausführungsbeispiel sind
Detektorelemente 20 in einer Vielzahl von parallelen Zeilen
angeordnet, so daß während einer Abtastung einer Vielzahl
von parallelen Schnitten entsprechende Projektionsdaten
gleichzeitig erfaßt werden können.
Die Drehung des Portals 12 und der Betrieb der
Röntgenquelle 14 werden durch eine Steuervorrichtung 26 des
CT-Systems 10 gesteuert. Die Steuervorrichtung 26 umfaßt
eine Röntgensteuereinrichtung 28, die Leistungs- und
Zeitablaufsteuersignale für die Röntgenquelle 14
bereitstellt, und eine Portalmotorsteuereinrichtung 30, die
die Drehgeschwindigkeit und -position des Portals 12
steuert. Ein Datenerfassungssystem (DAS) 32 in der
Steuervorrichtung 26 tastet analoge Daten von
Detektorelementen 20 ab und wandelt die Daten in digitale
Signale für eine nachfolgende Verarbeitung. Eine
Bildrekonstruktionseinrichtung 34 empfängt abgetastete und
digitalisierte Röntgendaten von dem DAS 32 und führt eine
Hochgeschwindigkeitsbildrekonstruktion durch. Das
rekonstruierte Bild wird einem Computer 36 als Eingabe
zugeführt, der das Bild in einer Massenspeichervorrichtung
38 speichert.
Der Computer 36 empfängt auch Befehle und Abtastparameter
von einer Bedienungsperson über eine Bedieneinheit 40, die
eine Tastatur aufweist. Eine zugehörige
Kathodenstrahlröhrenanzeige 42 ermöglicht es der
Bedienungsperson, das rekonstruierte Bild und andere Daten
von dem Computer 36 zu beobachten. Die von der
Bedienungsperson zugeführten Befehle und Parameter werden
von dem Computer 36 zur Bereitstellung von Steuersignalen
und Informationen für das DAS 32, die
Röntgensteuereinrichtung 28 und die
Portalmotorsteuereinrichtung 30 verwendet. Darüber hinaus
betätigt der Computer 36 eine Tischmotorsteuereinrichtung
44, die einen motorisierten Tisch 46 zur Positionierung des
Patienten 22 in dem Portal 12 steuert. Insbesondere bewegt
der Tisch 46 Abschnitte des Patienten 22 durch eine
Portalöffnung 48.
Wenn eine HQ-Teilung verwendet wird, heißt es, daß ein CT-
Abbildungssystem in einer Betriebsart "hoher Qualität" (HQ-
Betriebsart) arbeitet. Für eine gegebene Spiralenteilung in
einer HQ-Betriebsart ist ein Satz von Radon-Punkten
vorhanden, für den die konjugierten Messungen um annähernd
die Hälfte einer erfaßten Apertur z versetzt sind. Die
"Detektor-Apertur" bezieht sich auf eine Apertur in einem
mit einer "Makrodetektorzeile" oder einer Messung
verknüpften Isozentrum. Beispielsweise sind bei einer in
einer HQ-Betriebsart arbeitenden bekannten Vier-Schnitt-
Abtasteinrichtung mit einem Detektorwinkel ã gleich Null
und einer Teilung von 3 : 1 die konjugierten Messungen um
Äz/2, d. h. die Hälfte der Apertur, versetzt. In einer
derartigen Situation wird häufig eine Interpolationsbreite
mit minimalem z verwendet. Diese ideale Abtastsituation
tritt jedoch lediglich bei einer begrenzten Anzahl von
Radon-Punkten auf. Der Versatz weicht allmählich von dem
Optimum ab, je weiter die Entfernung von dem Isozentrum
ist. Je größer die Anzahl der Zeilen und je höher
entsprechend die Teilung ist, desto schneller findet die
Abweichung von dem Optimum statt.
Ein in einer HQ-Betriebsart arbeitendes bekanntes Acht-
Schnitt-CT-Abbildungssystem mit einer Spiralenteilung von
entweder 5 : 1 oder 7 : 1, das bekannte
Rekonstruktionsverfahren verwendet, erzeugt Bildartefakte
in einem sich zu einer äußeren Kante eines "Gesichtsfelds"
(FOV) in dem Bild erstreckenden periodischen Riffelmuster,
das die Bilder für eine klinische Verwendung ungeeignet
macht. Die Bildartefakte ergeben sich aus Änderungen der
spiralenförmigen Gewichte als einer Funktion eines
Detektorwinkels y. Die Änderung der Gewichte wird durch die
Kombination einer schnelleren Tischgeschwindigkeit bei
einer Acht-Schnitt-Abtasteinrichtung und eines durch die
konjugierten Abtastwerte zu der POR gebildeten großen
Neigungswinkels verursacht. Bei einem bekannten
Rekonstruktionsalgorithmus wird die POR als eine zu der
Tischtranslationsachse senkrechte Ebene ausgewählt, und die
POR kann durch die nachstehende lineare Beziehung
beschrieben werden:
b = βk, (1)
wobei βk den Projektionswinkel darstellt, bei dem die
Detektorzeile k die POR kreuzt. Auf der Grundlage der
Fächerstrahlenbündelabtastgeometrie befinden sich die
konjugierten Abtastwerte entlang der Linie:
β' = βk±π-2γ, (2)
wobei βk den Projektionswinkel darstellt, bei dem die
Detektorzeile k die POR kreuzt, und γ ein Detektorwinkel
ist. Bei der Gleichung (2) handelt es sich um die
konjugierte Abtastlinie β', die die konjugierten
Abtastwerte zu der POR für alle Detektorzeilen beschreibt,
die einen Satz von geneigten parallelen Linien bilden. Bei
einer spiralenförmigen Mehrfachschnitt-CT werden die
konjugierten Abtastwerte zu der POR der Detektorzeile k mit
anderen Detektorzeilen als k gesammelt.
Die Fig. 3 zeigt ein Flußdiagramm 50, das die durch das
(in der Fig. 1 gezeigte) CT-System 10 zur Verringerung von
durch nichtdeterministische Schnittpunkte erzeugten
Bildartefakten ausgeführten Schritte veranschaulicht. Das
durch das Flußdiagramm 50 in der Fig. 3 veranschaulichte
Verfahren kann durch das (in der Fig. 2 gezeigte) DAS 32,
die (in der Fig. 2 gezeigte)
Bildrekonstruktionseinrichtung 34 oder den (in der Fig. 2
gezeigten) Computer 36 praktiziert werden. Im allgemeinen
ist eine Verarbeitungseinrichtung in dem DAS 32, der
Rekonstruktionseinrichtung 34 und/oder dem Computer 36 dazu
programmiert, die nachstehend beschriebenen Prozeßschritte
auszuführen. Selbstverständlich ist das Verfahren nicht
darauf beschränkt, in dem CT-System 10 praktiziert zu
werden, und kann in Verbindung mit vielen anderen Arten und
Variationen von Abbildungssystemen verwendet werden.
Mit näherem Bezug auf die Fig. 3 wird bei der Durchführung
einer Abtastung ein Satz von Ursprungsprojektionsdaten
erfaßt 52, der zur Erzeugung eines Satzes von
Projektionsansichtsdaten vorverarbeitet wird. Es wird eine
Rekonstruktionsebene definiert 54. Es wird eine konjugierte
Abtastlinie definiert 56, die die Rekonstruktionsebene
nicht kreuzt. Wenn die POR und die konjugierte Abtastlinie
definiert werden, wird ein minimaler Wert von α ausgewählt,
um es sicherzustellen, daß die POR die konjugierte
Abtastlinie nicht kreuzt. Wenn die POR die konjugierte
Abtastlinie nicht kreuzt, werden Bildartefakte beseitigt.
Der Wert von α sollte eine glatte Gewichtungsfunktion
sicherstellen. Die Gewichtungsfunktion sollte keine
Unstetigkeiten aufweisen und wird auf alle konjugierten
Abtastwerte angewendet 58. Es wird ein Bild rekonstruiert
60, indem die gewichteten Daten gefiltert und
rückprojiziert werden.
Mit Bezug auf die Fig. 4 bilden eine POR-Linie 62 und eine
konjugierte Abtastlinie 64 ein Dreieck mit einer Basislänge
68, das sich an einem Schnittpunkt 66 kreuzt. Wie es
veranschaulicht ist, stellt die Basislänge 68 einen
Winkelabstand zwischen der POR-Linie 62 und der
konjugierten Abtastlinie 64 dar. Es ist für eine
Spiralenteilung von 7 : 1 bekannt, daß der Winkelabstand für
einen Isostrahl π/7 beträgt. Ähnlich sind für andere
Spiralenteilungen verschiedene Winkelabstände zwischen der
POR und der konjugierten Abtastlinie vorhanden. Zur
Vermeidung von durch den Schnittpunkt der POR-Linie 62 und
der konjugierten Abtastlinie 64 verursachten Bildartefakten
wird bei einem Ausführungsbeispiel eine "neue" POR-Linie 70
als eine gerade Linie definiert, die eine konjugierte
Abtastlinie 72 nicht kreuzt. Genauer ist die neue POR-Linie
70 mit Bezug auf eine horizontale Linie in dem Radon-Raum
eine geneigte Linie. Durch eine Unstetigkeit der Gewichte
während der Interpolation verursachte potentielle
Bildartefakte an dem Schnittpunkt 66 werden vermieden, da
die neue POR-Linie 70 und die konjugierte Abtastlinie 72
sich nicht kreuzen. Mathematisch werden die neue POR-Linie
70, β, und die konjugierte Abtastlinie 72, β', durch die
nachstehenden Gleichungen beschrieben:
β = βk-αγ, (3)
β' = βk±π-(2-α)γ, (4)
wobei βk den Projektionswinkel darstellt, bei dem die
Detektorzeile k die POR kreuzt, wenn γ = 0, γ einen
Detektorwinkel darstellt und α einen Wert darstellt, der es
sicherstellt, daß die neue POR-Linie 70 und die konjugierte
Abtastlinie 72 sich nicht kreuzen.
Der minimale Wert von α, der es sicherstellt, daß die zwei
Linien sich nicht kreuzen, ist α ≧1-π(2Qγm), wobei Q eine
Spiralenteilung ist. Bei einem Ausführungsbeispiel für eine
Spiralenteilung von 7 : 1 ist α ≧ 1-π/(14γm), wobei γm ein
maximaler Detektorwinkel ist. Für eine Spiralenteilung von
Q : 1 ist der minimale Wert von α
Der Wert von α
wird auf einem Minimum gehalten, da große Werte von α eine
große Abweichung der POR von einer idealen POR angeben. Die
ideale POR ist eine flache Ebene senkrecht zu der
Tischtranslationsachse. Ferner sollte der Wert von α einen
glatten Gewichtungsfunktionsübergang sicherstellen, um
potentielle Bildartefakte zu vermeiden. Bei einem
bevorzugten Ausführungsbeispiel ist α für eine
Spiralenteilung von 7 : 1 das 1,1-fache des minimalen Werts
von 1-π/(14γm), und für eine Spiralenteilung von Q : 1 ist α
das 1,1-fache des minimalen Werts von
Für andere
Spiralenteilungen wird eine ähnliche Berechnung
durchgeführt, um den minimalen Wert von α zu schätzen.
Mit näherem Bezug auf die Fig. 5 sind eine weitere POR-
Linie als eine stückweise POR-Linie 80 und eine konjugierte
Abtastlinie 82 definiert. Bei einem Ausführungsbeispiel
sind die stückweise POR-Linie 80 und die konjugierte
Abtastlinie 82 stückweise gerade Linien in dem Sinogramm-
Raum, die sich nicht kreuzen. Die stückweise POR-Linie 80
umfaßt einen Abschnitt 84 einer idealen POR-Linie 86. Für
eine Spiralenteilung von 7 : 1 ist ein Winkelabstand 88
zwischen der stückweisen POR-Linie 80 und der konjugierten
Abtastlinie 82 für einen Isostrahl π/7.
Mit näherem Bezug auf die Fig. 6 sind noch eine weiter
POR-Linie als eine kurvenförmige POR-Linie 90 und eine
konjugierte Abtastlinie 92 definiert. Die kurvenförmige
POR-Linie 90 und die konjugierte Abtastlinie 92 sind bei
einem Ausführungsbeispiel Kurven in dem Sinogramm-Raum, die
sich nicht kreuzen. Die kurvenförmige POR-Linie 90 umfaßt
einen Abschnitt 94 einer idealen POR-Linie 96. Für eine
Spiralenteilung von 7 : 1 ist ein Winkelabstand 98 zwischen
der kurvenförmigen POR-Linie 90 und der konjugierten
Abtastlinie 92 für einen Isostrahl π/7.
Bei der Auswahl einer POR besteht ein Ziel darin, die
Verwendung der idealen POR mit flacher Ebene zu minimieren.
Die Form der POR kann auf eine spezielle klinische
Anwendung zugeschnitten werden. Beispielsweise ist in dem
Fall einer Kopfabtastung das mittlere "Gesichtsfeld" (FOV)
von 25 cm der zu rekonstruierenden Bereich von Interesse.
Daher sollte eine POR mit einem Design einer stückweisen
geraden Linie ausgewählt werden, um die Abweichung der POR
für das mittlere FOV von 25 cm zu minimieren. Falls die
Anwendung ein ganzes FOV von 50 cm erfordert, sollte
demgegenüber eine POR mit einem Design einer einzelnen
geraden Linie ausgewählt werden.
Sobald die POR und die konjugierten Abtastwerte zu der POR
bestimmt sind, wird eine Vielzahl der spiralenförmigen
Gewichte bestimmt. Mit näherem Bezug auf die Fig. 7 sind
eine Vielzahl von konjugierten Abtastbereichen 102 für eine
Spiralenteilung von 7 : 1 und eine gegen eine Vielzahl von
Projektionswinkeln 106 aufgezeichnete Vielzahl von Zeilen
104 gezeigt. Eine Vielzahl von konjugierten Bereichen ist
durch identische Zahlen bezeichnet, z. B. ist der
konjugierte Bereich 108 mit der Zahl 5 bezeichnet. Es ist
unter Vernachlässigung des Kegelstrahlenbündeleffekts eine
Drehverzögerung von 2π zwischen einer ersten Detektorzeile
110 und einer letzten Detektorzeile 112 vorhanden. Den
Detektorzeilen entsprechend sind drei konjugierte
Abtastbereiche vorhanden.
Eine spiralenförmige Gewichtungsfunktion für jede Zeile
kann durch die nachstehende Beziehung beschrieben werden:
wobei
β0,k(γ,β) die POR für die Detektorzeile k ist, z. B. bei einem beispielhaften Ausführungsbeispiel
β0,k(γ,β) die POR für die Detektorzeile k ist, z. B. bei einem beispielhaften Ausführungsbeispiel
β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und
einen Wert auf der Detektorzeile k konjugiert zu der POR
einer weiteren Zeile darstellt, z. B. bei einem
beispielhaften Ausführungsbeispiel
β2,k(γ,β) für k = 0, . . . 7 ein Abtastwert auf der Detektorzeile
k konjugiert zu der POR einer weiteren Zeile ist, der nach
β0,k(γ,β) gesammelt wird, z. B. bei einem beispielhaften
Ausführungsbeispiel
und
ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist, die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1. Bei einem Ausführungsbeispiel ist keine Unstetigkeit bei der Gewichtungsfunktion vorhanden.
ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist, die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1. Bei einem Ausführungsbeispiel ist keine Unstetigkeit bei der Gewichtungsfunktion vorhanden.
Bei einem Ausführungsbeispiel wird eine einfache lineare
Interpolation für die Schätzung der spiralenförmigen
Gewichte verwendet. Bei einem alternativen
Ausführungsbeispiel können Interpolationsfunktionen einer
höheren Ordnung zur Stimmung spiralenförmiger Gewichte
ausgewählt werden. Beispielsweise wird eine auf wk(γ,β)
bezogene Gewichtungsfunktion θk(γ,β) durch die nachstehende
Beziehung beschrieben:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β) (6)
Bei der Auswahl einer Gewichtungsfunktion summiert sich der
gesamte Beitrag von allen konjugierten Abtastwerten zu
Eins, die Gewichtungsfunktion ist kontinuierlich, und die
Gewichtungsfunktion sollte in der γ-Richtung so glatt wie
möglich sein. Eine zusätzliche Glättung in der z-Richtung
verringert restliche Bildartefakte weiter.
Bei noch einem weiteren Ausführungsbeispiel umfaßt das CT-
System 10 ein sich auf einem von einem Computer lesbaren
Träger in dem Massenspeicher 38 befindendes
Computerprogramm zur Rekonstruktion des Bilds. Das Programm
umfaßt eine Vielzahl von Datensätzen von Projektionsdaten,
die zur Definition einer Rekonstruktionsebene verwendet
werden. Eine Vielzahl von Regeln definiert eine konjugierte
Abtastlinie auf der Grundlage eines Datensatzes von
konjugierten Abtastwerten. Das Programm verwendet eine
Vielzahl von Regeln, um es sicherzustellen, daß die
Rekonstruktionsebene die konjugierte Abtastlinie nicht
kreuzt. Darüber hinaus umfaßt das Programm eine Vielzahl
von Regeln zur Bestimmung einer Gewichtungsfunktion, die
zur Verringerung von Bildartefakten auf die konjugierten
Abtastwerte angewendet wird, und zur Rekonstruktion des
Bilds.
Während die Erfindung in Form von verschiedenen
spezifischen Ausführungsbeispielen beschrieben ist, erkennt
es der Fachmann, daß die Erfindung mit einer Modifikation
im Rahmen des Inhalts und Schutzbereichs der
Patentansprüche in die Praxis umgesetzt werden kann.
Es werden Verfahren und eine Vorrichtung zur Verringerung
von Bildartefakten bei der Rekonstruktion eines Bilds mit
einer Mehrfachschnitt-
Computertomographieabbildungsabtasteinrichtung
(Mehrfachschnitt-CT-Abbildungsabtasteinrichtung) (10)
bereitgestellt. Ein Objekt (22) wird spiralenförmig
abgetastet, um eine Vielzahl von Projektionsansichten eines
Objekts zu erhalten (52), es wird eine Rekonstruktionsebene
(70) definiert (54), es wird eine konjugierte Abtastlinie
(72) definiert (56), die die Rekonstruktionsebene nicht
kreuzt, es wird eine Gewichtungsfunktion auf die
konjugierten Abtastwerte angewendet (58), um Bildartefakte
zu verringern, und das Bild wird rekonstruiert (60),
nachdem die gewichteten Daten gefiltert und rückprojiziert
sind. Das Verfahren verringert potentielle Bildartefakte,
ohne daß zusätzliche Hardware verwendet oder ersetzt werden
muß.
Claims (31)
1. Verfahren zur Rekonstruktion zumindest eines ein Objekt
(22) darstellenden Bilds mit einem Abbildungssystem (10),
mit den Schritten:
Abtasten eines Objekts zur Erfassung (52) einer Vielzahl von Projektionsansichten des Objekts;
Definieren (54) einer Rekonstruktionsebene (70) derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) des Bilds nach einer Filterung und Rückprojektion der Daten.
Abtasten eines Objekts zur Erfassung (52) einer Vielzahl von Projektionsansichten des Objekts;
Definieren (54) einer Rekonstruktionsebene (70) derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) des Bilds nach einer Filterung und Rückprojektion der Daten.
2. Verfahren nach Anspruch 1, wobei der Schritt des
Definierens (54) einer Rekonstruktionsebene (70) ferner den
Schritt des derartigen Auswählens der Rekonstruktionsebene,
daß sie durch eine lineare Beziehung, eine stückweise
gerade Linie und/oder eine Kurve definiert ist, umfaßt.
3. Verfahren nach Anspruch 2, wobei die lineare
Rekonstruktionsebene (70) die Beziehung:
β = βk-αγ
umfaßt, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und eine konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β = βk-αγ
umfaßt, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und eine konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
4. Verfahren nach Anspruch 1, wobei die konjugierte
Abtastlinie (72) die Beziehung:
β' = βk ± π-(2-α)γ
umfaßt, wobei β' eine konjugierte Abtastlinie ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β' = βk ± π-(2-α)γ
umfaßt, wobei β' eine konjugierte Abtastlinie ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
5. Verfahren nach Anspruch 1, wobei der Schritt des
Anwendens (58) einer Gewichtungsfunktion den Schritt des
Definierens eines spiralenförmigen Gewichts gemäß der
Beziehung:
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
6. Verfahren nach Anspruch 5, wobei der Schritt des
Anwendens (58) einer Gewichtungsfunktion den Schritt des
Bestimmens einer Gewichtungsfunktion gemäß der Beziehung:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
7. Verfahren nach Anspruch 1, wobei das Abtasten eines
Objekts (22) den Schritt des spiralenförmigen Abtastens
eines Volumens des Objekts mit einer Teilung, die während
der Abtastung variiert, umfaßt.
8. Abbildungssystem (10) mit einem Computer (36), einem
Portal (12) mit einem Detektorarray (18), einer
Röntgenquelle (14) zur Ausstrahlung eines
Röntgenstrahlenbündels (16) entlang einer Abbildungsebene
zu einem Detektorarray einschließlich einer Vielzahl von
Detektorzellen (20) hin, wobei der Computer mit der
Röntgenquelle und dem Portal gekoppelt ist, wobei das
Abbildungssystem konfiguriert ist zum:
Abtasten eines Objekts (22) zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene (70) derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
Abtasten eines Objekts (22) zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene (70) derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
9. Abbildungssystem (10) nach Anspruch 8, wobei die
Rekonstruktionsebene (70) eine lineare Beziehung, eine
stückweise gerade Linie und/oder eine Kurve umfaßt.
10. Abbildungssystem (10) nach Anspruch 9, wobei das
Abbildungssystem ferner dazu konfiguriert ist, die lineare
Rekonstruktionsebene (70) gemäß der Beziehung:
β = βk-αγ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β = βk-αγ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
11. Abbildungssystem (10) nach Anspruch 8, wobei das
Abbildungssystem ferner dazu konfiguriert ist, die
konjugierte Abtastlinie (72) gemäß der Beziehung:
β' = βk±π-(2-α)γ
zu bestimmen, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β' = βk±π-(2-α)γ
zu bestimmen, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
12. Abbildungssystem (10) nach Anspruch 8, wobei die
Gewichtungsfunktion ein spiralenförmiges Gewicht gemäß der
Beziehung:
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist;
β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist;
β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
13. Abbildungssystem (10) nach Anspruch 8, wobei die
Gewichtungsfunktion die Beziehung:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
14. Abbildungssystem (10) nach Anspruch 8, ferner dazu
konfiguriert, ein Volumen des Objekts (22) mit einer
Teilung abzutasten, die während der Abtastung variiert.
15. Verarbeitungseinrichtung, die zur Verringerung von
Bildartefakten bei einem Computertomographiesystem (10)
programmiert ist, wobei die Verarbeitungseinrichtung
konfiguriert ist zum:
Abtasten eines Volumens eines Objekts (22) mit einer variierenden Teilung zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene mit einer mathematischen Beziehung derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene (70) nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
Abtasten eines Volumens eines Objekts (22) mit einer variierenden Teilung zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene mit einer mathematischen Beziehung derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene (70) nicht kreuzt;
Anwenden (58) einer Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
16. Verarbeitungseinrichtung nach Anspruch 15, wobei die
Verarbeitungseinrichtung ferner dazu konfiguriert ist, die
Rekonstruktionsebene (70) derart zu bestimmen, daß sie
einer stückweise geraden Linie und/oder einer Kurve
entspricht.
17. Verarbeitungseinrichtung nach Anspruch 15, wobei die
Verarbeitungseinrichtung ferner dazu konfiguriert ist, die
lineare Rekonstruktionsebene (70) gemäß der Beziehung:
β = βk-α,γ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β = βk-α,γ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
18. Verarbeitungseinrichtung nach Anspruch 15, wobei die
Verarbeitungseinrichtung ferner dazu konfiguriert ist, die
konjugierte Abtastlinie (72) gemäß der Beziehung:
β' = βk ±π-(2-α)γ
zu bestimmen, wobei β' die konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β' = βk ±π-(2-α)γ
zu bestimmen, wobei β' die konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
19. Verarbeitungseinrichtung nach Anspruch 15, wobei die
Verarbeitungseinrichtung ferner dazu konfiguriert ist, ein
Gewicht gemäß der Beziehung:
zu bestimmen, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
zu bestimmen, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
20. Verarbeitungseinrichtung nach Anspruch 19, wobei die
Verarbeitungseinrichtung ferner dazu konfiguriert ist, eine
Gewichtungsfunktion gemäß der Beziehung:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
anzuwenden (58), wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
anzuwenden (58), wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
21. Von einem Computer lesbarer Träger in einem
Abbildungssystem (10), wobei der von einem Computer lesbare
Träger umfaßt:
Datensätze von Projektionsdaten, die zur Definition (54) einer Rekonstruktionsebene verwendet werden;
eine Vielzahl von Regeln zur Definition einer Rekonstruktionsebene derart, daß eine konjugierte Abtastlinie auf der Grundlage einer Vielzahl von Datensätzen von konjugierten Abtastwerten die Rekonstruktionsebene nicht kreuzt;
eine Vielzahl von Regeln zur Bestimmung einer Gewichtungsfunktion, die auf die konjugierten Abtastwerte angewendet wird (58); und
eine Vielzahl von Regeln zur Rekonstruktion (60) eines Bilds.
Datensätze von Projektionsdaten, die zur Definition (54) einer Rekonstruktionsebene verwendet werden;
eine Vielzahl von Regeln zur Definition einer Rekonstruktionsebene derart, daß eine konjugierte Abtastlinie auf der Grundlage einer Vielzahl von Datensätzen von konjugierten Abtastwerten die Rekonstruktionsebene nicht kreuzt;
eine Vielzahl von Regeln zur Bestimmung einer Gewichtungsfunktion, die auf die konjugierten Abtastwerte angewendet wird (58); und
eine Vielzahl von Regeln zur Rekonstruktion (60) eines Bilds.
22. Von einem Computer lesbarer Träger nach Anspruch 21,
wobei die Rekonstruktionsebene (70) eine lineare Beziehung,
eine stückweise gerade Linie und/oder eine Kurve umfaßt.
23. Von einem Computer lesbarer Träger nach Anspruch 22,
wobei die lineare Rekonstruktionsebene (70) die Beziehung,
β = βk-αγ
umfaßt, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β = βk-αγ
umfaßt, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile k die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
24. Von einem Computer lesbarer Träger nach Anspruch 21,
wobei die konjugierte Abtastlinie (72) die Beziehung:
β' = βk±π-(2-α)γ
umfaßt, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β' = βk±π-(2-α)γ
umfaßt, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
25. Von einem Computer lesbarer Träger nach Anspruch 21,
wobei die Gewichtungsfunktion ein spiralenförmiges Gewicht
gemäß der Beziehung:
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist;
β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als b = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist;
β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als b = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
26. Von einem Computer lesbarer Träger nach Anspruch 21,
wobei die Gewichtungsfunktion die Beziehung:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
27. Abbildungssystem (10) mit einem Computer (36), einem
Portal (12) mit einem Detektorarray (18), einer
Röntgenquelle (14) zur Ausstrahlung eines
Röntgenstrahlenbündels (16) entlang einer Abbildungsebene
zu einem Detektorarray einschließlich einer Vielzahl von
Detektorzellen (20) hin, wobei der Computer mit der
Röntgenquelle und dem Portal gekoppelt ist, wobei das
Abbildungssystem konfiguriert ist zum:
Abtasten eines Objekts (22) zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene (70) mit einer linearen Beziehung, einer stückweise geraden Linie und/oder einer Kurve derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt, wenn eine mathematische Beziehung zwischen einem Projektionswinkel, Detektorzeilen und Detektorwinkeln verwendet wird;
Anwenden (58) einer spiralenförmigen Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
Abtasten eines Objekts (22) zur Erfassung (52) einer Vielzahl von Projektionsansichten eines Objekts;
Definieren (54) einer Rekonstruktionsebene (70) mit einer linearen Beziehung, einer stückweise geraden Linie und/oder einer Kurve derart, daß eine konjugierte Abtastlinie (72) die Rekonstruktionsebene nicht kreuzt, wenn eine mathematische Beziehung zwischen einem Projektionswinkel, Detektorzeilen und Detektorwinkeln verwendet wird;
Anwenden (58) einer spiralenförmigen Gewichtungsfunktion auf die konjugierten Abtastwerte zur Verringerung von Bildartefakten; und
Rekonstruieren (60) zumindest eines Bilds des Objekts.
28. Abbildungssystem (10) nach Anspruch 27, wobei das
Abbildungssystem ferner dazu konfiguriert ist, die lineare
Rekonstruktionsebene (70) gemäß der Beziehung:
β = βk-αγ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile K die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β = βk-αγ
zu bestimmen, wobei β die Rekonstruktionsebene ist, βk einen Projektionswinkel darstellt, bei dem eine Detektorzeile K die Rekonstruktionsebene kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie (72) sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
29. Abbildungssystem (10) nach Anspruch 27, wobei das
Abbildungssystem ferner dazu konfiguriert ist, die
konjugierte Abtastlinie (72) gemäß der Beziehung:
β' = βk±π-(2-α)γ
zu bestimmen, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
β' = βk±π-(2-α)γ
zu bestimmen, wobei β' eine konjugierte Abtastlinie ist, βk den Projektionswinkel darstellt, bei dem die Detektorzeile k die Rekonstruktionsebene (70) kreuzt, wenn γ = 0, α ein Wert ist, der es sicherstellt, daß die Rekonstruktionsebene und die konjugierte Abtastlinie sich nicht kreuzen, und α als α ≧ 1-π/(2Qγm) definiert ist, wobei Q eine Spiralenteilung ist, γm ein maximaler Detektorwinkel ist und γ ein Detektorwinkel ist.
30. Abbildungssystem (10) nach Anspruch 27, wobei die
Gewichtungsfunktion ein spiralenförmiges Gewicht gemäß der
Beziehung:
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
umfaßt, wobei wk(γ,β) das spiralenförmige Gewicht ist; β0,k(γ,β) die POR für eine Detektorzeile k ist; β1,k(γ,β) ein vor β0,k(γ,β) gesammelter Abtastwert ist und einen Wert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile darstellt; β2,k(γ,β) für k = 0, . . . 7 nach β0,k(γ,β) gesammelt wird und ein Abtastwert auf der Detektorzeile k konjugiert zu der POR einer weiteren Zeile ist; ηk ein zeilenabhängiger Skalierungsfaktor ist, der für k = 0 und k = 7 gleich 0,5 ist; die mittlere Ansicht der Rekonstruktion als β = 0 definiert ist und Q eine Spiralenteilung ist. Für andere Werte von k ist ηk gleich 1.
31. Abbildungssystem (10) nach Anspruch 27, wobei die
spiralenförmige Gewichtungsfunktion die Beziehung:
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
θk(γ,β) = 3w2 k(γ,β)-2w3 k(γ,β)
umfaßt, wobei θk(γ,β) die spiralenförmige Gewichtungsfunktion ist und wk(γ,β) das spiralenförmige Gewicht ist.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/810,925 US6452996B1 (en) | 2001-03-16 | 2001-03-16 | Methods and apparatus utilizing generalized helical interpolation algorithm |
Publications (1)
Publication Number | Publication Date |
---|---|
DE10211581A1 true DE10211581A1 (de) | 2002-10-24 |
Family
ID=25205054
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE10211581A Ceased DE10211581A1 (de) | 2001-03-16 | 2002-03-15 | Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen Interpolationsalgorithmus |
Country Status (4)
Country | Link |
---|---|
US (1) | US6452996B1 (de) |
JP (1) | JP4356863B2 (de) |
DE (1) | DE10211581A1 (de) |
IL (1) | IL148502A (de) |
Families Citing this family (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6778630B2 (en) * | 2001-03-23 | 2004-08-17 | Kabushiki Kaisha Toshiba | Method and system for reconstructing computed tomography images using redundant data |
US6529576B2 (en) * | 2001-06-29 | 2003-03-04 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for pitch-adaptive reconstruction in multislice helical computed tomographic systems |
US6873676B2 (en) * | 2003-03-05 | 2005-03-29 | Ge Medical Systems Global Technology Company, Llc | Convolution reconstruction algorithm for multi-slice CT |
US6977984B2 (en) * | 2003-10-07 | 2005-12-20 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for dynamical helical scanned image production |
US7415145B2 (en) * | 2003-12-30 | 2008-08-19 | General Electric Company | Methods and apparatus for artifact reduction |
US7050527B2 (en) * | 2004-02-09 | 2006-05-23 | G.E. Medical Systems Global Technology Company, Llc | Methods and apparatus for artifact reduction in cone beam CT image reconstruction |
DE102004020861B4 (de) * | 2004-04-28 | 2009-10-01 | Siemens Ag | Verfahren zur Rekonstruktion von Projektionsdatensätzen bei Dosis-reduzierter abschnittsweiser spiralförmiger Abtastung in der Computertomographie |
US7173996B2 (en) * | 2004-07-16 | 2007-02-06 | General Electric Company | Methods and apparatus for 3D reconstruction in helical cone beam volumetric CT |
US7424088B2 (en) * | 2004-09-29 | 2008-09-09 | Kabushiki Kaisha Toshiba | Image reconstruction method using Hilbert transform |
WO2006038142A1 (en) * | 2004-10-08 | 2006-04-13 | Koninklijke Philips Electronics N.V. | Image reconstruction with voxel dependent interpolation |
US7359478B2 (en) * | 2004-11-18 | 2008-04-15 | Toshiba Medical Systems Corporation | Method for restoring truncated helical cone-beam computed tomography data |
WO2006056942A1 (en) * | 2004-11-23 | 2006-06-01 | Philips Intellectual Property & Standards Gmbh | Image reconstruction device and method |
JP2007236662A (ja) * | 2006-03-09 | 2007-09-20 | Ge Medical Systems Global Technology Co Llc | X線ct装置およびそのx線ct画像再構成方法、x線ct画像撮影方法。 |
US7729890B2 (en) | 2006-08-22 | 2010-06-01 | Analog Devices, Inc. | Method for determining the change of a signal, and an apparatus including a circuit arranged to implement the method |
CN101620191B (zh) * | 2008-07-04 | 2013-01-02 | Ge医疗系统环球技术有限公司 | 具校正ct平扫图像失真的图像获取方法及其装置 |
US8938104B2 (en) * | 2008-08-29 | 2015-01-20 | Varian Medical Systems International Ag | Systems and methods for adaptive filtering |
SE533704C2 (sv) | 2008-12-05 | 2010-12-07 | Flatfrog Lab Ab | Pekkänslig apparat och förfarande för drivning av densamma |
TW201203052A (en) | 2010-05-03 | 2012-01-16 | Flatfrog Lab Ab | Touch determination by tomographic reconstruction |
CN102456228B (zh) * | 2010-10-29 | 2015-11-25 | Ge医疗系统环球技术有限公司 | 图像重建方法和装置及ct机 |
US10168835B2 (en) | 2012-05-23 | 2019-01-01 | Flatfrog Laboratories Ab | Spatial resolution in touch displays |
WO2014168567A1 (en) | 2013-04-11 | 2014-10-16 | Flatfrog Laboratories Ab | Tomographic processing for touch detection |
US9874978B2 (en) | 2013-07-12 | 2018-01-23 | Flatfrog Laboratories Ab | Partial detect mode |
US10146376B2 (en) | 2014-01-16 | 2018-12-04 | Flatfrog Laboratories Ab | Light coupling in TIR-based optical touch systems |
US10126882B2 (en) | 2014-01-16 | 2018-11-13 | Flatfrog Laboratories Ab | TIR-based optical touch systems of projection-type |
JP6321405B2 (ja) * | 2014-02-27 | 2018-05-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成装置、放射線断層撮影装置及び画像生成方法並びにプログラム |
WO2015199602A1 (en) | 2014-06-27 | 2015-12-30 | Flatfrog Laboratories Ab | Detection of surface contamination |
CN104657962B (zh) * | 2014-12-12 | 2017-08-25 | 西安电子科技大学 | 基于级联线性回归的图像超分辨重建方法 |
CN107209608A (zh) | 2015-01-28 | 2017-09-26 | 平蛙实验室股份公司 | 动态触摸隔离帧 |
US10318074B2 (en) | 2015-01-30 | 2019-06-11 | Flatfrog Laboratories Ab | Touch-sensing OLED display with tilted emitters |
EP3256936A4 (de) | 2015-02-09 | 2018-10-17 | FlatFrog Laboratories AB | Optisches berührungssystem mit mitteln zur projektion und detektion von lichtstrahlen über und in einer durchlässigen platte |
CN107250855A (zh) | 2015-03-02 | 2017-10-13 | 平蛙实验室股份公司 | 用于光耦合的光学部件 |
CN108369470B (zh) | 2015-12-09 | 2022-02-08 | 平蛙实验室股份公司 | 改进的触控笔识别 |
EP3545392A4 (de) | 2016-11-24 | 2020-07-29 | FlatFrog Laboratories AB | Automatische optimierung von berührungssignalen |
KR102495467B1 (ko) | 2016-12-07 | 2023-02-06 | 플라트프로그 라보라토리즈 에이비 | 개선된 터치 장치 |
US10963104B2 (en) | 2017-02-06 | 2021-03-30 | Flatfrog Laboratories Ab | Optical coupling in touch-sensing systems |
WO2018174786A1 (en) | 2017-03-22 | 2018-09-27 | Flatfrog Laboratories | Pen differentiation for touch displays |
EP4036697A1 (de) | 2017-03-28 | 2022-08-03 | FlatFrog Laboratories AB | Optische berührungserfassungsvorrichtung |
CN117311543A (zh) | 2017-09-01 | 2023-12-29 | 平蛙实验室股份公司 | 触摸感测设备 |
US11567610B2 (en) | 2018-03-05 | 2023-01-31 | Flatfrog Laboratories Ab | Detection line broadening |
US11943563B2 (en) | 2019-01-25 | 2024-03-26 | FlatFrog Laboratories, AB | Videoconferencing terminal and method of operating the same |
JP2023512682A (ja) | 2020-02-10 | 2023-03-28 | フラットフロッグ ラボラトリーズ アーベー | 改良型タッチ検知装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5270923A (en) | 1989-11-02 | 1993-12-14 | General Electric Company | Computed tomographic image reconstruction method for helical scanning using interpolation of partial scans for image construction |
US5233518A (en) | 1989-11-13 | 1993-08-03 | General Electric Company | Extrapolative reconstruction method for helical scanning |
US5216601A (en) | 1989-11-22 | 1993-06-01 | General Electric Company | Method for fan beam helical scanning using rebinning |
US5208746A (en) | 1989-11-22 | 1993-05-04 | General Electric Company | Method for helical scanning with a stationary detector using rebinning and splicing to create detector vertex projection sets |
US5546439A (en) | 1995-11-02 | 1996-08-13 | General Electric Company | Systems, methods and apparatus for incrementally reconstructing overlapped images in a CT system implementing a helical scan |
US6381297B1 (en) * | 2000-08-15 | 2002-04-30 | Ge Medical Systems Global Technology Company, Llc | High pitch reconstruction of multislice CT scans |
US6404842B1 (en) * | 2001-10-15 | 2002-06-11 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for image reconstruction in twin helical computed tomographic systems |
-
2001
- 2001-03-16 US US09/810,925 patent/US6452996B1/en not_active Expired - Lifetime
-
2002
- 2002-03-05 IL IL148502A patent/IL148502A/en active IP Right Grant
- 2002-03-15 JP JP2002071540A patent/JP4356863B2/ja not_active Expired - Fee Related
- 2002-03-15 DE DE10211581A patent/DE10211581A1/de not_active Ceased
Also Published As
Publication number | Publication date |
---|---|
IL148502A (en) | 2006-12-31 |
JP4356863B2 (ja) | 2009-11-04 |
IL148502A0 (en) | 2002-09-12 |
US6452996B1 (en) | 2002-09-17 |
JP2002330956A (ja) | 2002-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE10211581A1 (de) | Verfahren und Vorrichtung unter Verwendung eines verallgemeinerten spiralenförmigen Interpolationsalgorithmus | |
DE69830986T2 (de) | Rekonstruktion von on-line-bildern bei spiral-computertomographen | |
DE60026681T2 (de) | Verfahren und Gerät zur Mehrschichttomografie mit partieller Abtastung | |
DE60224770T2 (de) | Verfahren und Vorrichtung zur Rauschverminderung in Computertomographen | |
DE69931750T2 (de) | Verfahren und gerät zur verkalkungsmessung | |
DE60024907T2 (de) | Verfahren und Vorrichtung zur Kompensation von Artefakten mittels veränderlicher Winkelabtastung | |
DE69030555T2 (de) | Rechnergesteuerte Tomographie mit schraubenförmiger Abtastung | |
DE69932691T2 (de) | Bildrekonstruktion | |
DE60030507T2 (de) | Volumetrische Bildrekonstruktion | |
DE69218839T2 (de) | Tomographische Bilderzeugung mittels Schrägstrahlen | |
DE10164287A1 (de) | Computertomographiegewichtung (CT-Gewichtung) zur Bildrekonstruktion mit hoher Qualität | |
DE602004007930T2 (de) | TOMOGRAPHISCHES FOURIER-BILDREKONSTRUKTIONSVERFAHREN FÜR Fächerstrahl-DATEN | |
DE19945636A1 (de) | Verfahren und Vorrichtung zum dreidimensionalen CT-Abbilden eines interessierenden Bereichs in einem Objekt | |
DE102007039573A1 (de) | Verfahren zur analytischen Rekonstruktion für eine Mehrfachquellen-Inversgeometrie-CT | |
DE10164162A1 (de) | Rückprojektionsverfahren und -vorrichtung für Computertomographieabbildungssysteme | |
DE19904369A1 (de) | Wendelgewichtungsalgorithmen zur schnellen Rekonstruktion | |
DE10211578A1 (de) | Verfahren und Vorrichtung zur Verarbeitung von CT-Erkundungsbildern | |
DE19647435A1 (de) | Systeme, Verfahren und Vorrichtungen zur Rekonstruktion von Bildern in eine Wendelabtastung verwendenden Computer-Tomographie-Systemen | |
DE60019537T2 (de) | Halbabtastungsalgorithmus für eine schnelle mehrreihige Detektoranordnung | |
DE19743220A1 (de) | Verfahren und Vorrichtung zur dreidimensionalen Bildrekonstruktion mit maximaler Intensitätsprojektion bei einem Computer-Tomographie-System | |
DE102005005839A1 (de) | Verfahren und Vorrichtung zur Verminderung von Artefakten bei Kegelstrahl-CT-Bildrekonstruktionen | |
DE3689166T2 (de) | Substanzquantifizierung in Tierkörpern. | |
DE102005049628A1 (de) | Verfahren und Vorrichtung zur Wiederherstellung geneigter Konusstrahldaten | |
DE19650529A1 (de) | Verfahren und Vorrichtungen zur Mehrfach-Schnitt-Wendel-Bildrekonstruktion in einem Computer-Tomographie-System | |
DE69932390T2 (de) | Verfahren und Vorrichtung zur Bildrekonstruktion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8110 | Request for examination paragraph 44 | ||
R016 | Response to examination communication | ||
R002 | Refusal decision in examination/registration proceedings | ||
R003 | Refusal decision now final |
Effective date: 20131126 |