US 11,896,416 B2
Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images
Yunfei Huo, Suzhou (CN); Guangzhi Liu, Suzhou (CN); Xingyun Wu, Suzhou (CN); and Zhiyuan Wang, Suzhou (CN)
Assigned to SUZHOU RAINMED MEDICAL TECHNOLOGY CO., LTD., Suzhou (CN)
Filed by SUZHOU RAINMED MEDICAL TECHNOLOGY CO., LTD., Suzhou (CN)
Filed on May 24, 2021, as Appl. No. 17/328,550.
Application 17/328,550 is a continuation of application No. PCT/CN2019/071203, filed on Jan. 10, 2019.
Claims priority of application No. 201811438744.6 (CN), filed on Nov. 28, 2018.
Prior Publication US 2021/0275124 A1, Sep. 9, 2021
This patent is subject to a terminal disclaimer.
Int. Cl. G06T 7/155 (2017.01); A61B 6/00 (2006.01); G16H 50/30 (2018.01); G16H 30/40 (2018.01); A61B 6/03 (2006.01); G06T 11/00 (2006.01)
CPC A61B 6/507 (2013.01) [A61B 6/032 (2013.01); A61B 6/5217 (2013.01); A61B 6/5235 (2013.01); A61B 6/5247 (2013.01); G06T 7/155 (2017.01); G06T 11/008 (2013.01); G16H 30/40 (2018.01); G16H 50/30 (2018.01); G06T 2207/10081 (2013.01); G06T 2207/20156 (2013.01); G06T 2207/30104 (2013.01); G06T 2211/412 (2013.01)] 9 Claims
OG exemplary drawing
 
1. A method for calculating coronary artery fractional flow reserve based on myocardial blood flow and CT image, comprising the following steps:
S01: segmenting the CT image of heart, obtaining an image of heart via a morphological operation, subjecting the image of heart to a histogram analysis to obtain an image contains ventricular and atrial, obtaining a myocardial image by making a difference between the image of heart and the image contains ventricular and atrial, determining a myocardial volume by the myocardial image;
S02: obtaining a full aortic complementary image by processing an aortic image, obtaining the aortic image containing a coronary artery inlet through regional growth of the full aortic complementary image, and obtaining an image containing the coronary artery inlet according to the aortic image containing the coronary artery inlet and the full aortic complementary image to determine the coronary artery inlet by the image containing the coronary artery inlet;
S03: extracting a coronary artery through regional growth by taking the coronary artery inlet as a seed point on the myocardial image, calculating an average gray and average variance of the coronary artery, and along a direction of the coronary artery, extracting a coronary artery tree according to a gray distribution of the coronary artery;
S04: binarizing the coronary artery image, drawing an isosurface image to obtain a three-dimensional grid image of the coronary artery;
S05: calculating a total flow at the coronary artery inlet in a maximum hyperemia state, Qtotal=myocardial volume×myocardial blood flow×CFR, and CFR being the coronary flow reserve;
S06: calculating a blood flow velocity V1 in a hyperemia state;
S07: calculating a pressure drop ΔP from the coronary artery inlet to a distal end of a coronary artery stenosis using V1 as an inlet flow velocity of coronary artery stenosis blood vessel, using Pd=Pa−ΔP calculating a mean intracoronary pressure Pd at the distal end of the coronary artery stenosis, wherein Pa is a mean aortic pressure, and obtaining the fractional flow reserve using FFR=Pd/Pa.