DE69921154T2 - Anwendung von streu- und abschwächungs-korrekturen auf emissions-tomographie-bildern unter verwendung von einer atlas-abgeleiteten anatomie - Google Patents

Anwendung von streu- und abschwächungs-korrekturen auf emissions-tomographie-bildern unter verwendung von einer atlas-abgeleiteten anatomie Download PDF

Info

Publication number
DE69921154T2
DE69921154T2 DE69921154T DE69921154T DE69921154T2 DE 69921154 T2 DE69921154 T2 DE 69921154T2 DE 69921154 T DE69921154 T DE 69921154T DE 69921154 T DE69921154 T DE 69921154T DE 69921154 T2 DE69921154 T2 DE 69921154T2
Authority
DE
Germany
Prior art keywords
emission tomography
atlas
interest
images
computer model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
DE69921154T
Other languages
English (en)
Other versions
DE69921154D1 (de
Inventor
Frank S. London Prato
Robert Z. Etobicoke Stodilka
Brad J. London Kemp
Richard L. Arva Nicholson
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Application granted granted Critical
Publication of DE69921154D1 publication Critical patent/DE69921154D1/de
Publication of DE69921154T2 publication Critical patent/DE69921154T2/de
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/1611Applications in the field of nuclear medicine, e.g. in vivo counting using both transmission and emission sources sequentially
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/037Emission tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • G01T1/161Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/164Scintigraphy
    • G01T1/1641Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
    • G01T1/1647Processing of scintigraphic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Luminescent Compositions (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Description

  • Technisches Gebiet
  • Die vorliegende Erfindung bezieht sich auf Emissionstomographie und insbesondere auf ein Verfahren und eine Vorrichtung zum Anwenden einer Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder unter Verwendung einer aus einem Atlas abgeleiteten Anatomie.
  • Stand der Technik
  • Einzelphotonenemissions-Computertomographie (SPECT) und Positronenemissions-Tomographie (PET) sind nuklearmedizinische diagnostische bildgebende Techniken, die verwendet werden, um die dreidimensionale Verteilung eines Radiopharmazeutikums innerhalb des Körpers zu messen. Hirn-SPECT- und -PET-Bildtechniken werden primär verwendet, um den regionalen zerebralen Blutfluss in einem Patienten, dem ein Radiopharmazeutikum injiziert wurde, zu messen, um bei der Evaluierung eines Schlaganfalls und der Diagnose von Demenzien wie etwa Alzheimers Krankheit zu assistieren.
  • Obwohl SPECT und PET nützliche Bildtechniken sind, ist ihre geringe quantitative Genauigkeit ein Hindernis bei der Möglichkeit gewesen, eine verbesserte diagnostische Zuverlässigkeit zu erzielen. Die quantitative Genauigkeit von Gehirn-SPECT- und -PET-Bildgebung ist jedoch durch die Anwendung von Streuungs- und Abschwächungskorrektur bei SPECT- und PET-Gehirnbildern signifikant verbessert worden. Um hinreichend akkurat zu sein, muss die Anwendung von Streuungs- und Abschwächungskorrektur an SPECT- und PET-Gehirnbildern durch die Verteilung der Dichte innerhalb des Kopfes geleitet werden. Unglücklicherweise kann die Verteilung der Dichte innerhalb des Kopfes nicht aus den SPECT- und PET-Gehirnscans erhalten werden und daher sind getrennte Messungen notwendig.
  • Die Transmissionsbildgebung ist verwendet worden, um die Verteilung der Dichte innerhalb des Kopfes zu messen und zu gestatten, dass Zerstreuungs- und Abschwächungskorrektur auf SPECT- und PET-Gehirnbilder angewendet wird. Unglücklicherweise ist die für die Durchführung von Transmissionsmessungen notwendige Ausstattung komplex, unzuverlässig und erfordert ausgedehnte Wartung. Auch vergrößert die Notwendigkeit, transmissionsbildgebende Messungen zusätzlich zu den SPECT- oder PET-Gehirnbildern zu machen, die Zeit, die notwendig ist, um die gesamte bildgebende Prozedur abzuschließen. SPECT- und PET-Bildprozeduren sind selbst langwierig und erfordern, dass ein Patient bewegungslos bleibt, um genaue Gehirnbilder sicherzustellen. Für kranke und ältere Patienten ist dies eine schwierige Aufgabe. Ein Vergrößern der Länge der bildgebenden Prozedur vergrößert die Wahrscheinlichkeit, dass die Patienten nicht bewegungslos verharren. Die Bewegung eines Patienten während der Transmissionsbildgebungsprozedur führt zu ungenauen Messungen der Verteilung der Dichte innerhalb des Kopfs. Dies ergibt natürlich eine ungenaue Führung für die Anwendung der Streuungs- und Abschwächungskorrektur auf SPECT- und PET-Hirnbilder. Dementsprechend sind verbesserte Verfahren, die die diagnostische Zuverlässigkeit von Emissionstomographiebildern steigern, wünschenswert.
  • Es ist daher eine Aufgabe der vorliegenden Erfindung, ein neues Verfahren und eine neue Vorrichtung zum Anwenden von Zerstreuungs- und Abschwächungskorrektur auf Emissionstomographiebilder bereitzustellen.
  • Offenbarung der Erfindung
  • Allgemein gesagt, stellt die vorliegende Erfindung ein Verfahren und eine Vorrichtung zum Anwenden einer Zerstreuungs- und Abschwächungskorrektur auf Emissionstomographiebilder bereit, welche die Dichteverteilung innerhalb eines interessierenden Bereichs eines Patienten unter Beobachtung abschätzt oder "ableitet", unter Verwendung eines dreidimensionalen Computermodells des interessierenden Bereichs. Es ist herausgefunden worden, dass eine von einem Computermodell des interessierenden Bereichs unter Beobachtung geführte Streuungs- und Abschwächungskorrektur Ergebnisse ergibt, die ähnlich denen sind, wenn Transmissionsbilder des tatsächlich interessierenden Bereichs als Führung für das Anwenden der Streuungs- und Abschwächungskorrektur verwendet werden.
  • Gemäß einem Aspekt der vorliegenden Erfindung wird ein Verfahren des Anwendens von Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder eines interessierenden Bereichs eines beobachteten Subjekts bereitgestellt, welches die Schritte umfasst:
    Ausrichten eines dreidimensionalen Computermodells, welches die Dichteverteilung innerhalb des interessierenden Bereichs repräsentiert, zu den Emissionstomographiebildern; und
    Anwenden von Streuungs- und Abschwächungskorrektur auf die Emissionstomographiebilder unter Verwendung des ausgerichteten Computermodells als eine Führung.
  • In einer bevorzugten Ausführungsform liegt das Computermodell in Form eines Zweikomponentenatlasses vor. Während des Ausrichtschrittes wird zuerst eine funktionelle Komponente des Atlasses zu den Emissionstomographiebildern ausgerichtet, um einen Satz von räumlichen Transformationsparametern zu erzeugen. Danach wird eine anatomische Komponente des Atlasses mit den Emissionstomographiebildern unter Verwendung des Satzes von räumlichen Transformationsparametern ausgerichtet.
  • Der Atlas kann aus einer Datenbank von Atlanten ausgewählt werden, basierend auf dem Ausmaß an Deckung mit den Emissionstomographiebildern. Alternativ können mehrere Atlanten kombiniert werden, um einen Ergebnisatlas zu erzeugen, der besser zu den Emissionstomographiebildern ausgerichtet ist.
  • Der Atlas kann aus einer Datenbank von Atlanten ausgewählt werden, basierend auf dem Deckunggrad mit den Emissionstomographiebildern. Alternativ können mehrere Atlanten kombiniert werden, um einen Ergebnisatlas zu ergeben, der besser zur Emissionstomographiebildern ausgerichtet ist.
  • Gemäß einem anderen Aspekt der vorliegenden Erfindung wird ein Emissionstomographiebildgebungsverfahren bereitgestellt, bei dem Emissionstomographiebilder eines interessierenden Bereichs eines Subjektes für eine Analyse aufgenommen werden und hinsichtlich der Streuung und Abschwächung korrigiert werden, wobei das Verfahren weiterhin die Schritte umfasst:
    Verwenden eines dreidimensionalen Computermodells, welches die Dichteverteilung innerhalb des interessierenden Bereichs als eine Führung für die Anwendung von Streuungs- und Abschwächungskorrektur annähert.
  • Gemäß noch einem anderen Aspekt der vorliegenden Erfindung wird ein Emissionstomographiebildverarbeitungssystem bereitgestellt, das umfasst:
    einen Speicher, der Emissionstomographiebilder eines interessierenden Bereichs einer Person speichert;
    wobei der Speicher ebenfalls zumindest ein dreidimensionales Computermodell des interessierenden Bereichs speichert, wobei das Computermodell die Dichteverteilung innerhalb des interessierenden Bereichs repräsentiert; und
    einen Prozessor zum Indeckungbringen des Computermodells mit den Emissionstomographiebildern und zum Anwenden einer Streuungs- und Abschwächungskorrektur auf die Emissionstomographiebilder unter Verwendung des zur Deckung gebrachten Computermodells als eine Führung.
  • Gemäß noch einem anderen Aspekt der vorliegenden Erfindung wird ein Emissionstomographiebildgebungssystem bereitgestellt, das umfasst:
    Mittel zum Aufnehmen von Emissionstomographiebildern eines interessierenden Bereichs einer Person, um ein dreidimensionales Bild des interessierenden Bereiches zu bilden;
    einen Speicher zum Speichern der Emissionstomographiebilder, wobei der Speicher auch zumindest ein dreidimensionales Computermodell des interessierenden Bereichs speichert, wobei das Computermodell die Dichteverteilung innerhalb des interessierenden Bereichs repräsentiert; und
    einen Computerprogrammcode zum Ausrichten eines dreidimensionalen Computermodells, welches die Dichteverteilung innerhalb des interessierenden Bereiches repräsentiert, zu den Emissionstomographiebildern; und
    Computerprogrammcode zum Anwenden einer Streuungs- und Abschwächungskorrektur auf die Emissionstomographiebilder unter Verwendung des ausgerichteten Computermodells als eine Führung.
  • Die vorliegende Erfindung stellt Vorteile bereit, indem sie durch Verwenden eines dreidimensionales Computermodell des interessierenden Bereichs einer Person unter Beobachtung, welches sich seiner Dichte als eine Führung für die Anwendung einer Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder annähert, die Notwendigkeit für Transmissionsbildgebung beseitigt. Daher müssen im Fall von SPECT- und PET-Bilderzeugung die bildgebenden Verfahren nicht verlängert werden. Auch wird, da die Dichteverteilung im interessierenden Bereich unter Beobachtung durch ein dreidimensionales Computermodell angenähert wird, zusätzliche Hardware nicht erforderlich, um die Führung für die Anwendung der Streuungs- und Abschwächungskorrektur zu erzeugen. Dies macht die vorliegende Erfindung signifikant weniger teuer und flexibler als transmissionsbildgebende Systeme. Zusätzlich kann, da ein dreidimensionales Computermodell des interessierenden Bereichs unter Beobachtung als Führung für die Anwendung der Streuungs- und Abschwächungskorrektur verwendet wird, die Streuungs- und Abschwächungskorrektur retrospektiv auf existierende Datenbanken angewendet werden, die signifikante Anzahlen von Emissionstomographiebildern beinhalten, für welche keine transmissionsbildgebenden Messungen erhoben worden sind.
  • Kurze Beschreibung der Zeichnungen
  • Eine Ausführungsform der vorliegenden Erfindung wird nun vollständiger unter Bezugnahme auf die beigefügten Zeichnungen beschrieben werden, in denen:
  • 1a ein zweidimensionales Emissionstomographie-Hirnbild ohne die Anwendung einer Streuungs- und Abschwächungskorrektur zeigt;
  • 1b das zweidimensionale Bild von 1a unter Anwendung der Streuungs- und Abschwächungskorrektur zeigt;
  • 1c ein zweidimensionales Transmissionshirnbild ist;
  • 1d und 1e zweidimensionale Emissionstomographie-Hirnbilder mit einer aus Transmissionsbildern hergeleiteten, auf sie angewendeten zugehörigen Anatomie sind;
  • 2 ein Blockdiagramm ist, welches ein Verfahren zum Anwenden von Streuungs- und Abschwächungskorrektur auf Emissionstomographiebildern gemäß der vorliegenden Erfindung zeigt;
  • 3a ein zweidimensionales Emissionstomographie-Hirnbild eines Kopfphantoms mit Streuungs- und Abschwächungskorrektur zeigt;
  • 3b ein zweidimensionales Transmissionshirnbild ist;
  • 3c und 3d zweidimensionale Emissionstomographie-Hirnbilder mit einer aus Transmissionsbildern hergeleiteten, auf sie angewendeten zugehörigen Anatomie zeigen;
  • 3e und 3f zweidimensionale Emissionstomographie-Hirnbilder mit aus einem Atlas abgeleiteter, darauf angewendeter Anatomie sind;
  • 4 und 5 zweidimensionale Hirnbilder und ein Graph zeigen, die gleichförmige zugehörige Anatomie, gleichförmige abgeleitete Anatomie und Hirnkontur vergleichen;
  • 6 zweidimensionale Hirnbilder mit nicht-gleichförmiger zugehöriger und nicht-gleichförmiger abgeleiteter Anatomie, die darauf angewendet sind, zeigt;
  • 7 und 8 eine quantitative Evaluierung von Kopfphantom-SPECT-Rekonstruktionen zeigen;
  • 9 zweidimensionale SPECT-Hirnbilder zeigt;
  • 10 SPECT-Hirnbilder zeigt, die einer Transmissionskonstruktion und abgeleiteter Anatomie überlagert sind;
  • 11 Profile zeigt, die durch Rekonstruktionen eines SPECT-Hirnbilds aufgenommen sind;
  • 12 ein Graph ist, der einen Vergleich von regionalem Hirnblutfluss aus rekonstruierten SPECT-Bildern zeigt, die durch abgeleitete Anatomie und Transmissionsscans geführt werden; und
  • 13 ein Graph ist, der die Korrelation zwischen den rekonstruierten SPECT-Bildern von 12 zeigt.
  • Bestes Vorgehen zum Ausführen der Erfindung
  • Während der Emissionstomographiebildgebung, wie z.B. SPECT oder PET, wird einem Patienten ein Radiopharmazeutikum injiziert. Zweidimensionale Bilder oder Projektionen eines interessierenden Bereichs (ROI = region of interest) eines Patienten werden längs einer beliebigen Achse aufgenommen. Die Projektionen werden digitalisiert und an einen Computer zur Speicherung in einen dreidimensionalen Raster im Computerspeicher übertragen, um ein dreidimensionales Bild des interessierenden Bereichs zu rekonstruieren. Wenn das dreidimensionale Bild einmal rekonstruiert ist, können zweidimensionale Scheiben des interessierenden Bereiches längs praktisch jeder beliebigen Achse unter Verwendung konventioneller Software angezeigt werden. 1a zeigt ein zweidimensionales Emissionstomographie-Hirnbild.
  • Während des bildgebenden Verfahrens werden die Emissionen des Radiopharmazeutikums durch unterschiedlich dichte Gewebe, Lufthohlräume und/oder Knochen im interessierenden Bereich unter Beobachtung gestreut und/oder abgeschwächt. Wie aus 1a ersichtlich, beeinträchtigt die Streuung und Abschwächung von Radiopharmazeutika die quantitative Bildqualität. Im Ergebnis sind die gewonnenen Emissionstomographiebilder oft unzuverlässig. Das Anwenden einer Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder unter Verwendung von Transmissionsbildern desselben interessierenden Bereiches, die während derselben bildgebenden Prozedur aufgenommen worden sind, ist bekannt, leidet aber an den zuvor diskutierten Nachteilen. 1b zeigt das zweidimensionale Bild von 1a, bei dem Streuungs- und Abschwächungskorrektur unter Verwendung von, von derselben Anatomie aufgenommenen Transmissionsbildern, angewendet wurde. 1c zeigt ein Beispiel eines Transmissionsbildes, und die 1d und 1e zeigen zugehörige Anatomien (sagittal bzw. transversal), die von den auf Emissionstomographie-Hirnbildern angewendeten Transmissionsbildern abgeleitet sind.
  • Bei der vorliegenden Erfindung wird ein dreidimensionales Computermodell oder "Atlas" des interessierenden Bereichs, der eine genaue Dichteverteilung des interessierenden Bereichs bereitstellt, in einer Datenbank im Computerspeicher gespeichert und wird als Führung für die Anwendung der Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder verwendet. Bei der vorliegenden Ausführungsform enthält der Atlas zwei Komponenten, nämlich eine, einen SPECT- oder PET-Scan des interessierenden Bereichs simulierende funktionelle Komponente und eine, einen Transmissionsscan des interessierenden Bereichs simulierende anatomische Komponente. Der Atlas kann aus existierenden Transmissionsbildern oder Röntgen-CT-Scans von ähnlichen interessierenden Bereichen aus anderen Patienten erzeugt und gemittelt werden, um ein geeigneter Computermodell oder einen Atlas zu erzeugen. Mehrere Modelle für jeden interessierenden Bereich und Modelle für unterschiedliche interessierende Bereiche können in der Datenbank gespeichert werden und es kann während der Streuungs- und Abschwächungskorrektur von Emissionstomographiebildern individuell darauf zugegriffen werden. Da ein Computermodell des interessierenden Bereichs verwendet wird, ist keine zusätzliche Hardware und Verfahrenszeit notwendig, um eine Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder anzuwenden. Nunmehr wird ein Beispiel der Anwendung von Streuungs- und Abschwächungskorrektur auf Emissionstomographie-Hirnbilder gemäß der vorliegenden Erfindung beschrieben.
  • Nun unter Bezugnahme auf 2 wird ein Blockdiagramm, welches das hiesige Verfahren zum Anwenden von Streuungs- und Abschwächungskorrektur auf Emissionstomographie-Hirnbilder illustriert, gezeigt. Anfangs werden Emissionstomographie-Hirnbilder eines Patienten (Kasten 50) in bekannter Weise erfasst. Die erfassten Hirnbilder liegen in der Form zweidimensionaler Projektionen der radiopharmazeutischen Verteilung im Hirn vor.
  • Anfangs wird eine vorläufige Rekonstruktion der erfassten Hirnbilder durchgeführt, um die zweidimensionalen Projektionen zu digitalisieren und in eine dreidimensionale Anordnung im Computerspeicher zu assemblieren, um ein dreidimensionales Bild des Gehirns zu erzeugen (Kasten 52). Ein Atlas des Kopfs (in diesem Falldes interessierenden Bereichs) wird dann aus der Datenbank in eine Datendatei (Kasten 54) heruntergeladen und die funktionelle Komponente des Kopfatlasses wird identifiziert. Wenn sie einmal identifiziert ist, wird die funktionelle Komponente des Kopfatlasses, in dieser Ausführungsform die Hirnkomponente, in ein dreidimensionales Array kopiert. Daraufhin wird eine Ausrichtprozedur durchgeführt, um die Hirnkomponente des Kopfatlasses mit dem dreidimensionalen Hirnbild räumlich zur Deckung zu bringen (Block 56). Bei der vorliegenden Ausführungsform wird ein Simplexalgorithmus verwendet, um die Hirnkomponente des Kopfatlas mit dem Hirnbild räumlich zur Deckung zu bringen, wie in "Numerical Recipes In C" von Press et al., 2. Auflage, New York, New York, Cambridge University Press, 1992 beschrieben. Fachleute werden jedoch erkennen, dass andere Deckungsalgorithmen verwendet werden können, um die Hirnkomponente des Kopfatlasses am dreidimensionalen Hirnbild auszurichten. Während dieser Ausrichtprozedur wird ein Satz von räumlichen 3D-Transformationsparametern, welche die dreidimensionalen Transformationen repräsentieren, einschließlich, aber nicht beschränkt auf Rotation, Verschiebung und Skalierung, die notwendig sind, um die Hirnkomponente des Kopfatlasses mit dem dreidimensionalen Hirnbild zur Deckung zu bringen, berechnet und im Computerspeicher als eine Matrix gespeichert.
  • Wenn einmal der Satz von 3D-Transformationsparametern berechnet und gespeichert ist, werden die räumlichen 3D-Transformationsparameter im Satz auf die anatomische Komponente des Kopfatlasses angewendet, um ihn mit dem dreidimensionalen Hirnbild zur Deckung zu bringen (Kasten 58). Wenn die anatomische Komponente des Kopfatlasses zum dreidimensionalen Hirnbild ausgerichtet ist, kann der Atlas als eine Dichteverteilungsführung für die Anwendung der Streuungs- und Abschwächungskorrektur an dem dreidimensionalen Hirnbild verwendet werden.
  • Mit einer hergestellten genauen Dichteverteilungsführung wird die Streuungskorrektur auf das dreidimensionale Hirnbild angewendet, gefolgt von einer Abschwächungskorrektur in bekannter Weise (Kasten 60). Wenn die Streuungs- und Abschwächungskorrektur auf das dreidimensionale Hirnbild angewendet worden ist, wird das Hirnbild in ein dreidimensionales Array in einem Computerspeicher rekonstruiert, um den Bildkorrekturprozess abzuschließen (Kasten 62). Das korrigierte dreidimensionale Hirnbild kann dann analysiert werden (Kasten 64). Anhang A enthält Pseudocode, der den oben angegebenen Prozess repräsentiert.
  • 3a bis 3f zeigen einen Vergleich von zweidimensionalen Hirnbildern eines Kopfphantoms mit angewandter, aus Transmissionsbildern (zugehöriger Anatomie) hergeleiteter Anatomie und aus einem Kopfatlas (abgeleiteter Anatomie) hergeleiteter Anatomie.
  • Das vorliegende Verfahren und die vorliegende Vorrichtung wurden unter Verwendung eines anthropomorphen Kopfphantom getestet, das Weichgewebe, harte Gewebe und Lufthohlräume innerhalb eines Schädels modelliert und ein Zweikompartimentenhirnreservoir enthält. Die zwei Kompartimente des Reservoirs wurden getrennt mit zwei Wasserlösungen von Tc-99m mit einem spezifischen Aktivitätsverhältnis von 4 : 1 gefüllt. Ein Fächerstrahl-SPECT wurde erfasst, gefolgt von einem Tc-99m-Transmissionsscan 60 Stunden später. Das rekonstruierte Transmissionsbild wird als zugehörige Anatomie bezeichnet.
  • Fünf Streuungs- und Abschwächungskorrekturschemata wurden basierend auf nicht-gleichförmiger zugehöriger Anatomie, gleichförmiger zugehöriger Anatomie, nicht-gleichförmiger abgeleiteter Anatomie, gleichförmiger abgeleiteter Anatomie und gleichförmiger Hirnkontur evaluiert. Für eine gleichförmige Streuungs- und Abschwächungskorrektur wurden die abgeleiteten und zugehörigen Anatomien segmentiert und ihnen gleichförmigen Abschwächungskoeffizienten von Weichgewebe (0,15 cm–1 für Tc-99m) zugewiesen. Eine SPECT-Rekonstruktion des Kopfphantoms wurde in ähnlicher Weise bearbeitet, um eine durch die Hirnkontur geführte Streuungs- und Abschwächungskorrektur zu erleichtern. Die Streuungskorrektur basierte auf einer nicht-stationären Dekonvolutstreuungssubtraktion, wie in dem von Stodilka et al. verfassten Aufsatz mit dem Titel "The Relative Contributions of Scatter and Attenuation Correction Toward Brain SPECT Quantification", Phys. Med. Biol., 1998. Die Abschwächungskorrektur/Rekonstruktion wurde nachfolgend durch Erwartungsmaximierung geordneter Untermengen durchgeführt, wie in dem von HM Hudson et al. verfassten Aufsatz mit dem Titel "Accelerated Image Reconstruction using Ordered Subsets of Projection Data", IEEE Trans Med Imaging 13 601–609, 1994, dargelegt. Obwohl besondere Beispiele von Streuungskorrektur, Abschwächungskorrektur und Rekonstruktionsalgorithmen beschrieben werden, ist es für Fachleute ersichtlich, dass andere Streuungskorrektur-, Abschwächungskorrektur- und Rekonstruktionsverfahren verwendet werden können.
  • Eine gleichförmige Streuungs- und Abschwächungskorrektur sollte durch die Kontur des abschwächenden Mediums geführt werden, falls ein akzeptables Genauigkeitsniveau für eine objektive Diagnose erwartet wird. Die 4 und 5 zeigen, dass eine gleichförmige zugeordnete Anatomie besser durch gleichförmige abgeleitete Anatomie angenähert wird als durch Hirnkontur. Wie in 6 ersichtlich, sind nicht gleichförmige zugeordnete Anatomie und nicht-gleichförmige abgeleitete Anatomie ähnlich.
  • Kopfphantoms-SPECT-Rekonstruktionen wurden quantitativ durch vier Metriken evaluiert, nämlich Beeinflussung (Bias), Gleichförmigkeit, Kontrastwiedergewinnung und relative Quantifizierung (siehe 7 und 8). Wie ersichtlich, stellt eine von abgeleiteter Anatomie geführte Streuungs- und Abschwächungskorrektur eine quantitative Genauigkeit bereit, die ganz klar der von der Hirnkontur geführten Streuungs- und Abschwächungskorrektur überlegen ist.
  • Die Anwendung einer Streuungs- und Abschwächungskorrektur auf Emissionstomographiebilder einer anderen Anatomie als dem Kopf kann ebenfalls durchgeführt werden. Beispielsweise kann eine abgeleitete Anatomie ebenfalls verwendet werden, um eine Streuungs- und Abschwächungskorrektur auf Herzbilder anzuwenden.
  • Während der Konstruktion des kardialen Atlasses wird eine anatomische Komponente des kardialen Atlasses, welche die anatomischen Merkmale des Thorax repräsentiert, erzeugt und beinhaltet:
    Weichgewebe, wie etwa Herz, Leber, Muskel und Fett;
    Weichgewebe sehr niedriger Dichte, wie etwa die Lungen; und
    hochdichte Gewebe, wie etwa Knochen und Knorpel in den Rippen und der Wirbelsäule.
  • Es wird ebenfalls eine funktionelle Komponente des kardialen Atlasses erzeugt und sie simuliert ein SPECT- oder PET-Herzbild. Geeignete Daten zum Konstruieren des Herz-Atlasses können auf einer Reihe von Wegen erhalten werden, einschließlich, aber nicht beschränkt auf Abbildung eines Phantoms oder einer menschlichen Person durch Röntgen-CT, MRI oder gamma-Kamera-Transmissionscomputertomographie; oder Computersimulation. Der kardiale Atlas kann auch durch Verschmelzen einer Vielzahl von Patientenscans konstruiert werden.
  • Das Verfahren des Verwendens abgeleiteter kardialer Anatomie zum Anwenden von Streuungs- und Abschwächungskorrektur auf Herzbilder bleibt das gleiche wie bei der oben beschriebenen Hirnbildgebung. Zuerst wird eine vorläufige Rekonstruktion der kardialen SPECT-Daten des Patienten durchgeführt. Die funktionelle Komponente des kardialen Atlasses wird dann mit der vorläufigen Rekonstruktion zur Deckung gebracht. Die Deckung beinhaltet eine räumliche Transformation, die Verschieben, Rotieren, Skalieren und/oder nicht-lineare Operationen, wie Verbiegen, enthalten kann. Die Deckungs-Prozedur berechnet eine Matrix, welche die räumliche Transformation repräsentiert, die den Atlasraum in den patientenspezifischen Raum abbildet. Wenn die Matrix einmal berechnet worden ist, wird die Matrix auf die anatomische Komponente des kardialen Atlas angewendet, womit die Anatomie in der Brust abgeleitet wird.
  • Obwohl die obigen Verfahren die Verwendung eines einzelnen generischen Atlasses der Anatomie unter Beobachtung beschreiben, erkennen Fachleute, dass maßgeschneiderte Atlanten entwickelt und in der Datenbank gespeichert werden können. Beispielsweise können krankheitsspezifische Atlanten, wie etwa einen Atlas zur Alzheimer-Krankheit oder ein Schlaganfall-Atlas entwickelt werden und diese können verwendet werden, wenn Emissionstomographiebilder von an diesen Krankheiten leidenden Patienten korrigiert werden. Die krankheitsspezifischen Atlanten können tracer- oder läsionsspezifisch sein, um lokale Konkavitäten bei der Radiopharmazeutikumaufnahme zu gestatten. Patienten mit schwerer Alzheimer- oder Pick-Krankheit haben keinen normalen zerebralen Blutfluss. Bereiche mit Flussdefiziten können die Genauigkeit des Zurdeckungbringen der Emissionstomographiebilder mit der funktionellen Komponente des Kopfatlasses begrenzen, falls der Kopfatlas einen normalen Blutfluss und damit eine gleichförmige Radiopharmazeutikumaufnahme annimmt. Atlanten können auch entwickelt werden, um physische Merkmale, wie etwa beispielsweise einen außergewöhnlich großen nasalen Sinus zu berücksichtigen. Während der Streuungs- und Abschwächungskorrektur kann der Bediener den geeigneten, zu verwendenden Atlas auswählen. Alternativ kann der Atlas automatisch durch Computersoftware ausgewählt werden. In diesem Fall führt die Software eine vorläufige Rekonstruktion unter Verwendung jedes maßgeschneiderten Atlasses aus und bringt den Atlas der vorläufigen Rekonstruktion zur Deckung, um die Genauigkeit der Deckung zu messen. Der Atlas mit der höchsten Deckungsgenauigkeit wird dann ausgewählt. Falls gewünscht, kann die Software fuzzy-logic, theoretische Kalkulationen oder andere Kriterien verwenden, um zwei oder mehr Atlanten zu kombinieren, um einen Einzelergebnisatlas zu ergeben, der einen besseren Grad der Deckung bereitstellt.
  • Beispiel
  • Beispiel 1: Abgeleitete Anatomie und Hirnabbildung
  • Das folgende Beispiel wird für illustrative Zwecke beschrieben und soll nicht den Schutzumfang der vorliegenden Erfindung beschränken.
  • Abgeleitete Anatomie aus einem Kopfatlas
  • Ein Kopfatlas wurde wie folgt vorbereitet. Ein dreidimensionales digitalisiertes MRI-Kopfphantom nach Zubal [Zubal et al., 1994] wurde segmentiert, um zwei Datensätze herzustellen, nämlich einen SPECT-Atlas, der einen SPECT-Scan des Phantoms simuliert, und einen anatomischen Atlas, der einen Transmissionsscan des Phantoms simuliert. Der SPECT-Atlas bestand aus Voxeln, die graue Materie und weiße Materie enthielten, denen eine relative Aufnahme von 99mTc-HMPAO von 4 bzw. 1 zugewiesen wurden. Der anatomische Atlas bestand aus Voxeln, die Hartgewebe, Weichgewebe und nasalen Sinus enthielten, denen die entsprechenden 140 keV Schmalstrahl-Abschwächungskoeffizienten von 0,25, 0,15 und 0,075 m–1 zugewiesen wurden.
  • Ein Patienten-SPECT-Scan wurde dann ohne Streuungs- und Abschwächungskorrektur rekonstruiert. Die faciale Aktivität wurde entfernt, um einen Datensatz zu erzeugen, der als die vorläufige Patientenrekonstruktion bezeichnet wurde. Der SPECT-Atlas wurde dann mit der vorläufigen Patientenrekonstruktion zur Deckung gebracht, und die räumliche Transformation wurde aufgezeichnet. Diese Transformation wurde dann auf den anatomischen Atlas angewendet, womit die Lokation von Weichgewebe, Hartgewebe und Lufthohlräumen des Patienten abgeleitet wurde (siehe 9).
  • Ein radiologisches Mehrzweckanalyseprogramm (Hermes, Nuclear Diagnostics, Stockholm, Schweden) wurde verwendet, um die Unimodalitätsdeckung durchzuführen. Die große Variation bei der Kopforientierung erforderte, dass zuerst eine manuelle Zurdeckungbringung durchgeführt wurde. Dieser folgte eine automatische Verfeinerung. Die Kostenfunktion der automatischen Zurdeckungbringung wurde als die Summe der absoluten Zähldifferenzen definiert [Hoh et al. 1993]. Ein globales Minimum wurde durch Simplexsuche [Neider und Mead 1965] innerhalb eines Parameterraums gesucht, der aus Rotation, Verschieben und linearer Skalierung in x-, y- und z-Richtungen bestand [Holman et al. 1991, Slomka et al. 1995].
  • Sequenzielle Transmissions- und Emissionsabbildung
  • Es wurden 10 Demenzpatienten (5 weiblich und 5 männlich, mit einem mittleren Alter von 64,3 Jahre) analysiert. Für jeden Patienten wurde zuerst ein Transmissionsscan erfasst. Die Patienten wurden während des ruhigen Transmissionsscan sehr entspannt und ihnen wurden dann 740 MBq 99mTc-HMPAO injiziert. Die SPECT-Prozedur wurde ungefähr 5 Minuten post injektionem gestartet. Das SPECT-System, das Transmissionsfähigkeiten hat, ist im Detail beschrieben [Kemp et al. 1995]. Es besteht aus einer General Electric Gamm-Kamera 400AC/T (General Electric, Milwaukee, WI) mit einem kreisförmigen Sichtfeld mit 409,6 mm Durchmesser. Projektionen wurden durch einen Fächerstrahl-Kollimator (fan-beam collimator) (Nuclear Fields, Des Planes, IL) mit einer Brennweite von 600 mm und einer Ebenen-zu-Ebenen-hexagonalen Lochbreite von 1,5 mm erfasst. Die Transmissionskomponente beinhaltete einen auf dem Kamera-Kollimator montierten Rahmen, der eine Tantal-kollimierte 99mTc-Linienquelle längs der Fokuslinie des Fächerstrahl-Kollimators hält. Die Kollimation sowohl der Linienquelle als auch der Kamera minimiert die Streuung. Als Ergebnis misst das Transmissionssystem effektiv Schmaltrahlabschwächungskoeffizienten [Tsui et al. 1989, Kemp et al 1995]. SPECT- und Transmissionsscans wurden mit einem 20% Energiefenster erfasst, das auf den 99mTc-Photopeak von 140 keV zentriert war. Die Scans bestanden aus 128 Projektionen, die gleichmäßig über 360° beabstandet waren. Jede kreisförmige Projektion wurde in einer 128 × 128 Pixel-Quadratmatrix (1 Pixel = 3,2 mm) erfasst. Sowohl Transmissions- als auch SPECT-Scans betrugen 10 Sekunden pro Projektion und die Zählraten waren ungefähr 70 bzw. 1,5 Kilozählungen pro Sekunde. Alle Scans wurden bezüglich Gleichförmigkeit unter Verwendung von 100 Millionen Zählflutbildern korrigiert, und die Transmissionsscans wurden auf 50 Millionen Zählleerbilder normalisiert. Rotationsradii variierten unter den Patienten, wobei der kleinste 170 mm und der größte 205 mm war. Vor der Rekonstruktion wurden alle Scans auf Objektebenen-Parallellochgeometrie über eine zweidimensionale kubische Interpolation "rebinned".
  • Streuungs- und Abschwächungskorrektur und Rekonstruktion
  • Die SPECT-Daten wurden unter Verwendung eines Maximalwahrscheinlichkeits-Abschätzers mit einer ungeregelter 32-Ebenen-geordneten Untermengen- [Hudson and Larkin 1994] Implementierung des Erwartungsmaximierungsalgorithmus [Shepp und Vardi 1982, Lange und Carson 1984] rekonstruiert (OSEM). Die vier Projektionen, die pro Sub-Iteration verwendet wurden, wurden über 360° gleich beabstandet. Die Abschwächung wurde im passenden Projektor/Rückprojektorpaar moduliert und eine Streuungsschätzung [Stodilka et al. 1998b] wurde als ein a priori Hintergrund nach der Vorwärtsprojektion [Lange und Carson 1984, Bowsher et al. 1996, Kadrmas et al. 1998] hinzugefügt. Sowohl die Streuungs- als auch die Abschwächungsmodulierung inkorporierten die Schmalstrahl-Abschwächungskoeffizienten aus der Transmissionsbildgebung der abgeleiteten Anatomie. Die Detektorantwort wurde nicht eingeschlossen. Es wurden vier Iterationen von OSEM nach der Initialisierung mit einem gleichmäßigen [Nunez und Llacer 1990]-Träger, der zuvor aus der Transmissionsrekonstruktion oder der abgeleiteten Anatomie abgeleitet wurde, verwendet. Die Rekonstruktionen wurden dann unter Verwendung eines dreidimensionalen Butterworth-Filters mit einer Ordnung von 8 und einem Cutoff bei 0,42 cm–1 nachgefiltert [Nunez und Llazer 1990].
  • Die Transmissionsdaten wurden auch durch den Emissions-OSEM-Algorithmus nach einer Leer-Scan-Normalisierung und log-Transformation rekonstruiert.
  • Die Linienquellenkollimation, gekoppelt mit Beschränkungen der Detektorzählratenfähigkeit und der Patientenverträglichkeit führten zu weniger als idealen Transmissionsstatistiken. Um die Effekte der Transmissionsbildgebungs-Rauschfortpflanzung in die SPECT-Rekonstruktion [Xu et al 1991, Tung und Gullberg 1994] zu vermindern, wurden die Transmissionskonstruktionen wie folgt segmentiert. Das Weichgewebe in den rekonstruktuierten Transmissionsvolumina wurde auf gleichförmige Dichte gezwungen. Zuerst wurde ein großer interessierender Bereich (ROI) um Weichgeweberegionen gezeichnet, und es wurden Mittelwert und Varianzabschätzungen berechnet. Dann wurde allen Voxeln mit Zelldichten innerhalb ±2 Standardabweichungen dieses Mittelwerts 0,15 cm–1 zugewiesen. Somit waren die Transmissionsrekonstruktionen durch rauschfreies Weichgewebe charakterisiert, beinhalteten jedoch Hartgewebe und Lufthohlräume.
  • Vorlagen-basierte Quantifizierung
  • Frühe Arbeiten haben gezeigt, dass eine Hauptkonfusion für eine reproduzierbare Quantifizierung aus manueller und schwellenwertabhängiger Platzierung von anatomischen ROIs auf SPECT-Scans herstammt [Msaki et al. 1998]. Um diese subjektive Fehlerquelle zu vermindern, wurde eine normale Vorlage [Msaki et al. 1998], auf der zwölf bilateralen volumetrischen ROIs abgegrenzt [Karbe et al. 1994] sind, verwendet. Die Fähigkeit, die Vorlage und ihre ROIs zu speichern, sichert die Reproduzierbarkeit der Analyse. Diese Quantifizierungsprozedur führt auch die Standardisierung der Analyse ein, was den Austausch von Daten zwischen verschiedenen Institutionen erleichtert [Evans et al. 1988]. Alle rekonstruierten Scans wurden mit der normalen Vorlage zur Deckung gebracht, was hier als "räumliche Normalisierung" bezeichnet wird. Vor der räumlichen Normalisierung wurden zuvor als Facialaktivität identifizierte Voxel auf Null gesetzt. Nach der Überlagerung der Vorlagen ROIs auf jeden Scan, wurde der kortikale rCBF für jeden ROI auf den zerebellaren rCBF normalisiert [Karbe et al. 1994] und bezüglich eines Blutflussabhängigen Tracer-Reflux korrigiert [Lassen et al. 1988]. Die Analyse der absoluten Konzentration von Pharmazeutikum wurde nicht durchgeführt, da derzeit die absolute rCBF-Quantifizierung bei SPECT selten verwendet wird [Bakker und Pausels 1997].
  • Quantitative Fehleranalyse
  • Eine Probe einer abgeleiteten Anatomie und Transmissionsrekonstruktion ist in 10 illustriert. In 11 werden Profile durch die von Transmissionsscans und abgeleiteter Anatomie geführten SPECT-Rekonstruktionen präsentiert. Die Profile sind nicht auf die cerebellare Zähldichte normalisiert worden.
  • Die Mittelwerte und Standardfehler des regionalen zerebralen Blutflusses aus von abgeleiteter Anatomie und Transmissionsscans geführten SPECT-Scans wurden ROI-für-ROI vergleichen (siehe 12). Die statische Analyse wurde auch ROI-für-ROI unter Verwendung wiederholter Varianzanalyse durchgeführt, um zu bestimmen, wo es signifikante Differenzen zwischen der abgeleiteten Anatomie und den Transmissions-geführten Rekonstruktionsverfahren gab. Gepoolte ROI und ROI-abhängige Korrelationskoeffizienten wurden zwischen abgeleiteter Anatomierekonstruktion und Transmissions-geführten SPECT-Rekonstruktionen berechnet. 13 zeigt die Korrelation zwischen den zwei Rekonstruktions- und Quantifizierungsverfahren durch Poolen aller ROIs und Patienten. Der p-Wert für die Signifikanz wurde für alle Tests auf 0,05 eingestellt.
  • Abgeleitete Anatomiefehlerfortpflanzungssanalyse
  • Vier Fehlerquellen wurden identifiziert, die zu Diskrepanzen bei der ROI-Quantifizierung beitragen:
    • (1) Durch abgeleitete Anatomie geführte Streuungskorrektur;
    • (2) durch abgeleitete Anatomie geführte Abschwächungskorrektur;
    • (3) durch abgeleitete Anatomie geführte räumliche Normalisierung, die mit einer ROI-Fehlplatzierung äquivalent ist; und
    • (4) Patientenbewegung zwischen Transmissions- und SPECT-Scans.
  • Die ersten beiden sind inhärente Beschränkungen des Prinzips abgeleiteter Anatomie, während die beiden letzten Artefaktübertreibungen von Fehlern im Kontext dieses Beispiels repräsentieren. Die ersten drei Fehlerquellen, nämlich Streuung, Abschwächung und ROI-Fehlplatzierung wurden gemessen.
  • Die Fortpflanzung der ROI-Quantifizierung auf Unterschiede zwischen abgeleiteter Anatomie-Streuungskorrektur und Transmissions-geführter Streuungskorrektur wurden analysiert, indem durchgeführt wurde:
    • (1) durch abgeleitete Anatomie hergeleitete Streuungskorrektur;
    • (2) Transmissions-hergeleitete Abschwächungskorrektur; und
    • (3) Anwenden der räumlichen Normalisierung, die als Zurdeckungbringen der Transmissions-hergeleiteten SPECT-Rekonstruktion als optimal berechnet war, auf die Quantifizierungvorlage.
  • Eine ähnliche Analyse zum Evaluieren der Effekte von abgeleiteter Anatomieabgeleiteter Abschwächungskorrektur wurde ausgeführt, indem durchgeführt wurde:
    • (1) Transmissions-hergeleitete Streuungskorrektur;
    • (2) von abgeleiteter Anatomie hergeleitete Abschwächungskorrektur; und
    • (3) Anwenden der räumlichen Normalisierung, die als Zurdeckungbringen der Transmissions-hergeleiteten SPECT-Rekonstruktion als optimal berechnet war, auf die Quantifizierungsvorlage.
  • Die Effekte der ROI-Missplatzierung wurden durch Durchführen von einer Transmissions-hergeleiteten Streuungs- und Abschwächungskorrektur und dann Anwenden der räumlichen Transformation, die als zum Zurdeckungbringen der abgeleiteten Anatomie hergeleiteten SPECT-Rekonstruktion mit der Quantifizierungsvorlage als optimal berechnet waren, quantifiziert. Somit führte die volle Fortpflanzungsanalyse zu drei Rekonstruktionen, von denen Jede quantitativ mit dem "Gold-Standard" der Transmissions-hergeleiteten SPECT-Rekonstruktion und der räumlichen Normalisierung verglichen wurden (über die oben beschriebene Vorlagen-basierte Quantifizierungsprozedur). Diese Prozedur wurde für die 10 Patienten durchgeführt und die Ergebnisse wurden gemittelt. Wenn jedoch einmal die Fehler aufgrund von Streuung, Abschwächung und ROI-Fehlplatzierung getrennt quantifiziert wurden, wurde gefunden, dass ihre Gesamtheit weniger war als die Fehler, die vom vollständig abgeleiteten Anatomieprotokoll verursacht wurden. Diese zusätzliche Fehlerquelle wird in Tabelle 2 als "nicht verbuchbar" bezeichnet, und es wird angenommen, dass sie durch eine Patientenbewegung verursacht wurde.
  • Qualitative Analyse
  • Ein Probenvergleich von abgeleiteter Anatomie und Transmissionsrekonstruktion wird in 10 für mittel-sagittale und kortikal-axiale Schnitte gezeigt. Eine gute Reproduktion von Weichgewebe und Hartgewebe auf dem kortikalem Nideau wird angemerkt. In der Nähe des Vertex wird eine gewisse Diskrepanz sichtbar, jedoch ist dieser Bereich selten in einer quantitativen Analyse eingeschlossen. Diskrepanzen nahe am Nasensinus sind die markantesten. Glücklicherweise involvieren diese Strukturen meist niederdichte Bereiche wie etwa Luft und in einem viel kleineren Ausmaß Weichgewebe und Knorpel, die Photonen nicht so stark streuen und abschwächen wie hochdichte Strukturen.
  • Mittel-Hirnprofile, in 11 gezeigt, vergleichen eine SPECT-Rekonstruktion, die durch Transmissionsbildgebung geführt wird, mit derselben SPECT-Rekonstruktion, die durch abgeleitete Anatomie geführt wird. Das Profil wurde längs der längsten Achse des Kopfes genommen, die nach Streuungs- und Abschwächungskorrektur für falsch zur Deckung gebrachte Transmissionskarten am sensitivsten ist [Huang et al. 1979].
  • Quantitative Fehleranalyse
  • Tabelle 1 unten zeigt die Ergebnisse der wiederholten Analyse von Varianz- und Korrelationsanalyse, die Transmissionsrekonstruktion und quantitative und abgeleitete Anatomie-geführte Rekonstruktion und Quantifikation vergleicht. Die linken und zentralen Sulcus ROIs zeigten die höchste Wahrscheinlichkeit eines wirklichen Unterschieds und erreichten eine statistische Signifikanz mit p = 0,001 bzw. 0,002. Diese ROIs hatten auch die höchsten Korrelationskoeffizienten, die Transmissionssrekonstruktion und Quantifizierung und abgeleitete Anatomiegeführte SPECT-Rekonstruktion und Quantifizierung korrelierten. Diese gesteigerte Korrelation kann ein Artefakt der größeren Differenzen in den ROI-Mittelwerten sein. Interessanterweise wurde gefunden, dass die linken frontalen und zentralen Sulcus ROIs auch markante rCBF-Defizite aufwiesen, was suggeriert, dass eine abgeleitete Anatomie Schwierigkeiten in der Nähe von Bereichen mit maßgeblich verminderter Radiopharmazeutikaufnahme haben kann.
  • 12 zeigt die Mittelwerte und Standardfehler für Transmissionsrekonstruktion und abgeleitete Anatomie bei jedem der 12 ROIs, gemittelt über die gesamte Population. Die mittlere absolute Differenz für alle ROIs über die gesamte Population betrug 7,5 %. Die Korrelation für alle ROIs und alle Patienten wurde als hoch gefunden: r2 = 0,92, wie in 13 illustriert. Tabelle 1
    Figure 00210001
  • Fehlerfortpflanzungsanalyse
  • Die Ergebnisse der Fortpflanzungsanalyse werden in Tabelle 2 unten präsentiert. Diese Ergebnisse werden als Fehler relativ zum "Gold-Standard" der Transmissionsrekonstruktions-geführten Rekonstruktion und räumlichen Normalisierung präsentiert. Drei Fehlerquellen aufgrund der abgeleiteten Anatomie wurden analysiert, nämlich Fehler bei den Streuungsverteilungsschätzungen, Fehler aufgrund missgeleiteter Abschwächungskompensation und Fehler aufgrund von Fehlplatzierung interessierender Bereiche. Die Fehlerkomponenten, gemittelt über 10 Patienten, werden für jeden der 12 bilateralen interessierenden Bereiche gezeigt. Die Tabelle fasst die Prozentzahlen zusammen, die jede dieser Fehlerquellen zu den gesamten quantitativen Differenzen zwischen SPECT-Rekonstruktionen beitrugen, die durch Transmissionsscans oder abgeleiteter Anatomie geführt wurden. Die vierte nummerische Säule zeigt die Prozentzahlen der Gesamtdifferenzen an, die nicht bei abgeleiteter Anatomie verbucht werden konnten.
  • Im Durchschnitt wurde gefunden, dass Fehler bei Streuungsverteilungsabschätzungen zu ungefähr 10,0 % des quantitativen Gesamtfehlers führten; Abschwächungskorrektur: 36,6 % und ROI-Fehlplatzierung: 27,0 %. Die relativen Beiträge der abgeleiteten Anatomie-hergeleiteten Streuungs- und Abschwächungskorrektur für den Gesamtfehler erscheinen glaubwürdig. Eine Kompensation der Abschwächung hat beachtlich größere Konsequenzen als das Entfernen gestreuter Photonen [Rosenthal et al. 1995]. Ungefähr 26,5 % der Gesamtdiskrepanz zwischen abgeleiteten Anatomien und Transmissionsbildgebung konnte nicht bei der Fehlerfortpflanzungsanalyse von Streuung, Abschwächung und ROI-Fehlplatzierung verbucht werden, sondern es wird angenommen, dass sie ein Ergebnis der Patientenbewegung während der Datenerfassung ist. Tabelle 2
    Figure 00220001
    Figure 00230001
  • Qualitative und quantitative Vergleiche
  • Ein Vergleich der abgeleiteten Anatomie mit Transmissionsrekonstruktionen zeigte eine gute Reproduktion von Weichgewebe und Hartgewebe in kortikalen Bereichen für alle 10 Patienten an. Jedoch wurde bei vielen Scans eine Diskrepanz nahe der Sinushöhle angezeigt, wie in 10 gezeigt. Trotzdem ist die abgeleitete Anatomie robuster und akkurater beim Bereitstellen von Abschätzungen unterliegender Gewebeverteilung als das Passendmachen von Ellipsen an Photopeak-Emissionsdaten. Die Technik des Passendmachens von Ellipsen hängt von einer adäquaten Facialaktivität ab, da die Konturen des Gehirns und des Kopfes sich auf dem Niveau des Cerebellums so beachtlich unterscheiden. Die Facialaktivität sollte nicht die Basis zum Abschätzen unterliegender Gewebeverteilungen sein, da die Aufnahme mit dem Radiopharmazeutikum und der Zeit zwischen Injektion und Scannen variiert [Leveille et al. 1992, Van Dyck et al. 1996], was sie zu einer unzuverlässigen Abhängigkeit macht. Übrigens zeigten die qualitativen Ähnlichkeiten, die zwischen abgeleiteter Anatomie und Transmissionsrekonstruktion demonstriert wurden, das Zutrauen in das genaue Führen der Streuungs- und Abschwächungskorrektur. Jedoch ist es wichtig anzumerken, dass eine ähnliche Form eine hinreichende, aber nicht notwendige Voraussetzung für eine genaue Streuungs- und Abschwächungskorrektur ist [Welch et al. 1997, Natterer 1993].
  • Es wurden auf den Transmissionsbildern bei 3 der 10 Scans leichte Abflachungseffekte bemerkt. Eine Abflachung trat bei kyphotischen Patienten mit breiten Schultern oder mit kurzen Hälsen auf. Diese Beschränkung war allgemein auf transaxiale Schnitte unterhalb der Ebene des Cerebellums oder sehr nahe an seiner Basis beschränkt, wo das zirkuläre Sichtfeld der gamma-Kamera sich als unpraktisch erwies. Die Abflachung betraf nur Nasenknorpel, und es wird daher nicht erwartet, dass sie die Ergebnisse signifikant beeinflusst, wie auch durch die auf cerebellarer Ebene gezeigte quantitative Genauigkeit demonstriert (siehe 12).
  • Die abgeleitete Anatomie hatte Schwierigkeiten in Regionen mit markantem rCBF-Defizit, wie etwa den linken frontalen Lappen. Obwohl Frontallappen allgemein eine hohe Variabilität bei der HMPAO-Aufnahme zeigen [Deutsch et al. 1997] hatte nur der linke Frontallappen einen statistisch signifikanten Fehler. Dies suggeriert, dass regionale quantitative Fehler, die durch abgeleitete Anatomie auftraten, mit rCBF-Defiziten assoziiert sind. Jedoch zeigen vorherige Arbeiten, dass sie auch für den räumlichen Bereich sensitiv sein können. Das Erzielen einer guten quantitativen Genauigkeit in ausgedehnten Quellen ist oft schwierig, und dies ist insbesondere bei peripheren ROIs wahr, wo die rekonstruierte Aktivität für eine fehlerhafte Zurdeckungbringung des Abschwächungsmediums am sensitivsten ist [Huang et al. 1979]. Beispielsweise kann im Gehirn eine Fehlpassung zwischen Emissions- und Transmissionsdaten von nicht mehr als 5 mm einen 10%igen Fehler in einem 10 mm-dicken peripheren kortikalen ROI erzeugen [Andersson et al. 1995b, Huang et al. 1979]. Im Allgemeinen sind ausgedehnte Quellen, wie etwa das Gehirn, von einer fehlerhaften Zurdeckungbringung mehr betroffen als kompakte Quellen, wie etwa das Herz [Andersson et al. 1995b, McCord et al. 1992, Ter-Pogossian 1992, Matsunari et al. 1998]. Da das Gehirn elliptisch ist, wird erwartet, dass Bereiche längs der Peripherie der Längsachse des Kopfes gegenüber einer Abschwächungsplanfehldeckung sensitiver sind als jene längs der kurzen Achse.
  • Obwohl eine bevorzugte Ausführungsform der vorliegenden Erfindung beschrieben worden ist, werden Fachleute sehen, dass Variationen und Modifikationen an der vorliegenden Erfindung gemacht werden können, ohne von deren Schutzumfang, wie durch die anhängigen Ansprüche definiert, abzuweichen.
  • Anhang A
    Figure 00250001
  • Anhang A
    Figure 00260001
  • Anhang A
    Figure 00270001
  • Anhang A
    Figure 00280001
  • Referenzen
    • 1. Acton PD and Friston KJ (1998) Statistical parametric mapping in functional neuroimaging: beyond PET and fMRI activation studies Eur J Nucl Med 25 663–667
    • 2. Andersson JLR, Sundin A and Valind S (1995a) A Method for Coregistration of PET and MR Brain Images J Nucl Med 36 1307–1315
    • 3. Andersson JLR, Vagnhammer BE and Schneider H (1995b) Accurate Attenuation Conection Despite Movement during PET Imaging J Nucl Med 36 670–678
    • 4. Andersson JLR (1995e) A Rapid and Accurate Method to Realign PET Scans Utilizing Image Edge Information J Nucl Med 36 657–669
    • 5. Ardekani B, Braun M and Hutton B (1993) Improved quantification with the use of anatomical information in PET image reconstruction. In: Uemura K, Lassen NA, Jones T, Kanno I, eds. Quantification of Brain Function, tracer kinetics and image analysis in brain PET. International congress series no 1030 351–359
    • 6. Ashare AB and Chakraboriy DP (1994) Artificial Neural Networks: Better Than the Real Thing? J Nucl Med 35 2048–2049
    • 7. Bacharach SL and Buvat I (1995) Attenuation correction in cardiac positron emission tomography and single-photon emission computed tomography J Nucl Cardiol 2 246–255
    • 8. Bailey DL (1998) Transmission scanning in emission tomography Eur J Nucl Med 25 774–787
    • 9. Bakker D and Pauwels EKJ (1997) Stroke: the role of functional imaging Eur J Nucl Med 24 2–5
    • 10. Bergstrom M, Litton J, Eriksson L, Bohm C and Blomkvist G (1982) Determination of Object Contour from Projections for Attenuation Conection in Cranial Positron Emission Tomography J Comp Assist Tomo 6 365–372
    • 11. Bergstrom M, Eriksson L, Bohm C, Blomqvist G and Litton J (1983) Conection for Scatter Radiation in a Ring Detector Positron Camera by Integral Transformation of the Projection J Comput Assist Tomo 7 42–50
    • 12. Bailey DL, Hutton BF and Walker PJ (1987) Improved SPECT using simultaneous emission and transmission tomography. J Nucl Med 28 844–851
    • 13. Bowsher JE, Johnson VA, Turkington TG, Jaszczak RJ, Floyd CE and Coleman RE (1996) Bayesian Reconstruction and use of Anatomical A Priori Information for Emission Tomography IEEE Trans Med Imag 15 673–686
    • 14. Budinger TF and Gullberg GT (1977) Transverse section reconstruction of gamma-ray emitting radionuclides in patients. In: Ter-Pogossian MM, Phelps ME, Brownell GL, eds Reconstruction tomography in diagnostic radiology and nuclear medicine Baltimore: University Park Press 315–342
    • 15. Chabriat H, Lavasseur M, Vidailhet M, Loc'h C, Maziere B, Bourguignon MH, Bonnet AM, Zilbovicius M, Raynaud C, Agid Y, Syrota A and Samson Y (1992) In-Vivo SPECT Imaging of D2 Receptor with Iodine-Iodolisuride: Results in Supranuclear Palsy J Nucl Med 33 1481–1485
    • 16. Chan KH, Johnson KA, Becker JA, Satlin A, Mendelson J, Garada B, and Holman BL (1994) A Neural Network Classifier for Cerebral Perfusion Imaging J Nucl Med 35 771–774
    • 17. Claus JJ, van Harskamp F, Breteler MMB, Krenning EP, van der Cammen TJM, Hofman A and Hasan D (1994) Assessment of cerebral perfusion with single-photon emission tomography in normal subjects and in patients with Alzheimer's disease: effects of region of interest selection. Eur J Nucl Med 21 1044–1051
    • 18. DeFigueiredo RJP, Shankle WR, Maccato A, Dick MB, Mundkur P, Mena I and Cotman CW (1995) Neural-network-based classification of cognitively normal, demented, Alzheimer disease and vascular dementia from single photon emission with computed tomography image data from brain Proc Natl Acad Sci USA 92 5530–5534
    • 19. Deutsch G, Mountz JM, Katholi CR, Liu HG and Hanell LE (1997) Regional Stability of Cerebral Blood Flow Measured by Repeated Technetium-99m-HMPAO SPECT: Implications for the Study of State-Dependent Change J Nucl Med 38 6–13
    • 20. Eberl S, Kanno I, Fulton RR, Ryan A, Hutton BF and Fulham MJ (1996) Automated Interstudy Image Registration Technique for SPECT and PET J Nucl Med 37 137–145
    • 21. Evans AC, Beil C, Marrett S, Thompson CJ and Hakim A (1988) Anatomical-Functional Corretation Using an Adjustable MRI-Based Region of interest Atlas with Positron Emission Tomography J Cereb Blood Flow Metab 8 513–530
    • 22. Floyd CE, Jaszczak RJ, Greer KL and Coleman RE (1985) Deconvolution of Compton Scatter in SPECT J Nucl Med 26 406–408
    • 23. Frey EC, Tsui BMW and Perry JR (1992) Simultaneous acquisition of emission and transmission data for improved 71–201 cardiac SPECT imaging using a Tc-99m transmission source. J Nucl Med 33 2238–2245
    • 24. Hoh CK, Dahlbom M, Harris G, Choi Y, Hawkins RA, Phelps ME and Maddahi J (1993) Automated iterative three-dimensional registration of positron emission tomography images J Nucl Med 34 2009–2018
    • 25. Hollinger EF, Loncaric S, Yu D-C, Ali A and Chang W (1998) Using Fast Sequential Asymmetric Fanbeam Transmission CT for Attenuation Correction of Cardiac SPECT Imaging J Nucl Med 39 1335–1344
    • 26. Halman BL, Zimmerman RE, Johnson KA, Carvalho PA, Schwartz RB, Loeffler JS, Alexander E, Pelizzari CA and Chen GTY (1991) Computer-Assisted Superimposition of Magnetic Resonance and High-Resolution Technetium-99m-HMPAO and Thallium-201 SPECT Images of the Brain J Nucl Med 32 1418–1484
    • 27. Hasoba M, Wani H, Toyama H, Murata H and Tanaka E (1986) Automated Body Contour Detection in SPECT: Effects on Quantitative Studies J Nucl Med 27 1184–1191
    • 28. Huang S-C, Hoffman EJ, Phelps ME and Kuhl DE (1979) Quantitation in Positron Emission Computed Tomography: 2. Effects of Inaccurate Attenuation Correction J Comp Assist Tomo 3 804–814
    • 29. Huang S-C, Carson RE, Phelps ME, Hoffman EJ, Schetbert HR, and Kuhl DE (1981) A Boundary Method for Attenuation Correction in Positron Computed Tomography J Nucl Med 22 627–637
    • 30. Hudson HM and Larkin RS (1994) Accelerated Image Reconstruction using Ordered Subsets of Projection Data IEEE Trans Med Imag 13 601–604
    • 31. Iida H, Narita Y, Kado H, Kashikura A, Sugawara S, Shoji Y, Kinoshita T, Ogawa T and Eberl S (1998) Effects of Scatter and Attenuation Correction on Quantitative Assessment of Regional Cerebral Blood Flow with SPECT J Nucl Med 39 181–189
    • 32. Jaszczak RJ, Chang L-T, Stein NA and Moore FE (1979) Whole-body Single-photon Emission Computed Tomography using Dual, Large-field-of-view Scintillation Cameras Phys Med Biol 24 1123–1143
    • 33. Juni J (1994) Taking Brain SPECT Seriously: Reflections on Recent Clinical Reports in The Journal of Nuclear Medicine J Nucl Med 35 1891–1895
    • 34. Kadrmas DJ, Frey EC and Tsui BMW (1997) Analysis of the Reconstructibility and Noise Properties of Scatter Photons in 99mTc SPECT Phys Med Biol 42 2493–2516
    • 35. Kadrmas DJ, Frey EC, Karimi SS and Tsui BMW (1998) Fast implementations of reconstruction-based scatter compensation in Jully 3D SPECT image reconstruction Phys Med Biol 43 857–873
    • 36. Kaplan MS, Miyaoka RS, Kohlmyer SK, Haynor DR, Harrison RL and Lewellen TK (1996) Scatter and Attenuation Correction for In-111 Imaging Using Energy Spectrum Fitting Med Phys 23 1277–1285
    • 37. Karbe H, Kertesz A, Davis J, Kemp BJ, Prato FS and Nicholson RL (1994) Quantification of functional deficit in Alzheimer's disease using a computer-assisted mapping program for 99mTc-HMPAO SPECT Neuroradiology 361-6
    • 38. Kemp BJ, Prato FS, Nicholson RL and Reese L (1995) Transmission Computed Tomography Imaging of the Head with a SPECT System and a Collimated Line Source J Nucl Med 36 328–335
    • 39. Lange K and Carson R (1984) EM Reconstruction Algorithms for Emission and Transmission Tomography J Comp Assist Tomo 8 306–316
    • 40. Lassen NA, Anderson AR, Friberg L and Paulson OB (1988) The Retention of 99mTc-d,I-HM-PAO in the human brain after intracarotid bolus injection: a kinetic analysis J Cereb Blood Flow Metab 8 S13–S22
    • 41. Leveille 1, Demonceau G and Waiovitch RC (1992) Intrasubject comparison between technetium-99m-ECD and technetium-99m-HMPAO in healthy human subjects J Nucl Med 33 480–484
    • 42. Licho R, Glick SJ, Xia W, Pan T-S, Penney BC and King MA (1999) Attenuation Compensation in 99mTc SPECT Brain Imaging: A Comparison of the Use of Attenuation Maps Derived from Transmission Versus Emission Data in Normal Scans J Nucl Med 40 456–463
    • 43. Mangin JF, Frouin V, Bloch I, Bendriem B and Lopez-Krahe J (1994) Fast nonsupervised three-dimensional registration of PET and MR images of the brain J Cereb Blood Flow Metab 14 749–762
    • 44. Malko JA, Van Heertum RL, Gullberg GT and Kowalsky WP (1986) SPECT liver imaging using an iterative attenuatian correction and an external flood source J Nucl Med 27 701–705
    • 45. Matsunari I, Boning G, Ziegler SI, Kosa I, Nekolla SG, Ficaro EP and Schwaiger M (1998) Effects of Misalignment Between Transmission and Emission Scans on Attenuation-Corrected Cardiac SPECT J Nucl Med 39 411–416
    • 46. McIntosh AR, Bookstein FL, Haxby JV and Grady CL (1996) Spatial Pattern Analysis of Functional Brain Images Using Partial Least Squares Neuroimage 3 143–157
    • 47. McCord ME, Bacharach SL, Bonow RO, Dilsizian V, Cuocolo A and Freedman N (1992) Misalignment Between PET Transmission and Emission Scans: Its Effects on Myocardial Imaging J Nucl Med 33 1209–1214
    • 48. Meikle SR, Dahlbom M and Cherry SR (1993) Attenuation Correction Using Count-Limited Transmission Data in Positron Emission Tomography J Nucl Med 34 143–150
    • 49. Meikle SR, Hutton BF and Bailey DL (1944) A Transmission-Dependent Method for Scatter Correction in SPECT J Nucl Med 35 360–367
    • 50. Mielke R, Pietrzyk U, Jacobs A, Fink GR, Ichimiya A, Kessler J, Herholz K and Heiss WD (1994) HMPAO SPET and FDG PET in Alzheimer's disease and vascular dementia: comparison of perfusion and metabolic pattem Eur J Nucl Med 21 1052–1060
    • 51. Minoshima S, Berger KL, Lee KS and Mintun MA (1992) An Automated Method for Rotational Correction and Centering of Three-Dimensional Functional Brain Images J Nucl Med 33 1579–1585
    • 52. Minoshima S, Koeppe RA, Frey KA and Kuhl DE (1994) Anatomic Standardization: Linear Scaling and Nonlinear Warping of Functional Brain Images J Nucl Med 35 1528–1537
    • 53. Msaki P, Prato FS, Stodilka RZ, Davis J, Kemp BJ, Kertesz A and Nicholson RL (1998) The Merits of Fully-Automated rCBF Determination in Eliminating Subjective Errors in Quantitative SPECT [Abstract] J Nucl Med 39 106P
    • 54. Natterer F (1993) Determination of tissue attenuation in emission tomography of optically dense media Inverse Problems 9 731–736
    • 55. Nelder JA and Mead R (1965) A Simplex method for function minimization Comput J 7 308–313
    • 56. Nunez J and Llacer J (1990) A Fast Bayesian Reconstruction Algorithm for Emission Tomography with Entropy Prior Converging to Feasible Images IEEE Trans Med Imag 9 159–171
    • 57. Panin VY, Zeng GL and Gullberg GT (1998) Reconstructions of Truncated Projections Using an Optimal Basis Expansion Derived from the Cross-Correlation of a "Knowledge Set" of a priori Cross-Sections IEEE Trans Nucl Sci 45 2119–2125
    • 58. Panin VY, Zeng GL and Gullberg GT (1999) Registration Problem in Attenuation Map Estimation using only Emission Data with Data Consistency Conditions [Abstract] J Nucl Med 40 294P
    • 59. Pellizari CA, Chen GTY, Spelbring DR, Weichselbaum RR and Chen C-T (1989) Accurate Three-Dimensional Registration of CT, PET, and/or MR Images of the Brain J Comp Assisi Tomo 13 20–26
    • 60. Patterson JC, Early TS, Martin A, Walker MZ, Russell JM and Villanueva-Meyer H (1997) SPECT Image Analysis Using Statistical Parametric Mapping: Comparison of Technetium-99m-HMPAO and Technetium-99m-ECD J Nucl Med 38 1721–1725
    • 61. Phillips RL, London ED, Links JM and Cascella NG (1990) Program for PET Image Alignment: Effects on Calculated Differences in Cerebral Metabolic Rates for Glucose J Nucl Med 31 2052–2057
    • 62. Pietrzyk U, Herholz K, Fink G, Jacobs A, Mielke R, Slansky I, Wurker M and Heiss W-D (1994) An Interactive Technique for Three-Dimensional Image Registration: Validation for PET, SPECT, MRI and CT Brain Studies J Nucl Med 35 2011–2018
    • 63. Rosenthal MS, Cullom J, Hawkins W, Moore SC, Tsui BMW and Yester M (1995) Quantitative SPECT Imaging: A Review and Recommendations by the Focus Committee of the Society of Nuclear Medicine Computer and Instrumentation Council J Nucl Med 36 1489–1513
    • 64. Seibyl JP, Woods SW, Zoghbi SS, Baldwin RM, Dey HM, Goddard AW, Zea-Ponce Y, Zubal G, Germine M, Smith EO, Heninger GR, Charney DS, Kung HF, Alavi A, Hoffer PB and Innis RB (1992) Dynamic SPECT Imaging of Dopamine D2 Receptors in Human Subjects with Iodine-123-IBZM J Nucl Med 33 1964–1971
    • 65. Shepp LA and Vardi Y (1982) Maximum likelihood reconstruction for emission tomography IEEE Trans Med Imag 1 113–122
    • 66. Slomka PJ, Hurwitz G, Stephenson JA and Cradduck TD (1995) Automated aligrnment and sizing of myocardial stress and rest scans to three-dimensional normal templates using an image registration algorithm: a method for reproducible quantification J Nucl Med 36 1115–1122
    • 67. Stodilka RZ, Kemp BJ, Prato FS and Nicholson RL (1998a) Importance of Bone Attenuation in Brain SPECT Quantification J Nucl Med 39 190–197
    • 68. Stodilka RZ, Kemp BJ, Msaki P, Prato FS and Nicholson RL (1998b) The Relative Contributions of Scatter and Attenuation Conections toward improved brain SPECT Quantification Phys Med Biol 43 2991-3008
    • 69. Stodilka RZ, Kemp BJ, Msaki P, Prato FS and Nicholson RL (1998c) Retrospective Scatter and Attenuation Correction by Inferring Anatomy from a Head Atlas IEEE Medical Imaging Conference [Abstract]
    • 70. Stokking R, van Isselt JW, van Rijk PP, de Klerk JMH, Huiskes WLC, Mertens IJR, Buskens E and Viergever MA (1999) Integrated Visualization of Functional and Anatomic Brain Data: A Validation Study J Nucl Med 40 311–316
    • 71. Sugiura M, Kawashima R, Sadato N, Senda M, Kanno I, Oda K, Sato K, Yonekura Y and Fukuda H (1999) Anatomic Validation of Spatial Normalization Methods for PET J Nucl Med 40 317–322
    • 72. Talairach J, Tournoux P and Rayport M, eds (1988) Co-planar stereotaxic atlas of the human brain. Three-dimensional proportional system: an approach to cerebral imaging New York: Thieme Inc 1-122
    • 73. Tan P, Bailey DL, Meikle SR, Eberl S, Fulton RR and Hutton BF (1993) A scanning line source for simultaneous emission and transmission measurements in SPECT J Nucl Med 34 1752–1760
    • 74. Ter-Pogossian MM (1992) Misalignment Between PET Transmission and Emission Scans: Effect on Myocardial Imaging J Nucl Med 33 1214–1215
    • 15. Tsao J, Stundzia A and Ichise M (1998) Fully Automated Establishment of Stereotaxic Image Orientation in Six Degrees of Freedom for Technetium-99m-ECD Brain SPECT J Nucl Med 39 503–508
    • 76. Tsui BMW, Gullberg GT, Edgerton ER, Ballard JG, Perry JR, McCartney WH and Berg J (1989) Correction of Nonuniform Attenuation in Cardiac SPECT Imaging J Nucl Med 30 497–507
    • 17. Tung C-H and Gullberg GT (1994) A simulation of emission and transmission noise propagation in cardiac SPECT imaging with nonuniform attenuation correction Med Phys 21 1565–1576
    • 78. Tung C-H, Gullberg GT, Zeng GL, Christian PE, Datz FL and Morgan HT (1992) Nonuniform attenuation conection using simultaneous transmission and emission converging tomography IEEE Trans Nucl Sci 39 1134–1143
    • 79. Van Dyck CH, Lin H, Smith EO, Wisniewski G, Cellar J, Robinson R, Narayan M, Bennett A, Delaney RC, Bronen RA and Hoffer PB (1996) Comparison of Technetium-99m-HMPAO and Technetium-99m-ECD Cerebral SPECT Images in Alzheimer's Disease J Nucl Med 37 1749–1755
    • 80. Webb S, Flower MA, Ott RJ and Leach MO (1983) A comparison of attenuation correction methods for quantitative single photon emission computed tomography Phys Med Biol 28 1045–1056
    • 81. Welch A, Clack R, Natterer F, and Gullberg GT (1997) Toward accurate attenuation correction in SPECT without transmission measurements IEEE Trans Med Imag 16 532–541
    • 82. Woods K and Bowyer KW (1997) Generating ROC Curves for Artificial Neural Networks IEEE Trans Med Imag 16 329–337
    • 83. Woods RP, Cherry SR and Mazziotta JC (1992) Rapid Automated Algorithm for Aligning and Reslicing PET Images J Comp Assist Tomo 16 620–633
    • 84. Woods RP, Mazziotta JC and Cherry SR (1993) MRI-PET Registration with Automated Algorithm J Comp Assist Tomo 17 536–546
    • 85. Xu EZ, Mulani NA, Gould KL and Anderson WL (1991) A Segmented Attenuation Correction for PET J Nucl Med 32161-165
    • 86. Zubal IG, Spencer SS, Imam K, Seibyl J, Smith EO, Wisniewski G and Hoffer PB (1995) Difference Images Calculated from Ictal and Interictal Technetium-99m-HMPAO SPECT Scans of Epilepsy J Nucl Med 36 684–689
    • 87. Zubal IG, Hanell CR, Smith EO, Rattner Z, Gindi GR and Hoffer PB (2994) Computerized three-dimensional segmented human anatomy Med Phys 21 299–302

Claims (19)

  1. Verfahren zum Anwenden einer Streuungs- und Abschwächungskorrektur auf Emissionstomografiebilder eines interessierenden Bereichs einer beobachteten Person, umfassend die Schritte: Ausrichten eines dreidimensionalen Computermodells, welches die Dichteverteilung innerhalb des interessierenden Bereichs repräsentiert, mit den Emissionstomografiebildern; und Anwenden von Streuungs- und Abschwächungskorrektur auf die Emissionstomografie unter Verwendung des ausgerichteten Computermodells als Führung.
  2. Verfahren nach Anspruch 1, wobei das Computermodell in der Form eines Zweikomponentenatlasses vorliegt.
  3. Verfahren nach Anspruch 2, wobei während des Ausrichtschritts eine funktionelle Komponente des Atlasses zuerst mit den Emissionstomografiebildern ausgerichtet wird, um einen Satz von räumlichen Transformationsparametern zu erzeugen und danach die anatomische Komponente des Atlasses mit den Emissionstomografiebildern unter Verwendung des Satzes räumlicher Transformationsparameter ausgerichtet wird.
  4. Verfahren nach Anspruch 3, wobei die funktionelle Komponente einen SPECT- oder PET-Scan des interessierenden Bereichs simuliert und wobei die anatomische Komponente einen Transmissionsscan des interessierenden Bereichs simuliert.
  5. Verfahren nach Anspruch 4, wobei der interessierende Bereich der Kopf ist und wobei die funktionelle Komponente die Gehirnkomponente eines Kopfatlasses ist.
  6. Verfahren nach Anspruch 4, wobei der interessierende Bereich das Herz ist und die funlctionelle Komponente des Atlasses ein Herzbild simuliert und die anatomische Komponente des Atlasses anatomische Merkmale des Thorax repräsentiert.
  7. Verfahren nach Anspruch 6, wobei die anatomischen Merkmale des Thorax beinhalten: Weichgewebe, wie etwa Herz, Leber, Muskel und Fett; Weichgewebe sehr niedriger Dichte, wie etwa die Lungen; und hochdichte Gewebe, wie etwa Knochen und Knorpel in den Rippen und der Wirbelsäule.
  8. Verfahren nach Anspruch 3, weiterhin umfassend den Schritt des Auswählens eines Atlasses aus einer Datenbank von Atlanten vor dem Durchführen des Ausrichtschritts.
  9. Verfahren nach Anspruch 8, wobei der Auswählschritt manuell durchgeführt wird.
  10. Verfahren nach Anspruch 8, wobei der Auswählschritt automatisch durchgeführt wird, basierend auf dem Deckungsgrad jedes Atlasses in der Datenbank mit den Emissionstomografiebildern.
  11. Verfahren nach Anspruch 10, wobei der Deckungsgrad bestimmt ist durch: Durchführen einer vorläufigen Rekonstruierung jedes Atlasses; und zur Deckung Bringen des Atlasses mit der vorläufigen Rekonstruktion.
  12. Verfahren nach Anspruch 10, weiterhin umfassend den Schritt des Kombinierens mehrerer Atlanten, um einen sich ergebenden Atlas zu erhalten, der besser in Deckung mit den Emissionstomografiebildern ist.
  13. Verfahren nach Anspruch 8, wobei die Datenbank Krankheits-spezifische Atlanten, Körpermerkmal-spezifische Atlanten und/oder Indikator- oder Läsions-spezifische Atlanten beinhaltet.
  14. Emissionstomografiebildverarbeitungssystem, umfassend: einen Speicher, der Emissionstomografiebilder eines interessierenden Bereichs einer Person speichert; wobei der Speicher zumindest ein dreidimensionales Computermodell des interessierenden Bereichs speichert, wobei das Computermodell die Dichteverteilung innerhalb des interessierenden Bereichs repräsentiert; und einen Prozessor zum In-Deckung-Bringen des Computermodells mit den Emissionstomografiebildern und zum Anwenden einer Streuungs- und Abschwächungskorrektur auf die Emissionstomografiebilder unter Verwendung des zur Deckung gebrachten Computermodells als eine Führung.
  15. Emissionstomografiebildverarbeitungssystem gemäß Anspruch 14, wobei das zumindest eine Computermodell in der Form zumindest eines Zweikomponenten-Atlasses vorliegt.
  16. Emissionstomografiebildverarbeitungssystem gemäß dem Anspruch 15, wobei der Prozessor zuerst eine funktionelle Komponente des Atlasses mit den Emissionstomografiebildern zur Deckung bringt, um einen Satz von räumlichen Transformationsparameter zu erzeugen, und dann die anatomische Komponente des Atlasses mit den Emissionstomografiebildern zur Deckung bringt, unter Verwendung des Satzes von räumlichen Transformationsparametern.
  17. Computerlesbares Medium, das Computerprogrammcode beinhaltet, zum Anwenden einer Streuungs- und Abschwächungskorrektur auf Emissionstomografiebilder eines interessierenden Bereichs einer Person, wobei das computerlesbare Medium enthält: Computerprogrammcode zum Ausrichten eines dreidimensionalen Computermodells, welches die Dichteverteilung innerhalb des interessierenden Bereichs mit den Emissionstomografiebildern repräsentiert; und Computerprogrammcode zum Anwenden von Streuungs- und Abschwächungskorrekturen auf die Emissionstomografiebilder unter Verwendung des ausgerichteten Computermodells als eine Führung, wobei der Computerprogrammcode zum Ausrichten beinhaltet: Computerprogrammcode zum Ausrichten einer funktionellen Komponente des Computermodells, welche einen SPECT- oder PET-Scan des interessierenden Bereiches simuliert und zum Erzeugen eines Satzes von räumlichen Transformationsparametern; und Computerprogrammcode zum Ausrichten einer anatomischen Komponente des Computermodells, welches einen Transmissionsscan des interessierenden Bereiches simuliert, unter Verwendung des Satzes von räumlichen Transformationsparametern.
  18. Verfahren nach einem der Ansprüche 1 bis 13, wobei das Computermodell aus Transmissionsbildern oder Röntgen-CT-Scans des Bereichs anderer Personen erzeugt ist.
  19. Verfahren nach Anspruch 18, wobei die Transmissionsbilder oder Röntgen-CT-Scans von einer Vielzahl anderer Personen genommen werden und ein Durchschnitt gebildet wird, um das Computermodell zu erzeugen.
DE69921154T 1998-08-14 1999-08-13 Anwendung von streu- und abschwächungs-korrekturen auf emissions-tomographie-bildern unter verwendung von einer atlas-abgeleiteten anatomie Expired - Lifetime DE69921154T2 (de)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US9664998P 1998-08-14 1998-08-14
US96649P 1998-08-14
PCT/CA1999/000751 WO2000010034A1 (en) 1998-08-14 1999-08-13 Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas

Publications (2)

Publication Number Publication Date
DE69921154D1 DE69921154D1 (de) 2004-11-18
DE69921154T2 true DE69921154T2 (de) 2005-11-10

Family

ID=22258387

Family Applications (1)

Application Number Title Priority Date Filing Date
DE69921154T Expired - Lifetime DE69921154T2 (de) 1998-08-14 1999-08-13 Anwendung von streu- und abschwächungs-korrekturen auf emissions-tomographie-bildern unter verwendung von einer atlas-abgeleiteten anatomie

Country Status (6)

Country Link
EP (1) EP1105750B1 (de)
AT (1) ATE279736T1 (de)
AU (1) AU5274599A (de)
CA (1) CA2338199C (de)
DE (1) DE69921154T2 (de)
WO (1) WO2000010034A1 (de)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005119025A2 (en) 2004-06-01 2005-12-15 Spectrum Dynamics Llc Radioactive-emission-measurement optimization to specific body structures
US8909325B2 (en) 2000-08-21 2014-12-09 Biosensors International Group, Ltd. Radioactive emission detector equipped with a position tracking system and utilization thereof with medical systems and in medical procedures
US7225011B2 (en) 2001-04-02 2007-05-29 Koninklijke Philips Electronics, N.V. Heart modeling using a template
FI113615B (fi) 2002-10-17 2004-05-31 Nexstim Oy Kallonmuodon ja sisällön kolmiulotteinen mallinnusmenetelmä
US20050077459A1 (en) * 2003-02-19 2005-04-14 Engler Mark J. Radiation phantom
US9470801B2 (en) 2004-01-13 2016-10-18 Spectrum Dynamics Llc Gating with anatomically varying durations
US7968851B2 (en) 2004-01-13 2011-06-28 Spectrum Dynamics Llc Dynamic spect camera
WO2008010227A2 (en) 2006-07-19 2008-01-24 Spectrum Dynamics Llc Imaging protocols
SE526940C2 (sv) * 2004-02-03 2005-11-22 Sidec Technologies Ab Förfarande och anordning för att med hög upplösning återskapa ett observerat objekt
US9316743B2 (en) 2004-11-09 2016-04-19 Biosensors International Group, Ltd. System and method for radioactive emission measurement
US9943274B2 (en) 2004-11-09 2018-04-17 Spectrum Dynamics Medical Limited Radioimaging using low dose isotope
US8837793B2 (en) 2005-07-19 2014-09-16 Biosensors International Group, Ltd. Reconstruction stabilizer and active vision
EP1971257B1 (de) * 2005-12-28 2015-10-14 Biosensors International Group, Ltd. Gating mit anatomisch variierender dauer
US8894974B2 (en) 2006-05-11 2014-11-25 Spectrum Dynamics Llc Radiopharmaceuticals for diagnosis and therapy
US7601966B2 (en) 2006-06-28 2009-10-13 Spectrum Dynamics Llc Imaging techniques for reducing blind spots
DE102006033383A1 (de) 2006-07-12 2008-01-17 Eberhard-Karls-Universität Tübingen Universitätsklinikum Verfahren zum Bestimmen einer Eigenschaftskarte für einen Gegenstand, insbesondere für ein Lebewesen, basierend auf zumindest einem ersten Bild, insbesondere Kernspinresonanzbild
US8660635B2 (en) 2006-09-29 2014-02-25 Medtronic, Inc. Method and apparatus for optimizing a computer assisted surgical procedure
WO2008075362A2 (en) 2006-12-20 2008-06-26 Spectrum Dynamics Llc A method, a system, and an apparatus for using and processing multidimensional data
DE102008032996B4 (de) 2008-07-14 2014-08-21 Siemens Aktiengesellschaft Verfahren zur Bestimmung einer Schwächungskarte
EP3277169B1 (de) * 2015-04-02 2021-07-07 HeartFlow, Inc. Systeme und verfahren zur abschätzung von virtuellen durchblutungsbildern
CN108271346B (zh) * 2015-10-30 2022-04-29 皇家飞利浦有限公司 成像控制器和成像系统
US20230129584A1 (en) * 2021-09-16 2023-04-27 University Of Pittsburgh-Of The Commonwealth System Of Higher Education Real-time, artificial intelligence-enabled analysis device and method for use in nuclear medicine imaging
CN116306326B (zh) * 2023-05-25 2023-09-15 南方医科大学珠江医院 一种关节接触力学仿真模型建立方法、装置及电子设备

Also Published As

Publication number Publication date
CA2338199A1 (en) 2000-02-24
CA2338199C (en) 2008-10-28
DE69921154D1 (de) 2004-11-18
AU5274599A (en) 2000-03-06
ATE279736T1 (de) 2004-10-15
EP1105750A1 (de) 2001-06-13
EP1105750B1 (de) 2004-10-13
WO2000010034A1 (en) 2000-02-24

Similar Documents

Publication Publication Date Title
DE69921154T2 (de) Anwendung von streu- und abschwächungs-korrekturen auf emissions-tomographie-bildern unter verwendung von einer atlas-abgeleiteten anatomie
US6740883B1 (en) Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas
Zaidi et al. Magnetic resonance imaging‐guided attenuation and scatter corrections in three‐dimensional brain positron emission tomography
Dey et al. Estimation and correction of cardiac respiratory motion in SPECT in the presence of limited‐angle effects due to irregular respiration
DE102007040071A1 (de) Bereichsweise Rekonstruktion von räumlich verteilten Funktionen
Verhoeff et al. Estimation of dopamine D2 receptor binding potential in the striatum with iodine-123-IBZM SPECT: technical and interobserver variability
Lobaugh et al. Three brain SPECT region-of-interest templates in elderly people: normative values, hemispheric asymmetries, and a comparison of single-and multihead cameras
DE102010004384A1 (de) Verfahren zur Ermittlung von der Berechnung eines Bestrahlungsplans zugrunde zu legenden Informationen und kombinierte Magnetresonanz-PET-Vorrichtung
DE102012105331A1 (de) Artefaktbeseitigung in Nuklearbildern
DE112015002809T5 (de) Multiple Emissionsenergien bei Einzelphotonen-Emissionscomputertomographie
Farncombe et al. The incorporation of organ uptake into dynamic SPECT (dSPECT) image reconstruction
Stodilka et al. Scatter and attenuation correction for brain SPECT using attenuation distributions inferred from a head atlas
Itti et al. Assessment of ejection fraction with Tl-201 gated SPECT in myocardial infarction: precision in a rest-redistribution study and accuracy versus planar angiography
Kim et al. 3-D SPECT simulations of a complex 3-D mathematical brain model: effects of 3-D geometric detector response, attenuation, scatter, and statistical noise
Jonsson et al. Resting state rCBF mapping with single-photon emission tomography and positron emission tomography: magnitude and origin of differences
Muellauer et al. 18F, 11C and 68Ga in small animal PET imaging
Leslie et al. External reference markers for the correction of head rotation in brain single-photon emission tomography
DE112020006198T5 (de) Vorrichtung, System, Verfahren und Computerprogramm zum Bereitstellen eines Kernbildes einer interessiernden Region eines Patienten
DE69625546T2 (de) Verbessertes Gammakamera-Abbildungssystem
Ichise et al. Neuroanatomical localization for clinical SPECT perfusion brain imaging: a practical proportional grid method
Véra et al. Advantage of fan beam collimators for contrast recovery of hyperfixation in clinical SPECT
Sase et al. A neural network for fusing the MR information into PET images to improve spatial resolution
McGrath et al. Maximum entropy deconvolution of low count nuclear medicine images
Piez Jr et al. Single photon emission computed tomography
Aubert-Broche et al. Realistic simulated MRI and SPECT databases

Legal Events

Date Code Title Description
8364 No opposition during term of opposition