Nothing Special   »   [go: up one dir, main page]

DE102013201136A1 - Vorhersage eines voraussichtlichen Kontrastmittelverlaufs - Google Patents

Vorhersage eines voraussichtlichen Kontrastmittelverlaufs Download PDF

Info

Publication number
DE102013201136A1
DE102013201136A1 DE102013201136.9A DE102013201136A DE102013201136A1 DE 102013201136 A1 DE102013201136 A1 DE 102013201136A1 DE 102013201136 A DE102013201136 A DE 102013201136A DE 102013201136 A1 DE102013201136 A1 DE 102013201136A1
Authority
DE
Germany
Prior art keywords
heart rate
impulse response
response function
measurement
contrast agent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
DE102013201136.9A
Other languages
English (en)
Other versions
DE102013201136B4 (de
Inventor
Johannes Georg KORPORAAL
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.)
Siemens Healthineers Ag De
Original Assignee
Siemens AG
Siemens Corp
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 Siemens AG, Siemens Corp filed Critical Siemens AG
Priority to DE102013201136.9A priority Critical patent/DE102013201136B4/de
Priority to US14/149,919 priority patent/US9271656B2/en
Priority to CN201410022486.9A priority patent/CN103961117B/zh
Publication of DE102013201136A1 publication Critical patent/DE102013201136A1/de
Application granted granted Critical
Publication of DE102013201136B4 publication Critical patent/DE102013201136B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • A61B5/0275—Measuring blood flow using tracers, e.g. dye dilution
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/488—Diagnostic techniques involving pre-scan acquisition
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus 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/507—Apparatus 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024—Measuring pulse rate or heart rate
    • A61B5/0245—Measuring pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus 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/503—Apparatus 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 diagnosis of the heart
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus 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/504—Apparatus 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 diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • G16H20/17—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients delivered via infusion or injection

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Physiology (AREA)
  • Cardiology (AREA)
  • Hematology (AREA)
  • Signal Processing (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Dentistry (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Business, Economics & Management (AREA)
  • General Business, Economics & Management (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

Die Erfindung betrifft ein Verfahren zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs (CEP) für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts (2) mit folgenden Schritten: – Erfassung von patientenspezifischen Blutflussverlaufsdaten (35) in einer Vormessung unter Verwendung eines definierten Untersuchungsprotokolls (47) und Bestimmung einer individuellen Impulsantwortfunktion (37, 43) aus diesen Blutflussverlaufsdaten (35), – Erfassung eines der Vormessung zugeordneten Vormessungs-Herzratensignals (40) und – Vorhersage (39) des voraussichtlichen Kontrastmittelverlaufs (CEP) auf Basis der individuellen Impulsantwortfunktion (37, 43), des Vormessungs-Herzratensignals (40) sowie eines aktuellen Untersuchungsprotokolls (48) und eines aktuell erfassten Herzratensignals (44, 52). Außerdem betrifft die Erfindung ein Verfahren zur Ansteuerung eines medizintechnischen Bildgebungssystems (1). Weiterhin werden eine Vorrichtung (60) zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs (CEP), eine Steuereinrichtung (12) mit einer solchen Vorrichtung (60) für ein Bildgebungssystem (1) und ein Bildgebungssystem (1) umfassend eine solche Steuereinrichtung (12) beschrieben.

Description

  • Die Erfindung betrifft ein Verfahren und eine Vorrichtung zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts, insbesondere eines kontrastmittelunterstützten Bildgebungsverfahrens, bevorzugt im Rahmen eines angiographischen Untersuchungsverfahrens eines Patienten. Weiterhin betrifft die Erfindung ein Verfahren zur Ansteuerung eines medizintechnischen Bildgebungssystems, eine Steuereinrichtung für ein Bildgebungssystem und ein Bildgebungssystem mit einer solchen Steuereinrichtung.
  • In der älteren Patentanmeldung DE 10 2012 209 410.5 ist ein Verfahren zur Ermittlung einer individuellen, patientenspezifischen Kontrastmittel-Impulsantwortfunktion (im Folgenden auch kurz als "Impulsantwortfunktion" oder "Patientenfunktion" bezeichnet) anhand von Testbolus-Daten beschrieben, mit der später eine Vorhersage des Kontrastmittelverlaufs bei einer kontrastmittelunterstützten bildgebenden Messung erzeugt werden kann. Diese Patientenfunktion beschreibt die kardiovaskulären Eigenschaften des Patienten zum Zeitpunkt, als der Testbolus gemessen wurde. Im Prinzip könnte man damit den Kontrastmittelverlauf für jedes beliebige Injektionsprotokoll unter der Annahme vorhersagen, dass die individuelle Patientenfunktion zu dem späteren Zeitpunkt noch gültig ist. Ein Injektionsprotokoll ist dabei eine genaue zeitliche Vorgabe, gemäß der an den Patienten die Kontrastmittelmenge abgegeben wird, und umfasst z. B. den Startzeitpunkt, die Flussrate und den Endzeitpunkt der Kontrastmittelgabe.
  • Allgemein sind derartige sogenannte CEP-Algorithmen (engl. contrast enhancement prediction) aus verschiedenen Publikationen bekannt, wie beispielsweise:
  • Diese CEP-Algorithmen berechnen den voraussichtlichen Verlauf des Kontrastmittels im Patienten und sind meist aus zwei Subalgorithmen aufgebaut. Mit einem ersten Subalgorithmus wird eine patientenspezifische Kontrastmittel-Impulsantwortfunktion (engl. arterial impulse response, AIR) ermittelt, der sogenannte "AIR-Algorithmus". Der zweite Subalgorithmus berechnet eine Vorhersage (engl. prediction) des voraussichtlichen Kontrastmittelverlaufs im Patienten, der sogenannte "PRED-Algorithmus". Ein CEP-Algorithmus umfasst also meist eine Kombination aus einem AIR-Algorithmus und einem PRED-Algorithmus.
  • In allen diesen bekannten Verfahren wird der zu untersuchende Patient als ein lineares zeitinvariantes System (engl. LTIsystem) betrachtet. Dies bedeutet im Allgemeinen, dass man die in der Zeit gemessene Kontrastmittelanreicherung im Patienten C(t) mathematisch als eine Faltung (engl. convolution) des Injektionsprotokolls IF(t) mit der Impulsantwort- bzw. Patientenfunktion AIR(t) beschreiben kann: C(t) = IF(t) ⊗ AIR(t) [1].
  • Infolgedessen führen die CEP-Algorithmen alle Berechnungen in der Zeit- und/oder Fourier-Domäne durch, die von der Zeit-Domäne abgeleitet ist.
  • Dies setzt jedoch voraus, dass die treibende Kraft des Kontrastmittels, der Blutkreislauf, auch invariant und über die Zeit stabil ist. Dies ist jedoch oft nicht der Fall.
  • Ein typisches Anwendungsfeld für kontrastmittelunterstützte Bildgebungsverfahren ist z. B. die CT-Angiographie (CTA), d. h. die Bildgebung der Gefäße mittels Computertomographie. Ein CTA-Scan ist dabei ein Durchlauf eines Computertomographen für die Bildgebung der Gefäße, und eine Herz-CTA oder ein Herz-Scan ist die Bildgebung des Herzens und Darstellung der Herzgefäße.
  • Insbesondere bei solchen Herz-Scans in der CT-Angiographie (CTA) kann es sein, dass die Herzrate (Anzahl der Herzschläge pro Zeiteinheit) sich nach dem Testbolus dauerhaft erhöht oder gesenkt hat, beispielsweise unter Einfluss von Betablocker (z. B. Metoprolol), die typischerweise kurz vor dem Scan bei Patienten mit einer zu hohen Herzrate verabreicht werden. Da sich die Rezirkulationszeit des Blutes mit einer anderen Herzfrequenz ändert, hat sich in diesem Fall die Patientenfunktion auch geändert. Dies würde zu einer ungenauen Vorhersage des Kontrastmittelverlaufs führen.
  • Die Erfindung geht von der Aufgabe aus, ein Verfahren und eine Vorrichtung zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs für eine kontrastmittelunterstützte Untersuchung, insbesondere eines kontrastmittelunterstützten Bildgebungsverfahrens, bevorzugt im Rahmen eines angiographischen Untersuchungsverfahrens, der eingangs genannten Art derart auszubilden, so dass vor allem bei unregelmäßigen Herzraten eine bessere Vorhersage des Kontrastmittelverlaufs zu einem späteren Zeitpunkt ermöglicht wird, insbesondere falls sich die Herzrate zwischen der Durchführung des Testbolus und der eigentlichen späteren Untersuchung geändert haben sollte.
  • Diese Aufgabe wird zum einen durch ein Verfahren gemäß Patentanspruch 1 und zum anderen durch eine Vorrichtung gemäß Patentanspruch 12 gelöst.
  • Bei dem erfindungsgemäßen Verfahren werden zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts (im Folgenden ohne Beschränkung der Allgemeinheit als "Patient" bezeichnet, weshalb auch der Begriff "patientenspezifisch" stellvertretend für den Begriff "untersuchungsobjektspezifisch" verwendet wird) folgende Schritte durchgeführt:
    • – Zunächst werden in einer Vormessung patientenspezifische Blutflussverlaufsdaten unter Verwendung eines definierten Untersuchungsprotokolls erfasst und aus den Blutflussverlaufsdaten eine individuelle Impulsantwortfunktion für den Patienten bestimmt. Die Blutflussverlaufsdaten können dabei im Rahmen einer bildgebenden Messung, z. B. einer vorherigen CT-Messung, insbesondere CTA, unter Verwendung von entsprechendem Kontrastmittel wie in der späteren Untersuchung gewonnen werden. Dabei kann es sich um eine "echte" frühere kontrastmittelunterstützte Untersuchung handeln. In der Regel handelt es sich aber um eine spezielle Testmessung, die mit nur einer kleinen Menge Kontrastmittel (dem "Testbolus") durchgeführt wird. Die Blutflussverlaufsdaten können aber auch in anderer Weise bestimmt werden, beispielsweise optisch mit einem Testbolus in Form eines Mittels zur farblichen Markierung des Blutflussverlaufs, wie dies von F. Eisa et al in "Optical tracking of contrast medium bolus to optimize bolus shape and timing in dynamic computed tomography", Physics in Medicine and Biology 2012, 57, Seiten N173 bis N182, beschrieben ist. Das Untersuchungsprotokoll umfasst dabei z. B. das definierte Injektionsprotokoll zur Gabe des Kontrastmittels bzw. Testbolus.
    • – Außerdem wird ein der Vormessung zugeordnetes Vormessungs-Herzratensignal erfasst. Die Formulierung "der Vormessung zugeordnet" ist dabei so zu verstehen, dass das Vormessungs-Herzratensignal in ausreichender Weise, je nach vorhandenen technischen Möglichkeiten und genauer Ausführungsvariante der Erfindung, die zu einem Zeitpunkt der Vormessung bzw. während der Vormessung vorliegende Herzrate des Patienten repräsentiert. Zum Beispiel kann das Vormessungs-Herzratensignal im Rahmen der Vormessung, d. h. bei der Vormessung oder kurz vorher oder danach erfasst werden.
    • – Die Vorhersage des voraussichtlichen Kontrastmittelverlaufs erfolgt dann auf Basis der individuellen Impulsantwortfunktion, des Vormessungs-Herzratensignals (direkt oder – indem das Vormessungs-Herzratensignal bereits in die individuelle Impulsantwortfunktion eingegangen ist – indirekt) sowie eines aktuellen Untersuchungsprotokolls und eines aktuell erfassten Herzratensignals. Das aktuell erfasste Herzratensignal repräsentiert dabei die Herzrate bzw. den Herzschlag des Patienten zum Zeitpunkt der späteren Messung, für die der Kontrastmittelverlauf vorhergesagt werden soll. Es kann also beispielsweise in ausreichend kurzem Abstand vor oder bei der Untersuchung erfasst werden. Es wird an dieser Stelle ausdrücklich darauf hingewiesen, dass eine Kontrastmittelverlaufs-Vorhersage auch den Fall umfasst, dass während der späteren Untersuchung, beispielsweise der CTA, von Herzschlag zu Herzschlag in einer Art aktueller "Online-Vorausschau" vorhergesagt wird, wie sich das Kontrastmittel im Blutkreislauf ausbreitet. Dies wird später noch im Einzelnen ausgeführt.
  • Die Erfassung des Vormessungs-Herzratensignals und des Herzratensignals kann sich dabei sowohl auf jeden einzelnen Herzschlag (HB) beziehen, wenn z. B. die Herzrate mittels eines EKGs gemessen und aufgezeichnet worden ist, wie auch auf Informationen zu einer mittleren Herzrate (HR). Letztes kann beispielsweise bei Pulsmessern der Fall sein, die die Herzrate als Anzahl der Herzschläge pro Minute (engl. beats per minute, bpm) wiedergeben.
  • Mit den Herzraten- bzw. Herzschlag-Informationen bei der Vormessung kann die Patientenfunktion so gebildet werden, dass sie die Herzrate zum Zeitpunkt der Vormessung explizit berücksichtigt. Anschließend kann man dann im Falle einer Änderung der Herzrate den Kontrastmittelverlauf zu einem späteren Zeitpunkt besser vorhersagen. Durch die Erfassung und Aufzeichnung bzw. Speicherung des Vormessungs-Herzratensignals und durch die Erfassung des aktuellen Herzratensignals ergibt sich also in der Regel eine erheblich verbesserte Vorhersage eines voraussichtlichen Kontrastmittelverlaufs für nicht-konstante Herzraten.
  • Ein solches Verfahren kann dann wiederum bevorzugt in einem erfindungsgemäßen Verfahren zur automatischen Ansteuerung eines medizintechnischen Bildgebungssystems eingesetzt werden, wobei ein Referenzzeitpunkt – beispielsweise ein Start-Zeitpunkt zum Start einer Bildakquisition in einem bestimmten Bereich – für das Bildgebungssystem auf Basis eines Kontrastmittelgabe-Startzeitpunkts einer Kontrastmittelgabe und des vorhergesagten Kontrastmittelverlaufs bestimmt wird. Beispielsweise kann sich eine Verzögerungszeit zwischen dem Start der Kontrastmittelvergabe und einem Start der Abtastung durch das Bildgebungssystems aus der Lage eines lokalen Maximums eines vorhergesagten Kontrastmittelverlaufs ergeben.
  • Die Ermittlung eines solchen Referenzzeitpunkts kann aber auch für andere Verfahren sinnvoll sein, beispielsweise für eine spätere Auswertung von Bilddatensätzen einer Kontrastmittelaufnahme. So ist es in vielen Fällen wesentlich zu wissen, in welcher Phase bestimmte Bilder aufgenommen wurden, beispielsweise ob es sich um eine Phase handelt, in der sich das Kontrastmittel erst in einer bestimmten Struktur verteilt bzw. anreichert, oder ob die Bilder in einer Phase erzeugt wurden, in der das Kontrastmittel wieder aus der Struktur verschwindet, d. h. abgebaut wird.
  • Eine erfindungsgemäße Vorrichtung zur automatischen Vorhersage eines voraussichtlichen Kontrastmittelverlaufs für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts weist zumindest folgende Komponenten auf:
    • – Eine Eingangs-Schnittstelle zur Erfassung von patientenspezifischen Blutflussverlaufsdaten in einer Vormessung unter Verwendung eines definierten Untersuchungsprotokolls. Hierbei kann es sich beispielsweise um eine Schnittstelle zu einem Speicher handeln, in dem die betreffenden Daten hinterlegt sind, aber auch um eine Schnittstelle zu einer Auswerteeinheit, in der aktuell gerade die Messdaten aus einer Testbolus-Messung ausgewertet werden und von der die Daten direkt übernommen werden können.
    • – Eine Impulsantwortfunktions-Ermittlungseinheit zur Bestimmung einer individuellen Impulsantwortfunktion auf Basis der Blutflussverlaufsdaten.
    • – Eine Eingangs-Schnittstelle zur Erfassung eines Herzratensignals. Auch hier kann es sich beispielsweise um eine Schnittstelle zu einem Speicher handeln oder um eine Schnittstelle zu einem Messgerät zur Messung bzw. Aufzeichnung des Herzratensignals, beispielsweise eine Schnittstelle zu einem EKG. Grundsätzlich können die Eingangs-Schnittstellen zur Erfassung der patientenspezifischen Blutflussverlaufsdaten und zur Erfassung eines Herzratensignals auch als gemeinsame Schnittstelle ausgebildet sein, sofern die Daten bzw. Signale gemeinsam hinterlegt sind.
    • – Eine Berechnungseinheit zur Vorhersage des voraussichtlichen Kontrastmittelverlaufs auf Basis der individuellen Impulsantwortfunktion, eines der Vormessung zugeordneten Vormessungs-Herzratensignals (wobei das Vormessungs-Herzratensignal wieder direkt oder – indem es bereits in die individuelle Impulsantwortfunktion eingegangen ist – indirekt genutzt wird) sowie eines aktuellen Untersuchungsprotokolls und eines aktuell erfassten Herzratensignals.
  • Die Kontrastmittelverlaufs-Vorhersage bzw. eine entsprechende Funktion oder Tabelle, welche den zu erwartenden Kontrastmittelverlauf repräsentiert, kann dann über eine Ausgangs-Schnittstelle der erfindungsgemäßen Vorrichtung beispielsweise an andere Einheiten ausgegeben werden, die mit der Kontrastmittelverlaufs-Vorhersage weiter arbeiten, insbesondere weitere Berechnungen durchführen oder diese speichern.
  • Eine solche Vorrichtung zur automatischen Vorhersage des voraussichtlichen Kontrastmittelverlaufs ist besonders bevorzugt als Teil einer Steuereinrichtung für ein Bildgebungssystem ausgebildet. Das heißt, die Steuereinrichtung weist selbst diese Vorrichtung zur automatischen Vorhersage des voraussichtlichen Kontrastmittelverlaufs auf. In diesem Fall weist die Steuereinrichtung bevorzugt auch eine Schnittstelle zur Erfassung oder Ausgabe eines Kontrastmittelgabe-Start-Zeitpunkts sowie eine Referenzzeitpunkt-Ermittlungseinrichtung zur Ermittlung eines Referenzzeitpunkts für das Bildgebungssystem auf Basis des Kontrastmittelgabe-Startzeitpunkts und des voraussichtlichen Kontrastmittelverlaufs auf. Dies hat dann den Vorteil, dass die Steuereinrichtung vollautomatisch unmittelbar die Kontrastmittelverlaufs-Vorhersage auf Basis der Vormessung ermitteln und später entsprechend das Bildgebungssystem mit Hilfe dieser Steuereinrichtung so gesteuert werden kann, dass bei der eigentlichen Messung mit der Bilddatenakquisition zum richtigen Zeitpunkt nach einer Kontrastmittelgabe begonnen werden kann und/oder die Bilddaten zeitlich mit dem Referenzzeitpunkt verknüpft und gespeichert werden können.
  • Grundsätzlich kann eine erfindungsgemäße Vorrichtung zur automatischen Vorhersage eines voraussichtlichen Kontrastmittelverlaufs aber auch auf anderen Rechnern realisiert werden, die beispielsweise die benötigten Blutflussverlaufsdaten von anderen Komponenten übernehmen und beispielsweise mit dem Bildgebungssystem über ein Netzwerk verbunden sind. Dies ist insbesondere dann sinnvoll, wenn rechenintensivere Aufgaben aus der Steuereinrichtung verlagert werden sollen, um die Rechenkapazität der Steuereinrichtung voll zur Steuerung anderer Messungen zur Verfügung zu stellen.
  • Ein erfindungsgemäßes Bildgebungssystem zum Erzeugen von Bilddaten vom Inneren eines Untersuchungsobjekts umfasst vorzugsweise eine zuvor beschriebene Steuereinrichtung, welche auch eine erfindungsgemäße Vorrichtung zur automatischen Vorhersage eines voraussichtlichen Kontrastmittelverlaufs aufweist.
  • Insbesondere können die Impulsantwortfunktions-Ermittlungseinheit und die Berechnungseinheit der erfindungsgemäßen Vorrichtung zur automatischen Vorhersage eines voraussichtlichen Kontrastmittelverlaufs jeweils als Softwaremodule auf einem geeigneten Rechner mit entsprechenden Speichermöglichkeiten realisiert sein. Die Eingangs-Schnittstellen sowie die Ausgangs-Schnittstelle können ebenfalls in Form von reiner Software realisiert sein, sofern nur eine Übernahme der Blutflussverlaufsdaten bzw. eine Ausgabe der Kontrastmittelverlaufs-Vorhersage von anderen bzw. an andere – insbesondere auf der gleichen Rechnereinheit realisierte – weitere Einheiten oder Speicher erforderlich ist. Grundsätzlich können diese Schnittstellen aber auch als kombinierte Hardware-/Software-Schnittstellen realisiert sein, um eine externe Ein- und Ausgabe zu realisieren, beispielsweise mit Hilfe von Softwarekomponenten speziell konfigurierte Hardware-Schnittstellen. Eine weitgehend softwaremäßige Realisierung hat den Vorteil, dass beispielsweise auch schon bisher verwendete Steuereinrichtungen für Bildgebungssysteme auf einfache Weise durch ein Software-Update nachgerüstet werden können, um auf die erfindungsgemäße Weise zu arbeiten. Insofern wird die Aufgabe auch durch ein Computerprogrammprodukt gelöst, welches direkt in einen Speicher einer programmierbaren Steuereinrichtung ladbar ist, mit Programmabschnitten, um alle Schritte des erfindungsgemäßen Verfahrens auszuführen, wenn das Programm in der Steuereinrichtung ausgeführt wird.
  • Die abhängigen Ansprüche enthalten jeweils besonders vorteilhafte Weiterbildungen und Ausgestaltungen der Erfindung, wobei die Ansprüche einer Kategorie auch entsprechend den abhängigen Ansprüchen einer anderen Kategorie weitergebildet sein können.
  • Abhängig von den Herzraten-Informationen, d. h. der Art der Erfassung bzw. Aufzeichnung des Herzratensignals und/oder Implementation der CEP-Algorithmen, ergeben sich ein oder mehrere Möglichkeiten, die Vorhersage des Kontrastmittelverlaufs für schwankende Herzraten zu optimieren.
  • Bei einer bevorzugten Variante wird auf Basis des Vormessungs-Herzratensignals eine herzschlagbezogene Impulsantwortfunktion ermittelt. Diese individuelle "herzschlagbezogene Impulsantwortfunktion" definiert die individuelle Impulsantwortfunktion in Abhängigkeit von der Herzschlaganzahl. Mit anderen Worten: Während die oben erläuterte, bisher übliche individuelle Impulsantwortfunktion (im Folgenden auch zur Unterscheidung manchmal explizit "zeitbezogene Impulsantwortfunktion" genannt) in der Zeitdomäne betrachtet wird, d. h. eine Funktion über der Zeit t ist, wird eine "herzschlagbezogene" Impulsantwortfunktion in der "Herzratendomäne" bzw. "Herzschlagdomäne" betrachtet, d. h. es handelt sich um eine Funktion über der Herzschlaganzahl, die im Folgenden mit der Abkürzung HB (heart beat) benannt wird.
  • Hierzu ist es möglich, die patientenspezifischen Blutflussverlaufsdaten und das Untersuchungsprotokoll jeweils bereits in der Herzschlagdomäne zu ermitteln. Zum Beispiel können die Daten jeweils gleich in Abhängigkeit von einem parallel ermittelten Herzschlag aufgezeichnet werden. Das Untersuchungsprotokoll könnte auch in Abhängigkeit von der Herzschlaganzahl vorgegeben sein bzw. herzschlaggesteuert abgearbeitet werden. Es ist aber ebenso möglich, alle Daten zunächst in der Zeitdomäne zu erfassen und dann in die Herzschlagdomäne umzurechnen.
  • Analog zu der Gleichung [1] gilt auch in der Herzschlagdomäne C(HB) = IF(HB) ⊗ AIR(HB) [2] wobei (HB) hier wie gesagt die Abhängigkeit vom Herzschlag symbolisiert. Die herzschlagbezogene Impulsantwort- bzw. Patientenfunktion AIR(HB) lässt sich also auch bei bekannten patientenspezifischen Blutflussverlaufsdaten C(HB), d. h. der Kontrastmittelanreicherung im Patienten in Abhängigkeit vom Herzschlag, und bekanntem Injektionsprotokolls IF(HB) mathematisch durch eine Entfaltung ermitteln.
  • Da aber für die Entfaltung in der Zeitdömäne bzw. Frequenzdomäne schon eine Vielzahl von Standardprogrammen zur Verfügung stehen, kann es hinsichtlich des Gesamtaufwands vorteilhafter sein, zunächst in der üblichen Weise eine zeitbezogene individuelle Impulsantwortfunktion AIRTB(t) zu ermitteln und diese unter Nutzung des Vormessungs-Herzratensignals in eine herzschlagbezogene Impulsantwortfunktion AIRTB(HB) umzurechnen. Dies ist z. B. mittels folgender Gleichung möglich: AIRTB(HB) = AIRTB(t·HRTB) [3] HRTB ist hierbei die zu der Vormessung ermittelte Herzrate. Dabei und im Folgenden wird allgemein der Index "TB" (Testbolus) verwendet, um anzugeben, dass es sich um Daten aus der Vormessung handelt, unabhängig davon, ob hierzu ein Testbolus verwendet wurde. Das heißt, AIRTB ist die gesuchte individuelle Impulsantwortfunktion, die auf Basis der Daten der Vormessung ermittelt wurde.
  • Wie bereits oben erwähnt, kann die Herzrate unterschiedlich berechnet bzw. ermittelt werden. Häufig handelt es sich um eine mittlere Herzrate, die sich durch die Anzahl der Herzschläge in einem definierten Zeitraum ergibt. Dabei ist es möglich, den Zeitraum über alle bisherigen Herzschläge seit Messung der Herzrate zu legen und die Anzahl der gemessenen Herzschläge durch diesen Zeitraum zu teilen. Ebenso ist es auch möglich, eine Art gleitenden Mittelwert zu bilden, der sich über einen bestimmten zurückliegenden Zeitraum oder eine bestimmte zurückliegende Anzahl von Herzschlägen (beispielsweise die letzten zehn Herzschläge) erstreckt. Eine solche gleitende Mittelwertbildung wird beispielsweise bei einer üblichen Pulsmessung mit einem einfachen Pulsmesser implizit durchgeführt. Im Extremfall kann auch eine ganz aktuelle Herzrate nur auf dem letzten Herzschlag beruhen, indem als Herzrate einfach der Kehrwert des zeitlichen Abstands zwischen dem letzten und vorletzten Herzschlag verwendet wird.
  • Sofern zur Ermittlung der Herzrate HRTB nicht die gleiche Zeiteinheit verwendet wird wie zur Ermittlung der Zeit t, muss in Gleichung [3] zusätzlich ein Umrechnungsfaktor U berücksichtigt werden:
    Figure DE102013201136A1_0002
  • Wenn wie üblich die Zeit t in der Regel in der Einheit Sekunden (s) gemessen wird und die Herzrate in Herzschlägen pro Minute (bpm; beats per minute), gilt für den Umrechnungsfaktor U = 60, d. h. es gilt dann die Gleichung:
    Figure DE102013201136A1_0003
  • Ist die individuelle herzschlagbezogene Impulsantwortfunktion AIRTB(HB) bekannt und wird bei der späteren kontrastmittelunterstützten Messung für die eigentliche Untersuchung das Untersuchungsprotokoll, insbesondere Injektionsprotokoll IFCTA(HB) ebenfalls in der Herzschlagdomäne gesteuert oder in die Herzschlagdomäne umgerechnet (ein mögliches Verfahren hierzu wird später noch angegeben), so lässt sich auch die gesuchte Kontrastmittelverlaufs-Vorhersage CEP(HB) ohne weiteres mit der folgenden Gleichung (vgl. Gleichung [1]) direkt durch die übliche Faltungsoperation ermitteln: CEP(HB) = IFCTA(HB) ⊗ AIRTB(HB) [4]
  • Die Funktion CEP(HB) definiert hierbei letztlich die herzschlagbezogene Kontrastmittelverlaufs-Vorhersage, also eine erwartete Kontrastmittelanreicherung an einer bestimmten Position im Patienten in Abhängigkeit vom Herzschlag. Der Index "EP" steht für die Vorhersage der Anreicherung (Enhancement Prediction). Außerdem wird allgemein zur Veranschaulichung der Index "CTA" verwendet, um anzugeben, dass es sich um Daten der späteren kontrastmittelunterstützten Messung handelt, für die der Kontrastmittelverlauf vorherzusagen ist, und zwar unabhängig davon, ob es sich hierbei tatsächlich um eine CTA-Untersuchung oder um eine andere kontrastmittelunterstützte Untersuchung handelt.
  • Die Verwendung der Gleichung [4] ist insbesondere vorteilhaft, wenn während der späteren kontrastmittelunterstützten Messung das Herzratensignal in einer Weise erfasst und bearbeitet wird, dass Herzschlag für Herzschlag beobachtet wird, um die Vorhersage effektiv zu nutzen.
  • Sofern jedoch bei der späteren kontrastmittelunterstützten Messung nicht in der Herzschlagdomäne gearbeitet wird, bietet sich folgende alternativ bevorzugte Variante des erfindungsgemäßen Verfahrens an. Hierbei wird z. B. auf Basis eines (kurz vor oder während des Hauptscans, d. h. der eigentlichen Messung) aktuell erfassten Herzratensignals und der herzschlagbezogenen Impulsantwortfunktion eine zeitbezogene korrigierte individuelle Impulsantwortfunktion gebildet. Alternativ kann eine solche zeitbezogene korrigierte individuelle Impulsantwortfunktion auch direkt aus der zeitbezogenen individuellen Impulsantwortfunktion ermittelt werden. In beiden Fällen erfolgt dann die Vorhersage des voraussichtlichen Kontrastmittelverlaufs CEP(t) in der Zeitdomäne durch Kombination der korrigierten zeitbezogenen Impulsantwortfunktion AIRCORR(t) mit dem aktuellen (ebenfalls in üblicher Weise auf die Zeit bezogenen) Untersuchungsprotokoll bzw. Injektionsprotokoll IFCTA(t) gemäß CEP(t) = IFCTA(t) ⊗ AIRCORR(t) [5]
  • Bei diesen Varianten erfolgt also zur Untersuchung des Patienten unter Einsatz eines Kontrastmittels eine Korrektur einer patientenspezifischen arteriellen Kontrastmittel-Impulsantwortfunktion für nicht-konstante Herzraten zur verbesserten Vorhersage des voraussichtlichen Kontrastmittelverlaufs. Dabei können in vorteilhafter Weise in einer Variante Daten, die das Herzschlagsignal repräsentieren, und die Blutflussverlaufsdaten einem CEP-Algorithmus zugeführt werden, der aus diesen Werten die individuelle Impulsantwortfunktion berechnet. Aus dieser individuellen Impulsantwortfunktion und einer kurz vor der eigentlichen Untersuchung ermittelten CTA-Herzrate wird dann die korrigierte individuelle Impulsantwortfunktion abgeleitet. Oder es kann aus einer zunächst herkömmlich ermittelten zeitbezogenen Impulsantwortfunktion mittels der Daten, die das Herzratensignal bei der Vormessung repräsentieren, und der Daten, die ein kurz vor der eigentlichen kontrastmittelunterstützten Untersuchung ermitteltes Herzratesignal repräsentieren, als für die durchzuführende Untersuchung erforderliche Parameter die korrigierte individuelle Impulsantwortfunktion berechnet werden.
  • Sofern zunächst eine individuelle herzschlagbezogene Impulsantwortfunktion AIRTB(HB) ermittelt wurde, kann daraus die korrigierte individuelle Impulsantwortfunktion berechnet werden, indem die individuelle herzschlagbezogene Impulsantwortfunktion mit einem Zeitskalierungsfaktor skaliert wird, welcher durch den Kehrwert einer aus dem aktuell erfassten Herzratensignal ermittelten Herzrate HRCTA gegeben ist: AIRCORR(t) = AIRTB(HB/HRCTA) [6]
  • Wenn die Zeit t selbst und die Herzrate mit unterschiedlichen Zeiteinheiten ermittelt werden, kann wieder ein Umrechnungsfaktor U verwendet werden: AIRCORR(t) = AIRTB(HB/(HRCTA/U)) = AIRTB(HB/HRCTA·U) [6'] wobei bei den üblichen Einheiten (einerseits s und andererseits bpm) U = 60 gilt, wie dies oben in Zusammenhang mit Gleichung [3] schon erläutert wurde.
  • Wurde dagegen nur eine individuelle Impulsantwortfunktion AIRTB(t) wie üblich in der Zeitdomäne ermittelt, kann daraus die korrigierte individuelle Impulsantwortfunktion berechnet werden, indem die individuelle zeitbezogene Impulsantwortfunktion mit einem Zeitskalierungsfaktor skaliert wird, welcher durch eine aus dem Vormessungs-Herzratensignal ermittelten Herzrate HRTB im Verhältnis zu einer aus dem aktuell erfassten Herzratensignal ermittelten Herzrate HRCTA gegeben ist: AIRCORR(t) = AIRTB(t·HRTB/HRCTA) [7]
  • In allen Gleichungen [6] bis [7] kann die Herzrate wie oben erläutert wieder eine mittlere Herzrate oder eine nur auf einen einzelnen Herzschlag bezogene Herzrate sein.
  • In den obigen Gleichungen [6] bis [7] wurde jedoch noch nicht berücksichtigt, dass die Änderung de Herzrate letztlich auch dazu führt, dass bei einem bezüglich der Zeitdomäne unveränderten Injektionsprotokoll unterschiedliche Kontrastmittelmengen pro Herzschlag gegeben werden. Mit anderen Worten, eine Änderung der Herzrate führt zu einer Änderung des Injektionsprotokolls in der Herzschlagdomäne. Um dies zu berücksichtigen und die Vorhersage noch weiter zu verbessern, kann bei den obigen Gleichungen [6] bis [7] jeweils ein Korrekturfaktor hinzugenommen werden, der dies ausgleicht. Daher wird bevorzugt zur Bildung der korrigierten individuellen Impulsantwortfunktion die skalierte individuelle Impulsantwortfunktion mit dem Kehrwert des jeweiligen Zeitskalierungsfaktors multipliziert. Für Gleichung [6] erhält man dann AIRCORR(t) = (HRCTA)·AIRTB(HB/HRCTA) [6a] und für die Gleichungen [6'] und [7] entsprechend: AIRCORR(t) = (HRCTA/U)·AIRTB(HB/HRCTA·U) [6'a] AIRCORR(t) = (HRCTA/HRTB)·AIRTB(t·HRTB/HRCTA) [7a]
  • Anschaulich betrachtet, sorgt dieser Korrekturfaktor (HRCTA), (HRCTA/U) oder (HRCTA/HRTB) bzw. Vorfaktor in den zugehörigen Gleichungen jeweils dafür, dass bei einer Anpassung der AIR-Kurve in die Breite die Amplitude der individuellen Impulsantwortfunktion AIR sich reziprok ändert, damit die Fläche unter der AIR-Kurve konstant bleibt.
  • Setzt man beispielsweise Gleichung [7a] in Gleichung [5] ein, dann ergibt sich die Vorhersage des voraussichtlichen Kontrastmittelverlaufs CEP(t) auf Basis der üblichen zeitbezogenen individuellen Impulsantwortfunktion AIRTB(t) und der Herzrate HRTB bei der Vormessung sowie der aktuellen Herzrate HRCTA bei der späteren Messung wie folgt: CEP(t) = IFCTA(t) ⊗ (HRCTA/HRTB)·AIRTB(t·HRTB/HRCTA) [8]
  • Für die anderen Verfahren gemäß den Gleichungen [6] und [6'] erfolgt die Berechnung analog.
  • Alternativ kann der Korrekturfaktor (HRCTA/HRTB), welcher die Herzratenabhängigkeit des Untersuchungsprotokolls berücksichtigt, auch mit dem Untersuchungsprotokoll bzw. Injektionsprotokoll IFCTA(t) multipliziert werden: CEP(t) = (HRCTA/HRTB)·IFCTA(t) ⊗ AIRTB(t·HRTB/HRCTA) [9]
  • Da für die Faltungsoperation allgemein a·f ⊗ g = (a·f) ⊗ g = f ⊗ (a·g) [10] gilt, führt dies zum gleichen Ergebnis. Gleichung [9] gibt dabei ein Beispiel für die Berechnungsvariante mit Gleichung [7]. Dies gilt aber ebenso für die Varianten mit den Gleichungen [6] und [6'].
  • Wie bereits oben erwähnt, kann das Vormessungs-Herzratensignal auf verschiedene Weise erfasst werden. Es ist lediglich notwendig, dass es in ausreichender Weise die Herzrate während der Vormessung repräsentiert. Um eine möglichst genaue Korrektur zu erlauben, werden die Blutflussverlaufsdaten und das Vormessungs-Herzratensignal bevorzugt zeitlich korreliert, d. h. zeitlich parallel erfasst und aufgezeichnet. Besonders bevorzugt erfolgt dabei die Erfassung des Vormessungs-Herzratensignals durch die parallele Aufzeichnung von EKG-Daten während der Vormessung.
  • Implementationsbedingt kann dabei auch aus den EKG-Daten eine während der Vormessung vorliegende mittlere Herzrate abgeleitet werden, wenn beispielsweise bei der späteren Messung kein EKG vorhanden ist, sondern ohnehin nur mit der mittleren Herzrate gearbeitet werden kann. Dies erleichtert dann die Berechnungen, indem beispielsweise gemäß den obigen Gleichungen [7] oder [7a] eine Korrektur der individuellen Impulsantwortfunktion bzw. gemäß den Gleichungen [8] oder [9] eine Korrektur der Kontrastmittelverlaufs-Vorhersage erfolgt.
  • Die Erfindung wird nachfolgend anhand von in der Zeichnung dargestellten Ausführungsbeispielen näher erläutert. Es zeigen:
  • 1 in perspektivischer schematischer Darstellung ein erfindungsgemäßes Bildgebungssystem, welches zur Ausführung eines erfindungsgemäßen angiographischen Untersuchungsverfahrens zur Untersuchung eines Patienten unter Einsatz eines Kontrastmittels geeignet ist,
  • 2 eine während einer Voruntersuchung ermittelte Kontrastmittelkurve,
  • 3 ein bekanntes Verfahren zur Vorhersage des Kontrastmittelverlaufs in der Zeitdomäne gemäß dem Stand der Technik,
  • 4 ein Diagramm zur Veranschaulichung der erfindungsgemäßen Umrechnung der abgeleiteten Patientenfunktion AIRTB(t) als Funktion der Zeit in eine Funktion AIRTB(HB) der Anzahl der Herzschläge,
  • 5 ein Diagramm zur Erläuterung des Zusammenhangs zwischen Ereignissen in der Zeitdomäne und der Herzschlagdomäne,
  • 6 ein Diagramm mit einem Injektionsprotokoll eines Testbolus und EKG-Daten sowie einem auf eine Herzrate umgerechneten Injektionsprotokoll,
  • 7 eine Darstellung einer ersten Ausführungsform eines erfindungsgemäßen Verfahrens zur Vorhersage des Kontrastmittelverlaufs,
  • 8 eine Darstellung einer zweiten Ausführungsform eines erfindungsgemäßen Verfahren zur Vorhersage des Kontrastmittelverlaufs,
  • 9 ein Diagramm mit einem Beispiel der zweiten Ausführungsform zur Anpassung der abgeleiteten Patientenfunktion AIRTB(t) an veränderte Herzraten,
  • 10 ein Diagramm zur Erläuterung der Verbesserung der Vorhersage durch die erfindungsgemäße Korrektur gemäß der zweiten Ausführungsform und
  • 11 eine Darstellung einer dritten Ausführungsform eines erfindungsgemäßen Verfahrens zur Vorhersage des Kontrastmittelverlaufs in der Herzratendomäne.
  • In den folgenden Beispielen wird der Einfachheit halber und ohne Beschränkung der Allgemeinheit davon ausgegangen, dass es sich bei der Vormessung um eine Testbolus-Messung und bei der eigentlichen kontrastmittelunterstützten Untersuchung um eine CTA handelt.
  • In 1 ist ein medizintechnisches Bildgebungssystem (bzw. eine Diagnostikeinrichtung), in diesem Beispiel ein Computertomographiegerät 1, gezeigt, welches zur Ausführung des erfindungsgemäßen Untersuchungsverfahrens zur Untersuchung eines Patienten 2 unter Einsatz eines Kontrastmittels 3 geeignet ist. Anstelle des beschriebenen Computertomographiegeräts 1 kann z. B. auch ein Magnetresonanz- oder Ultraschalltomographiegerät eingesetzt werden.
  • Dem Computertomographiegerät 1 ist ein Patientenlagerungstisch 4 mit einer beweglichen Tischplatte 5 zugeordnet, auf der der Patient 2 lagerbar ist. Die Tischplatte 5 ist in Richtung einer Drehachse 6 verstellbar, so dass ein mit dem Patienten 2 verbundener Untersuchungsbereich 7 durch eine Öffnung einer Gantry 8 des Computertomographen 1 in den Messbereich eines Aufnahmesystems 9, 10 bewegt werden kann. Der Patient 2 und das Aufnahmesystem 9, 10 sind auf diese Weise in Richtung der Drehachse 6 relativ zueinander verstellbar, so dass unterschiedliche Abtastpositionen eingenommen werden können. Der Untersuchungsbereich 7 kann beispielsweise das Herz des Patienten sein.
  • Zur Erfassung von Projektionen weist das Aufnahmesystem 9, 10 einen Röntgenstrahler 9 in Form einer Röntgenröhre und einen diesem gegenüberliegend angeordneten Röntgendetektor 10 auf, wobei der Röntgendetektor 10 bogenförmig ausgebildet ist und mehrere zu Detektorzeilen aufgereihte Detektorelemente 11 umfasst. Der Röntgenstrahler 9 erzeugt eine Strahlung in Form eines fächerförmigen Röntgenstrahlenbündels, welches den Messbereich durchdringt und anschließend auf die Detektorelemente 11 des Röntgendetektors 10 auftrifft. Die Detektorelemente 11 erzeugen einen von der Schwächung der durch den Messbereich tretenden Röntgenstrahlung abhängigen Schwächungswert. Die Umwandlung der Röntgenstrahlung in einen Schwächungswert erfolgt beispielsweise jeweils mittels einer mit einem Szintillator optisch gekoppelten Photodiode oder mittels eines direkt konvertierenden Halbleiters. Der Röntgendetektor 11 erzeugt auf diese Weise einen Satz von Schwächungswerten, welcher auch als Projektion bezeichnet wird.
  • Das Aufnahmesystem 9, 10 ist an der Gantry 8 drehbar angeordnet, so dass Projektionen aus unterschiedlichen Projektionsrichtungen erfasst werden können. Je nach eingestellter Betriebsart des Computertomographiegerätes 1 erfolgt die Abtastung bei fest eingestellter oder veränderlicher Projektionsrichtung bei gleichzeitig fest eingestellter oder veränderlicher Abtastposition. Durch Rotation der Gantry 8 bei gleichzeitigem kontinuierlichem Vorschub des Patienten 2 in Richtung der Drehachse 6 werden beispielsweise Projektionen aus einer Vielzahl von unterschiedlichen Projektionsrichtungen an verschiedenen Positionen entlang der Drehachse 6 bzw. entlang des Patienten 2 erfasst. Die auf diese Weise durch eine spiralförmige Abtastung – Drehung des Aufnahmesystem 9, 10 und gleichzeitiger Vorschub der Tischplatte 5 – gewonnenen Projektionen des Aufnahmesystems 9, 10 werden an eine Recheneinheit bzw. Steuereinrichtung 12 des Bildgebungssystems 1 übermittelt und zu Bilddaten rekonstruiert, welche auf einem Display 13 darstellbar sind. Die Bilddaten können beispielsweise eine oder mehrere Schicht- oder Volumenbilder des Untersuchungsbereichs 7 umfassen.
  • Zur Untersuchung durchbluteter Organe, beispielsweise eines Herzens, einer Leber oder eines Gefäßes, kann dem Patienten 2 zur Erhöhung des sichtbaren Kontrastes gegenüber dem umliegenden Weichteilgewebe ein Kontrastmittel 3 mittels eines Kontrastmittel-Injektors 14 injiziert werden. Das Kontrastmittel 3 wird in automatisierter Form – vorgegeben durch einen Injektionsprotokoll – in der Regel zeitgesteuert aus einem Vorratsbehälter 15 über einen Kontrastmittelschlauch 16 in einer einstellbaren Menge und mit einer einstellbaren Flussgeschwindigkeit in eine Vene des Patienten 2 gepumpt. Durch eine elektrische Verbindung zwischen der Recheneinheit 12 und dem Kontrastmittel-Injektor 14 sind Steuerkommandos zur Kontrastmittelvergabe und Injektionsprotokolle übermittelbar.
  • Bei der Ausbreitung des Kontrastmittels 3 im Körperinneren durchläuft das eingebrachte Kontrastmittel 3 den Blutkreislauf des Patienten 2 und erreicht erst nach einer gewissen Zeit den Untersuchungsbereich 7. Um diese im Voraus unbekannte Verzögerung zwischen Start der Injektion und Start der Kontrastmittelanreicherung bestimmen zu können, kann vor dem Hauptscan (CTA) eine kleine Menge Kontrastmittel injiziert werden, der sogenannte Testbolus. Bei Erreichen des Untersuchungsbereichs 7 steigt die Konzentration des Kontrastmittels 3 des Testbolus im Blut zunächst an, durchläuft ein Maximum 18 und fällt anschließend wieder ab. Das zeitliche Verhalten der Konzentration ist in Form eines Kontrastmittelverlaufs 25 in der Zeit 20 darstellbar, der durch das Testbolussignal 19 im Wesentlichen repräsentiert wird. In der klinischen Routine werden auf der Grundlage dieses Testbolussignals 19 oder aufgrund einer daraus abgeleitete Kontrastmittel-Vorhersage geeignete Betriebsparameter des Computertomographiegerätes 1 ermittelt, so dass die Abtastung des Hauptscans zu einem Zeitpunkt und mit einer Geschwindigkeit erfolgt, zu dem die Konzentration des Kontrastmittels 3 in dem Untersuchungsbereich 7 möglichst hoch ist. Die Betriebsparameter umfassen beispielsweise den Startzeitpunkt der Abtastung und die Abtastgeschwindigkeit bzw. den Pitchwert, also das Verhältnis zwischen dem Vorschub der Tischplatte 5 pro Rotation der Gantry 8 und der Schichtdicke des Röntgendetektors 10.
  • Die Steuereinrichtung 12 ist in dem dargestellten Ausführungsbeispiel mit einer erfindungsgemäßen Vorrichtung 60 zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs CEP ausgestattet, die später noch näher erläutert wird. Zudem weist die Steuereinrichtung 12 eine Schnittstelle 70 zur Erfassung oder Ausgabe eines Kontrastmittelgabe-Startzeitpunkts und eine Referenzzeitpunkt-Ermittlungseinrichtung 80 auf, die dann unter Berücksichtigung des Kontrastmittelgabe-Startzeitpunkts sowie auf Basis des voraussichtlichen Kontrastmittelverlaufs CEP einen Referenzzeitpunkt für das Bildgebungssystem bestimmt, so dass beispielsweise automatisch kurz vor einem voraussichtlichen Zeitpunkt der maximalen Kontrastmittelanreicherung im untersuchten Bereich eine bildgebende Messung für diesen Bereich mit den passenden weiteren Parametern gestartet wird.
  • Ein typischer Kontrastmittelverlauf 25 eines Testbolus ist beispielhaft in 2 dargestellt. Entlang der x-Achse ist die Zeit t in Einheiten von Sekunden (s) aufgetragen. Die y-Achse entspricht den relativen Schwächungswerten in Einheiten von Hounsfield-Units (HU). Bei dem Beispiel wurden Schwächungswerte 19 alle zwei Sekunden in Abtastzeitpunkten 20 erfasst. Nicht sämtliche Schwächungswerte 19 und nicht sämtliche Abtastzeitpunkte 20 sind aus Gründen der Übersichtlichkeit mit einem Bezugszeichen versehen worden. Die Schwächungswerte 19 sind in dem Diagramm punktförmig eingetragen worden und können zur Berechnung einer patientenspezifischen Impulsantwortfunktion (AIR) dienen.
  • Zur Ermittlung des Kontrastmittelverlaufs 25 eines Testbolus im Zusammenhang mit einer Vormessung werden im Wesentlichen folgende Verfahrensschritte durchgeführt:
    • – Einstellen einer festen Abtastposition 21, vorzugsweise im Untersuchungsbereich (siehe 1),
    • – Durchführen einer Serie von Abtastungen zu vordefinierten Abtastzeitpunkten 20 zur Ermittlung des Kontrastmittelverlaufs, wobei zu jedem Abtastungszeitpunkt 20 ein Schwächungswert 19 erfasst wird, durch den eine Konzentration des Kontrastmittels 3 an der Abtastposition 21 repräsentiert wird, und
    • – Abspeichern der Schwächungswerte 19 des Kontrastmittelverlaufs und der Abtastzeitpunkte 20 zur Berechnung von Parametern und Vorhersagen für spätere Untersuchungen des Patienten 2.
  • In der 3 ist ein bekannter Rechenvorgang zur Vorhersage des Kontrastmittelverlaufs gemäß dem Stand der Technik dargestellt. Aus den Blutflussverlaufsdaten 35, im Folgenden auch als "Blutfluss-Signal" 35 (BCS(t)) (engl. blood circulation signal) bezeichnet, die beispielsweise in der Vormessung aus dem Testbolus TB oder optisch bestimmt werden, sowie einem definierten Untersuchungsprotokoll 47, nämlich dem bei der Vormessung verwendeten Injektionsprotokoll IFTB(t), wird mittels eines AIR-Algorithmus 36 die individuelle Patientenfunktion 37 (AIRTB(t)) berechnet. Mit einem PRED-Algorithmus 38 und dem Kontrastmittelinjektionsprotokoll 48 des Hauptscans IFCTA(t) wird eine Vorhersage 39 (CEP(t)) des Kontrastmittelverlaufs des Hauptscans berechnet. Die verschiedenen Berechnungen des CEP-Algorithmus 46, beispielsweise mit dem in der älteren Patentanmeldung DE 10 2012 209 410.5 beschriebenen Verfahren, erfolgen dabei immer in der Zeitdomäne und/oder der von der Zeitdomäne abgeleiteten Fourierdomäne.
  • Es sei an dieser Stelle erwähnt, dass es sich bei Blutflussverlaufsdaten BCS(t) (sofern diese mit einem Kontrastmittel-Testbolus gemessen wurden) und der Kontrastmittelverlaufsvorhersage CEP(t) im Prinzip um gleichartige Daten handelt, nämlich um eine Funktion der Schwächungswerte, welche die Kontrastmittelanreicherung an einem bestimmten Punkt innerhalb des Untersuchungsbereichs repräsentieren, in Abhängigkeit von der Zeit.
  • Bei allen bisherigen Ansätzen, aufgrund einer Testbolus-Messung eine Vorhersage über die Kontrastmittelanreicherung im Patienten zu machen, wurde die Herzrate des Patienten für konstant angenommen. Da aber der Blutkreislauf und damit die Strömungsgeschwindigkeit des Kontrastmittels von der Herzrate abhängen, ist die nachfolgend beschriebene Vorhersage, wobei Herzraten-Informationen berücksichtigt werden und in die Berechnungen mit einfließen, potentiell besser oder "sollte besser sein". Dies gilt vor allem bei Patienten mit zumindest über längeren Zeitabständen unregelmäßigen Herzraten.
  • Durch das Bestimmen der Herzrate mittels eines Pulsmessgeräts oder einer EKG-Vorrichtung 22 während einer Testbolus-Messung kann die individuelle Patientenfunktion (AIR – Arterial Impulse Response – patientenspezifische arterielle Kontrastmittel-Impulsantwortfunktion) zu einer Patientenfunktion umgestaltet werden, die direkt abhängig von der Herzschlaganzahl und somit quasi "unabhängig" von der Herzrate ist. Anders ausgedrückt, wird die Zeit-Achse zu einer Achse, die die Anzahl der Herzschläge wiedergibt. Graphisch veranschaulicht wird dies anhand eines Beispiels mit einer angenommenen regelmäßigen Herzrate von 90 Herzschlägen pro Minute (beats per minute (bpm)) in 4, in der im linken Diagramm die individuelle Patientenfunktion AIR(t) über der Zeit und im rechten Diagramm die individuelle Patientenfunktion AIR(HB) über der Herzrate aufgetragen ist.
  • Wichtig ist, dass diese Umstellung der Zeit-Achse prinzipiell nur dann möglich ist, wenn man die Herzrate für die Testbolus-Messung, bevorzugt während der Testbolus-Messung, gemessen hat. Man kann diese Umstellung der Zeit-Achse in erster Näherung mit der mittleren Herzrate während der Testbolus-Messung durchführen. Auch wäre es möglich, statt der mittleren Herzrate jeden einzelnen Herzschlag mit in Betracht zu ziehen, wenn z. B. ein EKG von Anfang der Testbolus-Injektion gemessen und registriert worden ist. Auch kann es z. B. implementationsbedingt vorkommen, dass nur die mittlere Herzrate in den CEP-Algorithmus mit einfließen kann, obwohl jeder einzelne Herzschlag während des Testbolus mittels EKG gemessen worden ist.
  • Vorzugsweise wird daher zeitgleich mit der Messung zur Ableitung einer AIR, d. h. der Testbolus-Messung, die Herzrate beispielsweise mit einem Pulsmessgerät oder einer EKG-Vorrichtung 22 gemessen (siehe 1). Eine solche Berücksichtigung der Herzraten-Informationen wird heutzutage nicht gemacht und ist für Testbolus-Scans an bekannten CT-Geräten derzeit auch nicht möglich. Mit diesen Informationen über die Herzrate kann die Patientenfunktion, welche aus dem Testbolus abgeschätzt worden ist, unabhängig von einer schwankenden oder wechselnden Herzrate gemacht werden, die zum Zeitpunkt des Testbolus vorhanden war. Anschließend kann man dann den Kontrastmittelverlauf zu einem späteren Zeitpunkt besser vorhersagen, falls sich die Herzrate geändert haben sollte.
  • Dies bringt folgende Vorteile mit sich:
    • 1. Eventuelle Änderungen in der Herzrate zwischen der Testbolus-Messung und dem CTA-Scan können in der Vorhersage des Kontrastmittelverlaufs im Patienten während des CTA-Scans berücksichtigt werden. Anhand dieser Vorhersage kann beispielsweise der Scan-Start optimiert werden.
    • 2. Insbesondere für Patienten mit einer unregelmäßigen Herzrate wird eine bessere Vorhersage des Kontrastmittelverlaufs erwartet.
    • 3. Eine Patientenfunktion, die unabhängig von der Herzrate ist, bietet neue Möglichkeiten, um den kardiovaskulären Status eines Patienten zu charakterisieren. Dieser neue Biomarker könnte vielleicht bei Patienten nützlich sein, die systemische Chemotherapie bekommen oder die an einer degenerativen kardiovaskulären Krankheit leiden.
    • 4. Da CTA-Scans des Herzens bereits EKG-getriggert akquiriert werden, verursacht das Messen der Herzrate während des Testbolus und der späteren Messung für diese Patienten keine extra Kosten und keinen extra Aufwand für den medizinisch-technischen Radiologieassistenten.
    • 5. Wenn eine herzschlagbezogene Impulsantwortfunktion AIR(HB) bestimmt wurde, kann diese auch abgespeichert und später nicht nur für einen relativ kurz Zeit später nachfolgenden Hauptscan, sondern für beliebige weitere Untersuchungen verwendet werden. Es kann dann auf einen weiteren Testbolus-Scan verzichtet werden, da die herzschlagbezogene Impulsantwortfunktion AIR(HB) ja immer an die aktuelle Herzrate angepasst werden kann.
    • 6. Durch die Messungen des Vormessungs-Herzratensignal und des aktuellen Herzratensignals kann sehr einfach auch die Ankunftszeit (Bolus arrival time) besser vorhergesagt werden, wann das Kontrastmittel nach dem Startpunkt der Injektion erstmals am Messort auftritt. Dies ist unabhängig von den weiteren Berechnungen bzw. Korrekturen des Kurvenverlaufs der Impulsantwortfunktion möglich, d. h die Impulsantwortfunktion wird hierdurch nur entlang der Zeitachse verschoben.
  • Die 5 dient zur Erläuterung der erfindungsgemäßen gleichzeitigen Akquirierung der Herzraten-Informationen und der Vormessung zur Ableitung der individuellen AIR und zeigt der Zusammenhang zwischen Ereignissen in der Zeitdomäne und in der Herzschlagdomäne. In der oberen Kurve der 5 ist das EKG eines Patienten über der Zeit t aufgetragen. Für jeden einzelnen Herzschlag HB ist bekannt, zu welchem Zeitpunkt t dieser stattgefunden hat. Weiterhin ist auf der unteren Kurve eine beispielhafte Kontrastmittel-Injektion CI[ml/s] eines Testbolus mit einer Gesamtmenge an Kontrastmittel von 15 ml bei einer Rate von 5 ml/s gezeigt (Bei dieser Kurve handelt es sich um ein typisches Injektionsprotokoll für einen Testbolus). Die Injektionszeit t[s] des Kontrastmittels von 3 Sek. ist äquivalent zu 4,3 Herzschlägen HB. Der CT-Scan beginnt nach 4 Sek., entsprechend 6,2 Herzschlägen HB. Analog zu diesen Beispielen kann man die Abtastzeitpunkte 20 der Testbolus-Messung 19 auch konvertieren von Sekunden nach Anzahl der Herzschläge HB[#]. Man erhält dann die Testbolus-Messung 19 als Funktion der Anzahl der Herzschläge HB[#]. Bei gleichzeitiger Aufzeichnung des Herzratensignals kann dieser Zusammenhang jederzeit im Nachhinein bestimmt werden.
  • Für eine korrekte Anwendung der Herzraten-Informationen in der Vorhersage sollte man nicht nur die Zeit-Informationen konvertieren, sondern auch ggf. implementationsbedingt die Informationen des Untersuchungsprotokolls, in diesem Beispiel das Injektionsprotokoll des Testbolus, ändern. Wenn man z. B. einen AIR-Algorithmus verwendet, der wie nachfolgend erläutert unmittelbar in der Herzschlagdomäne rechnen kann, kann zuvor das Injektionsprotokoll von der Zeitdomäne in die Herzschlagdomäne transferiert werden. Das Injektionsprotokoll IF (CI in der Figur) in [ml/s] kann dabei wie in 6 dargestellt in die Einheit [ml/HB] umgerechnet werden, wenn die Analyse mit jedem einzelnen Herzschlag erfolgt. Mit einer Kontrastmittel-Injektion CI[ml/s] eines Testbolus mit 15 ml Kontrastmittel bei einer Rate von 5 ml/s ergibt sich der gezeigte Verlauf der Kontrastmittel-Injektion CI[ml/HB]. Die jeweilige Flussrate in [ml/HB] ergibt sich einfach aus der Multiplikation des Zeitintervals zwischen zwei Herzschlägen [s] und der Flussrate in [ml/s].
  • 7 zeigt ein Flussdiagramm für ein erstes Beispiel eines erfindungsgemäßen Verfahrens, bei dem
    • – die Implementation des AIR-Algorithmus es erlaubt, die Ableitung der individuellen, patientenspezifischen Impulsantwortfunktion oder Patientenfunktion AIR in der Herzschlagdomäne (HB) statt in der Zeitdomäne (t) zu berechnen,
    • – das HR-Signal ein EKG-Signal ist, wobei jeder einzelne Herzschlag berücksichtigt werden kann, und
    • – die korrigierte Vorhersage des Kontrastmittelverlaufs während des Hauptscans auf einer kurz vor dem Hauptscan gemessenen mittleren Herzrate (HRCTA(t)) basiert.
  • Hierzu wird neben dem Blutfluss-Signal 35 (BCS) auch ein die Herzrate kennzeichnendes Herzratensignal 40 (HR-Signal(t)), in diesem Falle ein EKG, ermittelt und herangezogen, das, wie durch die gestrichelte Linie angedeutet – zeitlich korreliert mit dem Blutfluss-Signal 35 erfasst und aufgezeichnet wird.
  • Aus diesem Herzratensignal 40 wird jeder einzelne Herzschlag HBTB(t) während der Testbolusgabe als Funktion der Zeit bestimmt, wie dies in 5 durch den Verfahrensschritt 41 illustriert ist. Die Funktion HBTB(t) repräsentiert also auch die Anzahl der Herzschläge nach einer bestimmten Zeit t. Das Identifizieren 41 der einzelnen Herzschläge HBTB kann beispielsweise gemäß internationalen Standards erfolgen, wie sie z. B. von The CSE Working Party in "Recommendations for measurement standards in quantitative electrocardiography", European Heart Journal (1985), Vol. 6, Seiten 815 bis 825 und in "Common Standards for Quantitative Electrocardiography: Goals and Main Results" von J. L. Willems et al. for the CSE Working Party in Methods of Information in Medicine (1990) Bd. 29, Nr. 4, Seiten 263 bis 271, beschrieben sind.
  • Die so ermittelten Herzschlagdaten HBTB ermöglichen es, das Blutfluss-Signal 35 (BCS) und das Injektionsprotokoll 47 (IFTB(t)) in die Herzschlagdomäne zu konvertieren, wie dies bereits anhand der 5 und 6 gezeigt wurde. Das Blutfluss-Signal 35 (BCS) und das Injektionsprotokoll 47 (IFTB(t)) werden hierzu einem AIR-Algorithmus 42 zugeführt, der zunächst unter Berücksichtigung der ebenfalls zugeführten Herzschlagdaten 41 (HBTB) diese Umwandlungen durchführt und dann in der Herzschlag-Domäne aus diesen Werten unter Verwendung von Gleichung [2] die individuelle Patientenfunktion 43 (AIRTB(HB)) in Abhängigkeit vom Herzschlag berechnet (wobei die Funktion C(HB) durch die Blutflussverlaufsdaten BCS, umgerechnet in die Herzschlagdomäne, gegeben ist).
  • Mit der herzschlagbezogenen individuellen Patientenfunktion 43 (AIRTB(HB)) kann anschließend eine Vorhersage aufgrund einer kurz vor der eigentlichen Untersuchung ermittelten mittleren CTA-Herzrate 44 (HRCTA(t)) (in bpm) erfolgen.
  • Hierzu wird bei dem Ausführungsbeispiel gemäß 7 zunächst eine korrigierte Patientenfunktion 45 (AIRCORR(HB)) folgendermaßen berechnet: AIRCORR(t) = (HRCTA/60)·AIRTB(HB/HRCTA·60) [11].
  • Diese Gleichung gilt, wegen des Umrechnungsfaktors 60, wenn die Zeit wie üblich in Sekunden, aber die Herzrate in bpm gemessen wird. Da der Blutkreislauf des Patienten ein geschlossenes System ist, sollte bei einer Anpassung der AIR-Kurve in der Breite, die Amplitude der Patientenfunktion sich reziprok ändern, damit die Fläche unter der AIR-Kurve konstant bleibt. Dafür sorgt in Gleichung [11] der Vorfaktor (HRCTA/60). Im Grunde genommen handelt es sich bei Gleichung [11] um einen Spezialfall der Gleichung [6'a].
  • Mittels des PRED-Algorithmus 38 ergibt sich dann die Vorhersage 39 (CEP(t)) des Kontrastmittelverlaufs im Zeitbereich gemäß Gleichung [5].
  • 8 beschreibt ein zweites Beispiel eines erfindungsgemäßen Verfahrens, bei dem
    • – die Implementation des AIR-Algorithmus es nicht ermöglicht, die Ableitung der individuellen Impulsantwortfunktion AIR in der Herzschlagdomäne zu berechnen, wie es in 7 der Fall war,
    • – das HR-Signal entweder ein EKG-Signal ist, aus dem die mittlere Herzrate während des Testbolus HRTB(t) berechnet wird, oder die mittlere Herzrate HRTB(t) genommen wird, wie sie z. B. von einem Pulsmesser wiedergegeben wird, und
    • – die korrigierte Vorhersage des Kontrastmittelverlaufs während des Hauptscans auf einer kurz vor dem Hauptscan gemessenen mittleren Herzrate (HRCTA(t)) basiert.
  • In diesem Beispiel ist es implementationsbedingt unmöglich, die Berechnung der patientenspezifischen Impulsantwortfunktion AIR in der HB-Domäne auszuführen, weil beispielsweise der AIR-Algorithmus in der Fourier-Domäne rechnet und dafür eine homogene Abtastung der Signale erforderlich sein kann. Deshalb kann in diesem Beispiel nur die mittlere Herzrate HRTB(t) während des Testbolus mit in die Korrektur einfließen. Für Patienten, bei denen die Herzrate regelmäßig ist, aber zwischen dem Testbolus und dem Hauptscan Medikamente verabreicht worden sind, die die Herzrate dauerhaft gesenkt oder angehoben haben, wird dieses Ausführungsbeispiel trotzdem zu einer verbesserten Vorhersage führen, auch wenn nicht jeder einzelne Herzschlag berücksichtigt wird.
  • Erfindungsgemäß wird auch hier neben dem Blutfluss-Signal 35 (BCS(t)), vorzugsweise gleichzeitig, wieder ein Herzratensignal 40 ermittelt.
  • Die individuelle Patientenfunktion 37 (AIRTB(t)) eines CEP-Algorithmus 50 als Funktion der Zeit wird dann mit dem normalen, unveränderten, bereits vorhandenen AIR-Algorithmus 36 berechnet und nachfolgend einfach in eine korrigierte Patientenfunktion 45 (AIRCORR(t)) transformiert.
  • Hierzu könnte zum Beispiel zunächst gemäß den Gleichungen [3], [3'] oder [3''] in einem Zwischenschritt eine individuelle herzschlagbezogene Patientenfunktion AIRTB(HB) berechnet werden. Diese Umwandlung der abgeleiteten Patientenfunktion 37 AIRTB(t) als Funktion der Zeit in eine neue Patientenfunktion AIRTB(HB) als Funktion der Anzahl der Herzschläge ist graphisch anhand der 3 veranschaulicht. In diesem Falle sehen die Kurvenverläufe gleich aus, da die Herzrate über die Zeit konstant angenommen wird, es ist die mittlere Herzrate. Bei jedoch variierender Herzrate, und bei Berücksichtigung jedes einzelnen Herzschlags (7), würden die Kurvenverläufe unterschiedliche Formen aufweisen. Die Funktion der Anzahl der Herzschläge AIRTB(HB) ist aber unabhängig von der Herzrate.
  • Diese neue (zwischenzeitliche) Patientenfunktion AIRTB(HB) ermöglicht es, eine Änderung in der mittleren Herzrate 44 (HRCTA(t)) während des CTA-Scans zu berücksichtigen. Falls sich kurz vor dem CTA-Scan herausstellt, dass die mittlere Herzrate 44 (HRCTA(t)) signifikant anders ist als die mittlere Herzrate 51 während der Vormessung, dann kann diese Umrechnung z. B. gemäß Gleichung [6], [6'], [6a] oder [6'a] erfolgen, je nachdem, in welchen Einheiten gemessen wird und ob die korrigierte Funktion bereits den Vorfaktor erhalten soll oder der Vorfaktor gemäß Gleichung [9] später mit dem Injektionsprotokoll multipliziert wird.
  • Besonders bevorzugt kann aber auch ohne eine vorherige Umrechnung der zeitbezogenen Patientenfunktion AIRTB(t) in eine herzschlagbezogene Patientenfunktion AIRTB(HB) gemäß den Gleichungen [7] oder [7a] unmittelbar aus der zeitbezogenen Patientenfunktion AIRTB(t) eine zeitbezogene korrigierte Patientenfunktion AIRCORR(t) berechnet werden. Wie zu ersehen ist, ist es in diesem Ausführungsbeispiel nicht notwendig, einen vorhandenen CEP-Algorithmus weitreichend zu ändern, da die AIR- 36 und PRED-Algorithmen 38 unverändert bleiben können, und nur in einem neuen Zwischenschritt die individuelle Impulsantwortfunktion AIR korrigiert werden muss. Eine Berechnung in der Herzschlagdomäne ist in diesem Ausführungsbeispiel also nicht notwendig.
  • 9 zeigt an einem Beispiel eines 71-jährigen Patienten, bei dem während der Messung des Testbolus 30 die Herzrate 91 bpm betrug und bei dem sich nach Gabe eines Betablockers die Herzrate auf 74 bpm kurz vor dem CTA-Scan gesenkt hat, eine Anpassung der abgeleiteten Patientenfunktion AIRTB(t) 23 mit einem Verfahren nach 8, so dass man eine korrigierte Patientenfunktion 24 (AIRCORR(t)) erhält. Diese angepasste Patientenfunktion AIRCORR(t) führte bei diesem Patienten zu einer deutlich verbesserten Vorhersage des Kontrastmittelverlaufs, wie in 10 gezeigt wird.
  • Anhand der 10 wird dies nun näher erläutert. Dieses Ausführungsbeispiel bezieht sich auf die mittlere Herzrate während der Testbolus-Messung und kurz vor dem CTA-Scan. Die Punkte (o o o) zeigen die Kontrastmittel-Anreicherung 30 des gemessenen Testbolus (TB). Durch einen AIR-Algorithmus wurde eine Kurvenanpassung 31 ermittelt, die die Funktion IFTB(t) ⊗ AIRTB(t) aus der Gleichung [1] darstellt. Die gestrichelte Kurve gibt eine daraus abgeleitete unkorrigierte CTA Vorhersage 32 wieder. Die Punkte (
    Figure DE102013201136A1_0004
    ) zeigen die gemessene CTA Kontrastmittel-Anreicherung 33. Die durchgezogene Kurve gibt eine aus der Funktion AIRCORR(t) 24 abgeleitete korrigierte CTA Vorhersage 34 wieder.
  • Vergleicht man die beiden Vorhersagen, die unkorrigierte CTA Vorhersage 32 und die korrigierte CTA Vorhersage 34 mit der tatsächlich gemessenen CTA Kontrastmittel-Anreicherung 33, ist sofort ersichtlich, dass die korrigierte CTA Vorhersage 34 der gemessenen CTA Kontrastmittel-Anreicherung 33 näher kommt und somit besser bzw. genauer ist.
  • Die 11 beschreibt schließlich ein drittes Beispiel eines erfindungsgemäßen Verfahrens, bei dem
    • – die Implementation des AIR-Algorithmus es erlaubt, die Ableitung der individuellen Impulsantwortfunktion AIR in der Herzschlagdomäne (HB) statt in der Zeitdomäne (t) zu berechnen,
    • – das HR-Signal ein EKG-Signal ist, wobei jeder einzelne Herzschlag berücksichtigt werden kann, und
    • – die korrigierte Vorhersage des Kontrastmittelverlaufs während des Hauptscans auf einem während des Hauptscans gemessenen Herzratensignal (HBCTA) in Form eines EKG basiert.
  • Für Patienten mit einer sehr unregelmäßigen Herzrate wäre dieses Ausführungsbeispiel besonders vorteilhaft, da jeder einzelne Herzschlag, der während der Testbolus-Messung und während der späteren kontrastmittelunterstützten Messung gemessen worden ist, in die Berechnung mit einfließt.
  • Dieses Ausführungsbeispiel bezieht sich auf eine Aktualisierung der Vorhersage nach jedem einzelnen Herzschlag des aktuellen EKGs während der Hauptuntersuchung. Dabei erfolgt die Ableitung der individuellen Patientenfunktion AIRTB(HB) wie bei dem Ausführungsbeispiel gemäß 7. Alternativ könnte, sowohl hier als auch bei dem Ausführungsbeispiel nach 7, das Blutflusssignal BCS(HB) gleich in Abhängigkeit vom Herzschlag aufgezeichnet werden, und das Injektionsprotokoll IFTB(HB) könnte nicht zeitgesteuert, sondern herzschlaggesteuert sein. In diesem Fall erhält der AIR-Algorithmus 42 bereits alle Informationen in der Herzschlagdomäne.
  • Anders als bei dem Ausführungsbeispiel nach 7 erfolgt aber nun die Vorhersage des Kontrastmittelverlaufs kontinuierlich aktualisiert aufgrund der laufend während der kontrastmittelunterstützten Untersuchung gemessenen Herzrate bzw. dem aktuellen Herzschlag 52 (HBCTA(t)) des Patienten, welcher mithilfe des EKGs erfasst wird.
  • Hierzu muss lediglich das aktuelle Injektionsprotokoll IFCTA(t), sofern es zeitgesteuert ist, in ein herzschlagbezogenes Injektionsprotokoll 55 (IFCTA(HB)) umgerechnet werden. Außerdem müssen der CEP-Algorithmus 53 und damit der PRED-Algorithmus 54 die Vorhersage innerhalb der Herzschlagdomäne rechnen können, und es sollten genügend Computer-Ressourcen zur Verfügung stehen, um die Vorhersage in Echtzeit berechnen zu können.
  • Bei hier nicht dargestellten weiteren Alternativen wird gemäß den 7 und 8 gearbeitet, aber es wird permanent die aktuelle Herzrate während der kontrastmittelunterstützten Untersuchung aktualisiert (beispielsweise durch Neubildung eines Mittelwerts mit jedem weiteren Herzschlag) und dementsprechend die Kontrastmittelverlaufs-Vorhersage aktualisiert.
  • Um die verschiedenen Verfahrensvarianten durchzuführen, umfasst die in 1 innerhalb der Steuereinrichtung 12 dargestellte Vorrichtung 60 zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs CEP folgende Komponenten:
    • – eine Eingangs-Schnittstelle 61, mit der die patientenspezifischen Blutflussverlaufsdaten 35 erfasst werden,
    • – eine Impulsantwortfunktions-Ermittlungseinheit 62, welche so ausgebildet ist, dass sie eine individuellen Impulsantwortfunktion 37, 43 aus diesen Blutflussverlaufsdaten 35 bestimmen kann,
    • – eine Eingangs-Schnittstelle 63 zur Erfassung des Herzratensignals 40, 44, 52,
    • – eine Berechnungseinheit 64, welche ausgebildet ist, um z. B. nach einer der obigen Verfahrensvarianten den voraussichtlichen Kontrastmittelverlauf CEP vorherzusagen, und
    • – eine Ausgabeschnittstelle 65, um den voraussichtlichen Kontrastmittelverlauf CEP an andere Komponenten zu übergeben.
  • Falls eine herzschlagbezogene Impulsantwortfunktion AIR(HB) für eine spätere Verwendung gespeichert oder übergeben soll, könnte dies auch über die Ausgabeschnittstelle 65 erfolgen (die hierzu in geeigneter Weise ausgestaltet ist) oder über eine (in 1 nicht dargestellte) separate Ausgabeschnittstelle.
  • Voraussetzung für alle diese Ausführungsbeispiele ist immer die Erfassung eines geeigneten Herzratensignals (z. B. mittels eines EKGs) in Verbindung mit der Vormessung und kurz vor oder während der CT-Hauptuntersuchung. Dies ermöglicht es, die Vorhersage unabhängig von der Herzrate zu machen, die während der Vormessung vorhanden war. Schwankungen in der Herzrate sowie eine dauerhafte Erhöhung oder Senkung der Herzrate werden die Vorhersage dann nicht mehr beeinträchtigen.
  • Dabei gilt:
    • 1. Es gibt mehrere Möglichkeiten, die Herzrate aufzuzeichnen, beispielsweise über EKG, Pulsmesser oder Blutdruckmessgerät.
    • 2. Es gibt mehrere Möglichkeiten, die Messung zur Bestimmung der individuellen Impulsantwortfunktion AIR auszuführen, beispielsweise mit einem Testbolus oder optisch.
    • 3. Im Falle einer Testbolus-Messung gibt es mehrere bekannte Algorithmen, die eine patientenspezifische Impulsantwortfunktion aus dem Testbolus ableiten können, wie beispielsweise in der älteren Patentanmeldung DE 10 2012 209 410.5 beschrieben ist, oder eine Entfaltung im Fourierbereich oder eine Analyse mit einem pharmacokinetischen Modell.
    • 4. Abhängig vom CEP-, AIR- und/oder PRED-Algorithmus, der vorhanden ist, und von der Art des HR-Signals (vollständiges EKG oder nur Angaben zur mittleren Herzrate während der Akquisition), gibt es beliebig viele Ausführungsbeispiele, die Vorhersage für schwankende Herzraten zu korrigieren. Die Implementation der Herzraten-Korrektur sollte in jeden einzelnen Fall abgestimmt werden.
    • 5. Dieses erfindungsgemäße angiographische Untersuchungsverfahren beschränkt sich nicht nur auf Computertomographie, sondern kann potenziell auch für andere bildgebende Verfahren wie Magnetresonanz- oder Ultraschalltomographie angewendet werden.
  • ZITATE ENTHALTEN IN DER BESCHREIBUNG
  • Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
  • Zitierte Patentliteratur
    • DE 102012209410 [0002, 0071, 0107]
  • Zitierte Nicht-Patentliteratur
    • Dominik Fleischmann et al., "Mathematical Analysis of Arterial Enhancement and Optimization of Bolus Geometry for CT Angiography Using the Discrete Fourier Transform", Journal of Computer Assisted Tomography, 1999, Vol. 23, No. 3, Seiten 474 bis 484 [0003]
    • Andreas H. Mahnken et al, "Quantitative prediction of contrast enhancement from test bolus data in cardiac MSCT", Eur. Radiol. 2007, 17, Seiten 1310 bis 1319 [0003]
    • F. Eisa et al in "Optical tracking of contrast medium bolus to optimize bolus shape and timing in dynamic computed tomography", Physics in Medicine and Biology 2012, 57, Seiten N173 bis N182 [0012]
    • The CSE Working Party in "Recommendations for measurement standards in quantitative electrocardiography", European Heart Journal (1985), Vol. 6, Seiten 815 bis 825 [0082]
    • "Common Standards for Quantitative Electrocardiography: Goals and Main Results" von J. L. Willems et al. for the CSE Working Party in Methods of Information in Medicine (1990) Bd. 29, Nr. 4, Seiten 263 bis 271 [0082]

Claims (15)

  1. Verfahren zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs CEP für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts (2) mit folgenden Schritten: – Erfassung von patientenspezifischen Blutflussverlaufsdaten (35) in einer Vormessung unter Verwendung eines definierten Untersuchungsprotokolls (47) und Bestimmung einer individuellen Impulsantwortfunktion (37, 43) auf Basis der Blutflussverlaufsdaten (35), – Erfassung eines der Vormessung zugeordneten Vormessungs-Herzratensignals (40) und – Vorhersage (39) des voraussichtlichen Kontrastmittelverlaufs (CEP) auf Basis der individuellen Impulsantwortfunktion (37, 43), des Vormessungs-Herzratensignals (40) sowie eines aktuellen Untersuchungsprotokolls (48) und eines aktuell erfassten Herzratensignals (44, 52).
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass auf Basis des Vormessungs-Herzratensignals (40) eine herzschlagbezogene Impulsantwortfunktion (43) ermittelt wird, welche eine individuelle Impulsantwortfunktion in Abhängigkeit von einer Herzschlaganzahl definiert.
  3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass zunächst eine zeitbezogene individuelle Impulsantwortfunktion (37) ermittelt wird und dann unter Nutzung des Vormessungs-Herzratensignals (40) in eine herzschlagbezogene Impulsantwortfunktion (43) umgerechnet wird.
  4. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass auf Basis des aktuell erfassten Herzratensignals (44) und der herzschlagbezogenen Impulsantwortfunktion (43) oder der zeitbezogenen individuellen Impulsantwortfunktion (37) eine zeitbezogene korrigierte individuelle Impulsantwortfunktion (45) gebildet wird, und die Vorhersage (39) des voraussichtlichen Kontrastmittelverlaufs (CEP) durch Kombination der korrigierten Impulsantwortfunktion (45) mit dem aktuellen Untersuchungsprotokoll (48) erfolgt.
  5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass zur Bildung der korrigierten individuellen Impulsantwortfunktion (45) die individuelle herzschlagbezogene Impulsantwortfunktion (43) mit einem Zeitskalierungsfaktor skaliert wird, welcher durch den Kehrwert einer aus dem aktuell erfassten Herzratensignal (44) ermittelten Herzrate (HRCTA(t)) gegeben ist.
  6. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass zur Bildung der korrigierten individuellen Impulsantwortfunktion (45) die individuelle zeitbezogene Impulsantwortfunktion (37) mit einem Zeitskalierungsfaktor skaliert wird, welcher durch eine aus dem Vormessungs-Herzratensignal (37) ermittelten Herzrate (HRTB(t)) (HRTB(t)) im Verhältnis zu einer aus dem aktuell erfassten Herzratensignal (44) ermittelten Herzrate (HRCTA(t)) (HRCTA(t)) gegeben ist.
  7. Verfahren nach Anspruch 5 oder 6, dadurch gekennzeichnet, dass zur Bildung der korrigierten individuellen Impulsantwortfunktion (45) die skalierte individuelle Impulsantwortfunktion (37) mit dem Kehrwert des jeweiligen Zeitskalierungsfaktors multipliziert wird.
  8. Verfahren nach einem der Ansprüche 1 bis 7, dadurch gekennzeichnet, dass die Blutflussverlaufsdaten und das Vormessungs-Herzratensignal (40) zeitlich korreliert erfasst werden.
  9. Verfahren nach einem der Ansprüche 1 bis 8, dadurch gekennzeichnet, dass zur Erfassung des Vormessungs-Herzratensignals (40) eine Aufzeichnung von EKG-Daten erfolgt.
  10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, dass aus den EKG-Daten eine während der Vormessung vorliegende mittlere Herzrate (51) abgeleitet wird.
  11. Verfahren zur Ansteuerung eines medizintechnischen Bildgebungssystems (1), wobei ein Referenzzeitpunkt für das Bildgebungssystem (1) auf Basis eines Kontrastmittelgabe-Startzeitpunkts einer Kontrastmittelgabe und einem gemäß einem der Ansprüche 1 bis 10 vorhergesagten voraussichtlichen Kontrastmittelverlaufs (CEP) bestimmt wird.
  12. Vorrichtung (60) zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs (CEP) für eine kontrastmittelunterstützte Untersuchung eines Untersuchungsobjekts (2) mit – einer Eingangs-Schnittstelle (61) zur Erfassung von patientenspezifischen Blutflussverlaufsdaten (35) in einer Vormessung unter Verwendung eines definierten Untersuchungsprotokolls (47), – einer Impulsantwortfunktions-Ermittlungseinheit (62) zur Bestimmung einer individuellen Impulsantwortfunktion (37, 43) auf Basis der Blutflussverlaufsdaten (35), – einer Eingangs-Schnittstelle (63) zur Erfassung eines Herzratensignals (40, 44, 52), – einer Berechnungseinheit (64) zur Vorhersage (39) des voraussichtlichen Kontrastmittelverlaufs (CEP) auf Basis der individuellen Impulsantwortfunktion (37, 43), eines der Vormessung zugeordneten Vormessungs-Herzratensignals sowie eines aktuellen Untersuchungsprotokolls (48) und eines aktuell erfassten Herzratensignals (44, 52).
  13. Steuereinrichtung (12) für ein Bildgebungssystem (1) mit – einer Schnittstelle (70) zur Erfassung oder Ausgabe eines Kontrastmittelgabe-Startzeitpunkts, – einer Vorrichtung (60) zur Vorhersage eines voraussichtlichen Kontrastmittelverlaufs (CEP) nach Anspruch 11, – und einer Referenzzeitpunkt-Ermittlungseinrichtung (71) zur Ermittlung eines Referenzzeitpunkts für das Bildgebungssystem (1) auf Basis des voraussichtlichen Kontrastmittelverlaufs (CEP).
  14. Bildgebungssystem (1) zum Erzeugen von Bilddatensätzen vom Inneren eines Untersuchungsobjekts (2) umfassend eine Steuereinrichtung (12) nach Anspruch 13.
  15. Computerprogrammprodukt, welches direkt in einen Speicher einer Steuereinrichtung (12) eines Bildgebungssystems (1) ladbar ist, mit Programmcode-Mitteln, um alle Schritte eines Verfahrens nach einem der Ansprüche 1 bis 11 auszuführen, wenn das Programmprodukt auf der Steuereinrichtung (12) ausgeführt wird.
DE102013201136.9A 2013-01-24 2013-01-24 Vorhersage eines voraussichtlichen Kontrastmittelverlaufs Active DE102013201136B4 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102013201136.9A DE102013201136B4 (de) 2013-01-24 2013-01-24 Vorhersage eines voraussichtlichen Kontrastmittelverlaufs
US14/149,919 US9271656B2 (en) 2013-01-24 2014-01-08 Prediction of a likely contrast medium behavior
CN201410022486.9A CN103961117B (zh) 2013-01-24 2014-01-17 对预计的造影剂变化的预测

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102013201136.9A DE102013201136B4 (de) 2013-01-24 2013-01-24 Vorhersage eines voraussichtlichen Kontrastmittelverlaufs

Publications (2)

Publication Number Publication Date
DE102013201136A1 true DE102013201136A1 (de) 2014-07-24
DE102013201136B4 DE102013201136B4 (de) 2023-01-19

Family

ID=51064505

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102013201136.9A Active DE102013201136B4 (de) 2013-01-24 2013-01-24 Vorhersage eines voraussichtlichen Kontrastmittelverlaufs

Country Status (3)

Country Link
US (1) US9271656B2 (de)
CN (1) CN103961117B (de)
DE (1) DE102013201136B4 (de)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2990073B1 (de) 2004-11-24 2018-05-02 Bayer Healthcare LLC Vorrichtungen und systeme zur flüssigkeitsabgabe
US9302044B2 (en) 2006-12-29 2016-04-05 Bayer Healthcare Llc Patient-based parameter generation systems for medical injection procedures
AU2011270772C1 (en) 2010-06-24 2017-04-20 Bayer Healthcare Llc Modeling of pharmaceutical propagation and parameter generation for injection protocols
CN104428660A (zh) 2012-05-14 2015-03-18 拜耳医疗保健公司 用于基于x射线管电压确定药用流体注射协议的系统和方法
DE102012209410B4 (de) * 2012-06-04 2020-08-13 Siemens Healthcare Gmbh Ermittlung einer patientenspezifischen Kontrastmittel-Impulsantwortfunktion sowie darauf basierende Vorhersage eines voraussichtlichen Kontrastmittelverlaufs und Ansteuerung eines medizintechnischen Bildgebungssystems
DE102015105978B3 (de) * 2015-04-20 2016-09-15 Carl Mahr Holding Gmbh Haltevorrichtung für eine optische Messeinrichtung
EP3093678A1 (de) * 2015-05-13 2016-11-16 Bayer Pharma Aktiengesellschaft Verfahren zur optimierung der vorbestimmung des zeitlichen verlaufes einer kontrastmittelkonzentration bei der diagnostischen bildgebung mit einem magnetresonanzsystem
US10383590B2 (en) * 2015-09-28 2019-08-20 General Electric Company Methods and systems for adaptive scan control
JP6884367B2 (ja) * 2016-11-17 2021-06-09 株式会社根本杏林堂 注入プロトコルの生成装置、該生成装置を備える注入装置及び撮像システム、注入プロトコル生成方法、及び注入プロトコル生成プログラム
EP3570060A1 (de) * 2018-05-15 2019-11-20 Siemens Healthcare GmbH Verfahren zur rekonstruktion von magnetresonanztomographie-aufnahmen mit variabler zeitauflösung
CN110464326B (zh) 2019-08-19 2022-05-10 上海联影医疗科技股份有限公司 一种扫描参数推荐方法、系统、装置及存储介质
CN110852993B (zh) * 2019-10-12 2024-03-08 拜耳股份有限公司 一种造影剂作用下的成像方法与设备
CN110916704B (zh) * 2019-11-13 2023-09-01 上海联影医疗科技股份有限公司 校正方法、装置及存储介质
EP4362042A1 (de) 2022-10-31 2024-05-01 Siemens Healthineers AG Schätzung der pharmakokinetik eines kontrastmittels durch ein analytisches modell eines kardiovaskulären systems und/oder einen algorithmus eines neuronalen netzwerks

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060178836A1 (en) * 2004-08-31 2006-08-10 Er-Wei Bai System and method for adaptive bolus chasing computed tomography (CT) angiography
DE102005006657A1 (de) * 2005-02-14 2006-08-31 Siemens Ag Verfahren zur Vorhersage des Kontrastmittelflusses in einem lebenden Körper
DE102006055167A1 (de) * 2005-12-23 2007-07-05 Siemens Ag Betriebsverfahren für eine Röntgeneinrichtung zur Untersuchung eines Patienten unter Einsatz eines Kontrastmittels
US20080097197A1 (en) * 2004-11-16 2008-04-24 Kalafut John F Modeling of Pharmaceutical Propagation
US20100030073A1 (en) * 2006-12-29 2010-02-04 Medrad, Inc. Modeling of pharmaceutical propagation
US20100204572A1 (en) * 2007-07-17 2010-08-12 Medrad, Inc. Devices, Systems and Methods for Determination of Parameters for a Procedure, for Estimation of Cardiopulmonary Function and for Fluid Delivery
DE102012209410A1 (de) 2012-06-04 2013-12-05 Siemens Aktiengesellschaft Ermittlung einer patientenspezifischen Kontrastmittel-Impulsantwortfunktion

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060178836A1 (en) * 2004-08-31 2006-08-10 Er-Wei Bai System and method for adaptive bolus chasing computed tomography (CT) angiography
US20080097197A1 (en) * 2004-11-16 2008-04-24 Kalafut John F Modeling of Pharmaceutical Propagation
DE102005006657A1 (de) * 2005-02-14 2006-08-31 Siemens Ag Verfahren zur Vorhersage des Kontrastmittelflusses in einem lebenden Körper
DE102006055167A1 (de) * 2005-12-23 2007-07-05 Siemens Ag Betriebsverfahren für eine Röntgeneinrichtung zur Untersuchung eines Patienten unter Einsatz eines Kontrastmittels
US20100030073A1 (en) * 2006-12-29 2010-02-04 Medrad, Inc. Modeling of pharmaceutical propagation
US20100204572A1 (en) * 2007-07-17 2010-08-12 Medrad, Inc. Devices, Systems and Methods for Determination of Parameters for a Procedure, for Estimation of Cardiopulmonary Function and for Fluid Delivery
DE102012209410A1 (de) 2012-06-04 2013-12-05 Siemens Aktiengesellschaft Ermittlung einer patientenspezifischen Kontrastmittel-Impulsantwortfunktion

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
"Common Standards for Quantitative Electrocardiography: Goals and Main Results" von J. L. Willems et al. for the CSE Working Party in Methods of Information in Medicine (1990) Bd. 29, Nr. 4, Seiten 263 bis 271
Andreas H. Mahnken et al, "Quantitative prediction of contrast enhancement from test bolus data in cardiac MSCT", Eur. Radiol. 2007, 17, Seiten 1310 bis 1319
Dominik Fleischmann et al., "Mathematical Analysis of Arterial Enhancement and Optimization of Bolus Geometry for CT Angiography Using the Discrete Fourier Transform", Journal of Computer Assisted Tomography, 1999, Vol. 23, No. 3, Seiten 474 bis 484
F. Eisa et al in "Optical tracking of contrast medium bolus to optimize bolus shape and timing in dynamic computed tomography", Physics in Medicine and Biology 2012, 57, Seiten N173 bis N182
The CSE Working Party in "Recommendations for measurement standards in quantitative electrocardiography", European Heart Journal (1985), Vol. 6, Seiten 815 bis 825

Also Published As

Publication number Publication date
CN103961117B (zh) 2016-08-17
US20140206991A1 (en) 2014-07-24
CN103961117A (zh) 2014-08-06
DE102013201136B4 (de) 2023-01-19
US9271656B2 (en) 2016-03-01

Similar Documents

Publication Publication Date Title
DE102013201136B4 (de) Vorhersage eines voraussichtlichen Kontrastmittelverlaufs
DE102010027227B4 (de) Verfahren und Computertomographiegerät zur Durchführung einer angiographischen Untersuchung
DE102012209410B4 (de) Ermittlung einer patientenspezifischen Kontrastmittel-Impulsantwortfunktion sowie darauf basierende Vorhersage eines voraussichtlichen Kontrastmittelverlaufs und Ansteuerung eines medizintechnischen Bildgebungssystems
DE102012205351A1 (de) Darstellung von Blutgefäßen und Gewebe im Herzen
WO2014195077A1 (de) Verfahren und system zur ermittlung eines mess-startzeitpunktes
DE202017106016U1 (de) Medizinische Informationsverarbeitungsvorrichtung, Röntgen-CT-Vorrichtung und computerlesbares Speichermedium mit einem Programm für ein medizinisches Informationsverarbeitungsverfahren
EP3332710B1 (de) Charakterisierung von plaque
DE102011078278A1 (de) Verfahren zur Bilderzeugung und Bildauswertung
WO2013098074A2 (de) Kontrollverfahren und kontrollsystem
DE102015222853A1 (de) Verfahren zum automatischen Ermitteln eines Kontrastmittel-Injektionsprotokolls
DE102016215971A1 (de) Segmentierung einer Angiographie unter Verwendung einer bestehenden dreidimensionalen Rekonstruktion
DE102016219887A1 (de) Verfahren und System zur Nutzung von Messdaten
DE102009054990A1 (de) Messprotokoll-Optimierung für ein medizintechnisches Gerät
DE102008016891A1 (de) Betriebsverfahren für eine verschwenkbare Polyplan-Bildgebungsanlage zur zeitaufgelösten Abbildung eines Untersuchungsobjekts
DE102019201079A1 (de) Verfahren und bildgebendes Gerät zum Erzeugen eines bewegungskompensierten Bildes, Computerprogramm und Speichermedium
DE102008026035A1 (de) Betriebsverfahren für eine verschwenkbare Polyplan-Bildgebungsanlage zur Abbildung eines sich bewegenden Untersuchungsobjekts
DE102012205711B4 (de) Verfahren zum Betreiben eines bildgebenden Diagnosegerätes sowie medizinisches bildgebendes System
DE102011083629A1 (de) Bestimmung von potentiellen Perfusionsdefekten
DE102008010006A1 (de) Verfahren zur dreidimensionalen Darstellung einer bewegten Struktur durch ein tomographisches Verfahren
DE102008021835A1 (de) Verfahren und Tomographiegerät zur Normierung von Bilddaten hinsichtlich eines durch ein Kontrastmittel in den Bilddaten hervorgerufenen Kontrastes
DE102009043633A1 (de) Verbesserte Abtastung eines zyklisch bewegten Untersuchungsobjektes unter Einsatz eines Kontrastmittels im Rahmen einer Voruntersuchung mittels eines CT-Gerätes
DE102021208272A1 (de) Optimale Gewichtung von DSA-Maskenbildern
DE102006055167A1 (de) Betriebsverfahren für eine Röntgeneinrichtung zur Untersuchung eines Patienten unter Einsatz eines Kontrastmittels
DE102012215294A1 (de) Verfahren zum Erzeugen eines Volumenmodells von einem Herz und zugehörige C-Bogen-Röntgenanlage
DE102013202491A1 (de) Dosisreduzierte CT-Aufnahme mittels dynamischer Kollimierung

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R084 Declaration of willingness to licence

Effective date: 20140724

R079 Amendment of ipc main class

Free format text: PREVIOUS MAIN CLASS: A61B0019000000

Ipc: A61B0034000000

R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHCARE GMBH, DE

Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHINEERS AG, DE

Free format text: FORMER OWNER: SIEMENS HEALTHCARE GMBH, MUENCHEN, DE