CN104244827A - 灌注成像 - Google Patents

灌注成像 Download PDF

Info

Publication number
CN104244827A
CN104244827A CN201380020585.7A CN201380020585A CN104244827A CN 104244827 A CN104244827 A CN 104244827A CN 201380020585 A CN201380020585 A CN 201380020585A CN 104244827 A CN104244827 A CN 104244827A
Authority
CN
China
Prior art keywords
perfusion
radiography
tissue
interest
angiosomeses
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
CN201380020585.7A
Other languages
English (en)
Other versions
CN104244827B (zh
Inventor
R·卡尔米
G·S·卡夫里
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 CN104244827A publication Critical patent/CN104244827A/zh
Application granted granted Critical
Publication of CN104244827B publication Critical patent/CN104244827B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 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/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
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computerised tomographs
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. 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
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48—Diagnostic techniques
    • A61B6/486—Diagnostic techniques involving generating temporal series of image data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50—Clinical applications
    • A61B6/507—Clinical applications involving determination of haemodynamic parameters, e.g. perfusion CT
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus for radiation diagnosis, e.g. 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366—Perfusion imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • 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/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/0037—Performing a preliminary scan, e.g. a prescan for identifying a region of interest
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Hematology (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Pulmonology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Vascular Medicine (AREA)
  • Signal Processing (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Cardiology (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Nuclear Medicine (AREA)

Abstract

一种方法，包括获得被扫描对象的至少两个血管区域以及至少一种感兴趣组织的造影增强灌注成像数据，其中，所述至少两个血管区域一个在另一个的下游，所述至少一种感兴趣组织接收来自循环系统的血液。所述方法还包括基于图像数据来确定所述至少两个血管区域的造影剂峰值之间的血流时间差。所述方法还包括基于所述图像数据来确定所述感兴趣组织的绝对灌注。所述方法还包括基于所述绝对灌注和所述时间差来计算标准化灌注值。所述方法还包括显示所述标准化灌注值。

Description

灌注成像
技术领域
下文总体上涉及灌注成像并且结合对计算机断层摄影(CT)的具体应用加以描述；然而，下文也适用于其他成像模态，包括磁共振成像(MRI)、正电子发射断层摄影(PET)、单光子发射计算机断层摄影(SPECT)、超声(US)和/或其他成像模态。
背景技术
CT扫描器包括由旋转架支撑的X射线管。旋转架以及由此的X射线管绕检查区域旋转，并且X射线管发射贯穿检查区域和设置在其中的对象和/或目标的部分的辐射。辐射敏感探测器跨检查区域位于X射线管的对面，并且探测贯穿检查区域以及对象和/或目标的辐射。辐射敏感探测器生成指示所探测的辐射的投影数据。重建器重建投影数据并生成指示对象和/或目标的体积图像数据。图像处理器能够被用来处理体积图像数据并生成指示对象和/或目标的一幅或多幅图像。
CT灌注成像通过使用诸如碘造影药剂的造影剂对与器官和组织中的血流有关的特性进行成像来提供功能信息。典型的灌注研究要求在许多不同时间点对感兴趣体积进行重复成像，例如在相继的扫描之间具有数秒差异的3-30个重复扫描(例如，1到10秒每时间帧)。在常规灌注技术中，将一团造影药剂施予到患者的血管系统中，并且在涵盖该造影药剂团传输通过感兴趣区域中的组织的时段内收集感兴趣区域的图像。造影药剂的局部浓度随时间的变化(如能够从图像数据推测的)被用于分析生理参数。
灌注特别可用于研究具有例如由脑卒中、肿瘤、梗塞等引起的脑部、心脏或肝脏损伤的患者以及研究诸如肺和肾的额外器官的一般功能。近年来，已表明抗血管生成药剂不可以显著减小肿瘤体积，尤其在开始治疗不久以后，所以肿瘤大小变化的常规评估可能对处置响应不敏感或提供对处置响应的明显延迟的指示。微循环灌注成像提供额外的活体信息并且具有实际潜力成为用于评估多种类型的癌症和肿瘤的重要工具。在临床实践中，通常利用特殊分析算法对灌注图像序列进行定性检查或定量评估。
对于定量的绝对值，通常通过从主动脉组织时间衰减曲线(TAC)或从将血液引到组织的另一相关动脉导出的参数来将从TAC导出的参数归一化。该绝对灌注值(以1/时间或1/时间/密度为单位)可用于评估所检查组织的生理机能，但是其也可能受到诸如心输出量和全身血液体积的其他因素的影响从而可以导致临床诊断中的错误和不准确性。例如，具有特定微血管密度的组织内的灌注也能够通过简单地提高心输出量而得以增加。在大的患者中，还要求总的患者血液体积的校正，因为更高的心输出量将分布在整个更大血液体积中而不改变血液到给定组织的传递。文献已经建议使用标准化灌注值(SPV)作为度量来量化灌注。
一种公式化描述是：SPV＝(P×W)/(C/A)，其中，P是绝对灌注值，W是患者体重，C是总的造影剂剂量乘以到亨氏单位(HU)的校正转换，并且A是经再循环校正的动脉时间衰减曲线下的面积。比率C/A被假定为大约与患者心输出量成比例。遗憾的是，该方法具有限制和不准确性。例如，患者体重、注射剂量以及到HU的转换函数的外部信息可能不可得或不可靠。其次，比率W/(C/A)不是与全身灌注真的成比例(或相等)，因为患者体重不是适当的因子，这是由于骨骼和脂肪区域每周与灌注过程非常相关。另外，许多疾病和身体状况可以改变身体灌注。计算中的一些变化可能由于不准确的造影剂量化引起，例如，在由于变化的射束硬化而导致的不同的X射线管能量中能够获得不同的碘HU。
鉴于至少以上所述，存在对用于确定SPV(例如比以上讨论的公式化描述和/或其他公式化描述更准确可靠的SPV)的其他方法的未解决的需求。
发明内容
本文描述的各方面解决了以上提及的问题和其他问题。
在一个方面中，一种方法，包括获得被扫描对象的至少两个血管区域以及至少一种感兴趣组织的造影增强灌注成像数据，其中，所述至少两个血管区域一个(相对于造影药剂团流)在另一个的下游，所述至少一种感兴趣组织接收来自循环系统的血液。所述方法还包括基于图像数据来确定所述至少两个血管区域的造影剂峰值之间的血流时间差。所述方法还包括基于所述图像数据来确定所述感兴趣组织的绝对灌注。所述方法还包括基于所述时间差和所述绝对灌注来计算标准化灌注值。所述方法还包括显示所述标准化灌注值。
在另一方面中，一种造影数据处理器，包括造影剂水平确定器，所述造影剂水平确定器根据来自对象的扫描的造影增强灌注图像数据来确定至少两个血管区域和至少一种感兴趣组织中的造影增强的水平。所述造影数据处理器还包括参数确定器，所述参数确定器基于所述图像数据来确定所述至少两个血管区域的造影剂峰值之间的血流时间差并且基于所述图像数据来确定所述感兴趣组织的绝对灌注。所述造影数据处理器还包括度量确定器，所述度量确定器基于所述时间差和所述绝对灌注来计算标准化灌注值。
在另一方面中，一种成像系统，包括CT、MRI、US、PET或SPECT扫描器中的至少一个以及处理系统。所述处理系统包括：造影剂水平确定器，其根据来自对象的扫描的灌注图像数据来确定至少两个血管区域和至少一种感兴趣组织中的增强或活动水平；参数确定器，其基于所述图像数据来确定所述至少两个血管区域的增强或活动峰值之间的血流时间差并且基于所述图像数据来确定所述感兴趣组织的绝对灌注；以及度量确定器，其基于所述时间差和所述绝对灌注来计算标准化灌注值。
附图说明
本发明可以采取各种部件和部件的布置，以及各种步骤和步骤的安排的形式。附图仅出于图示优选实施例的目的，并且不得被解释为对本发明的限制。
图1示意性地图示了与造影数据处理器相连的范例成像系统；
图2图示了范例时间衰减曲线；
图3图示了能够在第一轴向扫描窗中一起扫描的范例血管；
图4图示了能够在第二轴向扫描窗中一起扫描的范例血管；
图5图示了针对血管和针对感兴趣组织的范例轴向扫描窗；
图6图示了用于确定不同时间衰减曲线中的造影剂峰值之间的时间差的范例；
图7图示了用于确定针对感兴趣组织的绝对灌注的范例；
图8图示了确定在针对血管区域和感兴趣组织的轴向扫描窗交叠的情况下的SPV的方法；
图9图示了确定在针对血管区域和感兴趣组织的轴向扫描窗不交叠的情况下的SPV的方法；
图10图示了使用实时图像数据定位对象以确定在针对血管区域和感兴趣组织的轴向扫描窗不交叠的情况下的SPV的方法。
具体实施方式
下文描述一种用于确定准确可靠的SPV的方法。来自CT、MRI、PET、SPET、US等的灌注成像数据能够被用来确定用于计算SPV的数据。然而，为简洁和清楚起见，下文描述一种使用CT灌注成像数据的方法。
首先参考图1，图示了一种诸如CT扫描器的成像系统100。成像系统100包括大体固定机架102和旋转机架104，旋转机架104可旋转地由固定机架102支撑并且关于z轴绕检查区域106旋转。诸如卧榻的对象支撑物108将目标或对象支撑在检查区域106中。
诸如X射线管的辐射源110可旋转地由旋转机架104支撑，随着旋转机架104旋转，并且发射贯穿检查区域106的多色辐射。辐射敏感探测器阵列112跨检查区域106在角弧度上对着辐射源110。辐射敏感探测器阵列112探测贯穿检查区域106的辐射并针对每个探测到的光子生成指示所述辐射的信号。
任选的注射器114被配置为例如针对诸如灌注扫描的造影增强成像流程来注射(一种或多种)造影剂。图示的注射器114由成像系统100控制，成像系统100可以触发或启动注射器114以施予造影剂。(一种或多种)造影剂或者能够由临床医师等手动施予。在手动施予(一种或多种)造影剂的情况下，能够省略注射器114。
重建器116重建投影，生成指示位于检查区域106中的对象或目标的被扫描部分的体积图像数据。通用计算系统或计算机服务器充当操作员控制台118。控制台118包括诸如监视器的人类可读输出设备和诸如键盘、鼠标等的输出设备。驻留在控制台118上的软件允许操作员经由图形用户接口(GUI)或以其他方式与扫描器100交互和/或操作扫描器100。
处理系统126包括造影剂水平确定器120、参数确定器122和/或度量确定器124。在另一实施例中，造影剂水平确定器120、参数确定器122和/或度量确定器124是控制台118的部分。在又一实施例中，造影剂水平确定器120、参数确定器122和/或度量确定器124是另一计算系统的部分，所述另一计算系统可以在系统100的本地或远离系统100(例如，位于另一房间、跨网络分布等)。
造影剂水平确定器120基于体积图像数据来确定在血管区域和/或感兴趣组织中摄取和洗出的造影剂的水平或浓度。参数确定器122基于摄取和洗出的造影剂的水平或浓度来确定一个或多个灌注参数。度量确定器124基于所确定的一个或多个灌注参数来确定一个或多个灌注度量。如以下更加详细地描述的，这包括基于两个不同血管区域中的血流之间的血流时间差和感兴趣组织中的绝对灌注来确定标准化灌注值(SPV)度量。
造影剂水平确定器120、参数确定器122和/或度量确定器124能够经由执行一个或多个计算机可读指令的一个或多个处理器来实施，所述一个或多个计算机可读指令被存储在诸如物理存储器和/或其他非暂态介质的计算机可读存储介质上。额外地或备选地，造影剂水平确定器120、参数确定器122和/或度量确定器124能够经由执行由载波、信号和/或暂态介质承载的一个或多个计算机可读指令的一个或多个处理器来实施。
下文更详细地描述了造影剂水平确定器120、参数确定器122和度量确定器124。一般地，为了确定SPV，使用造影增强灌注协议来扫描对象的第一血管区域和第二血管区域以及对象的第三区域(感兴趣组织)，所述第一血管区域和第二血管区域一个在另一个的下游，所述第三区域接收来自身体循环系统的任何分支的血液。能够经由针对三个不同区域的三个不同时间衰减曲线(TAC)202、204和206量化所得灌注图像数据，如图2所示。
在图2中，x轴210表示时间并且y轴208表示造影增强水平。TAC202表示根据时间的与第一血管区域相关的造影摄取和洗出，TAC 204表示根据时间的与第二血管区域相关的造影摄取和洗出，并且TAC 206表示根据时间的与第三区域或感兴趣组织相关的造影摄取和至少部分洗出。
在一个实例中，在相同扫描窗中捕获全部三个区域。在另一实例中，在相同扫描窗中捕获三个感兴趣区域中的两个，并且在不同扫描窗中捕获三个区域中的第三个。在又一实例中，全部三个区域在不同扫描窗中被捕获。图3示出了能够在相同轴向扫描窗306中捕获的两个血管区域肺动脉302和主动脉304的范例。图4示出了能够在相同轴向扫描窗406中捕获的两个血管区域主动脉402和肝门静脉404的范例。
图5示出了对应于针对灌注研究的感兴趣组织的非限制性范例的轴向扫描窗306和406以及其他轴向扫描窗502、504、506和508。
在图5中，轴向扫描窗502对应于感兴趣组织并在空间上与轴向扫描窗306和406分开。同样地，轴向扫描窗504对应于感兴趣组织并在空间上与轴向扫描窗306和406分开。
因此，肺动脉302和主动脉304(图3)能够在轴向扫描窗306中被捕获并且头部能够在轴向扫描窗502中被捕获(或者肾或前列腺能够在轴向扫描窗504中被捕获)。或者，主动脉402和肝门静脉404(图3)能够在扫描窗406中被捕获并且肾或前列腺能够在轴向扫描窗504中被捕获(或者头部能够在轴向扫描窗502中被捕获)。
轴向扫描窗506对应于感兴趣组织并与轴向扫描窗306交叠并在空间上与轴向扫描窗406分开，并且轴向扫描窗508对应于感兴趣组织并与轴向扫描窗406交叠并在空间上与轴向扫描窗306分开。
因此，肺动脉302、主动脉304和肺能够全部在轴向扫描窗506中被捕获，主动脉402、门静脉404和肝脏或胰腺能够在轴向扫描窗508中全部被捕获。在这些实例中，尽管可以根据是正在扫描第一血管区域或第二血管区域还是正在扫描感兴趣组织来使用不同扫描范围，但是对象支撑物108(图1)保持在相同位置处或者来回微动以覆盖轴向区域。
或者，肺动脉302和主动脉304能够在轴向扫描窗306中被捕获并且肺能够在轴向扫描窗506中被捕获，或者主动脉402和肝门静脉404两者能够在轴向扫描窗406中被捕获并且肝脏、前列腺、肾、胰腺等能够在轴向扫描窗504中被捕获。在这些实例中，对象支撑物108(图1)在轴向扫描窗之间移动对象。
在后一实例中，例如，对象支撑物108将对象定位以对应于扫描位置510以用于扫描第一血管和第二血管，并且扫描第一血管和第二血管直到达到特定时间点。之后，对象支撑物108将对象定位以对应于不与第一位置交叠的第二扫描位置512或514以用于扫描感兴趣组织，并且扫描感兴趣组织直到达到另一时间点。
血管区域和/或感兴趣组织的上述范例并非限制性的。例如，血管区域能够是心脏附近的肺动脉和主动脉区域，肝脏附近的主动脉和肝门静脉，肺动脉或主动脉之一中的主要的以及再循环的峰，主动脉以及上腔静脉，主动脉以及下腔静脉，颈区域中的颈动脉和邻近静脉，和/或血管区域的其他参考对。
参数确定器122基于TAC 202-206来确定一个或多个灌注参数。例如，参数确定器122能够针对TAC 202-206中的一个或多个确定从造影施予到峰值造影增强的时间(例如，达峰时间(TTP))。参数确定器122还能够基于TAC 202-206中的两个来确定时间差。所述时间差能够在两个TAC的平均渡越时间(即，“重心”或“第一动量”)之间，或者是表示最大梯度发生的两个时间点之间的差。任选地，伽马变量模型能够拟合至TAC。图6示出了确定TAC 202与204的峰值604与606之间的时间差600的范例。
参数确定器122还能够通过计算TAC 706的第一峰值702与再循环峰值704之间的时间差708的逆来根据TAC 706的第一峰值702和再循环峰值704确定全身灌注参数，其被示出在图7中。全身灌注值一般是在循环系统中全身血液体积循环一轮的速率。原则上，在能够测得第一通过曲线和第二通过(再循环)曲线的情况下能够计算出全身灌注值。
用于确定全身灌注的另一选项包括找到将真实全身灌注与测得的时间差相关的一般模型。该相关很可能被表示为针对所有患者的常数倍增因子(至少极近似)，因为例如肺系统血液体积与全身血液体积之间的比率可能几乎与体重、年龄、性别等无关。
用于确定全身灌注的另一选项是针对每个患者执行专属的校准测量。所述校准能够例如在开始治疗处置之前完成，所述治疗处置被规划为使用灌注成像随后进行。在校准流程期间，将测得全身灌注和诊断灌注协议的规划的时间差技术两者。通过计算全身灌注与在灌注协议中测得的时间差之间的相关，能够针对每个患者导出个体查找表(或甚至简单因子)。所述校准能够以若干不同心输出率完成以便实现更高准确度。
度量确定器124基于所确定的一个或多个灌注参数来确定一个或多个灌注度量。通过非限制性举例的方式，在一个实例中，度量确定器124基于等式1确定标准化灌注值(SPV)：
等式1：
SPV＝P×D，
其中，P表示绝对灌注(例如，在具体条件中测得的组织灌注等)并且D表示参考血流区域之间的流时间差(例如，差600、差708等)。
假设一个患者的心率在随后检查中更高并且绝对灌注仅因为更高心率而不是因为所检查的组织中的任何变化而变得更高(甚至在主动脉曲线的通常的归一化之后)，在这种情况下测得的时间差D将由于更快血流而变得更小。因此，标准化灌注值将很可能不会如其应该的那样相对于早先灌注研究发生变化，因为在检查组织中没有任何变化。
在等式2中示出了另一非限制性范例度量：
等式2：
SPV＝P×F(D)，
其中，F是从先验信息或预校准导出的模型函数或查找表，其将测得的时间差与确定全身灌注的另一方法相比较。
能够以类似于呈现常规灌注图的方式经由监视器视觉呈现所计算的SPV。另外，能够视觉呈现分析图和/或数值结果，例如在图2、6和7中示出的信息。
接下来预见变型。
作为一个选项，没有注射造影剂的一个或多个预扫描，或者在灌注成像之前采取的利用造影剂的测试团扫描，能够有助于实现更高准确度。在一个实例中，根据第一血管区域和第二血管区域确定的流时间差其本身可以在在检查范围的诊断灌注扫描之前或之后通过测试团灌注扫描获得。
在另一实例中，如果灌注扫描的注射的造影剂团时间短，可以甚至无需测量TAC的开始部分而准确地计算出组织灌注。在这种情况下，有足够的时间来测量第一血管区域和第二血管区域并且只有在那之后移动到检查区。如果规划的造影注射不够短，则能够使用(没有造影剂的)预扫描以用于测量基线组织增强。
在另一实例中，可以在诊断灌注扫描之前通过测试团扫描来确定用于将第一血管区域和第二血管区域移动到感兴趣组织的规划的时间点，所述测试团扫描使得能够估计具体患者血液循环的相关时间常数。
图8图示了基于针对血管区域和感兴趣组织的轴向扫描窗交叠的造影增强灌注扫描来确定SPV的方法。
应当理解，动作的顺序不是限制性的。因此，本文预见其他顺序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外动作。
在802，规划灌注扫描。这包括规划两个不同血管区域的一个或多个扫描，以及感兴趣组织的扫描，其中，所述两个不同血管区域一个在另一个的下游。
在804，任选地，执行无造影预扫描。所述无造影预扫描根据灌注计算方法是任选的。例如，如果使用最大斜率模型，则可以省略预扫描。
在806，对象支撑物定位对象以用于灌注扫描。
在808，基于计划来执行灌注扫描的第一部分，包括在第一持续时间内扫描对应于血管区域的轴向扫描窗。
在810，基于计划来执行灌注扫描的第二部分，包括在第二持续时间内扫描对应于感兴趣组织的轴向扫描窗。
在812，重建针对灌注扫描的第一部分和第二部分采集的数据。
在814，根据图像数据针对血管确定造影剂流峰值之间的时间差。
在816，基于来自图像数据的造影剂流来确定针对感兴趣组织的绝对灌注。
在818，基于所述差和所述绝对灌注例如如本文所描述的使用等式1或2和/或以其他方式来计算SPV。
图9图示了基于针对血管区域和感兴趣组织的轴向扫描窗不交叠的造影增强灌注扫描来确定SPV的方法。
应当理解，动作的顺序不是限制性的。因此，本文预见其他顺序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外动作。
在902，规划灌注扫描。这包括规划两个不同血管区域的一个或多个扫描，以及感兴趣组织的扫描，其中，所述两个不同血管区域一个在另一个的下游。
在904，任选地，执行无造影预扫描，如在图8的动作804中所描述的。
在906，对象支撑物将对象定位在对应于血管区域的轴向扫描范围的第一位置处。
在908，基于计划来执行灌注扫描的第一部分，包括在第一持续时间内扫描血管区域的轴向扫描范围。
在910，在第一持续时间后，对象支撑物将对象定位在对应于感兴趣组织的轴向扫描范围的第二不同位置处。
在912，基于计划来执行灌注扫描的第二部分，包括在第二持续时间内扫描感兴趣组织的轴向扫描范围。
在914，重建针对灌注扫描的第一部分和第二部分采集的数据。
在916，确定到达造影剂峰值的时间之间的差。
在918，基于对应图像数据针对感兴趣组织来确定绝对灌注。
在920，基于所述差和所述绝对灌注例如如本文所描述的使用等式1或2和/或以其他方式来计算SPV。
图10图示了使用实时图像数据定位对象，基于针对血管区域和感兴趣组织的轴向扫描窗不交叠的造影增强灌注扫描来确定SPV的方法。
应当理解，动作的顺序不是限制性的。因此，本文预见其他顺序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外动作。
在1002，选择灌注扫描协议。
在1004，执行低剂量预扫描。
在1006，从所述预扫描识别参考血管区域。
在1008，创建灌注扫描计划，识别要扫描的血管区域和感兴趣组织。
在1010，定位对象以用于扫描血管。
在1012，开始执行没有任何造影的灌注扫描并实时重建图像。
在1014，向对象施予造影。
在1016，利用造影执行灌注扫描。
在1018，实时重建图像。
在1020，分析所采集的数据以确定是否已经达到血管中的峰值造影增强。
如果还没有达到峰值造影增强，则重复动作1016-1020。
如果已经达到峰值造影增强，则在1022将对象移动到第二不同位置以用于扫描感兴趣组织。
在1024，针对感兴趣组织重复动作1016-1020。
在1026，基于所采集的数据例如如本文所描述的使用等式1或2和/或以其他方式来计算SPV。
可以通过编码于或嵌入于计算机可读存储介质上的计算机可读指令来实施以上描述的，所述计算机可读指令在由(一个或多个)计算机处理器执行时令所述(一个或多个)处理器执行所描述的动作。额外地或备选地，计算机可读指令中的至少一个由信号、载波或其他暂态介质承载。
尽管以上结合CT加以描述，但来自结合MRI、PET、SPECT等成像系统执行的造影灌注扫描的成像数据也能够被用来生成用于计算SPV的数据。这样的扫描的范例包括利用钆或氧化铁造影药剂的MRI、利用若干类型的辐射性示踪剂的PET和SPECT以及利用微泡造影药剂的超声。在动物临床前成像中，利用荧光药剂的光学断层摄影同样适用。
已经参考优选实施例描述了本发明。他人在阅读并理解前述详细说明后可以进行修改和变化。本发明旨在被解释为包括所有这样的修改和变化，只要它们落在权利要求书或其等价要件的范围内。

Claims (22)

1.一种方法，包括：
获得被扫描对象的至少两个血管区域以及至少一种感兴趣组织的造影增强灌注成像数据，其中，所述至少两个血管区域一个在另一个的下游，所述至少一种感兴趣组织接收来自循环系统的血液；
基于图像数据来确定所述至少两个血管区域的造影剂峰值之间的血流时间差；
基于所述图像数据来确定所述感兴趣组织的绝对灌注；
基于所述绝对灌注和所述时间差来计算标准化灌注值；并且
显示所述标准化灌注值。
2.如权利要求1所述的方法，还包括：
生成针对所述至少两个血管区域中的第一个的第一造影流行为；
基于所述第一行为来确定针对所述至少两个血管区域中的所述第一个的到达峰值造影增强的第一时间；
生成针对所述至少两个血管区域中的第二个的第二造影流行为；
基于所述第二行为来确定针对所述至少两个血管区域中的所述第二个的到达峰值造影增强的第二时间；
基于到达峰值造影增强的所述第一时间和到达峰值造影增强的所述第二时间来确定所述时间差。
3.如权利要求2所述的方法，还包括：
基于第一造影增强的平均渡越时间值和第二造影增强的平均渡越时间值来确定所述时间差。
4.如权利要求2所述的方法，还包括：
基于第一造影增强的最大梯度发生和第二造影增强的最大梯度发生来确定所述时间差。
5.如权利要求1到4中的任一项所述的方法，还包括：
生成针对所述感兴趣组织的造影流行为；
基于所述行为来确定第一峰值；
基于所述行为来确定最大梯度；
确定所述两个血管区域中的至少一个的曲线下面积和峰值增强；并且
基于所确定的参数来计算绝对组织灌注。
6.如权利要求1到5中的任一项所述的方法，还包括：
将所述标准化灌注值计算为所述绝对灌注和所述时间差的积。
7.如权利要求1到5中的任一项所述的方法，还包括：
将所述标准化灌注值计算为所述时间差的函数和所述绝对灌注的积。
8.如权利要求7所述的方法，其中，所述函数将所计算的各峰值之间的时间差与预定的全身灌注或归一化的血流指标中的至少一个相比较。
9.如权利要求1到8中的任一项所述的方法，还包括：
在相同轴向扫描窗内采集对应于所述至少两个血管区域的造影增强灌注成像数据。
10.如权利要求9所述的方法，还包括：
在所述相同轴向扫描窗内采集对应于所述感兴趣组织的造影增强灌注成像数据。
11.如权利要求9所述的方法，还包括：
在不同轴向扫描窗内采集对应于所述感兴趣组织的造影增强灌注成像数据。
12.如权利要求11所述的方法，还包括：
将所述对象在第一位置与第二不同位置之间移动，其中，在所述第一位置处采集对应于所述感兴趣组织的所述造影增强灌注成像数据，在所述第二不同位置处采集对应于所述感兴趣组织的所述造影增强灌注成像数据。
13.如权利要求12所述的方法，还包括：
在预定时间段过去之后将所述对象从所述第一位置移动到所述第二位置。
14.如权利要求12所述的方法，还包括：
基于实时重建的图像将所述对象从所述第一位置移动到所述第二位置，所述实时重建的图像视觉指示所述至少两个血管区域中的造影增强何时达到峰值。
15.一种处理系统(126)，包括：
造影剂水平确定器(120)，其根据来自对象的扫描的造影增强灌注图像数据来确定至少两个血管区域和至少一种感兴趣组织中的造影增强的水平；
参数确定器(122)，其基于所述图像数据来确定所述至少两个血管区域的造影剂峰值之间的血流时间差并且基于所述图像数据来确定所述感兴趣组织的绝对灌注；以及
度量确定器(124)，其基于所述时间差和所述绝对灌注来计算标准化灌注值。
16.如权利要求15所述的处理系统，其中，所述度量确定器将所述标准化灌注值计算为所述绝对灌注和所述时间差的积。
17.如权利要求15所述的处理系统，其中，所述度量确定器将所述标准化灌注值计算为所述时间差的函数和所述绝对灌注的积。
18.如权利要求17所述的处理系统，其中，所述函数将所计算的各峰值之间的时间差与预定的全身灌注或归一化的血流指标中的至少一个相比较。
19.一种成像系统，包括：
CT、MRI、US、PET或SPECT扫描器中的至少一个；以及
处理系统(126)，其包括：
造影剂水平确定器(120)，其根据来自对象的扫描的灌注图像数据来确定至少两个血管区域和至少一种感兴趣组织中的增强或活动水平；
参数确定器(122)，其基于所述图像数据来确定所述至少两个血管区域的增强或活动峰值之间的血流时间差并且基于所述图像数据来确定所述感兴趣组织的绝对灌注；以及
度量确定器(124)，其基于所述时间差和所述绝对灌注来计算标准化灌注值。
20.如权利要求19所述的成像系统，还包括：
对象支撑物(108)，其支撑对象以用于扫描；以及
控制台(118)，其控制所述对象支撑物的移动，其中，所述控制台将所述对象支撑物定位在第一位置处以用于扫描所述血管区域并定位在第二不同位置处以用于扫描所述感兴趣组织。
21.如权利要求20所述的成像系统，其中，所述对象支撑物在预定时间段过去之后将所述对象从所述第一位置移动到所述第二不同位置。
22.如权利要求20所述的成像系统，其中，所述对象支撑物基于实时重建的图像将所述对象从所述第一位置移动到所述第二不同位置，所述实时重建的图像视觉指示所述至少两个血管区域中的增强何时达到峰值。
CN201380020585.7A 2012-04-17 2013-04-10 灌注成像 Active CN104244827B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201261625143P 2012-04-17 2012-04-17
US61/625,143 2012-04-17
PCT/IB2013/052848 WO2013156901A1 (en) 2012-04-17 2013-04-10 Perfusion imaging

Publications (2)

Publication Number Publication Date
CN104244827A true CN104244827A (zh) 2014-12-24
CN104244827B CN104244827B (zh) 2017-11-14

Family

ID=48539316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380020585.7A Active CN104244827B (zh) 2012-04-17 2013-04-10 灌注成像

Country Status (6)

Country Link
US (1) US10617312B2 (zh)
EP (1) EP2838430B1 (zh)
JP (1) JP5820549B2 (zh)
CN (1) CN104244827B (zh)
RU (1) RU2014145862A (zh)
WO (1) WO2013156901A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110279417A (zh) * 2019-06-25 2019-09-27 沈阳东软智能医疗科技研究院有限公司 识别主动脉血管的方法、装置及设备
CN110461235A (zh) * 2017-03-24 2019-11-15 皇家飞利浦有限公司 心肌ct灌注图像合成

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR112015009354A2 (pt) * 2012-10-31 2017-07-04 Koninklijke Philips Nv método e sistema de processamento
JP6343004B2 (ja) 2014-07-15 2018-06-13 富士フイルムＲｉファーマ株式会社 コンピュータプログラム、画像処理装置及び方法
DE102015208202B3 (de) * 2015-05-04 2016-09-22 Bayer Pharma Aktiengesellschaft Verfahren zum automatischen Ermitteln eines Kontrastmittelprotokolls sowie Bildgebungsverfahren, Kontrastmittelprotokoll-Ermittlungseinrichtung, bildgebende medizinische Einrichtung, Computerprogrammprodukt und computerlesbares Medium
CN104873192A (zh) * 2015-06-10 2015-09-02 上海大学 微型脑功能监测与成像的装置和方法
RU2764500C2 (ru) * 2017-01-04 2022-01-17 Конинклейке Филипс Н.В. Система и способ параметрической ультразвуковой визуализации с контрастированием и привязкой ко времени

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6292683B1 (en) * 1999-05-18 2001-09-18 General Electric Company Method and apparatus for tracking motion in MR images
JP2004057812A (ja) * 2002-06-03 2004-02-26 Ge Medical Systems Global Technology Co Llc Ｃｂｆ定量解析方法および装置
CN101506842A (zh) * 2006-05-15 2009-08-12 卑尔根大学研究基金会 灌注成像的方法
CN101674774A (zh) * 2007-05-08 2010-03-17 皇家飞利浦电子股份有限公司 使用低剂量钙化积分扫描的冠状动脉选择性钙分配
CN102123665A (zh) * 2008-08-13 2011-07-13 皇家飞利浦电子股份有限公司 冠状血管和心肌灌注信息的动态可视化
CN102165455A (zh) * 2008-09-30 2011-08-24 皇家飞利浦电子股份有限公司 灌注成像

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6015384A (en) * 1998-08-31 2000-01-18 Acuson Corporation Ultrasonic system and method for tissue viability imaging
US6645147B1 (en) * 1998-11-25 2003-11-11 Acuson Corporation Diagnostic medical ultrasound image and system for contrast agent imaging
JP3920140B2 (ja) 2002-05-13 2007-05-30 株式会社東芝 Ｍｒｉ装置及びフロー定量化装置
US20040077941A1 (en) * 2002-10-21 2004-04-22 Ge Medical Systems Global Technology Company, Llc Method and system for image improvement with ECG gating and dose reduction in CT imaging
US20050283070A1 (en) * 2004-06-21 2005-12-22 Celina Imielinska Systems and methods for qualifying symmetry to evaluate medical images
US7394889B2 (en) * 2006-05-18 2008-07-01 Varian Medical Systems Technologies, Inc. Contrast-enhanced cone beam X-ray imaging, evaluation, monitoring and treatment delivery
US8099149B2 (en) 2006-11-20 2012-01-17 Northwestern University MRI method for quantification of cerebral perfusion
FR2913874A1 (fr) 2007-03-23 2008-09-26 Roussy Inst Gustave Determinations conjointes de parametres quantitatifs de perfusion et de permeabilite par irm d'un tissu, notamment d'une tumeur cerebrale.
DE102007028901B4 (de) * 2007-06-22 2010-07-22 Siemens Ag Verfahren und Vorrichtung zur automatischen Bestimmung von Perfusion mittels einer Magnetresonanzanlage
EP2189112A1 (en) 2008-11-24 2010-05-26 Bracco Research S.A. Real-time perfusion imaging and quantification
US20110245687A1 (en) * 2008-12-05 2011-10-06 Koninklijke Philips Electronics N.V. Optical detection method and device for optical detection of the condition of joints
FR2952452B1 (fr) * 2009-11-09 2016-01-08 Olea Medical Procede pour estimer des parametres hemodynamiques par estimation conjointe des parametres d'un modele global de perfusion contraint
DE102010062030A1 (de) * 2010-11-26 2012-05-31 Siemens Aktiengesellschaft Verfahren zur Berechnung von Perfusionsdaten aus 2-D-Angiographiedaten oder DSA-Sequenzen
JP5811491B2 (ja) * 2011-04-12 2015-11-11 株式会社ニコン 顕微鏡及びそのプログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6292683B1 (en) * 1999-05-18 2001-09-18 General Electric Company Method and apparatus for tracking motion in MR images
JP2004057812A (ja) * 2002-06-03 2004-02-26 Ge Medical Systems Global Technology Co Llc Ｃｂｆ定量解析方法および装置
CN101506842A (zh) * 2006-05-15 2009-08-12 卑尔根大学研究基金会 灌注成像的方法
CN101674774A (zh) * 2007-05-08 2010-03-17 皇家飞利浦电子股份有限公司 使用低剂量钙化积分扫描的冠状动脉选择性钙分配
CN102123665A (zh) * 2008-08-13 2011-07-13 皇家飞利浦电子股份有限公司 冠状血管和心肌灌注信息的动态可视化
CN102165455A (zh) * 2008-09-30 2011-08-24 皇家飞利浦电子股份有限公司 灌注成像

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KENNETH A. MILES等: "Standardized Perfusion Value: Universal CT Contrast Enhancement Scale that Correlates with FDG PET in Lung Nodules", 《RADIOLOGY》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110461235A (zh) * 2017-03-24 2019-11-15 皇家飞利浦有限公司 心肌ct灌注图像合成
CN110461235B (zh) * 2017-03-24 2023-12-01 皇家飞利浦有限公司 心肌ct灌注图像合成
CN110279417A (zh) * 2019-06-25 2019-09-27 沈阳东软智能医疗科技研究院有限公司 识别主动脉血管的方法、装置及设备

Also Published As

Publication number Publication date
JP5820549B2 (ja) 2015-11-24
JP2015519098A (ja) 2015-07-09
RU2014145862A (ru) 2016-06-10
US20150080707A1 (en) 2015-03-19
CN104244827B (zh) 2017-11-14
EP2838430A1 (en) 2015-02-25
US10617312B2 (en) 2020-04-14
EP2838430B1 (en) 2018-06-13
WO2013156901A1 (en) 2013-10-24

Similar Documents

Publication Publication Date Title
Miles Perfusion CT for the assessment of tumour vascularity: which protocol?
EP2399238B1 (en) Functional imaging
CN104244827A (zh) 灌注成像
JP6209520B2 (ja) 灌流イメージング
CN102740769B (zh) 功能成像
CN104519794B (zh) 基于动态对比增强成像的渗透率度量
Einstein Radiation risk from coronary artery disease imaging: how do different diagnostic tests compare?
CN105120738A (zh) 狭窄治疗规划
Hu et al. Design and implementation of automated clinical whole body parametric PET with continuous bed motion
Plathow et al. Quantification of lung tumor volume and rotation at 3D dynamic parallel MR imaging with view sharing: preliminary results
US9626759B2 (en) Method for the automatic recognition of anatomical structures in images obtained by positron emission tomography, system and computer program for performing said method
US20170014069A1 (en) Spectral imaging based fluid volume map
Ramirez-Giraldo et al. Evaluation of strategies to reduce radiation dose in perfusion CT imaging using a reproducible biologic phantom
US20140365189A1 (en) Simplified method for robust estimation of parameter values
Gremse et al. Qualitative and quantitative data analysis
Cao et al. Response to letter regarding article:“developing dce-ct to quantify intra-tumor heterogeneity in breast tumors with differing angiogenic phenotype”
EP2501287B1 (en) Functional imaging
WO2020242331A2 (en) Coincidence detection system for measuring arterial blood time-activity curves and methods of using same
Toufique et al. Simulation study of a coincidence detection system for non-invasive determination of arterial blood time-activity curve measurements
EP2601886B1 (en) Compartmental analyzer, compartmental analysis method, program, and recording medium
US20240074721A1 (en) Imaging method for diagnosing cardiovascular disease
JP2013510621A5 (zh)
Tognolini et al. Body tumor CT perfusion protocols: optimization of acquisition scan parameters in a rat tumor model
Caldeira et al. MR-based Corrections for Quantitative PET Image
Liang et al. Development of contrast-enhanced rodent imaging using functional CT

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant