CN103961117A - 对预计的造影剂变化的预测 - Google Patents

对预计的造影剂变化的预测 Download PDF

Info

Publication number
CN103961117A
CN103961117A CN201410022486.9A CN201410022486A CN103961117A CN 103961117 A CN103961117 A CN 103961117A CN 201410022486 A CN201410022486 A CN 201410022486A CN 103961117 A CN103961117 A CN 103961117A
Authority
CN
China
Prior art keywords
heart rate
contrast agent
impulse response
response function
independent
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
CN201410022486.9A
Other languages
English (en)
Other versions
CN103961117B (zh
Inventor
J.G.科波拉尔
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
Original Assignee
Siemens AG
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 filed Critical Siemens AG
Publication of CN103961117A publication Critical patent/CN103961117A/zh
Application granted granted Critical
Publication of CN103961117B publication Critical patent/CN103961117B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026Measuring blood flow
    • A61B5/0275Measuring blood flow using tracers, e.g. dye dilution
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/488Diagnostic techniques involving pre-scan acquisition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT 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/60ICT 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/63ICT 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504Apparatus 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices 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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10ICT 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/17ICT 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)
  • Epidemiology (AREA)
  • Primary Health Care (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)
  • General Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Business, Economics & Management (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明涉及一种对预计的造影剂变化(CEP)进行预测的方法,用于对检查对象(2)进行造影剂辅助的检查,具有以下步骤:在使用定义的检查协议(47)的条件下在预测量中采集特定于患者的血流变化数据(35),并且从该血流变化数据中确定单独的脉冲响应函数(37,43);采集与预测量对应的预测量-心率信号(40);并且基于单独的脉冲响应函数、预测量-心率信号以及当前的检查协议(48)和当前采集的心率信号(44,52)预测(39)预计的造影剂变化。本发明还涉及一种用于控制医学技术成像系统(1)的方法。另外描述了一种对预计的造影剂变化进行预测的设备(60)、一种具有这样的设备的控制装置(12)以及一种包括这样的控制装置的成像系统。

Description

对预计的造影剂变化的预测
技术领域
本发明涉及对预计的造影剂变化进行预测的一种方法和一种设备,用于对检查对象进行造影剂辅助的检查,特别是造影剂辅助的成像方法,优选是在对患者进行血管造影的检查方法的范围内。本发明还涉及一种用于控制医学技术成像系统的方法、一种用于成像系统的控制装置以及一种具有这种控制装置的成像系统。
背景技术
在早先的专利申请DE 10 2012 209 410.5中描述了一种根据测试推注数据(Testbolus-Daten)来确定单独的、特定于患者的造影剂脉冲响应函数(以下也简称为“脉冲响应函数”或“患者函数”)的方法,之后利用该函数可以在造影剂辅助的成像测量的情况下产生对造影剂变化的预测。该患者函数描述了在测量测试推注的时间点上的患者的心血管特征。原则上,在假设单独的患者函数在之后的时间点还有效的条件下,可以由此预测任意注射协议的造影剂变化。在此,注射协议是精确的时间预计,按照其给出了在患者上的造影剂量,并且例如包括给予造影剂的开始时间点、流速和终止时间点。
通常这类所谓的CEP算法(英语contrast enhancement prediction,对比度增强预测)在不同的出版物中公知,例如:
·Dominik Fleischmann等,“Mathematical Analysis of ArterialEnhancement and Optimization of Bolus Geometry for CT Angiography Usingthe Discrete Fourier Transform”,Journal of Computer Assisted Tomography,1999,Vol.23,No.3,第474至484页,和
·Andreas H.Mahnken等,“Quantitative prediction of contrastenhancement from test bolus data in cardiac MSCT”,Eur.Radiol.2007,17,第1310至1319页。
这种CEP算法计算患者内的造影剂的预计的变化,并且通常由两个子算法构成。利用第一子算法确定特定于患者的造影剂脉冲响应函数(英语arterial impulse response,AIR,动脉脉冲响应),即所谓的“AIR算法”。第二子算法计算对患者内的预计的造影剂变化进行预测(英语prediction),即所谓的“PRED算法”。因此,CEP算法通常包括AIR算法与PRED算法的组合。
在所有这些已知的方法中,待检测的患者被考虑为线性时不变系统(英语LTI-system)。这通常意味着,可以将时间上测量的患者内的造影剂积聚C(t)数学上描述为注射协议IF(t)与脉冲响应函数或患者函数AIR(t)的卷积(英语convolution):
C ( t ) = IF ( t ) ⊗ AIR ( t ) - - - [ 1 ] .
因此,CEP算法在时域和/或由时域导出的傅里叶域中执行全部的计算。
但是,这假定造影剂的驱动力(血液循环)也是不变的,并且是关于时间稳定的。但通常不是这种情况。
造影剂辅助的成像方法的典型应用领域例如是CT血管造影(CTA),即,借助计算机断层造影进行血管成像。在此,CTA扫描是用于血管成像的计算机断层造影的运行过程,而心脏CTA或心脏扫描是心脏成像以及显示心脏血管。
特别是在CT血管造影(CTA)的这种心脏扫描的情况下可能出现,例如在β受体阻滞剂(例如美托洛尔,Metoprolol)的影响下,心率(单位时间内的心跳次数)在测试推注后持续上升或下降,典型是在患者具有过高心率的情况下在扫描前的短期内给予该β受体阻滞剂。因为以其它的心率改变了血液的再循环时间,所以在这种情况下也改变了患者函数。这导致对造影剂变化的不精确预测。
发明内容
本发明要解决的技术问题是,构造对预计的造影剂变化进行预测的一种方法和一种设备,用于开头所述类型的造影剂辅助的检查,特别是造影剂辅助的成像方法,优选是在血管造影的检查方法的范围内,使得主要在不规则心率的情况下,特别是心率在执行测试推注与真正的之后的检查之间发生变化的情况下,能够对之后的时间点上的造影剂变化进行更好的预测。
在根据本发明的方法中,为了对预计的造影剂变化进行预测,用于对检查对象(以下在不限制普遍性的情况下称为“患者”,所以也使用术语“特定于患者”代替术语“特定于检查对象”)进行造影剂辅助的检查,执行以下步骤:
-首先,在使用定义的检查协议的条件下,在预测量中采集特定于患者的血流变化数据,并且从该血流变化数据中确定患者的单独的脉冲响应函数。在此,在成像测量(例如之前的CT测量,特别是CTA)的范围内,在使用合适的造影剂的条件下,可以如同在之后的检查中那样获得该血流变化数据。在此,可以是“真实的”先前的造影剂辅助的检查。但是,通常是只利用少量的造影剂(测试推注)执行的特殊的测试测量。但是,也可以按照其它方式确定血流变化数据,例如光学地利用测试推注按照用于对血流变化进行颜色标记的措施的方式来确定,如同F.Eisa等在“Optical tracking ofcontrast medium bolus to optimize bolus shape and timing in dynamic computedtomography”,Physics in Medicine and Biology2012,57,第N173至N182页中描述的那样。在此,检查协议例如包括用于给予造影剂或测试推注的定义的注射协议。
-另外,采集与预测量对应的预测量-心率信号。在此,应当将表述“与预测量对应的”理解为,预测量-心率信号以充分的方式,根据现有的技术可能性和本发明的精确的实施方案,表示在预测量的时间点或在预测量期间存在的患者的心率。例如可以在预测量范围内,即,在预测量时或者短暂之前或之后,采集预测量-心率信号。
-然后,基于单独的脉冲响应函数、预测量-心率信号(直接或间接通过将预测量-心率信号已经引入到单独的脉冲响应函数)、以及当前的检查协议和当前采集的心率信号,进行对预计的造影剂变化的预测。在此,当前采集的心率信号表示在之后的测量的时间点的患者的心率或心跳,对于之后的测量应当进行对造影剂变化的预测。因此,例如可以在检查前或检查时的足够短的间隔内进行采集。在此明确指出,造影剂变化预测也包括以下情况,在以后的检查期间(例如CTA),从心跳到心跳地按照当前的“在线预览”方式,预测造影剂如何在血液循环中扩散。这还将在以后详细地描述。
在此,对预测量-心率信号和心率信号的采集既可以涉及每个单独的心跳(HB)(如果例如借助EKG测量并记录了心率),也可以涉及有关平均的心率(HR)的信息。后者例如可以是脉搏测量计的情况,其将心率反映为每分钟的心跳次数(英语beats per minute,bpm)。
利用在预测量时的心率信息或心跳信息,可以这样构成患者函数,使得该患者函数明确地考虑在预测量的时间点的心率。随后,在心率发生变化的情况下,可以更好地预测在之后的时间点的造影剂变化。因此,通过采集和记录或存储预测量-心率信号,并通过采集当前的心率信号,通常得到显著改善的、针对非恒定心率的预计的造影剂变化的预测。
优选又可以将这种方法用于根据本发明的用于自动控制医学技术成像系统的方法,其中,基于给予造影剂的造影剂给予-开始时间点及预测的造影剂变化确定成像系统的参考时间点(例如用于开始在特定区域内的图像采集的开始时间点)。例如根据预测的造影剂变化的局部最大值的位置,可以得到开始给予造影剂与开始通过成像系统进行扫描之间的延迟时间。
然而确定这样的参考时间点可能对于其它方法也是有意义的,例如对于以后评估造影剂拍摄的图像数据组。因此,在许多情况下主要需要知道,特定的图像是在哪个相位上拍摄的,例如其是否是这样的相位,在该相位上,造影剂才分散或积聚在特定的结构中,或者图像是否是在这样的相位上产生的,在该相位上,造影剂又从结构中消失(即分解)。
一种根据本发明的自动对预计的造影剂变化进行预测的设备,用于对检查对象进行造影剂辅助的检查,至少具有以下部件:
-输入接口,用于在使用定义的检查协议的条件下,在预测量中采集特定于患者的血流变化数据。在此,例如可以是指向存储器的接口,在该存储器中储存了所涉及的数据,但是也可以是指向评估单元的接口,在该评估单元中,可以对当前的、刚从测试推注测量中得出的测量数据进行评估,并且可以从该评估单元中直接接收数据。
-脉冲响应函数确定单元,用于基于血流变化数据确定单独的脉冲响应函数。
-输入接口,用于采集心率信号。这里,例如也可以是指向存储器的接口,或者是指向用于测量或记录心率信号的测量设备的接口,例如指向EKG的接口。如果将数据或信号共同存储,原则上也可以将用于采集特定于患者的血流变化数据的输入接口与用于采集心率信号的输入接口构造成共同的接口。
-计算单元,用于基于单独的脉冲响应函数、与预测量对应的预测量-心率信号(其中,直接或间接通过将预测量-心率信号已经引入到单独的脉冲响应函数)、以及当前的检查协议和当前采集的心率信号,对预计的造影剂变化进行预测。
于是通过根据本发明的设备的输出接口,可以例如向其它单元输出造影剂变化预测或表示待预计的造影剂变化的相应的函数或表格,该其它单元利用造影剂变化预测进一步工作,特别是执行进一步的计算或将其存储。
特别优选地将这种用于对预计的造影剂变化进行预测的设备构造成用于成像系统的控制装置的组件。即,控制装置本身具有该用于自动对预计的造影剂变化进行预测的设备。在这种情况下,控制装置优选也包括:接口,用于采集或输出造影剂给予-开始时间点;以及参考时间点确定装置,用于基于造影剂给予-开始时间点和预计的造影剂变化确定用于成像系统的参考时间点。这具有以下优点,控制装置能够全自动地直接基于预测量确定造影剂变化预测,并且之后相应地借助控制装置这样控制成像系统,使得在实际测量时,能够在给予造影剂之后的正确时间点开始图像数据采集,和/或能够将图像数据在时间上与参考时间点相关联并存储。
原则上,也可以在其它计算机上实现根据本发明的用于自动对预计的造影剂变化进行预测的设备,该计算机例如从其它的部件接收所需的血流变化数据,并且例如通过网络与成像系统连接。应当将计算密集型的任务从控制装置中移出,以便将控制装置的计算容量全部提供用于控制其它测量,这是特别有意义的。
一种根据本发明的用于产生检查对象内部的图像数据的成像系统,有利地包括之前描述的控制装置,该控制装置又具有根据本发明的用于自动对预计的造影剂变化进行预测的设备。
特别地,根据本发明的用于自动对预计的造影剂变化进行预测的设备的脉冲响应函数确定单元和计算单元可以分别作为软件模块,在具有相应存储能力的合适的计算机上得以实现。如果只需要从其它(特别是在同一个计算单元上实现的)另外的单元或存储器接收血流变化数据或者向其它(特别是在同一个计算单元上实现的)另外的单元或存储器输出造影剂变化预测,则输入接口以及输出接口同样可以按照纯软件的形式得以实现。原则上,这些接口也可以作为组合的硬件软件接口得以实现,以便实现外部的输入和输出,例如借助由软件部件特别配置的硬件接口。大规模通过软件的实现具有以下优点,例如也能够以简单的方式通过软件升级来修正至今已经使用的成像系统的控制装置,以便按照根据本发明的方式进行工作。因此,也通过能直接加载到可编程的控制装置的存储器中的计算机程序产品来解决该技术问题,该计算机程序产品具有程序区段,以便在控制装置中执行该程序时,执行根据本发明的方法的全部步骤。
从属权利要求分别包括了本发明的特别优选的扩展和实施,其中,一种类型的权利要求也可以相应于另一种类型的从属权利要求进行扩展。
根据心率信息,即,根据采集或记录心率信号和/或实施CEP算法的类型,得到一个或多个可能性,来优化对不稳定的心率的造影剂变化的预测。
在优选的方案中,基于预测量-心率信号确定心跳相关的脉冲响应函数。该单独的“心跳相关的脉冲响应函数”定义了取决于心跳次数的单独的脉冲响应函数。换言之,在以上解释中,在时域中考虑迄今通常的单独的脉冲响应函数(以下为了区分有时也明确称为“时间相关的脉冲响应函数”),即,关于时间t的函数,而将“心跳相关的”脉冲响应函数在“心率域”或“心跳域”中考虑,即,它是关于心跳次数的函数,以下将其命名为缩写HB(heartbeat)。
为此,可以分别在心跳域中已经确定了特定于患者的血流变化数据和检查协议。例如,可以分别相同地根据并行确定的心跳来记录数据。也可以根据心跳次数来预先给定或者心跳控制地处理检查协议。同样也可以首先在时域中采集全部数据,然后换算到心跳域。
类似于公式[1],在心跳域中也成立:
C ( HB ) = IF ( HB ) ⊗ AIR ( HB ) - - - [ 2 ]
其中,这里(HB)如所述象征了心跳的依赖关系。因此,在特定于患者的血流变化数据C(HB)(即,取决于心跳的在患者内的造影剂积聚)与注射协议IF(HB)已知的情况下,能够数学上通过去卷积(Entfaltung)来确定心跳相关的脉冲响应函数或患者函数AIR(HB)。
因为,在时域或频域中已经为去卷积提供了大量的标准程序,关于总消耗更有利的可以是,首先以通常的方式确定时间相关的单独的脉冲响应函数AIRTB(t),并且在使用预测量-心率信号的条件下,将其换算成心跳相关的脉冲响应函数AIRTB(HB)。这可以例如借助以下公式:
AIRTB(HB)=AIRTB(t·HRTB) [3]
在此,HRTB是对应预测量确定的心率。在此并在以后,通常使用下标“TB”(Testbolus,测试推注),用于说明其是从预测量中得到的数据,而与为此是否使用了测试推注无关。这意味着,AIRTB是寻找的单独的脉冲响应函数,其基于预测量的数据而确定。
如上面已经提到的,可以不同地计算或确定心率。它经常是平均心率,其由在定义的时间段中的心跳次数得出。在此,可以从测量心率起定为关于全部迄今为止的心跳的时间段,并且将测量的心跳的次数除以该时间段。同样也可以构造一类浮动的平均值,其涉及特定过去的时间段或特定过去的心跳次数(例如最后十次心跳)。例如在利用简单的脉搏测量器的通常的脉搏测量中,隐含地执行了这样的浮动的平均值构造。在极端情况下,当前的心率也可以只以最后的心跳为基础,通过简单地使用最后的与次最后的心跳之间的时间间隔的倒数作为心率。
如果为确定心率HRTB没有使用与为确定时间t而使用的相同时间单位,则在公式[3]中必须附加考虑换算因子U:
AIR TB ( HB ) = AIR TB ( t · HR TB U ) - - - [ 3 ′ ]
如果如经常那样通常以单位秒(s)测量时间t,而以每分钟的心跳(bpm,beats per minute)测量心率,则对于换算因子U=60成立,即,成立以下公式:
AIR TB ( HB ) = AIR TB ( t · HR TB 60 ) - - - [ 3 ′ ′ ]
如果已知单独的心跳相关的脉冲响应函数AIRTB(HB),并且在以后的用于实际检查的造影剂辅助的测量中,检测协议,特别是注射协议IFCTA同样在心跳域中被控制或被换算到心跳域(一种为此可行的方法将在以后说明),则也可以无需其它,利用以下公式(参见公式[1])直接通过通常的卷积运算确定寻找的造影剂变化预测CEP(HB):
C EP ( HB ) = IF CTA ( HB ) ⊗ AIR TB ( HB ) - - - [ 4 ]
在此,函数CEP(HB)最终定义了心跳相关的造影剂变化预测,即,在患者内的特定位置上依据心跳的预计的造影剂积聚。下标“EP”代表对积聚的预测(Enhancement Prediction,增强预测)。另外,通常为了说明而使用下标“CTA”,用于说明其是以后的造影剂辅助的测量的数据,为该造影剂辅助的测量预测造影剂变化,更确切地说与它在此是否真正是CTA检查或者是其它的造影剂辅助的检查无关。
在以后的造影剂辅助的测量期间以一种方式采集和处理心率信号,使得观察每个心跳,以便有效地使用预测,使用公式[4]是特别有利的。
但是,如果在以后的造影剂辅助的测量中不在心跳域中工作,则提供了以下根据本发明的方法的替换的优选变换方案。在此,例如基于(在主扫描(即实际测量)的短暂之前或期间的)当前采集的心率信号和心跳相关的脉冲响应函数,构造时间相关的校正的单独的脉冲响应函数。替换地,也可以直接从时间相关的单独的脉冲响应函数中确定这样的时间相关的校正的单独的脉冲响应函数。在两种情况下,于是通过结合校正的时间相关的脉冲响应函数AIRCORR(t)与当前的(同样是以通常方式关于时间的)检查协议或注射协议IFCTA(t),按照
C EP ( t ) = IF CTA ( t ) ⊗ AIR TB ( t ) - - - [ 5 ]
对在时域中的预计的造影剂变化CEP(t)进行预测。
在这种变换方案中,为检查患者,在使用造影剂的条件下,实现了特定于患者的动脉的用于非恒定心率的造影剂脉冲响应函数的校正,以用于更好地对预计的造影剂变化进行预测。在此,可以按照有利的方式,在变换方案中将代表了心跳信号的数据和血流变化数据输送给CEP算法,该CEP算法从这些值中计算单独的脉冲响应函数。于是,从该单独的脉冲响应函数和在实际检查短暂之前确定的CTA心率中导出校正的单独的脉冲响应函数。或者可以从首先常规确定的时间相关的脉冲响应函数中,借助代表了在预测量中的心率信号的数据和代表了实际的造影剂辅助的测量短暂之前确定的心率信号的数据,作为待执行的检查所需要的参数,来计算校正的单独的脉冲响应函数。
如果首先确定了单独的心跳相关的脉冲响应函数AIRTB(HB),则通过利用时间定标因子将单独的心跳相关的脉冲响应函数定标,由此可以计算校正的单独的脉冲响应函数,通过从当前采集的心率信号中确定的心率HRCTA的倒数给出了该时间定标因子:
AIRCORR(t)=AIRTB(HB/HRCTA) [6]
如果以不同的时间单位确定了时间t本身和心率,又可以使用换算因子U:
AIRCORR(t)=AIRTB(HB/(HRCTA/U))=AIRTB(HB/HRCTA·U) [6′]
其中,在通常的单位(一侧是s,另一侧是bpm)下,成立U=60,如同以上结合公式[3]已经解释的那样。
如果相反地如通常那样只在时域中确定了单独的脉冲响应函数AIRTB(t),由此可以计算校正的单独的脉冲响应函数,其方式是,利用时间定标因子将单独的时间相关的脉冲响应函数定标,该时间定标因子通过从预测量-心率信号中确定的心率HRTB与从当前采集的心率信号中确定的心率HRCTA的比例来给出:
AIRCORR(t)=AIRTB(t·HRTB/HRCTA) [7]
在全部公式[6]至[7]中,如上所述,心率也可以是平均心率,或者是只针对一个单独的心跳的心率。
但是,在以上公式[6]至[7]中还未考虑,心率的变化最终也导致,在关于时域的不变的注射协议的情况下,每个心跳给出不同的造影剂剂量。换言之,心率的变化导致在心跳域中的注射协议的变化。为了考虑这些,并且更进一步改善预测,可以在以上公式[6]至[7]中分别添加补偿这些的校正因子。因此,优选为了构造校正的单独的脉冲响应函数,将定标的单独的脉冲响应函数与各自的时间定标因子的倒数相乘。对于公式[6]得到:
AIRCORR(t)=(HRCTA)·AIRTB(HB/HRCTA) [6a]
并且对于公式[6′]和[7]相应的:
AIRCORR(t)=(HRCTAU)·AIRTB(HB/HRCTA·U) [6′a]
AIRCORR(t)=(HRCTA/HRTB)·AIRTB(t·HRTB/HRCTA) [7a]
直观观察,校正因子(HRCTA)、(HRCTA/U)或(HRCTA/HRTB)或者在所属公式中的前置因子分别负责,在宽度上匹配AIR曲线时,幅度与单独的脉冲响应函数AIR是互逆变化的,以便保持AIR曲线下的面积恒定。
如果例如将公式[7a]代入到公式[5]中,则基于通常的时间相关的单独的脉冲响应函数AIRTB(t)和在预测量中的心率HRTB以及在以后的测量中的当前的心率HRCTA,如下得到对预计的造影剂变化CEP(t)的预测:
C EP ( t ) = IF CTA ( t ) ⊗ ( HR CTA / HR TB ) · AIR TB ( t · HR TB / HR CTA ) - - - [ 8 ]
对于其它的方法,按照公式[6]和[6′]类似地进行计算。
替换地,可以将考虑了检查协议的心率依赖关系的校正因子(HRCTA/HRTB)也与检查协议或注射协议IFCTA(t)相乘:
C EP ( t ) = ( HR CTA / HR TB ) · IF CTA ( t ) ⊗ AIR TB ( t · HR TB / HR CTA ) - - - [ 9 ]
因为对于卷积运算通常成立:
a · ( f ⊗ g ) = ( a · f ) ⊗ g = f ⊗ ( a · g ) - - - [ 10 ]
所以这导致相同的结果。在此,公式[9]给出了利用公式[7]的计算变换的例子。这对于利用公式[6]和[6′]的变换同样也成立。
如上已经提到的,可以按照不同的方式采集预测量-心率信号。只需要的是,以充足的方式代表在预测量期间的心率。为了允许尽可能精确的校正,优选在时间上校正血流变化数据和预测量-心率信号,即,时间上并行地采集和记录。特别优选地,在此通过并行地记录在预测量期间的EKG数据来进行对预测量-心率信号的采集。
在此,实施决定地,如果例如在以后的测量中不存在EKG,而总是只能利用平均心率工作,也可以从EKG数据中导出在预测量期间存在的平均心率。于是,这简化了计算,例如通过按照以上公式[7]或[7a]进行单独的脉冲响应函数的校正,或者按照公式[8]或[9]进行造影剂变化预测的校正。
附图说明
以下结合附图中示出的实施例详细解释本发明。附图中:
图1以透视示意图示出了根据本发明的成像系统,其适合用于实施根据本发明的血管造影的检查方法,以便在使用造影剂的条件下检查患者;
图2示出了在预检查期间确定的造影剂曲线;
图3示出了已知的用于按照现有技术在时域中对造影剂变化进行预测的方法;
图4示出了用于说明根据本发明的换算的图表,该换算将作为时间函数导出的患者函数AIRTB(t)换算成心跳次数的函数AIRTB(HB);
图5示出了用于解释在时域和心跳域中的事件之间的关系的图表;
图6示出了具有测试推注和EKG数据的注射协议以及换算成心率的注射协议的图表;
图7示出了根据本发明的用于预测造影剂变化的方法的第一实施方式的图示;
图8示出了根据本发明的用于预测造影剂变化的方法的第二实施方式的图示;
图9示出了具有第二实施方式的例子的图表,用于将导出的患者函数AIRTB(t)与变化的心率相匹配;
图10示出了用于解释通过根据本发明的校正按照第二实施方式改善预测的图表;
图11示出了根据本发明的用于在心率域中预测造影剂变化的方法的第三实施方式的图示。
具体实施方式
在以下的例子中,为简单起见并且不对普遍情况进行限定的条件下假定,预测量是测试推注测量,并且实际的造影剂辅助的检查是CTA。
在图1中示出了医学技术成像系统(或诊断装置),在该例子中是计算机断层造影设备1,其适合用于实施根据本发明的检查方法,以便在使用造影剂3的条件下检查患者2。代替所描述的计算机断层造影设备1,也可以例如使用磁共振断层造影设备或者超声波断层造影设备。
计算机断层造影设备1配备了具有可移动的台面5的患者安置台4,患者2可以被安置在该台面上。在旋转轴6的方向上可以调节台面5,从而使与患者2相连的检查区域7可以移动穿过计算机断层造影设备1的机架8的开口,进入拍摄系统9、10的测量区域。以这种方式可以在旋转轴6的方向上相对地互相调节患者2和拍摄系统9、10,从而可以采用不同的采样位置。检查区域7例如可以是患者的心脏。
为了采集投影,拍摄系统9、10具有以X射线管形式的X射线辐射器9和相对布置的X射线探测器10,其中,X射线探测器10被构造成弧形,并且包括多个排成探测器行的探测器元件11。X射线辐射器9产生以扇形的X射线束形式的辐射,该X射线束穿透测量区域,并且随后抵达X射线探测器10的探测器元件11。探测器元件11产生与穿过测量区域出现的X射线辐射的衰减有关的衰减值。例如分别借助与闪烁探测器光学耦合的光电二极管或者借助直接转换的半导体进行X射线辐射到衰减值的转变。X射线探测器11以这种方式产生一组衰减值,其也被称为投影。
拍摄系统9、10被可旋转地布置在机架8上,从而可以从不同的投影方向采集投影。根据设定的计算机断层造影设备1的运行方式,在固定设定的或可变的投影方向,在同时固定设定的或可变的采样位置的情况下,进行采样。通过在旋转轴6的方向上同时连续地进给患者2的情况下旋转机架8,例如从大量不同的投影方向中在沿着旋转轴6或沿着患者2的不同的位置处采集投影。以这种方式通过螺旋形采样(旋转拍摄系统9、10并同时进给台面5)获得的拍摄系统9、10的投影被传输到成像系统1的计算单元或控制装置12,并且被重建成可在显示器13上显示的图像数据。图像数据例如可以包括一个或多个检查区域7的层图像或体积图像。
为了检查血液流经的器官,例如心脏、肝脏或血管,可以借助造影剂注射器14给患者2注射造影剂3,以用于提高相对于周围的软组织部分的可见的对比度。造影剂3按照自动的形式(通过注射协议预先给定)通常时间控制地从储存容器15中通过造影剂软管16,以可调节的剂量并以可调节的流速被泵入患者2的静脉。通过计算单元12与造影剂注射器14之间的电连接,可以传输用于造影剂给予和注射协议的控制命令。
在造影剂3在身体内部扩散的情况下,注入的造影剂3流过患者2的血液循环,并且在一定时间后才到达检查区域7。为了能够确定预先未知的、在开始注射与开始造影剂积聚之间的延迟,可以在主扫描(CTA)之前注射少量的造影剂,即所谓的测试推注。在到达检查区域7时,测试推注的造影剂3的浓度在血液中首先增加,经过最大值18并随后再次下降。浓度的时间特性可以按照在时间20内的造影剂变化25的形式来表示,通过测试推注信号19基本上代表了该造影剂变化。按照临床惯例,在测试推注信号19的基础上,或者根据由此导出的造影剂预测,确定计算机断层造影设备1的合适的运行参数,从而在一个时间点并以一个速度进行主扫描的采样,对此,在检查区域7中的造影剂3的浓度尽可能高。运行参数例如包括采样的开始时间点和采样速度或节距值,即机架8每次旋转时台面5的进给与X射线探测器10的层厚之间的比例。
在所示的实施例中,控制装置12在所示的实施例中配备有根据本发明的用于对预计的造影剂变化CEP进行预测的设备60,以后还会对该设备进行详细解释。另外,控制装置12具有用于采集或输出造影剂给予-开始时间点的接口70和参考时间点确定装置80,于是在考虑造影剂给予-开始时间点以及基于预计的造影剂变化CEP的条件下,该参考时间点确定装置确定了用于成像系统的参考时间点,从而例如自动地在预计的检查区域内的最大造影剂积聚的时间点的短暂之前,利用合适的另外的参数开始对该区域的成像测量。
在图2中示例性示出了测试推注的典型的造影剂变化25。沿着X轴,以单位秒(s)标注了时间t。y轴对应于以单位豪斯菲尔德单元(HU)的相对的衰减值。在该例子中,每两秒钟在采样时间点20采集了衰减值19。出于清楚的原因,没有给所有的衰减值19并且没有给所有的采样时间点20设置附图标记。在图表中以点的形式标注了衰减值19,并且可以用于计算特定于患者的脉冲响应函数(AIR)。
为了结合预测量确定测试推注的造影剂变化25,主要执行以下方法步骤:
-设定固定的采样位置21,优选是在检查区域内(参见图1);
-执行一系列对预先定义的采样时间点20的采样,以用于确定造影剂变化,其中,对每个采样时间点20采集一个衰减值19,通过该衰减值表示在该采样位置21处的造影剂3的浓度;并且
-存储造影剂变化的衰减值19和采样时间点20,以用于计算参数并预测之后对患者2的检查。
在图3中示出了按照现有技术对造影剂变化进行预测的已知的计算过程。根据血流变化数据35,以下也称为“血流信号”35(BCS(t))(英语bloodcirculation signal),其例如在预测量中根据测试推注TB或光学地确定,以及根据定义检查协议47,即在预测量时使用的注射协议IFTB(t),借助AIR算法36计算单独的患者函数37(AIRTB(t))。利用PRED算法38和主扫描IFCTA(t)的造影剂注射协议48计算对主扫描的造影剂变化的预测39(CEP(t))。在此,例如利用在过去的专利申请DE102012209410.5中描述的方法,总是在时域和/或由时域导出的傅里叶域中进行CEP算法46的不同计算。
此处需要指出的是,血流数据BCS(t)(如果利用造影剂测试推注对其测量)和造影剂变化预测CEP(t)理论上是同类的数据,即衰减值关于时间的函数,该衰减值代表了检查区域内的特定的点处的造影剂积聚。
在迄今为止全部的方案中,根据测试推注测量进行关于患者内的造影剂积聚的预测,将患者的心率假定为恒定的。但是,因为血液循环和由此的造影剂的流动速度与心率相关,所以以下描述的预测可能是更好的,或者“应当是更好的”,其中,考虑了心率信息并且将其加入到了计算中。这尤其对于患者具有至少关于较长时间间隔的不规则的心跳下有效。
通过借助脉搏测量设备或EKG设备22在测试推注测量期间确定心率,可以将单独的患者函数(AIR,Arterial Impulse Resonse,动脉脉冲响应,特定于患者的动脉造影剂脉冲响应函数)转换成直接与心跳次数相关的、并因此某种程度上与心率“无关”的患者函数。换言之,时间轴变为了反映心跳次数的轴。这一点结合在图4中假定每分钟90次心跳(beats per minute(bpm))的规则的心率的例子图形地说明,其中,在左边图表中标出了关于时间的单独的患者函数AIR(t),并且在右边图表中标出了关于心率的单独的患者函数AIR(HB)。
重要的是,原则上只有在测量了对于测试推注测量的心率的情况下,优选在测试推注测量期间测量心率,才可以进行时间轴的这种转换。如果例如从测试推注注射开始时测量了EKG,并且进行了记录,在第一近似下,可以利用在测试推注测量期间的平均心率执行时间轴的转换。还可以代替平均心率而考虑每个单独的心跳。也可能出现例如由实施决定地,只能将平均心率代入CEP算法,尽管借助EKG测量了测试推注期间的每个单独的心跳。
因此,优选与用于导出AIR的测量(即,测试推注测量)时间同步地,例如利用脉搏测量设备或EKG设备22测量心率(参见图1)。目前不做这样地心率信息的考虑,并且现在为测试推注扫描而在已知的CT设备上考虑心率信息也是不可行的。利用这些关于心率的信息,可以使从测试推注中估计的患者函数独立于在测试推注的时间点存在的波动的或交变的心率。因此,如果心率应当发生改变,则可以更好地预测在以后的时间点的造影剂变化。
这带来以下优点:
1.在CTA扫描期间,在对患者内的造影剂变化的预测中,考虑测试推注测量与CTA扫描之间在心率上可能的改变。根据该预测例如可以优化扫描开始。
2.特别是对具有不规律心率的患者预计更好的造影剂变化的预测。
3.与心率无关的患者函数提供了新的可能性,用于表明患者的心血管的状态的特征。这个新的生物学标志也许可能对接受系统化疗或忍受退化的心血管疾病的患者有益。
4.因为已经EKG触发地采集了心脏的CTA扫描,在测试推注期间和在以后的测量期间的心率的测量不对患者造成额外的花费,并且不对医学技术的放射学工作人员造成额外的开支。
5.如果确定了心跳相关的脉冲响应函数AIR(HB),可以也将其存储,并且以后可以不仅用于相对短时间之后跟随的主扫描,而且用于任意的其它检查。于是可以免去另外的测试推注扫描,因为可以将心跳相关的脉冲响应函数AIR(HB)总是匹配于当前的心率。
6.通过测量预测量-心率信号和当前的心率信号,也可以非常简单地更好地预测到达时间(Bolus arrival time,推注到达时间),即造影剂在注射的开始时间点之后何时第一次出现在测量位置。这可以独立于其它的计算或者脉冲响应函数的曲线变化的校正,即,脉冲响应函数由此仅沿着时间轴延迟。
图5用于解释根据本发明的对心率信息的同时采集和用于导出单独的AIR的预测量,并且示出了在时域和在心跳域中的事件之间的关系。图5上面的曲线标注了关于时间t的患者的EKG。对于每个单独的心跳HB,已知其发生对应哪个时间点t。另外,在下面的曲线中示出了示例性的测试推注的造影剂注射CI[ml/s],具有在速率5ml/s下的造影剂总量15ml(该曲线是典型的用于测试推注的注射协议)。3秒的造影剂的注射时间t[s]等价于4.3次心跳HB。类似于该例子,也可以将测试推注测量19的采样时间点20从秒转换成心跳的次数HB[#]。于是得到作为心跳的次数HB[#]的函数的测试推注测量19。在同时记录了心率信号的情况下,可以在事后任意时间确定该关系。
对于心率信息在预测中的正确应用,应当不仅转换时间信息,而且也必要时由实施决定地改变检查协议的信息(在该例子中是测试推注的注射协议)。如果例如使用如以下解释那样可以直接在心跳域中计算的AIR算法,则可以首先将注射协议从时域传递到心跳域。在此,如果以每个单独的心跳进行分析,则可以将按照单位[ml/s]的注射协议IF(图中的CI)如图6所示那样换算到单位[ml/HB]。利用测试推注的造影剂注射CI[ml/s],其具有在速率5ml/s下的15ml造影剂,得到所示出的造影剂注射的变化CI[ml/HB]。简单地由两次心跳之间的时间间隔[s]与按照单位[ml/s]的流动速率相乘得到按照单位[ml/HB]的各自的流动速率。
图7示出了根据本发明的方法的第一例子的流程图,其中:
-AIR算法的实施允许,代替时域(t),而在心跳域(HB)中计算单独的、特定于患者的脉冲响应函数或者患者函数AIR的推导;
-HR信号是EKG信号,其中,可以考虑每个单独的心跳;并且
-对主扫描期间的造影剂变化的校正的预测基于主扫描短暂之前测量的平均心率(HRCTA(t))。
为此,除了血流信号35(BCS)外,还确定并引入表示心率的心率信号40(HR信号(t)),在该情况是EKG,如虚线显示地,与血流信号35时间相关地采集并记录该心率信号。
从该心率信号40中确定给予测试推注期间的每个单独的心跳HBTB(t),作为时间的函数,如在图5中通过方法步骤41图示的那样。因此,函数HBTB(t)也代表了按照一定时间t的心跳的次数。例如可以按照国际标准进行对单独的心跳HBTB的识别41,如其例如由The CSE Working Party在“Recommendations for measurement standards in quantitativeelectrocardiography”,European Heart Journal(1985),Vol.6,第815至825页和由J.L.Willems等在“Common Standards for QuantitativeElectrocardiography:Goals and Main Results”,the CSE Working Party inMethods of Information in Medicine(1990)Bd.29,Nr.4,第263到271页描述。
这样确定的心跳数据HBTB能够将血流信号35(BCS)和注射协议47(IFTB(t))转换到心跳域,如这一点已经结合图5和图6所示那样。为此,将血流信号35(BCS)和注射协议47(IFTB(t))传输到AIR算法42,该算法首先在考虑了同样被传输到的心跳数据41(HBTB)的条件下执行该转换,并且然后在心跳域中根据该值在使用公式[2]的条件下,计算关于心跳的单独的患者函数43(AIRTB(HB))(其中,通过换算到心跳域的血流变化数据BCS给出函数C(HB)。
利用心跳相关的单独的患者函数43(AIRTB(HB)),可以随后根据实际检查短暂之前确定的平均CTA心率44(HRCTA(t))(按照单位bpm)进行预测。
为此,在按照图7的实施例中,首先根据以下公式计算校正的患者函数45(AIRCORR(HB)):
AIRCORR(t)=(HRCTA/60)·AIRTB(HB/HRCTA·60) [11]。
由于换算因子60,该公式在以下情况下成立,即,如通常那样以秒测量时间,但是以bpm测量心率。因为患者的血液循环是封闭的系统,在匹配AIR曲线时,应当在患者函数的宽度、幅度上是互逆变化的,以便保持AIR曲线下的面积恒定。在公式[11]中由前置因子(HRCTA/60)对此负责。实质上,公式[11]是公式[6′a]的特殊情况。
然后,借助PRED算法38,按照公式[5]得到在时间范围内的造影剂变化的预测39(CEP(t))。
图8描述了根据本发明的方法的第二例子,其中:
-AIR算法的实施不可以如同在图7中的情况那样,在心跳域中计算单独的脉冲响应函数AIR的推导;
-HR信号是EKG信号,根据该EKG信号计算在测试推注HRTB(t)期间的平均心率,或者取平均心率HRTB(t),如其例如由脉搏测量器反映的那样;并且
-对主扫描期间的造影剂变化的校正的预测基于主扫描短暂之前测量的平均心率(HRCTA(t))。
在该例子中,由实施决定地,不能在HB域中实施特定于患者的脉冲响应函数AIR的计算,因为例如AIR算法在傅里叶域中进行计算,而为此可能需要对信号进行均匀的采样。因此在该例子中,只能够将测试推注期间的平均心率HRTB(t)代入校正。对于心率规则、但是在测试推注与主扫描之间给予了药物的患者,该药物使心率持续地下降或上升,即使没有考虑每个单独的心跳,该实施例仍然导致预测改善。
根据本发明,这里除了血流信号35(BCS(t))外,(优选同时地)还确定心率信号40。
然后,利用常规的、不改变的、已经存在的AIR算法36计算作为时间的函数的CEP算法50的单独的患者函数37(AIRTB(t)),并且随后简单地将其转换到校正的患者函数45(AIRCORR(t))。
为此,例如可以首先按照公式[3]、[3′]或[3′′]在中间步骤中计算单独的心跳相关的患者函数AIRTB(HB)。结合图3图解了从作为时间的函数的导出的患者函数37AIRTB(t)到作为心跳次数的函数的新的患者函数AIRTB(HB)的转换。在该情况下,曲线变化看起来相同,因为假定关于时间的心率恒定,它就是平均心率。但是在变化的心率的情况下,并且在考虑了每次单独的心跳的情况下(图7),曲线变化将具有不同的形式。但是,心跳次数的函数AIRTB(HB)与心率无关。
新的(中间时间的)患者函数AIRTB(HB)能够考虑在CTA扫描期间在平均心率44(HRCTA(t))中的变化。如果在CTA扫描短暂之前证实,平均心率44(HRCTA(t))明显不同于在预测量期间的平均心率51,则根据以哪个单位测量和校正的函数是否应当已经包含前置因子,或者该前置因子按照公式[9]随后与注射协议相乘,可以例如按照公式[6]、[6′]、[6a]或[6′a]进行换算。
特别优选地,但是也可以不进行之前的按照公式[7]或[7a]的从时间相关的患者函数AIRTB(t)到心跳相关的患者函数AIRTB(HB)的换算,直接根据时间相关的患者函数AIRTB(t)计算时间相关的校正的患者函数AIRCORR(t)。如可见那样,在该实施例中无需大规模地改变现有的CEP算法,因为AIR算法36和PRED算法38可以保持不变,并且仅须在新的中间步骤中校正单独的脉冲响应函数AIR。因此,在该实施例中不需要在心跳域中的计算。
图9以71岁患者为例示出了利用按照图8的方法匹配导出的患者函数AIRTB(t)23,其中,在测试推注30的测量期间,心率是91bpm,而在给予β受体阻滞剂之后,心率在CTA扫描短暂之前下降到74bpm,从而获得校正的患者函数24(AIRCORR(t))。对于该患者,匹配后的患者函数AIRCORR(t)导致明显改善的造影剂变化的预测,如图10中所示那样。
现在结合图10对其进行解释。该实施例涉及了测试推注测量期间的平均心率和CTA扫描短暂之前的平均心率。点(o o o)示出了测量的测试推注(TB)的造影剂积聚30。通过AIR算法确定曲线匹配31,其表示根据公式[1]的函数虚线的曲线描绘了由此导出的未校正的CTA预测32。点(◇◇◇)示出了测量的CTA造影剂积聚33。实线的曲线描绘了由函数AIRCORR(t)24导出的校正的CTA预测34。
如果将这两个预测、即未校正的CTA预测32和校正的CTA预测34与实际测量的CTA造影剂积聚33进行比较,可以马上看出,校正的CTA预测34更接近于测量的CTA造影剂积聚33,并且因此更好或更精确。
最后,图11描述了根据本发明的方法的第三例子,其中:
-AIR算法的实施允许,代替时域(t),而在心跳域(HB)中计算单独的脉冲响应函数AIR的推导;
-HR信号是EKG信号,其中,可以考虑每个单独的心跳;并且
-对主扫描期间的造影剂变化的校正的预测基于主扫描期间测量的以EKG形式的心率信号(HBCTA)。
该实施例特别有利于具有非常不规律的心率的患者,因为在计算中引入了在测试推注测量期间和在以后的造影剂辅助的测量期间测量了的每个单独的心跳。
该实施例涉及了在主检查期间在当前的EKG的每个单独的心跳之后对预测进行更新。在此,如同按照图7的实施例那样,进行单独的患者函数AIRTB(HB)的推导。替换地,既可以在这里也可以在按照图7的实施例中,同时根据心跳来记录血流信号BCS(HB),并且注射协议IFTB(HB)可以不是由时间控制的,而是由心跳控制的。在这种情况下,AIR算法42已经获得了心跳域里的全部信息。
与按照图7的实施例不同地,现在根据持续地在造影剂辅助的检查期间测量的患者的心率或借助EKG采集的当前的心跳52(HBCTA(t)),对造影剂变化的预测进行连续地更新。
为此,仅须将当前的注射协议IFCTA(t)(如果它是由时间控制的)换算成心跳相关的注射协议55(IFCTA(HB))。另外,CEP算法53和PRED算法54必须能够在心跳域中对预测进行计算,并且应当提供足够的计算机资源,以便能够实时地对预测进行计算。
在这里未示出的另外的替代方案中,按照图7和图8来工作,但是,它在造影剂辅助的检查期间永久性地更新当前的心率(例如通过利用每个新增的心跳重新构造平均值),并且相应地更新造影剂变化预测。
为了执行不同的方法方案,在图1中的控制装置12内示出的设备60包括用于对预计的造影剂变化CEP进行预测的以下部件:
-输入接口61,利用其采集特定于患者的血流变化数据35;
-脉冲响应函数确定单元62,将其构造成可以根据血流变化数据35确定单独的脉冲响应函数37、43;
-输入接口63,用于采集心率信号40、44、52;
-计算单元64,将其构造成用于例如按照以上方法方案中的任一个,对预计的造影剂变化CEP进行预测;以及
-输出接口65,用于将预计的造影剂变化CEP传输到其它部件。
如果为了以后的应用而应当存储或传输心跳相关的脉冲响应函数AIR(HB),则这也可以通过输出接口65进行(为此以合适的方式布置该输出接口),或者通过(在图1中未示出的)单独的输出接口进行。
全部这些实施例的前提条件是,总是结合预测量和在CT主检查短暂之前或期间采集合适的心率信号(例如借助EKG)。这能够使得预测与在预测量期间存在的心率无关。于是,心率的波动以及心率的持续上升或下降不再妨碍预测。
在此成立:
1.存在记录心率的多种可能性,例如通过EKG、脉搏测量器或血压测量设备。
2.存在实施用于确定单独的脉冲响应函数AIR的测量的多种可能性,例如利用测试推注或光学的。
3.在测试推注测量的情况下,存在多种已知的算法,这些算法可以从测试推注中导出特定于患者的脉冲响应函数,如同例如在过去的专利申请DE102012209410.5中描述的那样,或者在傅里叶范围内的去卷积或利用药物动力学模型的分析。
4.根据存在的CEP算法、AIR算法和/或PRED算法,以及根据HR信号的类型(完整的EKG或者只说明了获取期间的平均心率),存在任意多种实施例,其校正了对波动的心率的预测。应当在每个单独的情况中调谐心率校正的实施。
5.根据本发明的血管造影的检查方法不仅在计算机断层造影上应用,而且可以潜在性地也应用于其它的成像方法,诸如磁共振断层造影或超声波断层造影。

Claims (15)

1.一种对预计的造影剂变化(CEP)进行预测的方法,用于对检查对象(2)进行造影剂辅助的检查,包括以下步骤:
-在使用定义的检查协议(47)的条件下,在预测量中采集特定于患者的血流变化数据(35),并且基于该血流变化数据(35)确定单独的脉冲响应函数(37,43);
-采集与所述预测量对应的预测量-心率信号(40);并且
-基于单独的脉冲响应函数(37,43)、预测量-心率信号(40)以及当前的检查协议(48)和当前采集的心率信号(44,52),进行对预计的造影剂变化(CEP)的预测(39)。
2.按照权利要求1所述的方法,其特征在于,基于所述预测量-心率信号(40)确定心跳相关的脉冲响应函数(43),后者根据心跳次数定义了单独的脉冲响应函数。
3.按照权利要求2所述的方法,其特征在于,首先确定时间相关的单独的脉冲响应函数(37),然后在使用所述预测量-心率信号(40)的条件下,将其换算成心跳相关的脉冲响应函数(43)。
4.按照权利要求1至3中任一项所述的方法,其特征在于,基于当前采集的心率信号(44)和心跳相关的脉冲响应函数(43)或时间相关的单独的脉冲响应函数(37),构造时间相关的校正的单独的脉冲响应函数(45),并且通过结合所述校正的脉冲响应函数(45)与当前的检查协议(48)进行对预计的造影剂变化(CEP)的预测(39)。
5.按照权利要求4所述的方法,其特征在于,为了构造校正的单独的脉冲响应函数(45),利用时间定标因子对单独的心跳相关的脉冲响应函数(43)进行定标,所述时间定标因子通过从当前采集的心率信号(44)中确定的心率(HRCTA(t))的倒数给出。
6.按照权利要求4所述的方法,其特征在于,为了构造校正的单独的脉冲响应函数(45),利用时间定标因子对单独的时间相关的脉冲响应函数(37)进行定标,所述时间定标因子通过从预测量-心率信号(37)中确定的心率(HRTB(t))与从当前采集的心率信号(44)中确定的心率(HRCTA(t))的比例给出。
7.按照权利要求5或6所述的方法,其特征在于,为了构造校正的单独的脉冲响应函数(45),将定标的单独的脉冲响应函数(37)与各自的时间定标因子的倒数相乘。
8.按照权利要求1至7中任一项所述的方法,其特征在于,时间上校正地采集所述血流变化数据和所述预测量-心率信号(40)。
9.按照权利要求1至8中任一项所述的方法,其特征在于,为了采集所述预测量-心率信号(40),进行EKG数据的记录。
10.按照权利要求9所述的方法,其特征在于,从所述EKG数据中导出在预测量期间存在的平均心率(51)。
11.一种用于控制医学技术成像系统(1)的方法,其中,基于给予造影剂的造影剂给予-开始时间点及按照权利要求1至10中任一项所预测的预计的造影剂变化(CEP)确定所述成像系统(1)的参考时间点。
12.一种对预计的造影剂变化(CEP)进行预测的设备(60),用于对检查对象(2)进行造影剂辅助的检查,包括:
-输入接口(61),用于在使用定义的检查协议(47)的条件下在预测量中采集特定于患者的血流变化数据(35);
-脉冲响应函数确定单元(62),用于基于所述血流变化数据(35)确定单独的脉冲响应函数(37,43);
-输入接口(63),用于采集心率信号(40,44,52);
-计算单元(64),用于基于单独的脉冲响应函数(37,43)、与预测量对应的预测量-心率信号以及当前的检查协议(48)和当前采集的心率信号(44,52),对预计的造影剂变化(CEP)进行预测(39)。
13.一种用于成像系统(1)的控制装置(12),包括:
-接口(70),用于采集或输出造影剂给予-开始时间点;
-设备(60),用于按照权利要求11对预计的造影剂变化(CEP)进行预测;以及
-参考时间点确定装置(71),用于基于预计的造影剂变化(CEP)确定用于成像系统(1)的参考时间点。
14.一种用于产生检查对象(2)内部的图像数据组的成像系统(1),包括按照权利要求13所述的控制装置(12)。
15.一种计算机程序产品,其能够直接加载到成像系统(1)的控制装置(12)的存储器中,具有程序代码资源,以便在所述控制装置(12)中执行该程序产品时,执行按照权利要求1至11中任一项所述的方法的全部步骤。
CN201410022486.9A 2013-01-24 2014-01-17 对预计的造影剂变化的预测 Active CN103961117B (zh)

Applications Claiming Priority (2)

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

Publications (2)

Publication Number Publication Date
CN103961117A true CN103961117A (zh) 2014-08-06
CN103961117B CN103961117B (zh) 2016-08-17

Family

ID=51064505

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410022486.9A Active CN103961117B (zh) 2013-01-24 2014-01-17 对预计的造影剂变化的预测

Country Status (3)

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

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018079138A (ja) * 2016-11-17 2018-05-24 株式会社根本杏林堂 注入プロトコルの生成装置、該生成装置を備える注入装置及び撮像システム、注入プロトコル生成方法、及び注入プロトコル生成プログラム
CN108283017A (zh) * 2015-05-13 2018-07-13 拜耳医药股份公司 在使用磁共振系统的诊断成像中优化造影剂浓度的时间分布的预先确定的方法
CN110464326A (zh) * 2019-08-19 2019-11-19 上海联影医疗科技有限公司 一种扫描参数推荐方法、系统、装置及存储介质

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5232473B2 (ja) 2004-11-24 2013-07-10 メドラッド インコーポレーテッド 注入処理中に生成される圧力をモデル化するためのシステム及び装置
JP5486315B2 (ja) 2006-12-29 2014-05-07 メドラッド インコーポレーテッド 医療用注入手順用に患者ベースのパラメータを生成するシステム
BR112012033116A2 (pt) 2010-06-24 2016-10-25 Medrad Inc modelagem da propagação de fármaco e geração de parâmetro para protocolos de injeção
BR112014028442A2 (pt) 2012-05-14 2018-04-24 Bayer Medical Care Inc. sistemas e métodos de protocolos de injeção de fluido farmacêutico com base em voltagem de tubo de raios 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
US10383590B2 (en) * 2015-09-28 2019-08-20 General Electric Company Methods and systems for adaptive scan control
EP3570060A1 (de) * 2018-05-15 2019-11-20 Siemens Healthcare GmbH Verfahren zur rekonstruktion von magnetresonanztomographie-aufnahmen mit variabler zeitauflösung
CN110852993B (zh) * 2019-10-12 2024-03-08 拜耳股份有限公司 一种造影剂作用下的成像方法与设备
CN110916704B (zh) * 2019-11-13 2023-09-01 上海联影医疗科技股份有限公司 校正方法、装置及存储介质
EP4362042A1 (en) 2022-10-31 2024-05-01 Siemens Healthineers AG Estimating pharmacokinetics of a contrast medium through a cardiovascular system analytical model and/or neural network algorithm

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7522744B2 (en) 2004-08-31 2009-04-21 University Of Iowa Research Foundation System and method for adaptive bolus chasing computed tomography (CT) angiography
CN102488527B (zh) 2004-11-16 2015-06-10 拜耳医疗保健公司 从成像过程确定病人生理参数的系统和方法
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
JP2010514506A (ja) 2006-12-29 2010-05-06 メドラッド インコーポレーテッド 薬剤伝搬のモデリング
US8428694B2 (en) 2007-07-17 2013-04-23 Medrad, Inc. Methods for determination of parameters for a procedure, for estimation of cardiopulmonary function and for fluid delivery
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

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108283017A (zh) * 2015-05-13 2018-07-13 拜耳医药股份公司 在使用磁共振系统的诊断成像中优化造影剂浓度的时间分布的预先确定的方法
US10852378B2 (en) 2015-05-13 2020-12-01 Bayer Pharma Aktiengesellschaft Method for optimizing the predetermination of the time profile of a contrast agent concentration in diagnostic imaging using a magnetic resonance system
CN108283017B (zh) * 2015-05-13 2021-01-12 拜耳医药股份公司 在使用磁共振系统的诊断成像中优化造影剂浓度的时间分布的预先确定的方法和磁共振系统
JP2018079138A (ja) * 2016-11-17 2018-05-24 株式会社根本杏林堂 注入プロトコルの生成装置、該生成装置を備える注入装置及び撮像システム、注入プロトコル生成方法、及び注入プロトコル生成プログラム
CN110464326A (zh) * 2019-08-19 2019-11-19 上海联影医疗科技有限公司 一种扫描参数推荐方法、系统、装置及存储介质
CN110464326B (zh) * 2019-08-19 2022-05-10 上海联影医疗科技股份有限公司 一种扫描参数推荐方法、系统、装置及存储介质

Also Published As

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

Similar Documents

Publication Publication Date Title
CN103961117A (zh) 对预计的造影剂变化的预测
EP1970009B1 (en) Continuous x-ray image screening examination device, program, and recording medium
EP2633815B1 (en) Fractional flow reserve estimation
Kühl et al. Assessment of left atrial volume and function: a comparative study between echocardiography, magnetic resonance imaging and multi slice computed tomography
EP3403582A1 (en) Method and apparatus for determining blood velocity in x-ray angiography images
US10674986B2 (en) Methods for personalizing blood flow models
EP2062532B1 (en) X-ray CT apparatus and myocardial perfusion image generating system
EP2509505B1 (en) Method and apparatus for using time of flight information to detect and correct for motion in imaging scans
EP1967140B1 (en) X-Ray CT-Apparatus, Myocardial Perfusion Information Generating System, X-Ray Diagnostic Method and Myocardial Perfusion Information Generating Method
CN101277648B (zh) 用于在血管造影中预测增强的方法和设备
Macron et al. Single-beat versus multibeat real-time 3D echocardiography for assessing left ventricular volumes and ejection fraction: a comparison study with cardiac magnetic resonance
US20170325770A1 (en) Methods for personalizing blood flow models
JP6818492B2 (ja) 画像処理装置、画像処理方法、及びプログラム
EP3806742A1 (en) Method and apparatus for quantitative hemodynamic flow analysis
Jirik et al. Ultrasound perfusion analysis combining bolus-tracking and burst-replenishment
CN105283131A (zh) 用于确定测量开始时间的方法和系统
US10134127B2 (en) Method for post-processing flow-sensitive phase contrast magnetic resonance images
CN103356224B (zh) 用于成像诊断设备的运行的方法以及医学成像系统
Liu et al. Ultrasound‐guided identification of cardiac imaging windows
Thomas Tissue Doppler echocardiography–A case of right tool, wrong use
US10852378B2 (en) Method for optimizing the predetermination of the time profile of a contrast agent concentration in diagnostic imaging using a magnetic resonance system
US20140200461A1 (en) Cardiac analysis based on vessel characteristics
Holland et al. Disposable point‐of‐care portable perfusion phantom for quantitative DCE‐MRI
JP5591512B2 (ja) 血流動態解析装置およびその制御プログラム
EP4362042A1 (en) Estimating pharmacokinetics of a contrast medium through a cardiovascular system analytical model and/or neural network algorithm

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220128

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Address before: Munich, Germany

Patentee before: SIEMENS AG