[go: up one dir, main page]

CN105792738B - 用于改善功能狭窄分析的局部ffr估计和可视化 - Google Patents

用于改善功能狭窄分析的局部ffr估计和可视化 Download PDF

Info

Publication number
CN105792738B
CN105792738B CN201480066240.XA CN201480066240A CN105792738B CN 105792738 B CN105792738 B CN 105792738B CN 201480066240 A CN201480066240 A CN 201480066240A CN 105792738 B CN105792738 B CN 105792738B
Authority
CN
China
Prior art keywords
ffr
image
local
transfer functions
pressure
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.)
Active
Application number
CN201480066240.XA
Other languages
English (en)
Other versions
CN105792738A (zh
Inventor
H·尼克基施
M·格拉斯
H·施米特
J·蒂默
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 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 NV filed Critical Koninklijke Philips NV
Publication of CN105792738A publication Critical patent/CN105792738A/zh
Application granted granted Critical
Publication of CN105792738B publication Critical patent/CN105792738B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/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/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • 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/74—Details of notification to user or communication with user or patient; User input means
    • A61B5/742—Details of notification to user or communication with user or patient; User input means using visual displays
    • 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/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • 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/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • 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
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/026—Measuring blood flow
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/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/48—Diagnostic techniques
    • A61B6/481—Diagnostic techniques involving the use of contrast agents
    • 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
    • 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/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G06T2207/10121—Fluoroscopy
    • 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/20—Special algorithmic details
    • G06T2207/20081—Training; Learning
    • 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
    • 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/30172—Centreline of tubular or elongated structure

Landscapes

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

Abstract

一种用于流量储备分数FFR模仿的系统(IPS)和相关方法。针对脉管系统部分的FFR值范围的模仿基于复合转移函数，根据全局效应转移函数he的加权和来组合所述复合转移函数，每个全局效应转移函数表示引起压降的不同物理效应。相对与各自的脉管系统几何结构相关联的压力pi与流率fi样本测量，从先前的训练阶分割获得权重we。根据在各自的间隔内的压力值和流率值，所模拟的FFR值的范围在图形显示器(GD)中被可视化。

Description

用于改善功能狭窄分析的局部FFR估计和可视化
技术领域
本发明涉及一种图像处理方法、图像处理系统、计算机程序元件以及计算机可读介质。
背景技术
世界上的群体受到心血管疾病的严重影响。心血管疾病仍是世界范围内导致死亡的主要原因。缺乏锻炼、不健康的饮食和其它因素例如可能导致狭窄，即，一部分心血管系统收缩的状况。狭窄可能引起狭窄下游的组织供给不足。尽管可以对狭窄进行处置，但是这些处置大部分是介入性的且具有风险。例如，一种介入性过程是经由导管将支架引入到脉管系统中。然后支架将导航到狭窄部位，并在该位置处扩张以便放大，从而恢复丢失的血管横截面。但是这些类型的介入性处置具有风险，并且它们给已经缺乏现金的国民医疗制度带来压力(有时甚至是不必要的)。例如，如果狭窄发生在血管供应已经垂死的组织的部分，则介入是否将实际导致任何积极的健康效果是值得怀疑的。垂死的组织不能复活，并且因为狭窄受限的血流在一些情况下仍是刚好足够的。在此FFR(流量储备分数)技术确定是有帮助的，因为其允许依据分数(即，数字)来评估狭窄的严重性。FFR是对除冠状动脉几何结构之外的狭窄严重性的度量，因为其包括狭窄对血流的影响。在导管插入过程期间有创式FFR测量的临床关联已经在临床研究中得以证明。随后，医疗人员可以更好地判断在所述情况下的介入是否是有保证的。一种建立这种数字的方法是采取原位测量，其包括跨狭窄部位的压差和通过狭窄的血液流率。为此使用一种特殊的导管，其装备了合适的探头以获得所述测量结果。导管被放置在远端，然后接近狭窄部位以特别建立压差。但是，同样地，这是有风险的，因为导管需要强制通过狭窄部位。另外，介入性工具本身是复杂的且昂贵的，通常需要给予潜在的有害物质，例如腺苷，以增加心脏区域的充血。
因此，过去已经提出用于导管插入实验室应用的无创局部FFR估计方法，以避免使用昂贵的原位压力测量设备。这些方法的变型基于根据CT投影重建的体积数据。该方法的另一变型一般基于根据一组血管造影投影(3D冠状动脉造影(3DCA))生成的冠状动脉的体积模型。两种变型都基于计算流体力学(CFD)模拟。
这些已知的方法被认为有时是计算量大和/或倾向于招致针对患者相对高的剂量。
发明内容
因此需要有一种替代方式来更有效地确定流量储备分数值，并为患者减少麻烦。
本发明的目的是通过独立权利要求的主题来解决的，其中，其它的实施例被并入到从属权利要求中。应该注意的是，本发明以下描述的方面同样应用于图像处理方法、计算机程序元件和计算机可读介质上。
根据本发明的第一方面，提供一种图像处理系统，包括：
输入端口，其用于接收包括感兴趣对象的投影视图的至少一幅图像；
分割器，其被配置为分割所述图像，以依据所述投影视图获得针对所述对象的对象分割；
划分器，其被配置为将所述分割划分为部分；
适配单元，其被配置为将多个局部效应转移函数中的一个或多个适配到在各个所述部分中记录的所述对象的局部几何结构中，以获得多个全局效应转移函数；
组合器，其被配置为将所述多个全局效应转移函数组合为针对所述对象的复合转移函数；
评估器，其被配置为根据所述复合转移函数计算针对流经所述对象的流体的FFR估计。
例如，当使用提出的系统来模拟狭窄部位周围的冠状的相关部分的FFR时，局部效应转移函数描述基本单独归因于特定物理效应的局部压降，而全局效应转移函数描述由跨整个成像的血管树或者至少在狭窄周围的整个感兴趣区上的各自的效应所引起的压降。复合转移函数则描述跨对象的整体压降，同时解释在各个效应之间的相关性和相互作用。换言之，提出一种“集中”参数模型，在这里是由用于模拟考虑中的流体行为的几何导出转移函数定义的，例如，在一个实施例中，人或动物的血液循环。
根据一个实施例，评估器被配置为通过改变至少一个物理或生理参数来计算FFR估计的范围。
换言之，评估器对由在物理或生理参数之间的功能关系定义的FFR表面进行采样。
根据一个实施例，所述图像处理系统还包括可视化器，其被配置为在显示单元上绘制以显示根据至少一个物理或生理参数的FFR估计的范围。在一个实施例中，然后可以在显示单元上以FFR估计的函数依赖来显示FFT估计的范围。对血管中的生理状态的范围的FFR值进行可视化允许解释估计FFR计算的边界条件的解剖知识是有限的事实，从而使得计算更简单并因此更快。
根据一个实施例，所述参数是压力和流率中的至少一个，并且其中，根据压力和流率中的至少一个来显示FFR估计的范围。所显示的FFR估计还可以基于改变用于计算FFR的其它变量/参数。例如，在一些实施例中，替代或除了流率/压力变化外，可以改变以下中的任何一个(或以下的任何组合)：血比容水平、血液粘度或血液密度。在一个实施例中，用户可以选择改变哪些变量。
根据一个实施例，使用先前从压力与流率采样测量结果/模拟学习的权重，对局部效应转移函数进行线性组合。换言之，本发明管理先前原位压力/流率测量和/或冠状动脉CFD计算的数据语料库，以学习各种物理效应(例如，Bernoulli和/或(Poiseuille)摩擦等)对跨冠状动脉的整体或净压降的(由所述权重测量的)相对贡献。在将权重训练到足够大的数据语料库之后，提出的系统准备好以供使用。当使用提出的系统时，可以避免需要为将来FFR模拟进行原位测量和/或运行计算量大的CFD模拟。
根据一个实施例，通过C型臂类型的平面X射线装置采集所述图像。特别地，没有必要向提出的装置馈送重建的CT体积图像。尽管使用来自不同投影的投影越多，则模拟的保真度越好，投影视图的若干离散角度样本是足够的。在一个实施例中，所述适配基于由3D血管中心线和CSA(横截面积)总结的几何结构，这两者例如都可以根据来自沿着不同视角采集的一幅或多幅2D血管造影片的旋转X射线采集(圆形或其它)而生成。这里设想用于投影图像的采集的多个不同成像模式，例如CT、旋转C型臂、MRI及其它。在这里还可以设想相衬成像或光谱图像。在一个实施例中，相衬成像在没有造影剂注射的情况下使用。
尽管提出的方法可以用于管理先前的原位(即，介入性)压力测量以训练用于组合复合函数的权重，但是一旦获得了权重就不再依靠或没有必要为了压力测量目的而进行介入，尽管这种原位压力测量仍在此使用。换言之，不要求用于采取压力测量的冠状动脉介入来获得给定患者的FFR信息。
提出的系统允许在CathLab(导管实验室)中使用关于血管几何结构的3D信息来生成虚拟FFR信息，并将所述信息呈现给介入性心脏病学家。
另外，已经观察到，CT冠状动脉造影(CTCA)和CFD依赖于边界条件扩展集合的易于误差估计。例如，基于CTCA的虚拟FFR要求评估在血管系统的流入或流出处的边界条件。但是难以估计与向血管系统供应压力的心肌相关的边界条件。因此，在一个实施例中，这里提出的针对边界条件的范围的一组FFR值的相对快速计算是有利的。另外，不需要给予对充血最大化物质。
尽管在一个实施例中，提出的系统用于FFR模拟以调查在狭窄冠状动脉中的血流行为，但该系统的益处还在于在人(或动物)脉管系统中的其它部分中使用。另外，所述系统的用途不限于医疗环境。例如，在此还可以想到检查土木工程中不可访问的复杂管道系统，或检查生产设备或交通工具(例如飞机或轮船)的水利管网。
附图说明
现在将参考附图描述本发明的示例性实施例，其中：
图1示出了用于模拟FFR信息的图像处理系统的示意性框图；
图2A、2B示出了图示血管中心线和横截面积估计的人类脉管系统的投影视图；
图3示出了3D中心线点的对极几何估计；
图4示出了用于显示FFR值的范围的图形显示；
图5示出了图像处理方法的流程图；
图6示出了用于适配局部转移函数的几何相关参数的表格。
具体实施方式
参考图1，示出了允许FFR模拟的图像处理系统IPS。该系统允许根据图像Ai计算FFR值，所述图像Ai编码图像对比度足以识别患者PAT的血管树形式的冠状动脉的各自的投影视图(“足迹”)。
通过成像器IM采集图像Ai，例如平面X射线放射线照相装置或旋转X射线C型臂系统或者其它适当的成像模态。成像器IM包括X射线源XR和探测器D。患者在相关感兴趣区域(例如，心脏区域)处暴露于从X射线源XR发出的辐射。辐射与感兴趣区域(在该情况下是心脏区域)中的物质相互作用，然后在探测器D处被探测到。在探测器D处探测到的辐射被转换为所述数字图像Ai，所述数字图像由图像处理系统IPS在其输入端口IN处输出并接收。
在一个实施例中，供应图像Ai的成像器IM是C型臂X射线荧光透视系统。尽管原则上FFR模拟基于所采集的感兴趣冠状动脉的单个图像Ai，但是优选实施例使用至少两幅图像(在一些实施例中，恰好2幅图像)。在一个实施例中，使用双平面X射线C型臂系统，并且在不同的(例如，正交)投影方向上采集两幅图像A1、A2，并且已经发现这样的图像对A1、A2供应关于血管树的3D结构的足够的图像信息以达到该目的。由于脉管系统组织的低辐射吸收属性，在一个实施例中，所述图像是血管造影片，即，它们是在造影剂驻留在患者的脉管系统中的同时采集的X射线图像。然而，这只是在一个实施例中，使用具有启用相衬成像设备的成像器，在此也设想为这些类型的成像系统在一些情况下能够提供足够的对比，即使在脉管系统中没有驻留造影剂。在一个实施例中，相衬成像器包括被布置在X射线管与探测器D之间的一系列干涉光栅。
广义上，本系统IPS为特定患者提供计算FFR值，而在计算或模拟时不采取原位压力测量，并不运行扩展的CFD模拟。与这些方法相反，所述系统使用在过去已经获得的测量结果的数据语料库。更具体地，训练模块TR以血流和压降样本对<fi,Δpi>形式利用该知识。从该数据语料库中学习一组权重we。通过组合多个局部适配的局部效应转移函数或基础函数模板，权重we能够被用于“构成”复合转移函数。每个局部效应转移函数独立表示每个可以引起某个压降的不同的个体物理效应。那些基本函数的模板被保存在库或数据库DB中，并被适配到接收到的实际患者的图像A1、A2中的图像信息。现在将之前仅独立考虑的个体物理效应组合考虑，以便对那些效应之间的各种相互作用进行更好的建模。这允许达成复合转移函数，该函数已经被发现以合理的计算成本提供对狭窄区域周围流体动力学的非常现实的建模。已经发现，由于本发明的系统与计算量大的CFD相比相对简单，所提出的系统高度响应。当处理图像时，近乎实时地递送FFR估计以识别(例如，经由分割)血管树，这在忙碌的导管实验室环境中是有益的，在该环境中需要快速决定是否进行介入。在将复合函数适配或“裁剪”到即时患者的图像信息之后，复合函数用于计算FFR值的范围，然后由可视化器VS在屏幕MT上进行可视化，以便更好处理剩余的不确定性，并由此向用户提供狭窄部位的FFR情形的现实的高保真图片。然后，医生可以更好地评估狭窄的严重性。现在将更详细地解释图像处理系统的操作。下述描述的第一部分将集中于计算方面，而第二部分将集中于可视化方面。
下文总结了在不同实施例中使用不同的成像模式来采集一幅或多幅投影图像。如上所述，在一个实施例中，经由通过导管的局部造影剂注射而生成血管造影投影。
根据一个实施例，通过平面X射线成像器采集单个血管造影投影。
根据备选实施例，使用具有单个造影剂团注的双平面系统，以及时采集在相同点处的两个投影，这允许生成更准确的3D血管模型。
根据备选实施例，使用具有两次采集以及两个或更多造影剂注射的单个平面C型臂系统，来根据不同的投影生成更准确的模型。在一个实施例中，人们采集ECG(心电图)，同时组合相同心脏相位的投影。
人们能够使用具有单个造影剂注射的单个平面C型臂系统，并操作C型臂系统以获取旋转采集，然后从该序列中选择ECG相位对应的血管造影片以到达不同的采集角度，以生成更准确的模型树血管模型。
根据一个实施例，使用CT扫描器来提供投影图像。
最后，人们可以使用CT或MR扫描器，并对来自体积数据集的3D冠状动脉进行分割，从而确定中心线和横截面，尽管操作并不必然依赖于体积图像数据。
在优选实施例中，接收到恰好两幅血管造影图像A1、A2，优选但不是必要地在正交投影方向采集，并且将参考这两幅图像解释以下实施例，但是可以理解的是，在备选实施例中(如上所列出的)使用单个图像，或者在又一实施例中使用多于两幅图像，例如三幅或四幅图像。
以任何速率，不管成像模态，对于本发明目的，投影图像输入的数量应该在足够的程度上允许计算在感兴趣区域处血管树的局部横截面面积的近似。换言之，如果只使用单个图像，这里涉及的计算(在下文中将更详细地解释)是固有的近似，并将最终依赖于合理的假设，因为血管沿着未在单个图像中记录的维度延伸。
然而，当处理血管时，假设圆形或至少是椭圆形的横截面是合理的，因此两幅图像通常足以估计短轴和长轴，利用足够细节来估计横截面面积，以达到本发明的目的。
当接收仅在感兴趣区域周围的几个离散“采样”投影方向上采集的投影图像作为输入时，所提出的系统能够给出令人满意的结果。体积(即，重建的CT图像数据)不是必要的，但是这里也设想其在备选实施例中的使用。
操作
在一个非限制性实施例中，在输入端口IN处接收两幅图像A1、A2。
图像A1、A2然后被传递到分割器SEG。分割器SEG操作以针对在每个视图内的血管树足迹分割两幅图像。在一个实施例中，分割整个血管树，但是在其它实施例中，仅分割在实际狭窄周围由半径(ROI半径)限定的感兴趣区域，这通常是容易识别的。在一个实施例中，狭窄部位周围的ROI半径是可调整的，并延伸到狭窄每侧的大约2到5cm处，但这仅仅是非限制性示例，在其它使用场景中还可以使用有益的其它ROI半径定义。
随后在所述ROI中分割的血管足迹的图像坐标被传递给划分器PAR。然后划分器PAR沿着分割的血管树部分行进以按步长将其划分到部分(section)中，所述步长在实施例中被设置为大约5到10mm。在一些实施例中，所述步长是用户可调整的。由于在图像采集期间成像器的几何结构是已知的，所需的步长可以被转换为像素单位，划分器PAR然后可以在所述像素单位上操作以针对每幅输入图像A1、A2定义多个图像部分，后文称作“血管树分割”j。每个血管树分割j在各自的投影视图中记录特定的几何参数，所述几何参数定义对应的血管树部分的局部几何结构。在一个实施例中，几何参数包括以下中的任意一个(或其任意组合)：一部分血管中心线、血管的局部横截面面积A、局部中心线曲率κ、血管分割的长度l、局部血管周长P以及局部血管半径r。
然后将每幅图像A1、A2的血管树分割j传递给适配器ADP。适配器ADP被通信耦合到数据库DB上，在所述数据库中保存有“基本”函数he(f)的库。每个基本函数对应于局部转移函数的模板。每个转移函数定义针对将引起流体动力学行为改变的特定物理效应的专用流体动力学模型。更具体地，并且在一个实施例中，在给定流率f和给定感兴趣管状结构的特定局部几何结构的情况下，库中的每个局部效应转移函数描述由一个特定物理效应单独引起的特定压降Δp。
局部效应转移函数he(f)的个体模板被建模为具有次数d的奇次多项式：he(f)＝αesign(f)|f|d (1)
从包括表示为与各自的效应e相关联的局部几何系数αe的参数化的意义而言，该函数是模板。在一些情况下，血液的某些流体属性、局部几何系数αe单独取决于局部血管几何结构，即，取决于血管的分割轮廓或内腔以及由每个血管树分割j所捕获到的其中心线(根据其可以计算系数)。符号函数sign(f)表示流f的方向(“+”/“-”)。
根据一个实施例，依据局部几何系数αe，图6中的表格包括保存在库DB中的不同局部效应转移函数模板的示例。
在该表格中，局部几何参数A、κ、l、P和r如前所定义的，并且ρ和μ表示感兴趣流体的属性。在一个实施例中，ρ和μ表示生理参数，特别地，分别表示血液密度和血液粘度。
数据库保存集中参数模型，其能够被用于模拟人类血液循环，如本文更详细描述的。换言之，图6中的表格总结各个局部转移函数模板如何能被适配到血管的局部几何结构(如局部几何系数αe所表示的)，以便单独(即，独立地)为在最左列中命名的各个压力效应计数。例如，第一行描述了如何对仅由伯努利(Bernoulli)效应引起的压降进行建模。即，如果血管在血管部分处局部具有如最右列示出的“哑铃”形状的局部几何结构，以及如果只有伯努利效应是压力变化的唯一原因，则在对局部几何参数A输入、A输出适用了等式(1)之后，可以根据等式(1)对所述血管部分计算所述压力改变。
但是，局部压降通常是由多于一个效应引起的。因此，为了增强本文提出的FFR模拟的保真度，适配器ADP操作以为每个血管树分割j对(来自两幅图像A1、A2)不仅针对特定效应计算一个局部效应转移函数，而且为在所述分割j处的每个效应计算专用局部转移函数。可以理解的是，通常不是每个部分j针对每个效应e引起非零局部效应转移函数贡献。例如，相对直的血管分割将对博尔达-卡诺扩张效应返回基本为零的贡献。对在图像A1、A2的对中记录的每个血管部分j进行同样处理。对于每个血管树分割j，输出是多个局部效应转移函数，对不同效应一个函数。换言之，对于每个血管树分割j，计算不同的效应依赖系数αe,j。在以这种方式处理所有的血管树部分j之后，可以跨血管树部分j将属于同一效应的系数αe,j相加，以为各个效应计算整体的全局(即，为整个血管或血管ROI)系数αe＝∑jαe,j。使用通过该相加构造的全局系数，允许为整个血管或血管ROI的每个效应依据(1)用公式表达全局效应转移函数。结果是将来自每个血管部分j的关于个体、局部效应特定压降的信息依据等式(1)联合为各自的全局效应转移函数he(f)＝∑jαe,jsign(f)|f|d。
可以依据图6的表格，根据如在血管树分割j中记录的局部几何参数来计算局部系数。例如，适配器ADP可以被配置为使用冠状动脉建模技术来确定泊肃叶(Poiseuille)摩擦系数αP以及曲率效应系数αC。更具体地，可以基于根据两个A1、A2(或更多的投影)建模的3D血管中心线，使用如图3所示的对极采集几何结构来确定两个系数，其中
Figure BDA0001008897320000101
表示在采集两幅投影图像A1、A2时在成像器的X射线源XR中各个焦点位置。由于已知两幅投影A1、A2之间的几何关系，3D中心线点可以被确定为如图3所示。可以以对极几何方式逐点构造在每幅图像A1、A2中的中心线。可以从图像头部中的DICOM元数据检索关于焦点位置的空间信息，或者可以在采集期间从成像器的几何结构的记录中检索该信息。
可以基于在两幅血管造影投影图像A1、A2中对应的血管位置确定出的至少两个半径计算所有其它效应。换言之，可以通过沿着3D血管中心线所取的横截面面积CSA，来总结基底血管几何结构。因此，为了计算这些效应，并不要求生成血管的完全的体积数据集(体素元数据集)或3D表面模型。最小的要求是根据两幅或更多投影以及采集几何结构的知识来确定3D中心线。总之，使用3D中心线和采集几何结构，对于中心线上的一个3D点，可以基于相关的投影来确定对应的血管和血管半径，并且可以依据图6的表格来确定相关的几何参数，例如A、P或r。
图2A图示了A1或A2作为血管造影投影，以粗体示出了血管中心线。图2B图示了A1或A2作为血管造影投影，示出了跨中心线所取的局部血管半径。对几何参数的图像内测量能够基于灰度值阈值。如可以在图2B中看出的，对于中心线的例子，几何参数通常是从如在两幅投影A1、A2中的对应树部分对j中记录的集体图像信息导出的。
现在返回到在图1中提出的图像处理系统IPS中的过程流，然后将多个全局效应转移函数(即，每个效应e一个)传递到组合器∑，以通过各种全局效应转移函数的线性组合来形成复合转移函数Δp＝∑ewehe(f)。这允许用公式表达简化的模型来计算在狭窄周围的血管ROI中的流体行为。现在将在给定体积流率f处的狭窄血管分割上的压降Δp建模为不同全局压降效应he(f)的加权和，以说明由于以下中任一个或其组合引起的压力变化：i)伯努利原理，ii)由于血管椭圆度和/或曲率和/或横截面几何的变化(例如，扩张和收缩)引起的摩擦力造成的压力损失。
在形成复合转移函数时，从所述(或不同的)数据库DB中检索权重we。在一个实施例中，先前通过训练模块TR在学习阶分割通过统计训练学习过程来计算效应权重。具体地，给定一组例子，沿着血管树的几何结构的体积流fi和压降Δpi，在一个实施例中，权重的拟合实现最小二乘法拟合、正规化的最小二乘法拟合或者非负的最小二乘法拟合，其中通过类似于w＝arg minw∑i(Δpi-∑ewehe(fi))2最小化损失函数来发现权重向量w＝[we]e，其中i表示样本的索引，并且e是表示人们希望在训练阶段解释的各个效应e的索引。例如，使用依据图6的库，这里所表明的效应以及e＝5可以理解为图6库仅是一个实施例，并且可以在其它实施例中使用不同的效应。在这里假设训练样本不仅包括测量fi、Δpi，而且包括对应的至少描述各自的测量fi、Δpi的局部几何结构的几何参数，如A、r(如上所定义的)等。
可以通过CFD模拟或通过模拟压力/流或压力/速度测量结果来获得训练样本。如这里提出的，适应性选择权重来解释不同效应e之间的相互作用、相互依赖性和相关性。这与先前单独分离地检查效应的方法非常不同，而在此实现的方法目的在于从现有数据语料库推断各种压力损失效应之间的相互作用。
一旦得到复合转移函数，并将其适用于附近患者的血管树，评估器EVAL则能够为任何给定的对p0、f依照以下计算不同的FFR值：
Figure BDA0001008897320000111
换言之，在狭窄之后的远端压力p0除以在狭窄位置之前的近端压力。FFR值取决于近端压力和流经狭窄区域的量。然后在输出端口(未示出)输出这样计算出的FFT值，以进行存储和/或其他处理。
FFR等式(2)定义了2D表面，其中FFR取决于两个未知的p0、f。在一个实施例中，评估器EVAL操作以依据(2)对FFR表面进行采样，以通过改变在生理合理区间内的流量和/或近端压力而产生不同FFR值的范围。
在一个实施例中，随后将评估器EVAL输出的多个数据三元组<p0,f,FFR(p0,f)>传递给可视化器VS。在一个实施例中，可视化器VS操作以通过将采样的表面<p0,f,FFR(p0,f)>映射到平面上而形成图形显示，并且根据用户可定义的调色板来对FFR值的幅度进行颜色或灰度值编码。然后绘制图形显示器GD，以通过驱动可视化器VS输出所控制的系统视频卡而在监视器MT上进行显示。图4示出了这种图形显示的一个实施例。对经采样的FFR表面的图形绘制通过将FFR表示为近端压力和体积流率的函数而允许医生更好地评价狭窄的严重性。假设对近端血压的测量结果是可用的，则2D图形显示可以塌缩为1D曲线。可以从导管基于主动脉中的压力的测量结果或者根据外部血压测量结果(臂式电子血压计)而估计该近端血压值。颜色编码指示在特定工作制度中狭窄的严重性。如果图像采集是基于血管造影术，则造影剂通过或给予的导管可以用于获得近端压力的测量结果。没有必要(尽管仍可以)使用专用压力测量导管(压力导线)。
类似地，人们可以基于是计算的一部分但是尚未确切知道的其它量来将FFR依赖进行可视化，例如血比容、血液粘度、血液密度。
另外，由于根据多幅2D图像的CSA估计固有地是近似的，人们也能够基于FFR值可视化CSA估计的引入效应。用户随后能够在那些图之间切换，或者将其组合或聚集为3D图。换言之，通过在等式(1)、(2)中“冻结”或保持恒等的某些用户可选择值，可视化器允许用户检查附到FFR值的不确定性，然后让一个或多个剩余值在用户可定义误差区间内变化。例如，可以向用户呈现具有所有变量列表的用户界面。然后，用户点击或以其它方式指定他想要保持恒等的那些变量。剩余的变量随后在各自的误差容限内变化，并依赖于这些变量而显示得到的FFR值。例如，当改变CSA估计A时，在一个实施例中用户能够指定将变化应用于哪个血管树分割j处。然而，在更简单的实施例中，固定误差估计被应用于跨所有的血管树分割部分j的所有的CSA估计。
临床研究已经表明低于0.75或0.8的FFR值是临界值，这是因为它们是对狭窄严重性保证性介入的指示。然后在一个实施例中提出依据图1在图形显示器上在临界0.8FFR值处强加轮廓线。然而，0.8(或任何其它适当建议的阈值)轮廓线仅是一个实施例，并不应该解释为限制。在其它环境中，可以替代使用不同值的轮廓线。
在备选实施例中，不是FFR值本身显示在2D或3D图中，而是在计算出的图中显示针对给定的血管几何结构的高于或低于0.8阈值的FFR值的面积或体积。这允许将可用信息浓缩为单个数。
尽管在实施例中结合图1在上文解释了适配器ADP，将所有的效应都视为用于每个血管树部分j，但是可以理解的是，这仅是一个实施例。在备选实施例中，仅将在数据库中保存的局部效应转移函数的子选择应用到血管树分割部分。例如，在一个实施例中，用户可以决定他们希望考虑在各个部分j处的哪个效应。在一个实施例中，向用户提供了图形用户界面，其包括类似于图6的表格中最右列的图形符号，示出了针对各个局部几何结构的建议性符号。用户则可以点击或以其它方式指定在计算局部转移函数时由系统考虑哪些局部效应转移函数。
参考图5，示出了图像处理方法的流程图。
在步骤S505处，接收包括感兴趣对象的投影视图的至少一幅图像。在一个实施例中，感兴趣对象是人或动物心脏脉管系统的一部分。
在步骤S510处，针对在投影视图中捕获到的对象的足迹将图像进行分割，以获得对象分割，例如血管树分割。在一个实施例中，所述分割被限制在以狭窄部位为中心的感兴趣区域。
在步骤S515处，所述分割被划分到图像部分j。在一个实施例中，所述部分是沿着ROI内的血管树的部分。
在步骤S520处，多个局部效应转移函数的一个或多个(在一个实施例中是每一个)被适配到在各自的部件j中记录的对象的局部几何结构中，以获得多个全局效应转移函数he。
在步骤S525处，多个全局效应转移函数被组合为针对所述对象的复合转移函数。
在步骤S530处，根据复合转移函数计算针对流经对象的流体的FFR估计。在一个实施例中，计算步骤包括通过改变至少一个物理或生理参数来计算FFR估计的范围。
在步骤S535处，FFR估计的范围被显示为至少一个物理或生理参数的函数。
在一个实施例中，图1的图像处理系统IPS的部件全部运行在单个计算系统上。在备选实施例中，在此同样设想至少部分分布式架构，其中一个或多个部件远程定位，并在适当的通信网络中彼此连接和/或与图像处理器IPS连接。
在一个实施例中，图像处理IPS(或者至少其一些部件)被布置为专用FPGA或硬接线(独立)芯片。
在备选实施例中，图像处理IPS或至少其一部分部件驻留在成像器IM的工作站中。
图像处理系统IPS的部件可以在适当的科学计算平台(例如
Figure BDA0001008897320000141
)中编程，并被翻译成适合在计算系统(例如，成像器的工作站)上运行的C++或C例程。在本发明的另一示例性实施例中，提供计算机程序或计算机程序元件，其特征在于，适于在合适的系统上执行根据先前实施例中的一个所述的方法的方法步骤。
因此，计算机程序元件可以被存储于计算机单元中，所述计算机单元也可以是本发明实施例的一部分。计算单元可以适于执行或引起执行上述方法的步骤。此外，其可以适于操作上述装置的部件。计算单元能够适于自动操作和/或执行用户的命令。计算机程序可以加载到数据处理器的工作存储器中。数据处理器因此可以被装备以执行本发明的方法。
本发明的该示例性实施例覆盖从开始就使用本发明的计算机程序，以及通过更新将现有程序转为使用本发明的程序的计算机程序。
此外，计算机程序元件能够提供为完成上述方法的示例性实施例的过程的所有必需步骤。
根据本发明的另一示例性实施例，提出一种计算机可读介质，例如CD-ROM，其中，该计算机可读介质存储有计算机程序元件，所述计算机程序元件是通过先前段落描述的。
计算机程序可以存储于和/或分布于适当的介质上，例如光存储介质或与其它硬件一起供应或作为其它硬件一部分的固态介质。所述计算机程序还可以例如经由互联网或其它有线或无线电信系统以其它形式分布。
然而，计算机程序还可以通过诸如万维网的网络来呈现，并可以从这种网络下载到数据处理器的工作存储器中。根据本发明的另一示例性实施例，提供了使得计算机程序元件可用于下载的介质，所述计算机程序元件被布置为执行根据先前描述的本发明的实施例之一的方法。
必须注意的是，结合不同的主题描述了本发明的实施例。特别地，一些实施例通过方法类型权利要求进行描述，而其它实施例通过设备类型权利要求进行描述。然而，除非另有指出，否则本领域技术人员可以从上述和以下描述中获知，除了属于一种类型主题的特征的任意组合外，在属于不同主题的特征之间的任意组合也被视为被本申请所公开。然而，可以组合所有的特征来提供比简单相加特征更多的协同效应。
尽管已经在附图中以及前述描述详细图示和说明了本发明，但这种图示和描述被视为说明性或示例性而非限制性的。本发明并不限于所公开的实施例。在实践所要求保护的发明时，根据研究附图、公开物和随附权利要求，本领域技术技术人员可以理解和实现公开实施例的其它变型。
在权利要求中，词语“包括”并不排除其它元件或步骤，并且不定冠词“一(a)”或“一(an)”并不排除复数。单个处理器或其它单元可以完成在权利要求中记叙的若干项的功能。事实上，在相互不同的独立权利要求中记叙的某些测量并不表明不能有利地使用这些测量的组合。在权利要求中的任何附图标记不应被解释为限制范围。

Claims (10)

1.一种图像处理系统(IPS)，包括：
输入端口(IN)，其用于接收包括感兴趣对象的投影视图的至少一幅图像；
分割器(SEG)，其被配置为分割所述图像，以依据所述投影视图来获得针对所述对象的对象分割；
划分器(PAR)，其被配置为将所述分割划分为部分；
适配(ADP)单元，其被配置为将多个局部效应转移函数中的一个或多个适配到在各个所述部分中记录的所述对象的局部几何结构中，以获得多个全局效应转移函数；
组合器(Σ)，其被配置为将所述多个全局效应转移函数组合为所述对象的复合转移函数；
评估器(EVAL)，其被配置为根据所述复合转移函数来计算流经所述对象的流体的FFR估计。
2.根据权利要求1所述的系统，其中，评估器(EVAL)操作以通过改变至少一个物理或生理参数来计算FFR估计的范围，所述图像处理系统还包括：
可视化器(VS)，其被配置为在显示单元(MT)上绘制以显示根据所述至少一个物理或生理参数的FFR估计的所述范围；
所述显示单元(MT)，其操作从而以FFR估计的函数依赖来显示FFR估计的所述范围。
3.根据权利要求2所述的系统，所述参数是压力和流率中的至少一个，并且其中，根据压力和流率中的至少一个来显示FFR估计的所述范围。
4.根据权利要求2所述的系统，所述参数是以下中的至少一个：血比容、血液粘度、血液密度。
5.根据权利要求1-4中的任一项所述的系统，其中，所述局部效应转移函数是基于先前从压力相对于流率的采样测量结果学习的权重(we)而被线性组合的。
6.根据权利要求1-4中的任一项所述的系统，其中，所述图像是通过C型臂类型的平面X射线装置或旋转C型臂系统来采集的。
7.一种图像处理装置，包括：
用于接收(S505)包括感兴趣对象的投影视图的至少一幅图像的模块；
用于分割(S510)所述图像，以依据所述投影视图获得针对所述对象的对象分割的模块；
用于将所述分割划分(S515)为部分(j)的模块；
用于将多个局部效应转移函数中的一个或多个适配(S520)到在各个所述部分(j)中记录的所述对象的局部几何结构中，以获得多个全局效应转移函数(he)的模块；
用于将所述多个全局效应转移函数组合(S525)为针对所述对象的复合转移函数的模块；
用于根据所述复合转移函数，计算(S530)针对流经所述对象的流体的流量储备分数FFR估计的模块。
8.根据权利要求7所述的图像处理装置，其中，计算包括：通过改变至少一个物理或生理参数来计算FFR估计的范围，所述图像处理装置还包括：
用于显示(S535)根据所述至少一个物理或生理参数的FFR估计的所述范围的模块。
9.根据权利要求8所述的图像处理装置，所述参数是压力和流率中的至少一个，并且其中，根据压力和流率中的至少一个来显示FFR估计的所述范围。
10.一种计算机可读介质，其上存储有计算机程序单元，所述计算机程序单元用于控制根据权利要求1-6中的任一项所述的系统，当由处理单元执行时，所述计算机程序单元适用于：
接收(S505)包括感兴趣对象的投影视图的至少一幅图像；
分割(S510)所述图像，以依据所述投影视图获得针对所述对象的对象分割；
将所述分割划分(S515)为部分(j)；
将多个局部效应转移函数中的一个或多个适配(S520)到在各个所述部分(j)中记录的所述对象的局部几何结构中，以获得多个全局效应转移函数(he)；
将所述多个全局效应转移函数组合(S525)为针对所述对象的复合转移函数；
根据所述复合转移函数，计算(S530)针对流经所述对象的流体的流量储备分数FFR估计。
CN201480066240.XA 2013-12-04 2014-12-04 用于改善功能狭窄分析的局部ffr估计和可视化 Active CN105792738B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP13195675.7 2013-12-04
EP13195675 2013-12-04
PCT/EP2014/076477 WO2015082576A1 (en) 2013-12-04 2014-12-04 Local ffr estimation and visualisation for improved functional stenosis analysis

Publications (2)

Publication Number Publication Date
CN105792738A CN105792738A (zh) 2016-07-20
CN105792738B true CN105792738B (zh) 2020-03-13

Family

ID=49917425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480066240.XA Active CN105792738B (zh) 2013-12-04 2014-12-04 用于改善功能狭窄分析的局部ffr估计和可视化

Country Status (5)

Country Link
US (1) US10111633B2 (zh)
EP (1) EP3076854B1 (zh)
JP (1) JP6396468B2 (zh)
CN (1) CN105792738B (zh)
WO (1) WO2015082576A1 (zh)

Families Citing this family (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10210956B2 (en) 2012-10-24 2019-02-19 Cathworks Ltd. Diagnostically useful results in real time
EP3723041A1 (en) 2012-10-24 2020-10-14 CathWorks Ltd. Automated measurement system and method for coronary artery disease scoring
EP3061015A2 (en) 2013-10-24 2016-08-31 Cathworks Ltd. Vascular characteristic determination with correspondence modeling of a vascular tree
WO2016182508A1 (en) 2015-05-12 2016-11-17 Singapore Health Services Pte Ltd Medical image processing methods and systems
JP6738406B2 (ja) * 2015-08-05 2020-08-12 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 介入血行動態測定のための支援装置及び方法
US10898267B2 (en) * 2015-10-07 2021-01-26 Koninklijke Philips N.V. Mobile FFR simulation
IL263066B2 (en) 2016-05-16 2023-09-01 Cathworks Ltd Selecting blood vessels from images
EP4300419A3 (en) 2016-05-16 2024-04-03 Cathworks Ltd. System for vascular assessment
CN106168111B (zh) * 2016-08-01 2017-12-19 国网浙江省电力公司金华供电公司 等电位作业专用绝缘吊梯
JP7109429B2 (ja) * 2016-09-30 2022-07-29 コーニンクレッカ フィリップス エヌ ヴェ 狭窄評価用の機能的指標を決定する装置
US20180256038A1 (en) * 2017-03-08 2018-09-13 Medtronic Vascular, Inc. System, catheter, and method for calculating corrected fractional flow reserve
CN110612063B (zh) * 2017-04-06 2023-09-29 皇家飞利浦有限公司 标准化冠状动脉疾病度量
EP3644857B1 (en) * 2017-06-30 2025-02-26 Koninklijke Philips N.V. Machine learning spectral ffr-ct
EP3488774A1 (en) * 2017-11-23 2019-05-29 Koninklijke Philips N.V. Measurement guidance for coronary flow estimation from bernoulli´s principle
CN108186038B (zh) * 2018-02-11 2020-11-17 杭州脉流科技有限公司 基于动脉造影影像计算冠脉血流储备分数的系统
US10795054B2 (en) * 2018-03-20 2020-10-06 Mitsubishi Electric Research Laboratories, Inc. System and method for sensing wind flow passing over complex terrain
US20220125324A1 (en) * 2019-02-11 2022-04-28 University Of Louisville Research Foundation, Inc. System and method for determining a blood flow characteristic
CN109805949B (zh) * 2019-03-19 2020-05-22 苏州润迈德医疗科技有限公司 基于压力传感器和造影图像计算血流储备分数的方法
EP3948886A4 (en) 2019-04-01 2022-12-21 CathWorks Ltd. Methods and apparatus for angiographic image selection
JP7260218B2 (ja) * 2019-07-31 2023-04-18 ▲蘇▼州▲潤▼▲邁▼▲徳▼医▲療▼科技有限公司 造影画像間隔時間に基づいて血流速度を補正する方法及び装置
EP4027873A4 (en) * 2019-09-09 2023-09-20 Medhub Ltd AUTOMATED ANALYSIS OF IMAGE DATA TO DETERMINE A FRACTIONAL FLOW RESERVE
WO2021059165A1 (en) 2019-09-23 2021-04-01 Cathworks Ltd. Methods, apparatus, and system for synchronization between a three-dimensional vascular model and an imaging device
JP7369437B2 (ja) * 2019-09-27 2023-10-26 国立大学法人三重大学 評価システム、評価方法、学習方法、学習済みモデル、プログラム
CN111161342B (zh) * 2019-12-09 2023-08-29 杭州脉流科技有限公司 基于冠脉造影图像获取血流储备分数的方法、装置、设备、系统及可读存储介质
US12315076B1 (en) 2021-09-22 2025-05-27 Cathworks Ltd. Four-dimensional motion analysis of a patient's coronary arteries and myocardial wall

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3902765B2 (ja) * 2001-05-17 2007-04-11 シーメンス コーポレイト リサーチ インコーポレイテツド Ｍｒ心臓画像における左心室のセグメンテーションに対する様々なアプローチ
US8034352B2 (en) * 2005-04-06 2011-10-11 Ibc Pharmaceuticals, Inc. Tetrameric cytokines with improved biological activity
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
US10648918B2 (en) * 2011-08-03 2020-05-12 Lightlab Imaging, Inc. Systems, methods and apparatus for determining a fractional flow reserve (FFR) based on the minimum lumen area (MLA) and the constant
US9974508B2 (en) * 2011-09-01 2018-05-22 Ghassan S. Kassab Non-invasive systems and methods for determining fractional flow reserve
US10162932B2 (en) 2011-11-10 2018-12-25 Siemens Healthcare Gmbh Method and system for multi-scale anatomical and functional modeling of coronary circulation
US10034614B2 (en) 2012-02-29 2018-07-31 General Electric Company Fractional flow reserve estimation
US10373700B2 (en) * 2012-03-13 2019-08-06 Siemens Healthcare Gmbh Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance
CN108294735B (zh) * 2012-03-13 2021-09-07 西门子公司 用于冠状动脉狭窄的非侵入性功能评估的方法和系统
WO2013171644A1 (en) * 2012-05-14 2013-11-21 Koninklijke Philips N.V. Determination of a fractional flow reserve (ffr) value for a stenosis of a vessel
JP5953438B2 (ja) * 2012-12-11 2016-07-20 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 冠状動脈を通る血流量を決定する方法

Also Published As

Publication number Publication date
EP3076854B1 (en) 2022-04-20
US20160302750A1 (en) 2016-10-20
JP2017501775A (ja) 2017-01-19
JP6396468B2 (ja) 2018-09-26
WO2015082576A1 (en) 2015-06-11
EP3076854A1 (en) 2016-10-12
CN105792738A (zh) 2016-07-20
US10111633B2 (en) 2018-10-30

Similar Documents

Publication Publication Date Title
CN105792738B (zh) 用于改善功能狭窄分析的局部ffr估计和可视化
AU2020244586B2 (en) Method and system for sensitivity analysis in modeling blood flow characteristics
US10580526B2 (en) System and method for calculating vessel flow parameters based on angiography
JP7303260B2 (ja) 動脈網における流量および圧力勾配を患者特定コンピュータ断層撮影アルゴリズムに基づくコントラスト分布から判断するための方法
EP3160335B1 (en) Apparatus for determining a fractional flow reserve value
KR102172182B1 (ko) 시뮬레이션 정확도 및 성능을 위한 이미지 품질 평가
EP3127026B1 (en) Systems and methods for determining blood flow characteristics using flow ratio
CN108140430B (zh) 根据压力或流量测量结果及血管造影估计流量、阻力或压力
CN106537392A (zh) 用于冠状动脉中的血液动力学计算的方法和系统
JP2018537157A (ja) 非侵襲的血流予備量比（ｆｆｒ）に対する側副血流モデル化
EP3607527B1 (en) Quantitative evaluation of time-varying data
Tache Three Dimensional Reconstruction and Hemodynamic Information Extraction from Monoplane Angiography

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