DE102009021234B4

DE102009021234B4 - Verfahren zur Verarbeitung von Messdaten der Perfusions-Computertomographie

Info

Publication number: DE102009021234B4
Application number: DE102009021234A
Authority: DE
Inventors: Christian Canstein; Hendrik Ditt
Original assignee: Siemens AG
Current assignee: SIEMENS HEALTHINEERS AG, DE
Priority date: 2009-05-14
Filing date: 2009-05-14
Publication date: 2011-05-12
Anticipated expiration: 2029-05-15
Also published as: US20100290686A1; US9089308B2; DE102009021234A1

Abstract

Verfahren zur Verarbeitung von Messdaten der Perfusionscomputertomographie, bei dem
– aus den Messdaten sowohl Perfusionsparameter, für deren Berechnung arterielle und venöse Dichte-Zeit-Kurven (TACs) (8, 9) erzeugt werden, als auch dynamische CT-Angiographie-Bilder (5) berechnet werden und
– jeweils über mehrere dynamische CT-Angiographie-Bilder (5) zeitlich aufeinander folgender Zeitphasen gemittelt wird, um ein oder mehrere gemittelte CT-Angiographie-Bilder (6) zu erhalten,
– wobei die zeitlich aufeinander folgenden Zeitphasen für die Mittelung aus der oder den Dichte-Zeit-Kurven (8, 9) bestimmt werden,
dadurch gekennzeichnet,
dass für eine Segmentierung von Knochen aus den dynamischen (5) oder gemittelten CT-Angiographie-Bildern (6) zur Unterscheidung von Voxeln, die Knochen darstellen, und Voxeln, die Gefäße darstellen, eine zeitliche Änderung des HU-Wertes des jeweiligen Voxels über mehrere Zeitphasen ermittelt wird, wobei Voxel mit keiner oder nur einer innerhalb eines vorgegebenen Intervalls liegenden Variation des HU-Wertes den Knochen zugeordnet werden,
und/oder
dass für eine Segmentierung von...

Description

Die vorliegende Erfindung betrifft ein Verfahren zur Verarbeitung von Messdaten der Perfusions-Computertomographie, bei dem aus den Messdaten sowohl Perfusionsparameter, für deren Berechnung arterielle und venöse TACs (TAC: time attenuation curve) erzeugt werden, als auch dynamische CTA-Bilder berechnet werden.
Bei Untersuchungen mittels Computertomographie (CT) besteht grundsätzlich ein Interesse daran, sowohl die Menge an injiziertem Kontrastmittel als auch die Strahlendosis so gering wie möglich zu halten, um den Patienten keiner unnötigen Belastung auszusetzen. Dies ist jedoch nicht immer ohne weiteres möglich. So werden bspw. bei Untersuchung von Schlaganfallpatienten mittels Computertomographie nach einer Nativ-CT in der Regel zwei weitere CT-Scans unter Kontrastmittelgabe durchgeführt. Hierbei handelt es sich um eine dynamische Perfusionsmessung sowie eine daran anschließende CT-Angiographie (CTA). Für die dynamische Perfusionsmessung wird ein Injektionsprotokoll mit hoher Flussrate von mindestens 5 ml/s an Kontrastmittel verwendet, um einen kurzen, scharfen Kontrastmittelbolus zu produzieren, der Voraussetzung für die Berechnung von Perfusionsparametern, wie z. B. Blutfluss und Blutvolumen, ist. Bei der CT-Angiographie wird ein Kontrastmittel-Injektionsprotokoll verwendet, das einen längeren, breiten Kontrastmittelbolus produziert. Hier wird eine größere Gesamtmenge an Kontrastmittel verwendet und die Flussrate ist geringer. Das führt dazu, dass die Anreicherung des Kontrastmittels in den Gefäßen regelmäßiger erfolgt und längere Gefäßabschnitte in der Angiographie dargestellt werden können. Die CT-Angiographie dient dazu, den Ursprung einer Perfusionsstörung ermitteln zu können. Die beiden letztgenannten Untersuchungen führen zu einer nicht unerheblichen Menge an injiziertem Kontrastmittel und Strahlendosis, der der Patient ausgesetzt wird.
Neuere Techniken der Computertomographie, wie bspw. die sog. adaptive 4D Spirale (A4DS) der Firma Siemens ermöglichen Perfusionsmessungen des Schädels, die nicht mehr durch die Detektorbreite auf nur wenige Schichten begrenzt sind, sondern mit einer Abdeckung von bis zu 10 cm im gesamten Schädelvolumen durchgeführt werden können. Die Messdaten einer derartigen Perfusions-Computertomographie können nicht nur für die Berechnung von Perfusionsparametern, sondern auch für die Erzeugung von dynamischen CT-Angiographien eingesetzt werden, die das Ein- und Ausströmen des Kontrastmittels in den Gefäßen darstellen. Eine derartige dynamische Angiographie kann jedoch nicht die bisherige statische CT-Angiographie ersetzen, da in jeder durch die dynamische Angiographie dargestellten Zeitphase aufgrund des kurzen Kontrastmittelbolus nur ein kurzer Gefäßabschnitt mit Kontrastmittel angereichert ist.
Aus der US 2009/0028409 A1 ist ein Verfahren zur Verarbeitung von Messdaten der Perfusionscomputertomographie gemäß dem Oberbegriff von Patentanspruch 1 bekannt. Eine Extrahierung von Blutgefäßen erfolgt bei dieser Druckschrift mittels einer Schwellwerttechnik.
Die US 2007/0016016 A1 beschreibt unter anderem ein Verfahren zur Verarbeitung von Messdaten der Perfusionscomputertomographie, bei dem aus den Datensätzen einer dynamischen CT-Angiographie mittels MIP über den gesamten Zeitraum ein Bild oder durch Mittelung ein gemitteltes Bild erhalten wird. In P. Montes und G. Lauritsch: „Noise Reduction by Temporal Estimation in Perfusion Computed Tomography”, in: IEEE Nuclear Science Symposium Conference Record, 2005, M11-372, S. 2747–2751, sind Techniken zur Rauschreduzierung bei der Perfusionscomputertomographie beschrieben, wobei in einer Ausgestaltung eine zeitliche Glättung zur Verbesserung des Signal-Rauschverhältnisses durchgeführt wird. Die Glättung erfolgt beispielsweise durch Einsatz eines geeigneten Tiefpassfilters.
Die Aufgabe der vorliegenden Erfindung besteht darin, ein Verfahren anzugeben, durch das bei Durchführung einer Perfusionscomputertomographie in vielen Fällen auf eine anschließende statische CT-Angiographie für eine Gefäßdarstellung verzichtet werden kann und das eine eindeutige Identifizierung von Knochen- und/oder Gefäß-Voxeln für eine Segmentierung ermöglicht.
Die Aufgabe wird mit dem Verfahren zur Verarbeitung von Messdaten der Perfusions Computertomographie gemäß Patentanspruch 1 gelöst. Vorteilhafte Ausgestaltungen des Verfahrens sind Gegenstand der abhängigen Patentansprüche oder lassen sich der nachfolgenden Beschreibung sowie dem Ausführungsbeispiel entnehmen.
Bei dem Verfahren werden in bekannter Weise aus den Messdaten arterielle und venöse TACs erzeugt und daraus Perfusionsparameter eines Gewebes, insbesondere des Hirngewebes, berechnet. Die Messdaten werden auch zur Erzeugung von dynamischen CTA-Bildern genutzt. Das Verfahren zeichnet sich dadurch aus, dass jeweils über mehrere CTA-Bilder zeitlich aufeinander folgender Zeitphasen gemittelt wird, um ein oder mehrere gemittelte CTA-Bilder zu erhalten. Die zeitlich aufeinander folgenden Zeitphasen für die Mittelung werden dabei aus der oder den arteriellen und/oder venösen TACs bestimmt. Unter der Zeitphase wird der Zeitraum des jeweils zugrunde liegenden einzelnen Kontrastmittel-Scans verstanden, der einen Zustand der Kontrastmittelanflutung darstellt. Zusammen mit den vorhergehenden und nachfolgenden Kontrastmittel-Scans kann die An- und Abflutung des Kontrastmittels in den Gefäßen dynamisch dargestellt werden, wenn man sie hintereinander als Film zeigt. Man bezeichnet dies als dynamische CT-Angiographie.
Selbstverständlich werden in der Regel neben den arteriellen und venösen TACs auch Perfusionsparameter entsprechend der eingesetzten Perfusionssoftware generiert und dem Anwender dargestellt. Auch das oder die Bemittelten CTA-Bilder werden in geeigneter Weise an einem Bildanzeigegerät visualisiert.
Die Verarbeitung der Messdaten zur Erzeugung der Perfusionsparameter ist aus dem Stand der Technik derartiger Perfusionsmessungen bekannt. Hierzu stehen bereits zahlreiche kommerziell erhältliche Software-Pakete zur Verfügung. So geben bspw. K. A. Miles et al., „Perfusion CT: A worthwhile enhancement?”, The British Journal of Radiology, 76 (2003), 220 bis 231, einen Überblick über verfügbare Perfusionsverfahren, die auch in Softwarepaketen kommerziell erhältlich sind. TACs können dabei Voxel für Voxel oder gemittelt über vorab festgelegte Bereiche berechnet werden. In gleicher Weise ist die Berechnung von CTA-Bildern aus Computertomographie-Messdaten, die im Rahmen einer CT-Angiographie aufgezeichnet wurden, bekannt, so dass hierauf nicht näher eingegangen werden muss.
Bei dem vorliegenden Verfahren erfolgt nun eine Mittelung über mehrere CTA-Bilder zeitlich aufeinander folgender Zeitphasen, um dadurch den mit Kontrastmittel angereicherten Bereich in dem dargestellten Bild oder den dargestellten Bildern zu vergrößern. Dadurch wird ein Bildeindruck ähnlich dem CTA-Bild einer statischen CTA erreicht. Hierbei wird die Information aus den arteriellen oder venösen TACs genutzt, aus denen jeweils das Maximum der arteriellen oder venösen Kontrastmittelanreicherung ersichtlich ist. Für die Erzeugung eines arteriell gewichteten CTA-Bildes wird dann über Zeitphasen gemittelt, die um das Maximum der arteriellen TAC liegen. Für die Erzeugung eines venös gewichteten CTA-Bildes wird über Zeitphasen gemittelt, die um das Maximum der venösen TAC liegen. Für ein arteriell-venös gewichtetes CTA-Bild wird über Zeitphasen gemittelt, die um und zwischen den beiden Maxima liegen. Auf diese Weise werden jeweils optimale Bildergebnisse für die Mittelung erzielt. Unter Mittelung ist hierbei nicht ausschließlich die Bildung eines Mittelwertes aus mehreren Mess- oder Bilddaten sondern auch eine reine Addition dieser Daten zu verstehen.
In einer vorteilhaften Ausgestaltung kann die Anzahl der Zeitphasen, über die gemittelt wird, in Abhängigkeit von der Steilheit des Anstiegs oder Abfalls der jeweiligen arteriellen oder venösen TAC gewählt werden. Bei einem steileren Anstieg oder Abfall wird über eine geringere Anzahl an Zeitphasen gemittelt als bei einem flacheren Anstieg oder Abfall. Diese adaptive Wahl der Anzahl der zu mittelnden Zeitphasen berücksichtigt somit die zeitliche Auflösung, so dass bei einer hohen zeitlichen Dynamik des Kontrastmittelflusses auch eine höhere zeitliche Auflösung trotz Mittelung erreicht wird. Die Wahl der Zeitphasen, über die gemittelt wird, kann ebenso wie die Wahl ihrer Anzahl für die Mittelung automatisiert durch ein Verarbeitungsprogramm erfolgen, das die jeweiligen TACs geeignet auswertet. Für die Wahl der Anzahl der Zeitphasen muss dieses Programm Zugriff auf eine Tabelle oder eine Funktion haben, in der diese Anzahl mit der Steilheit der jeweiligen TAC verknüpft ist.
In einer weiteren Ausgestaltung erfolgt nicht nur eine Erzeugung eines gemittelten CTA-Bildes für Zeitphasen im Bereich des Maximums der jeweiligen TAC. Vielmehr werden durch Mittelung über Zeitphasen vor und/oder nach dem jeweiligen Maximum zusätzliche gemittelte CTA-Bilder erzeugt. In Zeitphasen, in denen die TAC einen flachen Verlauf zeigt, etwa deutlich nach dem Maximum, kann bspw. auch über eine größere Anzahl an Zeitphasen als im Bereich des Maximums gemittelt werden.
Die Mittelung über die dynamischen CTA-Bilder mehrerer Zeitphasen hat den weiteren Vorteil einer Verbesserung der Bildqualität auch bei einer zeitaufgelösten Darstellung. Wenn bspw. anstelle der Darstellung der dynamischen CTA-Bilder einzelner Zeitphasen jeweils ein gemitteltes Bild von bspw. drei aufeinander folgenden Zeitphasen dargestellt wird, erreicht man dadurch eine weichere Darstellung und ein besseres Signal-Rausch-Verhältnis. Derartige gemittelte CTA-Bilder können bspw. in zeitlicher Abfolge an einem Bildanzeigegerät visualisiert werden. Vorzugsweise ist dabei die Anzahl der Zeitphasen, über die gemittelt wird, durch einen Anwender interaktiv einstellbar. Dadurch kann der Anwender je nach klinischer Fragestellung und vorliegenden Daten jeweils das optimale Verhältnis aus gewünschter Zeitauflösung und Bild qualität einstellen.
Mit dem vorgeschlagenen Verfahren lassen sich aus den Volumenperfusionsdaten, d. h. aus den Messdaten der Perfusionscomputertomographie, arterielle, venöse und arteriell-venöse Gefäßdarstellungen erzeugen, die der Darstellung einer herkömmlichen statischen CTA entsprechen oder zumindest sehr nahe kommen. Das Verfahren nutzt dabei geschickt Ergebnisse aus, insbesondere die TACs, die bei einer Volumenperfusionsmessung mittels Computertomographie ohnehin erzeugt werden. Mit dem Verfahren kann in Fällen von Schlaganfalluntersuchungen, bei denen feststeht, dass keine Gefäßpathologie außerhalb des mit der Perfusionscomputertomographie abgetasteten Schädelbereiches Grund für eine Perfusionsstörung ist, auf eine separat durchgeführte statische CTA verzichtet werden. Damit wird sowohl die Menge an injiziertem Kontrastmittel und die Strahlendosis gegenüber einer herkömmlichen Untersuchung verringert, als auch die Untersuchungsdauer im CT verkürzt, was speziell bei Schlaganfalluntersuchungen von großer Bedeutung ist. Anders als bei der Darstellung der einzelnen Zeitphasen einer dynamischen CTA zeigen die optimierten Gefäßdarstellungen der gemäß dem Verfahren gemittelten CIA-Bilder die Gefäße in der Art und Weise, wie es der Anwender von herkömmlichen statischen CTA-Bildern gewohnt ist. Der kurze Kontrastmittelbolus der Perfusionscomputertomographie wird quasi gestreckt. Gegenüber den einzelnen Zeitphasen der dynamischen CTA wird durch die Mittelung über mehrere Phasen außerdem das Signal-Rausch-Verhältnis verbessert, sowohl für die so erzeugten arteriellen, venösen und arteriell-venösen Gefäßdarstellungen als auch für die durch Mittelung optimierten Gefäßdarstellungen.
Bei dem vorgeschlagenen Verfahren zur Verarbeitung von Messdaten der Perfusionscomputertomographie werden für die Segmentierung bzw. Entfernung der Knochen aus dem gemittelten oder dynamischen CTA-Datensatz zeitliche Informationen genutzt. Da Knochen und kontrastierte Gefäße oft sehr ähnliche HU-Werte haben, lässt sich die für die Entfernung der Knochen aus einem CTA-Datensatz erforderliche Unterscheidung zwischen Knochen und Gefäßen bei herkömmlichen CT Angiographien oft nicht eindeutig vornehmen. Bei der vorgeschlagenen Ausgestaltung wird für einen dynamischen, registrierten CTA-Datensatz die HU-Information jedes einzelnen Voxels über der Zeit genutzt, um eine eindeutige Unterscheidung zwischen Voxeln, die Knochen darstellen, und Voxeln, die Gefäße darstellen, zu erreichen. Verändert sich der HU-Wert des Voxels über mehrere Zeitphasen deutlich, so liegt dies am Ein- und Ausströmen des Kontrastmittels und das Voxel gehört zu einem Gefäß. Verändern sich die HU-Werte des Voxels dagegen nicht oder nur kaum, so handelt es sich eindeutig um ein Knochen-Voxel. Der maximale Variationsbereich der HU-Werte für eine Zuordnung des entsprechenden Voxel zu einem Knochen wird dabei durch den Anwender vorgegeben. Diese Technik lässt sich nicht nur dafür verwenden, Knochen-Voxel eindeutig zu erkennen und sie zu entfernen, sondern auch zur eindeutigen Identifizierung von Gefäß-Voxeln, bspw. für eine Segmentierung des Gefäßbaums.
In einer vorteilhaften Weiterbildung des vorgeschlagenen Verfahrens zur Verarbeitung von Messdaten der Perfusionscomputertomographie werden die zeitaufgelösten CTA-Daten auch für die Berechnung funktionaler Parameter für die Gefäße eingesetzt. Damit lassen sich zusätzliche Informationen über Gefäßpathologien ermitteln, insbesondere dynamische Parameter ähnlich den durch Dopplersonographie oder Phasenkontrast-Magnetresonanztomographie erhaltenen Parametern. Hierzu wird entweder wie bei Perfusionsmessungen für jedes Pixel innerhalb des Gefäßes jeweils eine TAC bestimmt, oder eine TAC z. B. für mehrere Gefäß-Pixel in einer Schicht bestimmt, also gemittelt. Aus diesen TACs lassen sich dann weitere Parameter ermitteln, wie bspw. die Höhe eines Peaks (Maximum), die Breite des Peaks (FWHM: full width half maximum), der zeitliche Abstand der arteriellen Kurve zur venösen Kurve in 2 verschiedenen Gefäßen, der Abstand des arteriellen zum venösen Maximums, die Steigung der jeweiligen Kurven, das Integral über die arterielle oder venöse Kurve, die Anzahl der Peaks usw.. Diese Parameter lassen sich ähnlich den Perfusionsparameterdarstellungen farbig auf den Gefäßen darstellen oder alternativ in Diagrammform, bspw. die Höhe des arteriellen Peaks in Abhängigkeit von der Position entlang der Gefäßmittellinie, visualisieren und liefern so eine potentiell klinisch relevante Information zusätzlich zur reinen Gefäßgeometrie.
Bei allen Ausgestaltungen des vorgeschlagenen Verfahrens ist eine Registrierung der Messdaten bzw. Bilder der einzelnen Zeitphasen erforderlich. Eine derartige Registrierung wird von der kommerziell verfügbaren Perfusionssoftware bereits automatisch durchgeführt, weil dies auch Vorraussetzung für die Berechnung der Perfusionsparameter ist. Im Schädel sind die Gefäße relativ statisch und bewegen sich während des Flusszykluses wenig, so dass hier eine gute Registrierung der Gefäße in den verschiedenen Zeitphasen aufeinander möglich ist. In der Regel wird hierzu eine rigide Registrierungstechnik genutzt.
Mit der bei dem Verfahren eingesetzten Technik der Unterscheidung von Knochen-Voxeln und Gefäß-Voxeln auf Basis der zeitlichen Entwicklung der HU-Werte wird eine eindeutige Trennung von Knochen- und Gefäß-Voxeln auch in den problematischen Bereichen der Schädelbasis erreicht. Diese Technik kann auch als Grundlage für die Gefäßsegmentierung dienen. Durch die Berechnung der Gefäßparameter basierend auf den dynamischen CTA-Bildern wird zusätzlich zur reinen Gefäßgeometrie auch Aufschluss über die Gefäßfunktion ermöglicht. So kann basierend auf den Messdaten der Perfusionscomputertomographie nicht nur die Ursache der Pathologie, bspw. eine Stenose, sondern potentiell auch die Auswirkung auf die Flussfunktion, bspw. ein veränderter Anstieg der TAC-Kurve, dargestellt werden.
Das vorgeschlagene Verfahren wird in einem Datenverarbeitungssystem durchgeführt, das in bekannter Weise mindestens einen Rechenprozessor, einen Speicher und entsprechende Schnittstellen für die Datenein- und Datenausgabe aufweist. Das Datenverarbeitungssystem umfasst hierbei ein oder mehrere Programme, die die einzelnen Verfahrensschritte zur Verarbeitung der eingelesenen oder gespeicherten Messdaten gemäß dem vorgeschlagenen Verfahren ausführen. Ein derartiges Datenverarbeitungssystem kann direkt mit einer Computertomographie-Anlage verbunden sein, um die Ergebnisse im unmittelbaren Anschluss an die Perfusionscomputertomographie an einem Bildanzeigegerät darzustellen.
Das vorgeschlagene Verfahren wird nachfolgend in Verbindung mit den Zeichnungen nochmals kurz erläutert. Hierbei zeigen:
1 eine schematische Darstellung eines Workflows für die Durchführung des vorgeschlagenen Verfahrens und
2 ein Beispiel für TACs, die bei dem Verfahren berechnet werden.
1 zeigt ein Beispiel für den Workflow bei der Untersuchung eines Schlaganfall-Patienten mittels Computertomographie. Hierbei wird zunächst ein nativer CT-Scan 1 durchgeführt, um Blutungen im Schädelbereich zu erkennen (Blutungsausschluss 12). Im Anschluss erfolgt eine Volumenperfusions-Computertomographie 2 gemäß dem vorgeschlagenen Verfahren. Aus den Messdaten dieser Perfusions-CT 2 werden zunächst nach einer Registrierung 3 der Einzelphasen arterielle und venöse TACs 4 erzeugt. Selbstverständlich werden hierbei auch Perfusionsparameter berechnet. Im vorliegenden Beispiel werden aus den Messdaten der Perfusions-CT 2 auch zeitaufgelöste (dynamische) CTA-Bilder 5, arteriell und venös gewichtete CTA-Bilder 6 und funktionelle Parameter der Gefäße berechnet (funktionelle CTA 7). Auf diese Weise kann auf den Einsatz einer anschließenden statischen CTA verzichtet werden, solange die Ursache für die Perfusionsstörung innerhalb des von der Perfusions-CT 2 abgedeckten Volumenbereiches liegt.
2 zeigt ein Beispiel für eine arterielle 8 und venöse TAC 9, wie sie bei der Auswertung der Messdaten der Perfusionscomputertomographie automatisch erzeugt werden. Weiterhin sind in dieser Abbildung die ebenfalls automatisch erzeugten TACs für das Gewebe der linken 10 und rechten Hirnhälfte 11 dargestellt. Diese Information wird bei dem vorliegenden Verfahren genutzt, um die Zeitphasen für eine Mittelung der CTA-Bilder auszuwählen. Diese TACs bieten die Information darüber, welche Zeitphasen idealerweise für eine arterielle, welche für eine venöse und welche für eine arteriell-venöse Gefäßdarstellung aufaddiert werden müssen. So hat bspw. die arterielle TAC 8 ihr Maximum bei Zeitphase 12. Für eine optimierte arterielle Gefäßdarstellung können dann bspw. die CTA-Bilder der Zeitphasen 10 bis 14 addiert werden.

Claims

Verfahren zur Verarbeitung von Messdaten der Perfusionscomputertomographie, bei dem – aus den Messdaten sowohl Perfusionsparameter, für deren Berechnung arterielle und venöse Dichte-Zeit-Kurven (TACs) (8, 9) erzeugt werden, als auch dynamische CT-Angiographie-Bilder (5) berechnet werden und – jeweils über mehrere dynamische CT-Angiographie-Bilder (5) zeitlich aufeinander folgender Zeitphasen gemittelt wird, um ein oder mehrere gemittelte CT-Angiographie-Bilder (6) zu erhalten, – wobei die zeitlich aufeinander folgenden Zeitphasen für die Mittelung aus der oder den Dichte-Zeit-Kurven (8, 9) bestimmt werden, dadurch gekennzeichnet, dass für eine Segmentierung von Knochen aus den dynamischen (5) oder gemittelten CT-Angiographie-Bildern (6) zur Unterscheidung von Voxeln, die Knochen darstellen, und Voxeln, die Gefäße darstellen, eine zeitliche Änderung des HU-Wertes des jeweiligen Voxels über mehrere Zeitphasen ermittelt wird, wobei Voxel mit keiner oder nur einer innerhalb eines vorgegebenen Intervalls liegenden Variation des HU-Wertes den Knochen zugeordnet werden, und/oder dass für eine Segmentierung von Gefäßen aus den dynamischen (5) oder gemittelten CT-Angiographie-Bildern (6) zur Unterscheidung von Voxeln, die Gefäße darstellen, und Voxeln, die andere Bereiche darstellen, eine zeitliche Änderung des HU-Wertes des jeweiligen Voxels über mehrere Zeitphasen ermittelt wird, wobei Bereiche mit einer oberhalb eines vorgegebenen Intervalls liegenden Variation des HU-Wertes den Gefäßen zugeordnet werden.
Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass für die Erzeugung eines arteriell gewichteten CT-Angiographie-Bildes über Zeitphasen gemittelt wird, die um ein Maximum der arteriellen Dichte-Zeit-Kurve (8) liegen.
Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass für die Erzeugung eines venös gewichteten CT-Angiographie-Bildes über Zeitphasen gemittelt wird, die um ein Maximum der venösen Dichte-Zeit-Kurve (9) liegen.
Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass die Anzahl der Zeitphasen, über die gemittelt wird, in Abhängigkeit von einer Steilheit des Anstiegs oder Abfalls der jeweiligen Dichte-Zeit-Kurve (8, 9) gewählt wird, wobei bei einem steileren Anstieg oder Abfall über eine geringere Anzahl an Zeitphasen gemittelt wird als bei einem flacheren Anstieg oder Abfall.
Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass die Wahl der Anzahl der Zeitphasen automatisiert durch ein Verarbeitungsprogramm erfolgt.
Verfahren nach einem der Ansprüche 1 bis 5, dadurch gekennzeichnet, dass für die Erzeugung von CT-Angiographie-Bildern mit verbesserter Bildqualität über Zeitphasen gemittelt wird, die um ein Maximum der venösen und/oder arteriellen Dichte-Zeit-Kurven (8, 9) liegen, wobei die Anzahl der Zeitphasen, über die gemittelt wird, interaktiv einstellbar ist.
Verfahren nach Anspruch einem der Ansprüche 1 bis 6, dadurch gekennzeichnet, dass für die Erzeugung weiterer gemittelter CT-Angiographie-Bilder auch über jeweils aufeinanderfolgende Zeitphasen vor und nach dem Maximum der Dichte-Zeit-Kurven (8, 9) gemittelt wird, wobei die erhaltenen Bilder in zeitlicher Abfolge dargestellt werden.
Verfahren einem der Ansprüche 1 bis 7, dadurch gekennzeichnet, dass aus den dynamischen CT-Angiographie-Bildern (5) Dichte-Zeit-Kurven von Gefäßen ermittelt werden, aus denen funktionelle Parameter für die Gefäße berechnet werden.
Verfahren nach Anspruch 8, dadurch gekennzeichnet, dass ein oder mehrere der Parameter in einer Bilddarstellung der dynamischen CT-Angiographie-Bilder (5) als farbige Darstellung der Gefäße oder in einer separaten Diagrammdarstellung visualisiert werden.