CN107580470A - 瞬时流量储备‑计算机断层摄影 - Google Patents

瞬时流量储备‑计算机断层摄影 Download PDF

Info

Publication number
CN107580470A
CN107580470A CN201680005857.XA CN201680005857A CN107580470A CN 107580470 A CN107580470 A CN 107580470A CN 201680005857 A CN201680005857 A CN 201680005857A CN 107580470 A CN107580470 A CN 107580470A
Authority
CN
China
Prior art keywords
resistance
stenosis
vessel
branch
metric
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
CN201680005857.XA
Other languages
English (en)
Other versions
CN107580470B (zh
Inventor
L·戈申
Y·拉马什
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 CN107580470A publication Critical patent/CN107580470A/zh
Application granted granted Critical
Publication of CN107580470B publication Critical patent/CN107580470B/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
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • 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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording 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/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • 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/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • 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
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion
    • 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/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders

Landscapes

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

Abstract

一种方法包括无创地确定瞬时无波比率度量，包括接收由成像系统生成的电子格式图像数据。图像数据包括具有表示具有狭窄的血管的强度的体素。所述方法还包括根据图像数据计算血管的出口的外围阻力。所述方法还包括基于边界条件的集合和计算流体动力学算法来计算血管入口的入口和血管的出口之间的狭窄的狭窄阻力。所述方法还包括计算瞬时无波比率度量。所述度量是基于狭窄阻力的数值，并且生成指示计算的瞬时无波比率度量的信号。

Description

瞬时流量储备-计算机断层摄影
技术领域
以下总体上涉及基于图像数据估计血管中的狭窄的瞬时无波比率或瞬时流量储备(iFR)度量，并且尤其应用于计算机断层摄影(CT)来描述。
背景技术
流量储备分数(FFR)是用于指示冠状动脉狭窄的功能意义的度量。其已在介入心脏病学流程的情况下被用于测量跨冠状动脉狭窄的压力差，以确定狭窄是否阻碍血液流动，并且因此阻碍针对心肌的氧气递送，这可能导致心肌缺血。FFR是狭窄后的压力Pd相对于狭窄前的压力Pin的比率，或Pd/Pin，其是狭窄病变远端的最大血流量与相同血管中的正常最大流量的比率。在冠状动脉导管插入期间执行FFR。为此，使用护套和导丝将导管插入到股动脉或桡动脉中。FFR使用线的尖端上的小传感器来测量压力、温度和流量，以确定病变的确切严重性。这是在最大血流(充血)期间完成的。执行压力线的回拉，并且记录跨血管的压力。
瞬时无波比率或瞬时流量储备(iFR)是能够用于指示冠状动脉狭窄的功能意义的另一种测量。iFR被定义为舒张期的特定时段(称为“无波期”)上的Pd与Pin的比率。在该无波期期间，影响冠状动脉血流的竞争力(波)是静止的，意味着压力和流量以类似的方式变化，这使得它们的比率(即阻力)与心脏周期的其余部分相比几乎恒定。iFR也使用侵入式冠状动脉压力线在心脏导管插入期间执行，侵入式冠状动脉压力线被放置在要评估的冠状动脉中。FFR和iFR两者为零(0)至1(1)的范围中的值，其中，较高值指示非显著性或较不显著性狭窄，并且较低值指示更显著的狭窄。遗憾的是，FFR和iFR两者都是侵入式流程，并且因此易受范围从感染到死亡的并发症影响。
FFR-CT是一种无创流程，其用于使用从无创心脏CT研究中提取的冠状动脉几何结构上的计算流体动力学(CFD)模拟来提供FFR的估计。遗憾的是，基于CFD的FFR估计需要具有在心脏收缩和松弛期间发生的许多过渡效应的复杂的动态模拟以及具有可能导致高估计误差的许多不确定性的复杂模型。
发明内容
本申请的方面解决上述问题和其他问题。
根据一个方面，一种方法包括用于无创地确定瞬时无波比率度量，包括接收由成像系统生成的电子格式图像数据。图像数据包括具有表示具有狭窄的血管的强度的体素。所述方法还包括根据图像数据计算血管的出口的外围阻力。所述方法还包括基于边界条件的集合和计算流体动力学算法来计算血管入口的入口与血管的出口之间的狭窄的狭窄阻力。所述方法还包括计算瞬时无波比率度量。所述度量是数值，基于狭窄阻力，并且生成指示计算的瞬时无波比率度量的信号。
在另一方面中，一种用于无创地确定瞬时无波比率度量的系统，参数确定器被配置为根据图像数据计算血管的出口的外围阻力。系统还包括狭窄阻力确定器，所述狭窄阻力确定器被配置为基于边界条件的集合和计算流体动力学算法的结果来计算血管入口的入口与血管的出口之间的血管的狭窄的狭窄阻力。所述系统还包括iFR确定器，所述iFR确定器被配置为基于所述狭窄阻力计算瞬时无波比率度量，其中，所述度量是数值，并且生成指示所计算的瞬时无波比率度量的信号。
在另一方面中，一种计算机可读存储介质被编码有计算机可读指令，所述计算机可读指令在由计算系统的处理器运行时，令所述处理器：仅根据无创地确定的参数计算iFR度量，所述无创地确定的参数是基于计算机断层摄影图像数据和计算流体动力学算法无创地确定的。
在阅读和理解了以下详细描述后，本领域普通技术人员将会理解本发明的另外的方面。
附图说明
本发明可以采取各种部件和各部件布置的形式，并且可以采取各种步骤和各步骤安排的形式。附图仅用于说明优选实施例，而不应被解释为对本发明的限制。
图1示意性地图示了与数据分析器连接的成像系统。
图2图示了图1的数据分析器的范例。
图3图示了用于确定iFR度量的范例方法。
具体实施方式
以下描述了基于根据心脏CT扫描估计的特定冠状动脉几何结构估计iFR度量的无创方法。图1示意性地示出了诸如CT扫描器的成像系统100。成像系统100包括大致固定的机架102和旋转机架104，旋转机架104由固定机架102可旋转地支撑并关于z轴围绕检查区域106旋转。
诸如X射线管的辐射源110由旋转机架104可旋转地支撑，与旋转机架104一起旋转，并且发射穿过检查区域106的辐射。辐射敏感探测器阵列112跨检查区域106与辐射源110相对围成一角弧。辐射敏感探测器阵列112探测穿过检查区域106的辐射，并且生成投影数据，或者针对每个探测到的光子指示其的信号。
重建器114重建投影，生成指示位于检查区域106中的对象或物体的被扫描部分的体积图像数据。操作者控制台116包括人类可读输出设备(例如监测器)和输入设备(诸如键盘、鼠标等)。驻留在控制台116上的软件允许操作者经由图形用户界面(GUI)或其它方式与扫描器100交互和/或操作扫描器100。
诸如卧榻的对象支撑体119支撑检查区域106中的物体或对象。数据存储库118存储包括体积图像数据和/或投影数据的电子格式数据。数据存储库的范例包括影像归档和通信系统(PACS)、放射信息系统(RIS)、医院信息系统(HIS)、电子病历(EMR)、数据库、服务器等。
数据分析器120被配置为处理表示包括狭窄的(一个或多个)感兴趣血管(例如，冠状动脉)的图像数据，其中，在舒张期的无波期期间采集图像数据。如下面更详细描述的，这样的处理包括使用CFD模拟估计血管出口处的外围阻力并且基于此来确定iFR度量。在舒张期中，冠状动脉压力和速度两者都单调减小。舒张无波期处的时间压力和速度可以使用线性函数来近似。由压力和速度之间的比率定义的每个血管出口处的无波外围阻力近似恒定并且可以被估计。在舒张期的无波期期间的CFD模拟提供了狭窄的功能意义的快速估计。所得到的基于CT的iFR(iFR-CT)度量导致狭窄的功能意义的更鲁棒和无创的估计。
可以利用一个或多个计算系统的一个或多个处理器(例如，微处理器、中央处理单元等)来实施数据分析器120，所述计算系统运行被存储在一个或多个非暂态计算机可读存储介质(诸如物理存储器)中的一个或多个计算机可读指令。由一个或多个处理器处理的至少一个指令可以额外地或备选地由载波、信号和/或其他暂态介质承载。数据分析器120可以是控制台116的部分、成像系统100外部的计算系统，跨计算系统分布，为其组合等。
图2图示了数据分析器120的范例。
感兴趣组织(TOI)识别器202获得包括具有带表示感兴趣组织(TOI)204的强度的体素的图像数据作为输入，并且识别图像数据中的感兴趣组织204。感兴趣组织204可以通过指示用户选择的感兴趣组织、缺省的感兴趣组织等的信号来预先确定或识别。TOI识别器202可采用自动和/或手动方法来识别感兴趣组织。感兴趣组织的范例是管状组织，例如具有狭窄的血管。然而，感兴趣组织可以是其他组织。
TOI几何结构提取器206从所识别的感兴趣组织提取几何信息。TOI识别器202可采用自动和/或手动方法来提取几何结构信息。通过范例，提取可以包括采用具有主动轮廓和针对冠状动脉血管调节的水平集的分割(其中，感兴趣组织是冠状动脉血管)，之后是任选的额外的手动编辑以产生高质量分割。根据该提取和/或其他提取，可以确定冠状动脉血管的感兴趣几何结构(例如，血管直径)。
参数确定器208基于所提取的感兴趣几何结构来确定至少一个参数。图示的参数确定器208被配置为估计血管出口处的外围阻力。为此，参数确定器208可以利用在将血液分布到组织中的生理规则，并考虑心脏CT扫描的特性和限制。例如，由于冠状动脉的远端部分可能处于切断的分辨率之下，所以可以发现代表性的切割点，使得出口面积与供应到由该动脉馈送的区域的血液的分数成比例。
切割点选择的范例规则包括：狭窄后；最后分支不久之后；表示血管的面积；并且针对不同的分割器是不变的。本文还考虑了其他和/或不同的、更多或更少等的规则。可以关于主要冠状动脉分支基于比例定律估计特定出口的无波外围阻力，如等式1所示：
其中，Ri表示感兴趣出口的无波外围阻力，R主要表示主要冠状动脉分支的阻力，Di表示感兴趣出口处的直径，并且D主要表示主要冠状动脉分支处的直径。可以通过机器学习和/或其他方法来估计R主要。
狭窄阻力确定器210确定图像数据中的狭窄的阻力。在一个范例中，狭窄阻力确定器210基于边界条件212和等式2确定每个入口和每个出口之间的狭窄的阻力r：
其中，P入口是边界条件并且表示入口处的压力，v出口是边界条件并且表示出口处的速度，P出口表示出口处的压力，其基于CFD()来估计，CFD()是计算流体动力学算法，该算法采取P入口和v出口作为输入。在一个实例中，针对出口速度和入口压力的初始边界条件分别被设置为正常平均速度和压力值。(例如，出口速度＝40cm/s，并且入口压力＝100mmHg)，出口压力的边界条件可以被设置为零梯度，并且通过CFD获得值。
CFD处理器214例如使用偏微分方程来执行计算流体动力学(CFD)模拟。通常，CFD是使用数值方法和/或算法来求解和分析涉及流体流动的问题的流体力学方法。CFD处理器214利用由边界条件定义的表面来执行计算。在一个实例中，输出包括所有点处的压力的完整体积信息，例如出口处的压力。在舒张无波期期间的CFD模拟不如其他时段期间复杂。
使用等式2，可以迭代地估计狭窄的阻力r。下面图示了n(n＝1、…N)次迭代上的范例迭代方法。边界条件P入口和v出口被初始化。在每次迭代期间，利用当前边界条件执行CFD模拟，生成计算经更新的估计阻力r(n+1)，如等式3中所示：
逻辑216检查是否已经满足停止标准。在一个实例中，停止标准基于最小化当前阻力和先前阻力之间的平方差的和(sum of squares difference，SSD)，如等式4所示：
其中，Rp是目标阻力，其使用数据库上的统计数据和机器学习技术来学习，Rc表示出口处的实际(当前)阻力(优化期间)并且使用出口的速度和出口的压力在每次CFD迭代处被计算。针对给定的出口：Rc＝P出口/V出口，其中，P出口从CFD迭代获得。等式4中的停止标准是Rc随着Rp的收敛。
如果停止标准不满足，则边界条件更新器218基于所计算的狭窄阻力r(n)更新边界条件v出口，如等式5所示：
并再次执行等式2、3和4。在等式5中，Ri指代出口i处的外围阻力，r指代入口和每个出口之间的狭窄的阻力。如果没有狭窄，则r～＝0。
如果满足停止标准，则iFR确定器220基于所计算的狭窄阻力来确定iFR。通常，在阻力Ri达到其目标值之后估计iFR度量。可以如等式6所示估计iFR度量：
其中，Pa是入口压力(或主动脉压力)，并且Pd是所讨论的狭窄的远端的模拟压力。
上述允许无创、鲁棒的、准确和快速的模拟。此外，在对FFR-CT建模中的不确定性水平远远大于估计iFR-CT中，并且模拟iFR(在安静阶段处执行)比基于CFD的FFR估计简单得多。这样一来，iFR-CT导致狭窄的功能意义的更鲁棒的无创估计。此外，使用本文描述的接近，iFR估计被按数量级加速，例如从五(5)分钟到十(10)秒。
图3图示了用于确定iFR度量的范例方法。
应当意识到，以上动作的顺序不是限制性的。这样一来，本文预期了其他顺序。此外，可以省略一个或多个动作，和/或可以包括一个或多个额外的动作。
在302处，扫描对象的区域。
在304处，在来自扫描的图像数据中分割诸如具有狭窄的冠状动脉的感兴趣组织。
在306处，从感兴趣组织提取诸如血管直径、血管半径等的几何信息。
在308处，基于几何信息确定冠状动脉血管出口处的外围阻力。
在310处，入口压力和出口速度边界条件被初始化。
在312处，使用入口压力和出口速度边界条件执行CFD，从而估计出口压力。
在314处，基于入口压力、出口速度和估计的出口压力计算狭窄阻力。
在316处，检查停止标准。
如果不满足停止标准，则在318处，基于所计算的狭窄阻力来更新出口速度，并且重复动作312-316。
如果满足停止标准，则在320处估计iFR度量。
iFR度量可以被保存到非暂态存储器，被视觉呈现，被处理等。如本文所讨论的，iFR表示冠状动脉狭窄的功能意义，例如狭窄将阻碍针对心肌的氧递送的可能性。
以上可以通过编码或嵌入在计算机可读存储介质上的计算机可读指令来实施，所述计算机可读指令在由(一个或多个)计算机处理器执行时令(一个或多个)处理器执行所描述的动作。额外地或备选地，计算机可读指令中的至少一个由信号、载波或其他暂态介质承载。
本文已经参考各种实施例描述了本发明。在阅读本文的描述后，其他人能够进行修改和变化。本发明旨在被解释为包括所有这些修改和变化，只要其在权利要求或其等价方案的范围内。

Claims (20)

1.一种用于无创地确定瞬时无波比率度量的方法，包括：
接收由成像系统生成的电子格式图像数据，其中，所述图像数据包括具有表示具有狭窄的血管的强度的体素；
根据所述图像数据计算所述血管的出口的外围阻力；
基于边界条件的集合和计算流体动力学算法来计算血管入口的入口与所述血管的所述出口之间的所述狭窄的狭窄阻力；并且
基于所述狭窄阻力来计算所述瞬时无波比率度量，其中，所述度量是数值，并生成指示所计算的瞬时无波比率度量的信号。
2.根据权利要求1所述的方法，还包括：
从所述图像数据分割所述血管；
根据所分割的图像数据确定所述血管的出口的血管直径；
根据所分割的图像数据确定主要冠状动脉分支的分支直径；
使用机器学习来估计主要冠状动脉分支的分支阻力；并且
其中，所述外围阻力是基于所述血管直径、所述分支直径和所述分支阻力计算的。
3.根据权利要求2所述的方法，其中，所述外围阻力被计算为由所述血管直径与所述分支直径的比率缩放的所述分支阻力。
4.根据权利要求3所述的方法，其中，所述外围阻力被计算为由所述血管直径与所述分支直径的比率的立方根缩放的所述分支阻力。
5.根据权利要求1至4中的任一项所述的方法，其中，所述边界条件的集合包括所述血管的入口处的压力和所述血管的输出处的速度。
6.根据权利要求5所述的方法，其中，所述计算流体动力学算法使用所述入口处的所述压力和所述输出处的所述速度作为输入参数来计算所述血管的所述出口处的估计的压力。
7.根据权利要求6所述的方法，其中，计算所述狭窄阻力包括通过以下步骤迭代地计算所述狭窄阻力：
使所述边界条件初始化；
使用所述边界条件来执行所述计算流体动力学算法以计算所述出口的所述估计的压力；
估计经更新的狭窄阻力；并且
响应于所述经更新的狭窄阻力不满足预定停止标准而执行第二次迭代。
8.根据权利要求7所述的方法，其中，所述预定停止标准包括将先前狭窄阻力和当前狭窄阻力之间的平方差的和最小化。
9.根据权利要求7至8中的任一项所述的方法，其中，所述第二次迭代包括：
基于所述经更新的狭窄阻力来计算所述出口处的经更新的速度；
使用所述入口处的所述压力和所述出口处的所述经更新的速度来执行所述计算流体动力学算法；
估计第二经更新的狭窄阻力；并且
响应于所述第二经更新的狭窄阻力不满足所述预定停止标准而执行另一次迭代。
10.根据权利要求7至9中的任一项所述的方法，还包括：
响应于所述经更新的狭窄阻力满足所述预定停止标准而基于所述狭窄阻力来计算所述瞬时无波比率度量。
11.根据权利要求1至10中的任一项所述的方法，还包括：
生成指示所述度量是否大于或小于预定阈值的分类信号，其中，高于所述阈值的度量值指示第一水平的狭窄，并且低于所述阈值的度量值指示第二水平的狭窄，其中，所述第一水平的所述狭窄比所述第二水平的所述狭窄更加严重。
12.一种用于无创地确定瞬时无波比率度量的系统(120)，包括：
狭窄阻力确定器(210)，其被配置为基于边界条件(212)的集合和计算流体动力学算法的结果来计算血管入口的入口与血管的出口之间的所述血管的狭窄的狭窄阻力；以及
iFR确定器(220)，其被配置为基于所述狭窄阻力来计算所述瞬时无波比率度量，其中，所述度量是数值，并且所述iFR确定器被配置为生成指示所计算的瞬时无波比率度量的信号。
13.根据权利要求12所述的系统，还包括：
感兴趣组织识别器(202)，其被配置为识别图像数据中的所述血管；
感兴趣组织几何结构提取器(206)，其被配置为根据所述图像数据来确定所述血管的出口的血管直径和主要冠状动脉分支的分支直径，并且使用机器学习来确定主要冠状动脉分支的分支阻力；以及
参数确定器(208)，其被配置为基于所述血管直径、所述分支直径和所述分支阻力根据图像数据来计算血管的出口的外围阻力。
14.根据权利要求13所述的系统，其中，所述外围阻力被计算为由所述血管直径与所述分支直径的比率缩放的所述分支阻力。
15.根据权利要求13所述的系统，其中，所述外围阻力被计算为由所述血管直径与所述分支直径的比率的立方根缩放的所述分支阻力。
16.根据权利要求12至15中的任一项所述的系统，其中，所述边界条件的集合包括所述血管的入口处的压力和所述血管的输出处的速度。
17.根据权利要求16所述的系统，其中，所述狭窄阻力被计算为减去所述出口处的压力的所述入口处的压力与所述出口处的所述速度的比率。
18.根据权利要求12至17中的任一项所述的系统，其中，所述狭窄阻力被迭代地确定。
19.根据权利要求18所述的系统，其中，每次迭代，经更新的狭窄阻力被计算，并且仅在所述经更新的狭窄阻力与先前狭窄阻力之间的平方和不满足停止标准时执行另一次迭代。
20.一种编码有计算机可读指令的计算机可读存储介质，所述计算机可读指令在由计算系统的处理器运行时令所述处理器：
仅根据无创地确定的参数计算iFR度量，所述无创地确定的参数是基于计算机断层摄影图像数据和计算流体动力学算法无创地确定的。
CN201680005857.XA 2015-01-15 2016-01-06 瞬时流量储备-计算机断层摄影 Active CN107580470B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201562103721P 2015-01-15 2015-01-15
US62/103,721 2015-01-15
PCT/IB2016/050043 WO2016113646A1 (en) 2015-01-15 2016-01-06 Ifr-ct

Publications (2)

Publication Number Publication Date
CN107580470A true CN107580470A (zh) 2018-01-12
CN107580470B CN107580470B (zh) 2021-06-22

Family

ID=55236828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201680005857.XA Active CN107580470B (zh) 2015-01-15 2016-01-06 瞬时流量储备-计算机断层摄影

Country Status (5)

Country Link
US (1) US10646185B2 (zh)
EP (1) EP3244790B1 (zh)
JP (1) JP6749917B2 (zh)
CN (1) CN107580470B (zh)
WO (1) WO2016113646A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111166315A (zh) * 2018-11-13 2020-05-19 苏州润迈德医疗科技有限公司 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法
CN111166316A (zh) * 2018-11-13 2020-05-19 苏州润迈德医疗科技有限公司 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法
WO2020186611A1 (zh) * 2019-03-19 2020-09-24 苏州润迈德医疗科技有限公司 基于压力传感器和造影图像计算瞬时无波形比率的方法
WO2021031355A1 (zh) * 2019-08-19 2021-02-25 苏州润迈德医疗科技有限公司 测量无波形期压力、比率方法、装置、系统及存储介质

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109843161B (zh) * 2016-09-30 2022-07-12 皇家飞利浦有限公司 用于确定针对狭窄评估的功能指数的装置
CN109616200A (zh) * 2018-11-06 2019-04-12 北京三普威盛科技有限公司 用于冠脉狭窄评估的方法，装置，存储介质及电子设备
CA3125211A1 (en) * 2019-01-06 2020-07-09 Covanos, Inc. Noninvasive determination of resting state diastole hemodynamic information
US10861157B2 (en) 2019-04-04 2020-12-08 Medtronic Vascular, Inc. System and methods for determining modified fractional flow reserve values

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1468305A (zh) * 2000-10-06 2004-01-14 亚文塔克有限公司 促炎基因产物转录之调制
US20090099755A1 (en) * 2007-10-15 2009-04-16 Harbert Richard H Even fire 90a°v12 ic engines, fueling and firing sequence controllers, and methods of operation by ps/p technology and ifr compensation by fuel feed control
US20100149920A1 (en) * 2008-12-17 2010-06-17 Medison Co., Ltd. Clutter Signal Filtering For Doppler Signal
CN101847208A (zh) * 2010-06-11 2010-09-29 哈尔滨工程大学 指纹与指静脉双模态识别的二级分类融合识别方法
CN101999002A (zh) * 2008-02-04 2011-03-30 彼帕科学公司 诊断和治疗parp-介导的疾病的方法
US20120053918A1 (en) * 2010-08-12 2012-03-01 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
WO2014002095A2 (en) * 2012-06-26 2014-01-03 Sync-Rx, Ltd. Flow-related image processing in luminal organs
US20140073977A1 (en) * 2012-09-12 2014-03-13 Heartflow, Inc. Systems and methods for estimating blood flow characteristics from vessel geometry and physiology
CN204089830U (zh) * 2014-02-25 2015-01-07 捷荣航材(广州)有限公司 一种vhf收发机综合试验装置

Family Cites Families (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000072037A1 (en) 1999-05-21 2000-11-30 Nycomed Imaging As Method of magnetic resonance imaging
GB0221434D0 (en) 2002-09-16 2002-10-23 Houston John G A method of analysing fluid flow in a conduit
IL165636A0 (en) 2004-12-08 2006-01-15 Paieon Inc Method and apparatus for finding the coronary velocity and flow and related parameters
WO2006061815A1 (en) 2004-12-08 2006-06-15 Paieon Inc. Method and apparatus for blood vessel parameter determinations
US8052611B2 (en) 2007-03-14 2011-11-08 Cardiac Pacemakers, Inc. Method and apparatus for management of heart failure hospitalization
DE102008014792B3 (de) 2008-03-18 2009-06-18 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Simulation eines Blutflusses in einem Gefäßabschnitt
US8200466B2 (en) 2008-07-21 2012-06-12 The Board Of Trustees Of The Leland Stanford Junior University Method for tuning patient-specific cardiovascular simulations
WO2010022762A2 (en) 2008-08-25 2010-03-04 ETH Zürich Method, system and device for enhancing flow field data
CN102165455B (zh) 2008-09-30 2015-08-26 皇家飞利浦电子股份有限公司 灌注成像
US20100125197A1 (en) 2008-11-18 2010-05-20 Fishel Robert S Method and apparatus for addressing vascular stenotic lesions
US20100130878A1 (en) 2008-11-24 2010-05-27 General Electric Company Systems, apparatus and processes for automated blood flow assessment of vasculature
US9405886B2 (en) 2009-03-17 2016-08-02 The Board Of Trustees Of The Leland Stanford Junior University Method for determining cardiovascular information
EP2742858B1 (en) * 2009-09-23 2024-06-05 Light-Lab Imaging Inc. Lumen morphology and vascular resistance measurements data collection systems, apparatus and methods
US20110307231A1 (en) 2010-06-09 2011-12-15 Jens Kirchner Method and arrangement for creating an individualized, computer-aided model of a system, and a corresponding computer program and a corresponding machine-readable storage medium
US8682626B2 (en) 2010-07-21 2014-03-25 Siemens Aktiengesellschaft Method and system for comprehensive patient-specific modeling of the heart
US8315812B2 (en) 2010-08-12 2012-11-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US9119540B2 (en) 2010-09-16 2015-09-01 Siemens Aktiengesellschaft Method and system for non-invasive assessment of coronary artery disease
DE102010043849B3 (de) 2010-11-12 2012-02-16 Siemens Aktiengesellschaft Vorrichtung und Computertomograph zur Bestimmung und Darstellung der Durchblutung des Herzmuskels
US10186056B2 (en) 2011-03-21 2019-01-22 General Electric Company System and method for estimating vascular flow using CT imaging
US20120296199A1 (en) 2011-03-21 2012-11-22 New York University Apparatus and Method of Non-Contrast Magnetic Resonance Angiography of Abdominal and Pelvic Arteries
US9135699B2 (en) * 2012-03-15 2015-09-15 Siemens Aktiengesellschaft Method and system for hemodynamic assessment of aortic coarctation from medical image data
BR112014028143A2 (pt) 2012-05-14 2017-06-27 Koninklijke Philips Nv sistema e meio de armazenamento legível por computador
US9675301B2 (en) 2012-10-19 2017-06-13 Heartflow, Inc. Systems and methods for numerically evaluating vasculature
CN104768465B (zh) * 2012-11-06 2018-06-19 皇家飞利浦有限公司 血流储备分数(ffr)指标
US9424395B2 (en) 2013-03-04 2016-08-23 Heartflow, Inc. Method and system for sensitivity analysis in modeling blood flow characteristics
US20140276137A1 (en) 2013-03-14 2014-09-18 Volcano Corporation Systems and methods for determining coronary flow reserve
US9595089B2 (en) 2014-05-09 2017-03-14 Siemens Healthcare Gmbh Method and system for non-invasive computation of hemodynamic indices for coronary artery stenosis
US9747525B2 (en) * 2014-06-16 2017-08-29 Siemens Healthcare Gmbh Method and system for improved hemodynamic computation in coronary arteries
US10130266B2 (en) * 2014-06-30 2018-11-20 Siemens Healthcare Gmbh Method and system for prediction of post-stenting hemodynamic metrics for treatment planning of arterial stenosis
CN106659400B (zh) * 2014-06-30 2021-01-05 皇家飞利浦有限公司 用于确定血流储备分数值的装置
US9888968B2 (en) * 2014-07-22 2018-02-13 Siemens Healthcare Gmbh Method and system for automated therapy planning for arterial stenosis

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1468305A (zh) * 2000-10-06 2004-01-14 亚文塔克有限公司 促炎基因产物转录之调制
US20090099755A1 (en) * 2007-10-15 2009-04-16 Harbert Richard H Even fire 90a°v12 ic engines, fueling and firing sequence controllers, and methods of operation by ps/p technology and ifr compensation by fuel feed control
CN101999002A (zh) * 2008-02-04 2011-03-30 彼帕科学公司 诊断和治疗parp-介导的疾病的方法
US20100149920A1 (en) * 2008-12-17 2010-06-17 Medison Co., Ltd. Clutter Signal Filtering For Doppler Signal
CN101847208A (zh) * 2010-06-11 2010-09-29 哈尔滨工程大学 指纹与指静脉双模态识别的二级分类融合识别方法
US20120053918A1 (en) * 2010-08-12 2012-03-01 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
WO2014002095A2 (en) * 2012-06-26 2014-01-03 Sync-Rx, Ltd. Flow-related image processing in luminal organs
US20140073977A1 (en) * 2012-09-12 2014-03-13 Heartflow, Inc. Systems and methods for estimating blood flow characteristics from vessel geometry and physiology
CN204089830U (zh) * 2014-02-25 2015-01-07 捷荣航材(广州)有限公司 一种vhf收发机综合试验装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SAYAN SEN: "《Development and Validation of a New Adenosine-Independent Index of Stenosis Severity From Coronary Wave-Intensity Analysis》", 《JOURNAL OF THE AMERICAN COLLEGE OF CARDIOLOGY》 *
TAYLOR C A: "《Computational Fluid Dynamics Applied to Cardiac Computed Tomography for Noninvasive Quantification of Fractional Flow Reserve: scientific basis》", 《JOURNAL OF THE AMERICAN COLLEGE OF CARDIOLOGY》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111166315A (zh) * 2018-11-13 2020-05-19 苏州润迈德医疗科技有限公司 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法
CN111166316A (zh) * 2018-11-13 2020-05-19 苏州润迈德医疗科技有限公司 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法
WO2020098140A1 (zh) * 2018-11-13 2020-05-22 苏州润迈德医疗科技有限公司 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法
CN111166316B (zh) * 2018-11-13 2023-03-21 苏州润迈德医疗科技有限公司 基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法
CN111166315B (zh) * 2018-11-13 2023-03-28 苏州润迈德医疗科技有限公司 基于造影图像计算瞬时无波型比率和静息态舒张期压力比率的方法
WO2020186611A1 (zh) * 2019-03-19 2020-09-24 苏州润迈德医疗科技有限公司 基于压力传感器和造影图像计算瞬时无波形比率的方法
WO2021031355A1 (zh) * 2019-08-19 2021-02-25 苏州润迈德医疗科技有限公司 测量无波形期压力、比率方法、装置、系统及存储介质

Also Published As

Publication number Publication date
JP2018501894A (ja) 2018-01-25
EP3244790A1 (en) 2017-11-22
JP6749917B2 (ja) 2020-09-02
WO2016113646A1 (en) 2016-07-21
US10646185B2 (en) 2020-05-12
EP3244790B1 (en) 2020-09-02
CN107580470B (zh) 2021-06-22
US20180271468A1 (en) 2018-09-27

Similar Documents

Publication Publication Date Title
CN107580470B (zh) 瞬时流量储备-计算机断层摄影
US10595806B2 (en) Fractional flow reserve (FFR) index with adaptive boundary condition parameters
JP6484760B2 (ja) 非侵襲的血流予備量比（ｆｆｒ）に対する側副血流モデル化
US11195278B2 (en) Fractional flow reserve simulation parameter customization, calibration and/or training
WO2016001017A1 (en) Apparatus for determining a fractional flow reserve value
WO2015164086A1 (en) Method and system for hemodynamic computation in coronary arteries
EP3606433B1 (en) Standardized coronary artery disease metric
US20190076196A1 (en) Vessel geometry and additional boundary conditions for hemodynamic ffr/ifr simulations from intravascular imaging
WO2018095791A1 (en) Vascular tree standardization for biophysical simulation and/or an extension simulation for pruned portions
EP3602485B1 (en) Interaction monitoring of non-invasive imaging based ffr

Legal Events

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