CN103619257A - 动态灌注成像 - Google Patents

动态灌注成像 Download PDF

Info

Publication number
CN103619257A
CN103619257A CN201280029809.6A CN201280029809A CN103619257A CN 103619257 A CN103619257 A CN 103619257A CN 201280029809 A CN201280029809 A CN 201280029809A CN 103619257 A CN103619257 A CN 103619257A
Authority
CN
China
Prior art keywords
contrast agent
time frame
interest
frame data
area
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
CN201280029809.6A
Other languages
English (en)
Other versions
CN103619257B (zh
Inventor
M·格拉斯
H·施米特
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN103619257A publication Critical patent/CN103619257A/zh
Application granted granted Critical
Publication of CN103619257B publication Critical patent/CN103619257B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • 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/486Diagnostic techniques involving generating temporal series of image data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • 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
    • 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/54Control of apparatus or devices for radiation diagnosis
    • A61B6/542Control of apparatus or devices for radiation diagnosis involving control of exposure
    • 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/54Control of apparatus or devices for radiation diagnosis
    • A61B6/545Control of apparatus or devices for radiation diagnosis involving automatic set-up of acquisition parameters
    • 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
    • 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
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04CROTARY-PISTON, OR OSCILLATING-PISTON, POSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; ROTARY-PISTON, OR OSCILLATING-PISTON, POSITIVE-DISPLACEMENT PUMPS
    • F04C2270/00Control; Monitoring or safety arrangements
    • F04C2270/04Force
    • F04C2270/042Force radial
    • F04C2270/0421Controlled or regulated

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Veterinary Medicine (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

一种方法,包括:在基于造影剂的灌注扫描期间,在感兴趣区域中造影剂摄取期间以预定时间采样率扫描所述感兴趣区域;并且及生成指示所扫描的感兴趣区域的时间帧数据。所述方法还包括:从所述时间帧数据,识别所述感兴趣区域中所述造影剂的量的预定改变。所述方法还包括:响应于识别出所述感兴趣区域中所述造影剂的量的所述预定改变,以低的时间采样率扫描所述感兴趣区域,所述低时间采样率低于所述造影剂摄取期间的所述时间采样率。

Description

动态灌注成像
技术领域
以下大体涉及动态灌注成像,并且结合对计算机断层摄影(CT)的应用进行描述;然而,以下也适用于其他成像模态。
背景技术
灌注成像这样一种成像技术:其捕获施予的造影剂的通过诸如血管的感兴趣脉管组织和/或如心脏的器官的运输。一般地,针对灌注成像,将造影剂团注剂施予给患者,并且扫描所述患者的感兴趣区域,包括感兴趣脉管组织。所述造影剂引起所述感兴趣脉管组织中的X射线密度因所述造影剂流动通过所述脉管组织而暂时增加。典型的灌注扫描包括在多个时间间隔上采集同一区域的数据,所述多时间间隔覆盖造影剂到达、摄取和洗出。针对心脏应用,所述扫描包括采集同一心脏相位的数据。
对所采集数据的分析可以被用于,例如,基于在扫描视场中对造影剂动力学的观察,确定所述感兴趣脉管组织的灌注状态。对于心脏应用,这可以包括量化心肌中随时间的造影剂分布。这种分析可以包括确定针对所述感兴趣脉管组织的各种灌注相关信息,例如时间衰减曲线、血流量、血体积、平均通过时间、最大向上斜率、到达峰值时间,等等。该信息能被用于识别缺血组织和/或用于在不可逆损伤(或坏死)组织与潜在可逆损伤(或风险)组织之间进行区分。
传统的灌注成像包括从造影剂到达之前直到造影剂洗出,连续扫描感兴趣区域。最近的灌注成像包括了时间间歇,以相等时间距离从造影剂到达直到造影剂洗出对所述感兴趣区域进行扫描。一般而言,时间间歇采样的频率是基于获得适于准确推导灌注参数的数据的所必需的时间采样,所述灌注参数例如到达峰值时间、最大向上斜率和/或其他相关灌注参数。对于心脏应用,这已包括了在每个或每隔一个心脏周期的一个或多个特定感兴趣心脏运动相位(例如平静相位)期间进行扫描。
遗憾的是,计算机断层摄影灌注成像使患者暴露于电离辐射,其能杀死或损伤细胞并且其可以增加癌症的风险,并且认为以连续和时间间隔成像所沉积的剂量是高的,并且一般不将这种成像用于筛查和/或常规临床实践中。更近一步地,经历这种成像的患者通常经历几次后续成像程序,这增加了累积辐射剂量。此外,简单地减少时间间隔采样可能在灌注参数中引入误差。因此,存在着对于进一步减少灌注成像的患者剂量的其他途径的待解决需要。
发明内容
本申请的各方面解决以上提到的问题及其他问题。
根据一个方面,一种方法包括:在基于造影剂的灌注扫描期间,在感兴趣区域中的造影剂摄取区间以预定时间采样率扫描所述感兴趣区域;并且生成指示所扫描的感兴趣区域的时间帧数据。所述方法还包括:从所述时间帧数据,识别所述感兴趣区域中造影剂的量的预定改变。所述方法还包括:响应于识别出所述感兴趣区域中所述造影剂的所述量的所述预定改变,以低的时间采样率扫描所述感兴趣区域,所述低的时间采样率低于所述造影剂摄取期间的所述时间采样率。
在另一方面中,一种系统包括具有造影剂识别器的造影剂检测器,所述造影剂识别器识别时间帧灌注扫描采集数据中造影剂的存在并且生成指示所述时间帧采集数据中所述造影剂的所述存在的造影剂信号。所述造影剂信号被用于:在所述造影剂不存在于所述时间帧数据中时的低时间采样率与所述造影剂存在于所述时间帧数据中时的较高时间采样率之间,变化时间采样率。
在另一方面中,一种被编码于计算机可读介质上的计算机可读指令,其在被计算系统的处理器运行时,令所述处理器:基于在得到的采集时间帧数据中造影剂的状态,在灌注扫描期间变化数据采集的时间采样率,其中,基于使用对所述时间帧数据的基于模型的分割确定的造影剂水平,至少在造影剂摄取与造影剂洗出之间变化所述时间采样率。
本领域技术人员在时阅读和理解以下详细描述后,将认识到本发明再进一步的方面。
附图说明
本发明可以采取各种部件和部件的布置,以及各种步骤和步骤的安排的形式。附图仅出于图示优选的实施例的目的,并且不应被解释为限制本发明。
图1示意性地图示了示范性成像系统,其中,基于在时间帧采集数据中的造影剂水平,随时间选择性地并且动态地变化用于灌注成像的时间采样。
图2示意性地图示了用于检测时间帧采集数据中的造影存在和/或水平的示范性造影检测器。
图3以图形的方式图示了用于检测造影存在和/或造影剂水平的示范性建模方法。
图4以图形的方式图示了用于心脏灌注扫描的示范性时间采样,其中,基于造影剂水平曲线,在造影摄取前、造影摄取和造影洗出期间,分别使用不同但恒定的时间采样。
图5以图形的方式图示了针对心脏灌注扫描的示范性时间采样,其中,随时间变化造影摄取期间的时间采样。
图6以图形的方式图示了针对心脏灌注扫描的示范性时间采样,其中,随时间变化在造影洗出期间时间采样。
图7以图形的方式图示了针对心脏灌注扫描的示范性时间采样,其中时间采样开始与造影摄取。
图8以图形的方式图示了用于后期增强扫描的时间采样的范例。
图9图示了用于基于时间帧采集数据中的造影剂水平来随时间变化灌注成像的时间采样的流程图。
具体实施方式
图1图示了成像系统100,例如计算机断层摄影(CT)扫描器。
成像系统100包括固定机架102和旋转机架104,旋转机架104由固定机架102以能够旋转的方式支撑。旋转机架104关于纵轴或z轴110,绕检查区域106以及检查区域106中的对象或受试者108的部分旋转。
辐射源112,例如X射线管,由旋转机架104支撑,并与旋转机架104一起绕检查区域106旋转。辐射源112发射辐射,所述辐射被准直以形成穿过检查区域106的大体为扇形、楔形或锥形的辐射束。
源控制器114控制辐射源112。该控制包括激活源112以在时间序列灌注成像期间发射辐射以采集时间帧数据。这种激活可以针对源112的一次或多次旋转(包括针对源112的每次旋转的整体或子部分),在成像检查期间为时间连续的或时间间歇的。
辐射敏感探测器阵列116探测由辐射源112发射的穿过检查区域106的光子,并且生成指示所探测的辐射的投影数据。辐射敏感探测器阵列116可以包括一维或二维探测器阵列。
重建器118重建所述投影数据,并且生成指示检查区域106的时间序列体积图像数据。可以采用各种重建算法,例如滤波后向投影、统计、迭代、稀疏采样和/或其他重建算法。
注射器120被配置为将诸如一种或多种造影剂的材料注射或施予给要被扫描的受试者108或对象(例如体模)。可以额外地或备选地由临床医师等手动施予造影剂。在向受试者108手动施予造影剂时,可以省略注射器120。
运动传感器122检测受试者108的移动结构的运动,并且生成指示所述运动的运动信号。移动结构的范例包括心、肺和/或其他移动结构。对于心脏应用,运动传感器122可以包括心电图仪(ECG);并且对于呼吸应用,运动传感器122可以包括呼吸带、发光标志等。在扫描期间不考虑运动时,可以省略运动传感器122。
运动相位识别器124识别所述运动信号中的运动状态(一个或多个),并且生成指示所述运动状态的运动相位信号。对于心脏应用,这可以包括识别心脏周期的持续时间,识别心脏周期中的标志(例如“R”波),并且然后基于所述心脏周期和所述“R”波的预定百分比(例如40%、70%等等),来识别相位。所述相位可以为平静状态,其中心脏运动相对于其他状态或另一相位更低。
造影剂检测器126生成指示在所述采集时间数据帧中造影剂的存在和/或造影剂水平的信号。下文中将更加详细地描述,在一种非限制性情况中,该生成包括基于针对所述结构的模型或以其他方式,确定在所述采集时间帧数据中结构的像素强度、所述结构的子集或所述结构的子集的子区域。
受试者支撑体128(例如躺椅)支撑检查区域106中的受试者108或对象,并且能够协同旋转机架104的旋转沿x轴、y轴和z轴方向移动,以方便螺旋、轴向或其他期望扫描轨迹。
通用计算系统充当操作者控制台130。驻留于控制台130上并由处理器运行的软件(例如)通过允许所述操作者采用来自扫描协议存储器132的扫描协议(例如(基于运动相位或非运动相位的)自适应时间采样动态灌注扫描协议),来允许所述操作者控制系统100的操作。
下文中将更加详细地描述,一种这样的协议包括基于在采集时间帧数据中的造影剂,在灌注扫描期间动态地调节时间采样。这允许在造影剂摄取前和/或造影剂洗出期间以低时间采样,并且在造影剂摄取期间以较高时间采样(适于确定相关灌注参数)来采集数据。这可以有利于相对于其中贯穿造影摄取前、造影摄取和造影洗出均使用单一时间采样的构造,减小剂量。
对于基于运动相位的灌注扫描,例如心脏灌注扫描,数据采集也可以包括针对心脏周期的整体或针对一个或多个相位(例如收缩和/或舒张,和/或它们的部分)采集数据,用于时间帧采集。来自运动相位识别器124的运动相位信号可以被用于基于心脏相位来触发心脏周期内的数据采集。相对于在整个心脏周期期间扫描,仅在心脏周期的子部分(或相位)期间扫描,可以进一步减小剂量。
可以通过减少空间采样和/或管电流来进一步减小剂量。对于空间采样,这包括将所述空间采样减小为稀疏采样,并且使用任意已知的或其他的稀疏采样重建算法来重建所述数据。在一种情况中,所述稀疏空间采样可以为约全空间采样的十分之一。可以在造影摄取前和/或造影洗出期间减少管电流。在一种情况中,所述管电流可以被减小到在造影摄取期间使用的管电流的三分之一。
图2结合源112、源控制器114、重建期118和控制台130示意性地图示了造影剂检测器126的范例。
所图示的造影剂检测器126包括造影剂识别器202,其被配置为例如基于来自造影算法存储器204的造影识别算法,来评价采集时间帧数据,并且生成指示所述数据中造影的存在(和/或不存在)和/或程度的造影剂信号。
通过非限制性范例的方式,造影剂识别器202利一种造影识别算法确定采集时间帧数据中结构的平均像素强度。造影剂识别器202然后比较所确定的强度与基线阈值强度。在该情况中,如果所确定的强度满足所述阈值,则造影剂识别器202生成造影存在的造影剂信号。否则,造影剂识别器202可以不生成造影剂信号,或者可以生成指示不存在造影(或造影不存在)的造影剂信号。
以另一种造影识别算法,造影剂识别器202采用基于解剖模型的分割方法,以将所述采集时间帧数据分割成子结构(例如心、肺,等等)。造影剂识别器202然后确定所述子结构的子集的平均像素强度,并且比较所述强度与所述阈值。类似地,造影剂识别器202生成指示存在造影的造影剂信号,或者生成指示造影不存在的造影剂信号或不生成信号。
在再另一算法中,造影剂识别器202还基于所述模型,将子结构(例如心)分割成子区域(例如左心室、右心室、左心房、右心房,等等)。图3中例示了这样的范例,其以图形的方式图示了被分割为子区域302、304、306和308的心的图像300。造影剂识别器202确定一个或多个所述子区域(例如所述造影最后达到的子区域)的平均像素强度,比较所述强度与预定阈值,并且基于其生成如上文讨论的造影剂信号。
本文也预期用于识别造影的存在和/或水平的其他方法。
峰值造影剂识别器206被配置为评价所述造影剂信号,并且响应于从所述造影剂信号识别出峰值造影而生成峰值造影剂信号。可以使用各种方法来识别所述峰值。通过非限制性范例的方式,峰值造影剂识别器206能基于对相继确定的造影剂水平的比较、造影剂水平改变速率和/或其他方法来识别所述峰值。
造影剂检测器126独自地或者组合地将所述造影剂信号和/或所述峰值造影剂信号作为造影剂检测器126的信号输出而传达到控制台130,控制台130采用所述信号,以控制源控制器114来激活源112,从而在灌注成像期间适应所述时间采样。
图4结合图2和图3以图形的方式图示了在灌注扫描期间变化时间采样的范例。
y轴包括表示造影的量(密度)的上部分402和表示图像采集的下部分404,x轴406表示时间,曲线408表示作为时间的函数的造影剂摄取和洗出,并且曲线409表示时间采样(或者作为时间的函数的数据采集)。
区域410表示的区域中造影剂已被施予但尚未到达正被扫描的感兴趣区域并且没有检测到造影剂。在该范例中,在该区域期间,源控制器114控制源112,以第一时间采样或造影剂摄取前时间采样412发射辐射。
相对于摄取期间的时间采样,该区域中的时间采样412可以低,因为该区域主要被用于建立基线强度,以在造影剂到达被扫描的感兴趣区域时进行识别。相对于使用与在造影剂摄取期间使用的相同的时间采样,低采样减少了该区域中的剂量。也可以减小针对该区域管电流。
区域414表示的区域中造影剂到达被扫描的感兴趣区域,并且来自造影剂检测器126的信号(例如来自造影剂识别器202的造影剂信号)指示在被扫描的感兴趣区域中造影剂的存在。在该范例中,在该区域期间,源控制器114控制源112,以第二时间采样或造影剂摄取时间采样416发射辐射。
时间采样416被用于灌注扫描,以用可以针对区域414确定的相关灌注信息(例如峰值时间、最大上斜率,等等)时的时间采样来采集数据。在心脏灌注的情况中,时间采样416可以基于正被扫面的患者的平均心脏周期。例如,时间采样416可以在每次心脏周期至每隔一次心脏周期的量级上(例如1至2秒)。
区域418表示这样的区域,其中由造影剂检测器126识别峰值造影(例如来自峰值造影剂识别器206的峰值造影剂信号)。
区域420表示这样的区域,其中被扫描感兴趣区域的造影剂洗出和来自造影剂检测器126的信号(例如来自造影剂识别器202的造影剂信号)指示在被扫描感兴趣区域中存在的造影剂减少。在该范例中,该区域期间,源控制器114控制源112,以等于或大于第一时间采样412并且低于第二时间采样416的第三时间采样或造影剂洗出时间采样422来发射辐射。
通常,以此方法,针对图像采集的时间采样根据灌注时间曲线(图4)进行调整,其中所述时间采样在造影摄取期间更密集(图4中的向上斜率)并且造影摄取前和造影洗出与造影摄取期间相比较不密集(图4中的向下斜率),并且对合适使用那种采样模式的决策是基于使用某些结构的像素强度的基于模型的图像分析(图2和图3)。
预期各种变型。
图5以图形的方式图示了在灌注扫描期间变化时间采样的另一范例。图5基本上类似于图4;但是,在图5中,在造影剂摄取期间使用变化时间采样502。
图6以图形的方式图示了在灌注扫描期间变化时间采样的另一范例。图6也基本上类似于图4;但是,在图6中,在造影剂洗出期间使用变化时间采样602。
图7基本上类似于图4,但是在造影剂摄取前区域410期间不执行扫描。而是,数据采集在从造影剂的施予的一段预定时间延迟之后开始,可以基于测试团注、所述患者的先前灌注扫描、其他患者的历史灌注扫描和/或以其他方式,确定所述一段预定时间延迟。
在另一种变型中,所述时间采样为图4、图5、图6和/或图7和或其他时间采样的组合。
在图示的实施例中,造影剂识别器126评价重建的时间帧数据。在变型中,造影剂识别器126可以单独地和/或结合所重建的时间帧数据,评价投影数据,并生成造影剂信号。
在图示的实施例中,运动相位识别器124评价运动信号,以识别其中的一个或多个运动状态。在变型中,运动相位识别器124单独地和/或结合所重建的时间帧数据来评价投影数据,以识别所述一个或多个运动状态。这可以包括比较连续采集的时间帧(例如,确定差别图像),以识别平静、类似和/或其他运动相位。
在另一种变型中,所述造影剂摄取时间采样可以被延长超过峰值造影,例如,越过预定时间延迟,例如,为了在所述造影剂摄取时间采样中扫描后期增强。图8示出了以距峰值418固定时间距离804时的后期增强扫描802的范例。在其他实施例中,固定时间距离804可以是从另一参数,例如最大向上斜率、到达最大向上斜率时间或任意其他参数。针对所述后期增强的扫描数据可以被用于识别不可逆损伤组织。
在图2和图3中,峰值造影剂识别器206被用于识别峰值造影剂水平,并然后改变时间采样。在另一种情况中,分析所述数据,以在图像300和/或子区域302至308中定位何时到达预定造影剂水平。在该情况中,响应于所述预定造影剂水平到达预定造影剂水平而改变时间采样。在再另一种情况中,分析所述数据,以确定所述预定的一个或多个子区域302至308中的造影剂水平何时到达或超过(或降至以下)一个或多个其他子区域302至308中的造影剂水平。在该情况中,响应于子区域之间的造影剂水平的关系满足条件而改变所述时间采样。
图9图示了范例灌注成像流程图。
要认识到,动作的顺序并非限制性的。这样,在其他实施例中,顺序可以不同,包括并行。更近一步地,可以省略一个或多个动作和/或可以增加一个或多个动作。
在902,将造影剂施予给受试者。
在904,基于造影剂摄取前时间采样扫描受试者的感兴趣区域。通过非限制性范例的方式,对于心脏灌注成像,可以扫描所述受试者以至少在每第五次心脏周期的子部分期间采集数据。所述子部分可以对应于所述感兴趣区域的感兴趣结构的一个或多个特定运动相位。也可以将管电流设定为造影剂摄取前水品。
在906,确定采集时间帧扫描数据的造影剂水平。如本文中所描述,可以基于所述数据中所述结构的平均像素强度、所述数据中所述结构的子集的平均像素强度、所述结构的子集的一个或多个子区域的平均像素强度,和/或以其他方式,确定所述造影剂水平。初始数据的强度被用于确定基线非造影剂水平。
在908,基于所述造影剂水平和基线非造影剂水平,确定在所述采集时间帧扫描数据中是否存在造影。在一种情况中,这通过比较所述造影剂水平与所述基线非造影剂水平,得以实现。
如果所施予造影在所述采集时间帧数据中不存在(例如,所述造影剂水平小于所述基线非造影剂水平,或在预定容许内),则重复动作904至908。
如果所施予的造影剂存在于所述采集时间帧数据中(例如,所述造影剂水平以所述预定容许大于所述基线非造影剂水平),则在910,基于大于所述造影摄取前时间采样的预定造影摄取时间,采样扫描所述受试者的所述感兴趣区域。通过非限制性范例的方式,以心脏灌注成像,所述造影摄取时间采样可以为每个心脏周期或每隔一个心脏周期。
在912,确定所述采集时间帧数据的所述造影剂水平,例如,如结合动作906所讨论的。
在914,基于所述采集时间帧数据,确定是否到达峰值造影剂摄取。这可以通过比较当前造影剂水平与先前确定的造影剂水平,并探测减少的造影剂水平,和/或以其他方式,得以实现。
如果在所述采集时间帧数据中未达到所述峰值造影剂水平,则重复动作910至914。
如果在所述采集时间帧数据中达到了所述峰值造影,则在916,基于小于所述造影剂摄取时间采样的预定造影剂洗出时间采样,扫描所述对象的所述感兴趣区域。通过非限制性范例的方式,对于心脏灌注成像,所述造影剂洗出时间采样可以为每第五个心脏周期。
在918,确定所述采集时间帧数据的所述造影剂水平,例如,如结合动作906和912所讨论的。
在920,确定造影剂是否洗出所述感兴趣区域。这可以通过比较所述造影剂水平与所述基线非造影剂水平,得以实现。
如果造影剂仍存在与所述采集时间帧数据中,则重复动作916至920。
如果在所述采集时间帧数据中不存在造影剂,则在922,停止扫描。
以上动作可以用计算机可读指令的方式来实现,所述计算机可读执行在被(一个或多个)计算机处理器运行时,引起(一个或多个)所述处理器执行所描述的动作。在这样的情况中,所述指令被存储在与相关计算机相关联或以其他方式由所述相关计算机可存取的计算机可读存储介质中。不需要与数据采集并行地执行所述行为。
尽管上文已结合常规扫描器进行了描述,但是应该理解系统100可以包括光谱CT扫描器。这样的扫描器可以包括以下中的一个或多个:被配置为发射具有不同发射光谱的辐射的多X射线管、被配置为在至少两个不同发射光谱之间切换系统的X射线管的管电压切换电路和/或能量解析探测器(例如,光子计数探测器、带有多个闪烁体/光电二极管对的探测器,每个对被配置为探测不同的预定能量范围内的光子)。在该情况中,所述灌注数据表示每组织体积的造影剂浓度,并且所述浓度可以通过多能量成像、材料分析,和/或k边缘成像来确定。
已参考多个实施例描述了本发明。他人在阅读了本文的描述之后可以进行修改和变型。目的是本发明被解释为包括所有这些修改和变型,只要它们落在权利要求或其等价方案的范围之内。

Claims (23)

1.一种方法,包括:
在基于造影剂的灌注扫描期间,在感兴趣区域的造影剂摄取期间以预定时间采样率扫描所述感兴趣区域,并且生成指示所扫描的感兴趣区域的时间帧数据;
从所述时间帧数据,识别所述感兴趣区域中造影剂的量的预定改变;并且
响应于识别出所述感兴趣区域中所述造影剂的所述量的所述预定改变,以低时间采样率扫描所述感兴趣区域,所述低时间采样率低于所述造影剂摄取期间的所述时间采样率。
2.如权利要求1所述的方法,其中,所述预定改变为所述造影剂的所述量的减少,其中,所述减少指示造影剂洗出。
3.如权利要求1至2中任一项所述的方法,还包括:
确定所述时间帧数据中的结构的平均像素强度;
将所述平均像素强度与前面的时间帧数据的平均像素强度进行比较;并且
基于所述比较的结果,识别所述预定改变。
4.如权利要求1所述的方法,还包括:
将所述时间帧数据分割成子结构;
确定所述子结构的子集的平均像素强度;
将所述平均像素强度与相应的前面的时间帧数据的平均像素强度进行比较;并且
基于所述比较的结果,识别所述预定改变。
5.如权利要求1所述的方法,还包括:
将所述时间帧数据分割成子结构;
将至少一个子结构分割成子区域;
确定所述子区域的子集的平均像素强度;
将所述平均像素强度与相应的前面的时间帧数据的平均像素强度进行比较;并且
基于所述比较的结果,识别所述预定改变。
6.如权利要求3至5中任一项所述的方法,还包括:
响应于前面的时间帧数据的所述平均像素强度超过所述时间帧数据的所述平均像素强度,识别所述感兴趣区域中所述造影剂的所述量的所述预定改变。
7.如权利要求1所述的方法,还包括:
将所述时间帧数据分割成子结构;
将至少一个子结构分割成子区域;
确定所述子区域中的至少一个的平均像素强度;
确定所述子区域中的至少另一个的平均像素强度;
比较所述的平均像素强度;并且
基于所述比较的结果,识别所述预定改变。
8.如权利要求1至7中任一项所述的方法,还包括:
在距所述预定改变的预定固定时间距离执行数据采集,以捕获后期造影剂增强。
9.如权利要求1至8中任一项所述的方法,还包括:
在所述感兴趣区域的所述造影剂摄取之前,以第二低时间采样率扫描所述感兴趣区域,并且生成指示所扫描的感兴趣区域的时间帧数据;
从所述时间帧数据,识别所述感兴趣区域中所述造影剂的存在;并且
响应于识别出所述感兴趣区域中所述造影剂的所述存在,以所述预定时间采样率扫描所述感兴趣区域。
10.如权利要求1至9中任一项所述的方法,还包括:
感测所述感兴趣区域的运动周期;
识别所述运动周期的子部分;并且
仅在所述运动周期的所述子部分期间,扫描所述感兴趣区域。
11.如权利要求10所述的方法,其中,所述子部分对应于所述运动周期的低运动相位。
12.如权利要求10至11中任一项所述的方法,还包括:
在所述子部分中的至少两个期间,扫描所述感兴趣区域。
13.如权利要求1至12中任一项所述的方法,还包括:
使用稀疏空间采样进行扫描;并且
使用稀疏重建算法,重建所述时间帧数据。
14.如权利要求1至13中任一项所述的方法,还包括:
在造影剂摄取前扫描或造影剂洗出扫描中的至少一个期间使用低管电流并且在造影剂摄取扫描期间使用高管电流,其中,所述高管电流大于所述低管电流。
15.一种系统,包括:
造影剂检测器(126),其包括:
造影剂识别器(202),其识别时间帧灌注扫描采集数据中造影剂
的存在,并且生成指示所述时间帧采集数据中所述造影剂的所述存在
的造影剂信号,
其中,所述造影剂信号被用于在所述造影剂不存在于所述时间帧数据中时的低时间采样率与所述造影剂存在于所述时间帧数据中时的高时间采样率之间改变时间采样率。
16.如权利要求15所述的系统,其中,所述造影剂识别器使用基于模型的分割方法,来识别在所述时间帧采集数据中造影剂的所述存在。
17.如权利要求16所述的系统,其中,所述基于模型的分割方法将感兴趣结构分割成子区域,并且所述造影剂识别器识别所述子区域的子集中所述造影剂的所述存在。
18.如权利要求16至17中任一项所述的系统,其中,所述造影剂识别器基于像素强度,识别所述造影剂的所述存在。
19.如权利要求15所述的系统,其中,造影剂检测器,还包括:
峰值造影剂识别器(206),其在时间帧采集数据中识别峰值造影剂,并且生成指示所述时间帧采集数据中的所述峰值造影剂的峰值造影剂信号,其中,所述峰值造影剂信号被用于响应于指示所述时间帧采集数据中的峰值造影剂的所述峰值造影剂信号,在所述较高时间采样率和第二低时间采样率之间变化时间采样率。
20.如权利要求15至19中任一项所述的系统,还包括:
成像系统(100),其执行生成所述时间帧灌注扫描采集数据的灌注扫描。
21.如权利要求20所述的系统,其中,所述成像系统为光谱计算机断层摄影扫描器。
22.如权利要求21所述的系统,其中,通过材料分离和/或k边成像来确定造影剂水平的浓度。
23.一种编码有计算机可读指令的计算机可读存储介质,所述计算机可读指令在由处理器运行时,令所述处理器:
基于得到的采集时间帧数据中的造影剂的状态,在灌注扫描期间变化数据采集的时间采样率,其中,基于使用对所述时间帧数据的基于模型的分割确定的造影剂水平,至少在造影剂摄取与造影剂洗出之间变化所述时间采样率。
CN201280029809.6A 2011-04-18 2012-04-17 动态灌注成像 Active CN103619257B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161476462P 2011-04-18 2011-04-18
US61/476,462 2011-04-18
PCT/IB2012/051914 WO2012143852A1 (en) 2011-04-18 2012-04-17 Dynamic perfusion imaging

Publications (2)

Publication Number Publication Date
CN103619257A true CN103619257A (zh) 2014-03-05
CN103619257B CN103619257B (zh) 2016-06-29

Family

ID=46149697

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280029809.6A Active CN103619257B (zh) 2011-04-18 2012-04-17 动态灌注成像

Country Status (5)

Country Link
US (1) US9848843B2 (zh)
EP (1) EP2699163B1 (zh)
CN (1) CN103619257B (zh)
RU (1) RU2013150999A (zh)
WO (1) WO2012143852A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105877842A (zh) * 2015-02-12 2016-08-24 西门子公司 分析动态的造影剂分布
CN108095751A (zh) * 2017-12-18 2018-06-01 上海联影医疗科技有限公司 Ct灌注扫描方法、系统及存储介质
CN110547822A (zh) * 2018-05-30 2019-12-10 通用电气公司 造影成像系统和方法
CN111429456A (zh) * 2020-05-06 2020-07-17 上海联影医疗科技有限公司 一种图像数据处理方法、装置、电子设备及存储介质

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102011003606A1 (de) * 2011-02-03 2012-08-09 Siemens Aktiengesellschaft Verfahren zum Entfernen unerwünschter perfusionsbedingter Signale aus MR-Bildern, Magnetresonanzanlage, Computerprogrammprodukt und elektronisch lesbarer Datenträger
DE102011076855B4 (de) * 2011-06-01 2017-12-07 Siemens Healthcare Gmbh Verfahren zur funktionalen Darstellung und Lokalisierung einer arteriovenösen Fehlbildung, rotierbares Bildgebungssystem und Kombination aus einem rotierbaren Bildgebungssystem und einer Bestrahlungseinheit
US9283317B2 (en) * 2013-03-15 2016-03-15 General Electric Company Functional brown adipose tissue imaging technique
DE102014202513B4 (de) * 2014-02-12 2016-09-01 Siemens Healthcare Gmbh Verfahren und Magnetresonanzanlage zur Erstellung von MR-Bildern
CN106061390B (zh) * 2014-02-21 2021-10-26 西门子医疗有限公司 用于记录医疗图像的方法和设备
JP2019202078A (ja) * 2018-05-25 2019-11-28 キヤノン株式会社 放射線撮影装置、放射線撮影システム、放射線撮影方法及びプログラム
US10638991B2 (en) * 2018-09-20 2020-05-05 General Electric Company System and method for imaging a subject with a contrast agent
US11276176B2 (en) 2019-09-04 2022-03-15 International Business Machines Corporation Intelligent boundary delineation of regions of interest of an organism from multispectral video streams using perfusion models
US11526703B2 (en) 2020-07-28 2022-12-13 International Business Machines Corporation GPU accelerated perfusion estimation from multispectral videos

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997040395A2 (en) * 1996-04-24 1997-10-30 Philips Electronics N.V. Method of and device for the imaging of objects by means of magnetic resonance
US20040023406A1 (en) * 2002-07-09 2004-02-05 Schomacker Kevin T. Optimal windows for obtaining optical data for characterization of tissue samples
US20040082846A1 (en) * 2002-10-25 2004-04-29 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
CN101001569A (zh) * 2004-02-06 2007-07-18 韦克福里斯特大学健康科学院 利用非侵入性成像的整体组织特性的组织评估和用于确定影像整体组织特性的系统
CN101291622A (zh) * 2005-10-17 2008-10-22 皇家飞利浦电子股份有限公司 利用生理模型来控制功能获取
CN101317766A (zh) * 2007-06-08 2008-12-10 通用电气公司 用ct投影数据进行高时间分辨率药丸检测的方法和系统
CN101502422A (zh) * 2008-02-07 2009-08-12 株式会社东芝 X射线ct装置以及断层摄影方法
CN101652100A (zh) * 2007-03-29 2010-02-17 皇家飞利浦电子股份有限公司 用于获取融合x射线图像的方法和装置
CN101686825A (zh) * 2007-06-21 2010-03-31 皇家飞利浦电子股份有限公司 使用动态模型调整用于动态医学成像的采集协议
CN101978375A (zh) * 2008-03-17 2011-02-16 皇家飞利浦电子股份有限公司 灌注成像

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6983182B2 (en) 1998-04-10 2006-01-03 Wisconsin Alumni Research Foundation Time resolved computed tomography angiography
US6201986B1 (en) 1998-11-24 2001-03-13 Mayo Foundation For Medical Education And Research Synchronized K-space sampling in magnetic resonance angiography
US6556856B1 (en) 1999-01-08 2003-04-29 Wisconsin Alumni Research Foundation Dual resolution acquisition of magnetic resonance angiography data with vessel segmentation
US6597937B2 (en) 2001-02-05 2003-07-22 Koninklijke Philips Electronics N.V. Self-adaptive tracking and phase encoding during data collection for contrast-enhanced MRA and dynamic agent uptake studies
FI116750B (fi) * 2002-08-28 2006-02-15 Instrumentarium Corp Lääketieteellisen röntgenkuvauksen menetelmä ja järjestely
JP5433299B2 (ja) * 2009-05-18 2014-03-05 株式会社東芝 医用画像診断装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1997040395A2 (en) * 1996-04-24 1997-10-30 Philips Electronics N.V. Method of and device for the imaging of objects by means of magnetic resonance
WO1997040395A3 (en) * 1996-04-24 1997-12-04 Philips Electronics Nv Method of and device for the imaging of objects by means of magnetic resonance
US20040023406A1 (en) * 2002-07-09 2004-02-05 Schomacker Kevin T. Optimal windows for obtaining optical data for characterization of tissue samples
US20040082846A1 (en) * 2002-10-25 2004-04-29 Koninklijke Philips Electronics N.V. Method and apparatus for volumetric cardiac computed tomography imaging
CN101001569A (zh) * 2004-02-06 2007-07-18 韦克福里斯特大学健康科学院 利用非侵入性成像的整体组织特性的组织评估和用于确定影像整体组织特性的系统
CN101291622A (zh) * 2005-10-17 2008-10-22 皇家飞利浦电子股份有限公司 利用生理模型来控制功能获取
CN101652100A (zh) * 2007-03-29 2010-02-17 皇家飞利浦电子股份有限公司 用于获取融合x射线图像的方法和装置
CN101317766A (zh) * 2007-06-08 2008-12-10 通用电气公司 用ct投影数据进行高时间分辨率药丸检测的方法和系统
CN101686825A (zh) * 2007-06-21 2010-03-31 皇家飞利浦电子股份有限公司 使用动态模型调整用于动态医学成像的采集协议
CN101502422A (zh) * 2008-02-07 2009-08-12 株式会社东芝 X射线ct装置以及断层摄影方法
CN101978375A (zh) * 2008-03-17 2011-02-16 皇家飞利浦电子股份有限公司 灌注成像

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105877842A (zh) * 2015-02-12 2016-08-24 西门子公司 分析动态的造影剂分布
CN108095751A (zh) * 2017-12-18 2018-06-01 上海联影医疗科技有限公司 Ct灌注扫描方法、系统及存储介质
CN108095751B (zh) * 2017-12-18 2021-02-19 上海联影医疗科技股份有限公司 Ct灌注扫描方法、系统及存储介质
CN110547822A (zh) * 2018-05-30 2019-12-10 通用电气公司 造影成像系统和方法
CN110547822B (zh) * 2018-05-30 2024-03-08 通用电气公司 造影成像系统和方法
US11972563B2 (en) 2018-05-30 2024-04-30 General Electric Company Contrast imaging system and method
CN111429456A (zh) * 2020-05-06 2020-07-17 上海联影医疗科技有限公司 一种图像数据处理方法、装置、电子设备及存储介质
CN111429456B (zh) * 2020-05-06 2023-09-05 上海联影医疗科技股份有限公司 一种图像数据处理方法、装置、电子设备及存储介质

Also Published As

Publication number Publication date
US20140037046A1 (en) 2014-02-06
EP2699163A1 (en) 2014-02-26
EP2699163B1 (en) 2019-06-26
CN103619257B (zh) 2016-06-29
WO2012143852A1 (en) 2012-10-26
RU2013150999A (ru) 2015-05-27
US9848843B2 (en) 2017-12-26

Similar Documents

Publication Publication Date Title
CN103619257A (zh) 动态灌注成像
JP5498787B2 (ja) エネルギー感受性コンピュータ断層撮影における動き補償
US8315353B1 (en) System and method of prior image constrained image reconstruction using short scan image data and objective function minimization
US7507968B2 (en) Systems and methods for correcting a positron emission tomography emission image
CN106691487B (zh) 成像方法和成像系统
US9684973B2 (en) Systems and methods for selecting imaging data for principle components analysis
JP2013513792A (ja) 飛行時間情報を用いて画像化スキャンにおける動き検出及び補正をする方法と装置
US11918390B2 (en) Methods and systems for motion detection in positron emission tomography
US11179128B2 (en) Methods and systems for motion detection in positron emission tomography
US9466132B2 (en) Systems and methods for motion mitigation determinations
CN110536640B (zh) 从pet列表数据中的呼吸运动信号的噪声鲁棒的实时提取
US20190133542A1 (en) Systems and methods for data-driven respiratory gating in positron emission tomography
CN103080977B (zh) 在发射图像中的心脏定位及验证
Li et al. Multienergy cone-beam computed tomography reconstruction with a spatial spectral nonlocal means algorithm
CN102933142A (zh) 基于幅度/斜率的运动时相映射
CN104379063A (zh) 基于谱成像的决策支持、处置规划和/或介入引导
JP2021511881A (ja) 定量分析が向上された低放射線量コンピュータ断層撮影灌流(ctp)
CN107072627A (zh) 用于针对心脏和肺图像的光子计数ct的心脏重建
US9760988B2 (en) Myocardial clusterification and orientation
US8867810B2 (en) Automatic identification of disruptive events in imaging scans
Daouk et al. Improved attenuation correction via appropriate selection of respiratory-correlated PET data
Breuilly et al. Amplitude-based data selection for optimal retrospective reconstruction in micro-SPECT
Holbrook System development and preclinical applications of photon counting micro-CT
Sanaat et al. A novel convolutional neural network for predicting full dose from low dose PET scans
Breuilly et al. Image-based motion detection in 4D images and application to respiratory motion suppression

Legal Events

Date Code Title Description
PB01 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