CN113693579A - 归一化冠状动脉微循环阻力指数计算方法 - Google Patents

归一化冠状动脉微循环阻力指数计算方法 Download PDF

Info

Publication number
CN113693579A
CN113693579A CN202110837710.XA CN202110837710A CN113693579A CN 113693579 A CN113693579 A CN 113693579A CN 202110837710 A CN202110837710 A CN 202110837710A CN 113693579 A CN113693579 A CN 113693579A
Authority
CN
China
Prior art keywords
normalized
calculating
coronary artery
coronary
blood flow
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.)
Pending
Application number
CN202110837710.XA
Other languages
English (en)
Inventor
谢辛舟
廉坤
陈焱
李成祥
张瑞晨
赵帅
郭晓雄
李玮昊
谢松云
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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110837710.XA priority Critical patent/CN113693579A/zh
Publication of CN113693579A publication Critical patent/CN113693579A/zh
Pending legal-status Critical Current

Links

Images

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
    • 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/021—Measuring pressure in heart or blood vessels
    • A61B5/0215—Measuring pressure in heart or blood vessels by means inserted into the body
    • 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/021—Measuring pressure in heart or blood vessels
    • A61B5/0215—Measuring pressure in heart or blood vessels by means inserted into the body
    • A61B5/02158—Measuring pressure in heart or blood vessels by means inserted into the body provided with two or more sensor elements
    • 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/7235—Details of waveform analysis
    • A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Abstract

本发明提出了归一化冠状动脉微循环阻力指数计算方法，包括：测量感兴趣血管最大充血状态下的远、近端压力值；基于造影图像重建感兴趣血管三维模型，同时提取血管远端截面面积；仿真计算感兴趣血管血流阻力模型参数；结合前述压力值和血流阻力模型参数计算对应血流量；基于前述压力值、血流量和血管远端截面面积计算归一化冠状动脉微循环阻力指数。本发明基于造影图像和侵入式测量的压力值，提出归一化冠状动脉微循环阻力指数计算方法，以使得可在冠状动脉不同部位、不同大小的血管间建立统一的微循环阻力计算指标。

Description

归一化冠状动脉微循环阻力指数计算方法
技术领域
本发明涉及计算冠状动脉生理学领域，尤其涉及应用血流动力学仿真计算方法计算冠状动脉微循环阻力的方法。
背景技术
本发明是基于本人已授权专利(“基于冠状动脉CT血管造影的狭窄病灶血流阻力计算方法”，ZL201710505979.1)的专利发明的发明。
冠状动脉微循环阻力(microcirculatory resistance,MR)是反映冠状动脉微循环功能的有效指标，其定义为冠脉远端动脉压力(Pd)除以最大充血状态下的冠脉血流量(Qmax)。目前可通过以下方法计算：
(1)第一类方法通过集成了压力传感器和血流量传感器的导丝置入冠状动脉远端，在药物诱导最大充血状态的条件下同时测量压力和血流量，进而计算MR(Williams,R.P.,et al.(2018)."Doppler Versus Thermodilution-Derived CoronaryMicrovascular Resistance to Predict Coronary Microvascular Dysfunction inPatients With Acute Myocardial Infarction or Stable Angina Pectoris."Am JCardiol 121(1):1-8.)；
(2)通过序列造影图像估算Qmax，基于血流动力学模型仿真计算最大充血状态下远、近端压力差(ΔP)，基于测量的静息态冠脉入口压力和深度学习模型估算最大充血状态下冠脉入口压力(Pa)(“快速计算微循环阻力的方法与系统”，中国201711258493.9[P])或直接测量最大充血状态下Pa(“基于造影图像和流体力学模型的微循环阻力指数计算方法”，中国201810413391.8[P])；
(3)结合侵入式压力导丝测量的充血态冠状动脉压力值和基于多角度冠脉造影图像重建的冠状动脉三维模型，采用计算流体力学方法仿真计算最大充血状态下的冠脉血流量(Qmax)(Morris PD,Gosling R,Zwierzak I,Evans H,Aubiniere-Robb L,Czechowicz K,et al.A novel method for measuring absolute coronary blood flow andmicrovascular resistance in patients with ischaemic heart disease.CardiovascRes.2021；117:1567-77.)。
但上述各种方法计算的微循环阻力均为绝对微循环阻力值，其在不同部位、不同大小的血管间数值分布差异巨大，难以直接为临床诊断冠状动脉微循环障碍提供精确的参考。
发明内容
有鉴于此，本发明提供一种归一化冠状动脉微循环阻力指数计算方法，该方法将侵入式压力导丝测量的远、近端压力值作为已知条件，基于多角度冠脉造影图像重建冠状动脉三维模型，利用血流阻力模型仿真计算冠脉血流量，并进一步基于压力值和重建三维模型中获取的血管远端截面面积，计算归一化微循环阻力指数。所采取的技术方案如下：
1.侵入式压力导丝测量感兴趣血管的远、近端压力值；
2.基于造影图像重建感兴趣血管三维模型，提取血管远端截面面积A；
3.通过数值方法计算感兴趣血管段血流阻力(压降-血流曲线)模型参数；
4.结合测量压力值和血流阻力模型参数计算冠脉血流量(Qmax)；
5.基于血管远端截面面积A和计算得到的Qmax，计算归一化微循环阻力指数MRnomalized。
附图说明
图1.归一化冠状动脉微循环阻力指数计算方法的整体流程图。
具体实施方式
下面结合附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
归一化冠状动脉微循环阻力指数计算方法的整体流程图如图1所示。以下将结合图1对具体实施方式进行详细说明。
1.测量感兴趣血管最大充血状态下的远、近端压力值：通过药物诱导最大充血状态，利用压力导丝测量感兴趣血管远、近端压力值(Pa、Pd)；
2.基于造影图像重建感兴趣血管三维模型：多角度拍摄的冠脉造影图像进行分割，分别得到冠脉中心线及直径；通过三维空间投影计算生成冠脉三维模型，同时提取血管远端截面面积A；
3.基于重建的冠脉三维模型，采用本人已授权的专利发明方法(“基于冠状动脉CT血管造影的狭窄病灶血流阻力计算方法”，ZL201710505979.1)计算血流阻力(压降-血流曲线)模型参数f和s；
4.基于f、s和前述测量压力值Pa、Pd，求解如下一元二次方程，得到Qmax：
Pa-Pd＝f·Qmax+s·Qmax 2 (1)
5.基于如下公式计算归一化冠状动脉微循环阻力指数：
Figure BDA0003177768640000021
其中MRnomalized为归一化冠状动脉微循环阻力指数，Pd为前述1步骤中测量的远端压力Pd，Qmax为上述4步骤中计算的最大充血血流量，k为固定的模型参数，典型取值在2到3之间，requ为等效血管半径，其基于步骤2中提取的血管远端截面面积A按如下公式计算得到：
Figure BDA0003177768640000022
其中π为圆周率。
本发明的有益效果在于：
基于造影图像和侵入式测量的压力值，提供了归一化的微循环阻力指数计算方法，以使得可在冠状动脉不同部位、不同大小的血管间建立统一的微循环阻力计算指标。

Claims (3)

1.归一化冠状动脉微循环阻力指数计算方法，包括如下步骤：
(1)测量感兴趣血管最大充血状态下的远、近端压力值：通过药物诱导最大充血状态，利用压力导丝测量感兴趣血管远、近端压力值(Pa、Pd)；
(2)基于造影图像重建感兴趣血管三维模型：多角度拍摄的冠脉造影图像进行分割，分别得到冠脉中心线及直径；通过三维空间投影计算生成冠脉三维模型，同时提取血管远端截面面积A；
(3)基于重建的冠脉三维模型，采用本人已授权的专利发明方法(“基于冠状动脉CT血管造影的狭窄病灶血流阻力计算方法”，ZL201710505979.1)计算血流阻力(压降-血流曲线)模型参数f和s；
(4)基于f、s和前述测量压力值Pa、Pd，求解如下一元二次方程，得到Qmax：
Pa-Pd＝f·Qmax+s·Qmax 2 (1)
(5)基于如下公式计算归一化冠状动脉微循环阻力指数：
Figure FDA0003177768630000011
其中MRnomalized为归一化冠状动脉微循环阻力指数，Pd为前述1步骤中测量的远端压力Pd，Qmax为上述4步骤中计算的最大充血血流量，k为固定的模型参数，典型取值在2到3之间，requ为等效血管半径，其基于步骤2中提取的血管远端截面面积A按如下公式计算得到：
Figure FDA0003177768630000012
其中π为圆周率。
2.如权利要求1的归一化冠状动脉微循环阻力指数计算方法，其特征在于：基于造影图像重建的三维模型，将侵入式压力导丝测量的远、近端压力值作为已知条件(Pa、Pd)，采用血流阻力模型仿真计算对应血流量Qmax。
3.如权利要求1的归一化冠状动脉微循环阻力指数计算方法，其特征在于：基于血管远端尺寸(等效血管半径)，对绝对微循环阻力进行归一化，得到归一化冠状动脉微循环阻力指数。
CN202110837710.XA 2021-07-23 2021-07-23 归一化冠状动脉微循环阻力指数计算方法 Pending CN113693579A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110837710.XA CN113693579A (zh) 2021-07-23 2021-07-23 归一化冠状动脉微循环阻力指数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110837710.XA CN113693579A (zh) 2021-07-23 2021-07-23 归一化冠状动脉微循环阻力指数计算方法

Publications (1)

Publication Number Publication Date
CN113693579A true CN113693579A (zh) 2021-11-26

Family

ID=78650348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110837710.XA Pending CN113693579A (zh) 2021-07-23 2021-07-23 归一化冠状动脉微循环阻力指数计算方法

Country Status (1)

Country Link
CN (1) CN113693579A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117481806A (zh) * 2022-12-30 2024-02-02 杭州脉流科技有限公司 基于ct图像的冠状动脉微循环阻力指数计算方法和装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107411767A (zh) * 2017-06-28 2017-12-01 西北工业大学 一种基于冠状动脉ct血管造影评估狭窄病灶血流阻力的非侵入式方法
CN107978371A (zh) * 2017-11-30 2018-05-01 博动医学影像科技（上海）有限公司 快速计算微循环阻力的方法及系统
WO2019210553A1 (zh) * 2018-05-03 2019-11-07 苏州润迈德医疗科技有限公司 基于造影图像和流体力学模型的微循环阻力指数计算方法
CN110520037A (zh) * 2017-02-28 2019-11-29 皇家飞利浦有限公司 使用差压原理和扩展流量传感器进行血管内血流测量
CN111067494A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110520037A (zh) * 2017-02-28 2019-11-29 皇家飞利浦有限公司 使用差压原理和扩展流量传感器进行血管内血流测量
CN107411767A (zh) * 2017-06-28 2017-12-01 西北工业大学 一种基于冠状动脉ct血管造影评估狭窄病灶血流阻力的非侵入式方法
CN107978371A (zh) * 2017-11-30 2018-05-01 博动医学影像科技（上海）有限公司 快速计算微循环阻力的方法及系统
WO2019210553A1 (zh) * 2018-05-03 2019-11-07 苏州润迈德医疗科技有限公司 基于造影图像和流体力学模型的微循环阻力指数计算方法
CN111067494A (zh) * 2019-12-27 2020-04-28 西北工业大学 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
骆秉铨, 陈兴新, 张义勤, 夏项, 王人鹏, 石怀岭, 李春梅, 路雯: "冠脉造影脉冲响应分析评价心肌微循环", 中国微循环, no. 02, 28 April 2003 (2003-04-28) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117481806A (zh) * 2022-12-30 2024-02-02 杭州脉流科技有限公司 基于ct图像的冠状动脉微循环阻力指数计算方法和装置
CN117481806B (zh) * 2022-12-30 2024-04-02 杭州脉流科技有限公司 基于ct图像的冠状动脉微循环阻力指数计算方法和装置

Similar Documents

Publication Publication Date Title
CN107978371B (zh) 快速计算微循环阻力的方法及系统
WO2017097073A1 (zh) 血管压力差与血流储备分数的计算方法及系统
Wentland et al. Aortic pulse wave velocity measurements with undersampled 4D flow‐sensitive MRI: comparison with 2D and algorithm determination
US10354744B2 (en) Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance
Bensalah et al. Geometry is a major determinant of flow reversal in proximal aorta
US20160000397A1 (en) Method for assessing stenosis severity through stenosis mapping
CN111067494B (zh) 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法
US11626211B2 (en) Modelling blood vessels and blood flow
CN112384137B (zh) 基于造影图像获取静息态下血管评定参数的方法及装置
CN109152563A (zh) 流体流分析
Vijayan et al. Assessing coronary blood flow physiology in the cardiac catheterisation laboratory
WO2013110929A1 (en) Aortic pulse wave velocity measurement
WO2015089222A2 (en) Method for post-processing flow-sensitive phase contrast magnetic resonance images
EP3836146A1 (en) Method and apparatus for simulating blood flow of subject-specific blood vessel
Liu et al. Physiologically personalized coronary blood flow model to improve the estimation of noninvasive fractional flow reserve
CN111091913A (zh) 基于血流储备分数和冠状动脉ct造影图像的微循环阻力计算方法
CN113693579A (zh) 归一化冠状动脉微循环阻力指数计算方法
Khalifé et al. Estimating absolute aortic pressure using MRI and a one-dimensional model
Kim et al. Quantification of wave reflection using peripheral blood pressure waveforms
Bianchini et al. Magnetic resonance imaging and computed tomography for the noninvasive assessment of arterial aging: a review by the VascAgeNet COST action
CN112384138B (zh) 获取心表大动脉的血流量的方法、装置、系统及存储介质
Sonnabend et al. Cerebrospinal fluid pulse wave velocity measurements: In vitro and in vivo evaluation of a novel multiband cine phase‐contrast MRI sequence
CN112704505B (zh) 一种利用cta和dsa测量冠状动脉血流储备分数的方法
US10332255B2 (en) Method for assessing stenosis severity in a lesion tree through stenosis mapping
KR101623186B1 (ko) 생리학적 압력-유량 관계를 이용한 관상동맥 유량과 압력을 구하는 방법

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