JP7236768B2 - Calculation method of coronary flow reserve ratio based on myocardial blood flow and CT image - Google Patents

Description

The present invention relates to the field of coronary artery imaging assessment, and more particularly to a method for calculating coronary flow reserve ratio based on myocardial blood flow and CT images.

Both coronary angiography and intravascular ultrasound are considered “golden standards” for diagnosing coronary heart disease. I don't know how big an impact it will have. Fractional coronary flow reserve (FFR) has already become an official index for assessing coronary artery stenosis functionality, and its most important function is to accurately assess the functional consequences of coronary artery stenosis, whose effects are unknown. is to do

Coronary flow reserve ratio (FFR) is the theoretical maximum blood flow obtained in the same region as the maximum blood flow obtained for the myocardial region of the target measurement vessel when there is a stenotic lesion in the coronary artery. refers to the ratio of FFR is primarily obtained by calculating the ratio of distal coronary stenosis pressure to aortic root pressure. Distal stenotic pressure is obtained by pressure guidewire measurement at maximal perfusion blood flow (intracoronary or intravenous injection of papaverine or adenosine or ATP). It can be simplified as the ratio of the stenotic distal intracoronary mean pressure (P _d ) to the coronary ostial aortic mean pressure (P _a ) at maximal myocardial hyperemia, ie, FFR=P _d /P _a .

Coronary artery CTA is a non-invasive, easy-to-operate coronary artery lesion diagnostic examination method that can accurately evaluate the degree of coronary artery stenosis and distinguish the characteristics of blood vessel wall plaque, and screen high-risk individuals. be the preferred method. Therefore, when intervening in a patient with coronary heart disease, CTA evaluation must be performed on the patient's coronary artery at an early stage.

When calculating the FFR obtained non-invasively by coronary CTA (CTFFR), no additional imaging tests or drugs are required, essentially avoiding unnecessary coronary angiography and revascularization treatment. As clearly demonstrated in the DeFacto study results, in coronary CT, analysis of CTFFR results provides physiologic information of true blood flow-limiting lesions that increase patient risk. CTFFR combines the advantages of coronary CTA and FFR to assess coronary artery stenosis both structurally and functionally, providing a novel non-invasive detection system for coronary lesion anatomic and functional information. . However, CTA cannot measure the coronary artery flow velocity in hyperemic conditions and can only be predicted by relying on numerical methods, which greatly limits the clinical application of CTFFR.

In order to solve the above technical problems, the present invention identifies resting-state myocardial blood flow and coronary flow reserve (CFR) by non-invasive measurement, and measures the maximal hyperemia of different vessel linings in the coronary artery tree. By specifying the flow rate and specifying the maximum hyperemic flow velocity V ₁ , the coronary flow reserve ratio FFR can be obtained quickly, accurately, and fully automatically. It is an object of the present invention to provide a method for calculating coronary flow reserve ratio.

The technical solution of the present invention is as follows:
A method for calculating coronary coronary flow reserve ratio based on myocardial blood flow and CT images includes the following steps:
S01: Divide a cardiac CT image, obtain a cardiac image by morphological manipulation, histogram-analyze the cardiac image to obtain a ventricular-atrial image, obtain a myocardial image by difference between the cardiac image and the ventricular-atrial image, and obtain a myocardial image. identify the myocardial volume based on the image;
S02: Processing the aorta image to obtain a total aorta complementary image, performing region growing to obtain an aorta image including the coronary ostia, and based on the aorta image including the coronary ostia and the total aorta complementary image, the coronary ostia obtaining an image containing the ostium of the coronary artery, identifying the ostium of the coronary artery based on the image containing the ostium of the coronary artery;
S03: Using the coronary artery ostium as a seed point in the myocardial image, extracting the coronary artery by region growing, calculating the mean grayscale and mean variance of the coronary artery, and constructing the coronary artery tree along the coronary artery direction based on the coronary artery grayscale distribution. and extract
S04: Binarize the coronary artery image, draw an isosurface image, obtain a three-dimensional grid image of the coronary artery,
S05: Calculate the total flow Q _total at the coronary artery ostium in the state of maximal hyperemia = myocardial volume × myocardial blood flow × CFR, where CFR is the coronary flow reserve (CFR);
S06: Calculate the blood flow velocity _V1 in the congested state,
S07: Let V1 be the inlet flow velocity of the coronary artery stenosis, calculate the pressure loss ΔP from the coronary artery inlet to the distal coronary artery stenosis, and the mean intracoronary artery pressure distal to the stenosis P _d = P _a - ΔP, where P _a is the aorta is the mean pressure, yielding the coronary flow reserve ratio FFR=P _d /P _a .

In the preferred technical solution, in step S02, after obtaining the image containing the ostia of the coronary arteries, perform connected region analysis on the image containing the ostia of the coronary arteries, mark each connected region with a different grayscale label, and identify the coronary artery Identify mouth.

In the preferred technical solution, in the above step S02, the ascending aorta and the central line are extracted by using the characteristic that the cross section of the aorta is circular in the heart image to obtain the aorta image.

In the preferred technical solution, the coronary artery image binarization in step S04 traverses voxels in the coronary artery image V1, if a voxel pixel is equal to 0, the pixel value remains unchanged; , set the pixel value to 1 and obtain one new data V2.

In a preferred technical solution, in step S05, by myocardial contrast echocardiography (MCE) or single photon emission tomography (SPECT) or positron emission tomography (PET) or cardiac nuclear magnetic resonance imaging (MRI) or CT perfusion , resting-state myocardial blood flow and coronary flow reserve (CFR).

In the preferred technical solution, the step S06 includes:
S61: Based on the flow-volume scaling law and the cardiac surface coronary artery tree of cardiac CT three-dimensional reconstruction, the blood flow in any one vessel in the tree Q = Q _total × (V/V _total ) ^3/4 where V _total is the sum of vessel volumes of all superficial coronary arteries in the cardiac CT three-dimensional reconstruction, and V is the volume of any one vessel in the superficial coronary artery tree and its downstream vessels is the sum of vessel volumes ,
S62: Based on the flow-volume scaling law and the cardiac surface coronary artery tree of cardiac CT three-dimensional reconstruction, identify the blood flow velocity V1=Q/D in any one vessel in the tree, where D is the is the average cross-sectional area of blood vessels.

但し、

Ｐ、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ₁であり、
三次元計算流体力学により各冠動脈狭窄の圧力損失ΔＰ₁、ΔＰ₂、ΔＰ₃・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰ_i（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐ_d＝Ｐ_a－ΔＰを計算する、但し、Ｐ_aは大動脈平均であることを含む。 In the preferred technical solution, the step S07 is specifically:
Solve for the vessel 3D grid and solve the continuity and Navier-Stokes equations using numerical methods:

however,

P, ρ, μ are flow velocity, pressure, blood density, blood viscosity, respectively,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia,
By three-dimensional computational fluid dynamics, the pressure losses ΔP ₁ , ΔP ₂ , _{ΔP 3} , . ), calculate the mean stenotic distal coronary pressure P _d =P _a -ΔP, where P _a is the aortic mean.

但し、ρは血液の密度を表し、ｕ_z、ｕ_rはそれぞれｚ方向、ｒ方向の流速を表し、μは血液の動粘度を表し、ｐは血液の圧の強さを表し、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ₁であり、
二次元流体力学により各冠動脈狭窄の圧力損失ΔＰ₁、ΔＰ₂、ΔＰ₃・・・、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰ_i（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐ_d＝Ｐ_a－ΔＰを計算する、但し、Ｐ_aは大動脈平均圧である、ことを含む。 In the preferred technical solution, the step S07 includes:
Straighten stenosed vessels, build a two-dimensional axisymmetric model, divide the two-dimensional grid, and solve the continuity and Navier-Stokes equations using numerical methods based on the geometric structure of the CT reconstruction. :

where ρ represents the density of blood, u _z and u _r represent the flow velocities in the z and r directions, respectively, μ represents the kinematic viscosity of blood, p represents the strength of blood pressure,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia,
Pressure loss ΔP ₁ , ΔP ₂ , ΔP _{3 ,} ΔP 1 , ΔP 2 , ΔP ₃ , . , to calculate the stenotic distal intracoronary mean pressure P _d =P _a -ΔP, where P _a is the aortic mean pressure.

In the preferred technical solution, the step S07 further includes calculating the pressure difference from the inlet to the outlet in the three-dimensional model for different types of curvature of the blood vessel, calculating it in the light of the two-dimensional axisymmetric model, and calculating each type of curvature. Build a database for storing correction factors for the two-dimensional axisymmetric result of the curve, obtain the pressure, then look up the correction factors in the database to obtain the corrected inlet-to-exit pressure difference, and then FFR including computing

The present invention can obtain the coronary flow reserve ratio FFR quickly, accurately, and fully automatically from the myocardial blood flow and cardiac CT images, and can greatly improve the accuracy of the conventional CTFFR (or FFRCT). Non-invasive measurement is easy to operate, greatly reduces the difficulty and risk of surgery, and can be applied clinically on a large scale.

The present invention will be further described in conjunction with the drawings and examples below.

FIG. 1 is a method flowchart of the present invention. FIG. 2 shows the myocardial segmentation result of a cardiac CT image. FIG. 3 is an aortic segmentation result with coronary ostia. FIG. 4 shows coronary artery ostial segmentation results. FIG. 5 shows coronary artery segmentation results. FIG. 6 is a grid model of coronary artery segmentation results. FIG. 7 is a schematic diagram of heart and coronary artery blood flow.

In order to more clearly describe the purpose, technical solution and advantages of the present invention, the present invention will be described in more detail below together with the detailed description of the invention with reference to the accompanying drawings. It should be understood that these descriptions are examples only and do not limit the scope of the invention. In addition, in the following description, descriptions of known structures and techniques are omitted to avoid unnecessary confusion with the concepts of the present invention.

For a predetermined cardiac CT image, the heart is extracted by the reverse method, and the non-target areas such as the descending aorta, spine, and ribs are processed, and non-target areas such as the chest wall, lungs, vertebrae, and descending aorta are extracted one by one. A heart image is obtained by removing the heart tissue. In the obtained cardiac image, the ascending aorta and the central line are extracted using the feature that the cross section of the aorta is circular to obtain the aorta image.

As shown in FIG. 1, the calculation method of the coronary artery flow reserve ratio (FFR) based on the myocardial blood flow and CT image of the present invention extracts the myocardial image, extracts the coronary ostium, and extracts the coronary artery to generate a coronary artery grid model, determine the resting-state myocardial flow and the coronary flow reserve (CFR), calculate the total flow at the coronary ostium at maximal hyperemia, Calculating the flow velocity V ₁ and identifying the coronary artery FFR. It specifically includes the following steps:

1: Extract myocardial images:
A cardiac CT image is segmented, a cardiac image is obtained by morphological manipulation, the cardiac image is histogram-analyzed to obtain a ventricular-atrial image, and a difference between the cardiac image and the ventricular-atrial image is obtained to obtain a myocardial image. See Figure 2.

2: Extract the coronary ostia:
Morphological dilation is performed on the binarized image of the aorta to obtain a binary image of the whole aorta, and pixel reversal to obtain a complemented whole aorta image.
Region growing is performed based on the mean grayscale of points on the aortic centerline to obtain an aortic image that includes the coronary ostia. See FIG.
Perform image multiplication by the aorta image including the coronary ostia and the total aorta complementary image to obtain the image including the coronary ostia, perform connected region analysis on the image including the coronary ostia, and label each connection with a different grayscale label. Mark the area and identify the coronary ostia. See FIG.

3: Extract the coronary arteries:
Using the coronary artery ostia as seed points in the myocardial image, extracting the coronary arteries by region growing, calculating the mean grayscale and mean variance of the coronary arteries, and extracting the coronary artery tree along the coronary artery direction based on the coronary artery grayscale distribution. do. See FIG.

4: Generate a coronary artery grid model:
In step 3, the coronary artery image data V1 is obtained, the voxels in the data constitute one cube in space, and the voxel pixel values belonging to the coronary artery are not 0 (pixel values of about -3000 to 3000). between), the remaining voxel pixel values are all zero.
In this step, it is necessary to transform the data into spatial three-dimensional grid data V3 for the FFR calculation in step 5.

(1) Coronary artery data binarization Traversal the voxels in the coronary artery image data V1 and perform a simple pixel value judgment, if pixel A1 is equal to 0, the pixel value is unchanged; If not, let the pixel value of A1 be 1.
Finally, one new image data V2 is obtained, in which the voxel pixel value belonging to the coronary artery portion is 1 and the rest is 0.

(2) Isosurface generation A voxel is defined as a minimal hexahedron, and four pixels between adjacent layers form eight vertices on the cube. An isosurface is the set of all points in space that have the same attribute value. This can be expressed as:
{(x,y,z)|f(x,y,z)=c}, where c is a constant c in this method is the pixel value 1 given in the 3D reconstruction process.

The flow for extracting isosurfaces is as follows:
(1) After preprocessing the raw data, it is read into a specific array.
(2) Extracting one unit body from the lattice data body to make it the current unit body, and simultaneously obtaining all the information of the unit body.
(3) Compare the function values of the eight vertices of the current unit body with the predetermined isosurface value C to obtain the state table of the unit body.
(4) From the state table index of the current unit body, find the unit body edge that intersects the isosurface, and calculate the position coordinates of each intersection by the method of linear interpolation.
(5) Find the normal vectors of the eight vertices of the current unit body by the central difference method, then obtain the normal of each vertex of the triangular patch by the method of linear interpolation.
(6) An isosurface image is created based on the vertex coordinates and vertex normal vectors of each triangular patch.
Finally, three-dimensional lattice image data V3 of the coronary arteries are obtained. See FIG.

5: Calculate the blood flow velocity V ₁ in the hyperemic state:
Resting-state myocardium by non-invasive measurements such as myocardial contrast echocardiography (MCE) or single photon emission tomography (SPECT) or positron emission tomography (PET) or magnetic resonance imaging (MRI) or CT perfusion Blood flow and coronary coronary flow reserve (CFR) are determined, and myocardial volume, myocardial blood flow, and CFR are the total flow of the coronary ostium (including the sum of the left and right coronary trees) at maximal hyperemia. Calculate Q _total = myocardial volume x myocardial blood flow x CFR.

Based on the flow-volume scaling law and the cardiac surface coronary artery tree of cardiac CT three-dimensional reconstruction, the blood flow Q in any one vessel in the tree: Q = Q _total × (V/V _total ) ^3/4 where V _total is the sum of the vascular volumes of all superficial coronary arteries (including the sum of the left and right coronary trees) in the cardiac CT three-dimensional reconstruction, and V is is the sum of vessel volumes in any one vessel of and its downstream vessels. See FIG. Based on the flow volume scaling law and the cardiac surface coronary artery tree of cardiac CT three-dimensional reconstruction, the blood flow velocity in any one vessel in a particular tree V ₁ : V ₁ =Q/D, where D is the It is the mean cross-sectional area of a vessel (the vessel's vascular volume divided by the vessel's length).

6: Calculation of coronary artery FFR:
Letting V ₁ be the flow velocity at the entrance of the coronary artery stenosis, the thickness loss ΔP ₁ , ΔP ₂ , ΔP ₃ etc. of each coronary artery stenosis is calculated by computational fluid dynamics (CFD) method, and the pressure loss ΔP from the coronary artery entrance to the distal coronary artery stenosis ΔP = ΣΔP _i ( _i = 1, 2, 3...), calculate the mean pressure in the distal stenotic coronary artery P _d = P _{a -} ΔP, where P _a is the mean aortic pressure, and finally the formula FFR = P Calculate the coronary flow reserve ratio as _d / _Pa .

但し、

Ｐ、ρ、μはそれぞれ流速、圧力、血液密度、血液粘度であり、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ₁であり、
式［Ａ１］及び［Ａ２］に基づき、三次元ＣＦＤを実行して各冠動脈狭窄の圧力損失ΔＰ₁、ΔＰ₂、ΔＰ₃等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰ_i（_i＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐ_d＝Ｐ_a－ΔＰを計算し、但し、Ｐ_aは大動脈平均圧である。 Processing steps for 3D models include:
Based on the geometrical structure of the CT reconstruction, the 3D grid is partitioned and the continuity and Navier-Stokes equations are solved by numerical methods (e.g., finite difference, finite element method, finite volume method, etc.). :

however,

P, ρ, μ are flow velocity, pressure, blood density, blood viscosity, respectively,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia,
Based on equations [A1] and [A2], three-dimensional CFD was performed to determine the pressure drop ΔP ₁ , ΔP ₂ , ΔP _{3 ,} etc. of each coronary artery stenosis, and the pressure loss ΔP from the coronary ostium to the distal coronary artery stenosis ΔP = ΣΔP _i ( _i = 1, 2, 3...), calculate the stenotic distal intracoronary mean pressure P _d = P _a -ΔP, where P _a is the aortic mean pressure.

但し、ρは血液の密度を表し、ｕ_z、ｕ_rはそれぞれｚ方向、ｒ方向の流速を表し、μは血液の動粘度を表し、ｐは血液の圧の強さを表し、
入口境界条件は、最大充血状態における冠動脈狭窄血管の入口流速Ｖ₁であり、
式［Ａ３］～［Ａ５］に基づき、二次元ＣＦＤを実行し、各冠動脈狭窄の圧力損失ΔＰ₁、ΔＰ₂、ΔＰ₃等、冠動脈入口から冠動脈狭窄遠位までの圧力損失ΔＰ＝ΣΔＰ_i（ｉ＝１、２、３・・・）、狭窄遠位冠状動脈内平均圧Ｐ_d＝Ｐ_a－ΔＰを計算し、但し、Ｐ_aは大動脈平均圧である。 For 2D models, it includes the following steps:
Based on the geometric structure of the CT reconstruction, straighten the stenosed vessels (2D axisymmetric model), divide the 2D grid, and use numerical methods (e.g., finite difference, finite element method, finite volume method, etc.) Find continuity and solutions to the Navier-Stokes equation:

where ρ represents the density of blood, u _z and u _r represent the flow velocities in the z and r directions, respectively, μ represents the kinematic viscosity of blood, p represents the strength of blood pressure,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia,
Based on the equations [A3] to [A5], two-dimensional CFD is performed and the pressure loss ΔP ₁ , ΔP ₂ , ΔP ₃ etc. of each coronary artery stenosis, and the pressure loss ΔP from the coronary ostium to the distal coronary artery stenosis ΔP = ΣΔP _i ( i ₌ 1 _, 2 _, 3, .

For different types of curvature of the vessel, the pressure difference from the inlet to the outlet is calculated in the 3D model, calculated against the 2D axisymmetric model, and the correction factor for the 2D axisymmetric result of each type of curvature is calculated. A database is constructed for storage, and after calculating the pressure, the correction factor in the database is used to obtain the corrected inlet-to-outlet pressure difference, and finally the calculated FFR is calculated by the formula.

It is to be understood that the above detailed description of the invention is merely illustrative for explaining or interpreting the principles of the invention and is not intended to be limiting of the invention. Therefore, any modification, equivalent replacement, improvement, etc., without departing from the spirit and scope of the present invention, shall fall within the protection scope of the present invention. In addition, the claims appended hereto are intended to cover all changes and modifications in the scope and boundaries of the claims appended hereto or in equivalent forms of those scopes and boundaries.

Claims (7)

S01: Segment a CT image including the heart, obtain a cardiac image including only the heart by morphological manipulation, histogram-analyze the cardiac image to obtain a ventricular- atrial image including only the ventricle-atrium, and combine the cardiac image with the above - mentioned obtaining an image containing the myocardium by subtracting the ventricular-atrial images, and identifying a myocardial volume based on the image containing the myocardium;
S02: Obtain a cardiac aorta image from the cardiac image, perform morphological dilation on the binarized image of the aortic image, obtain a binary image of the entire cardiac aorta, and perform pixel reversal to obtain the entire cardiac aorta image. An aorta complementary image is obtained, region growing is performed on the aorta image based on the average gray scale of points on the centerline of the aorta , an aorta image including the coronary ostium is obtained, and the aorta image including the coronary ostium and the total aorta image are obtained. obtaining an image including the coronary ostium based on the aorta complementary image, and identifying the location of the coronary ostium in the aorta based on the image including the coronary ostium;
S03: In the image including the myocardium, the ostium of the coronary artery is used as a seed point, the coronary artery is extracted by region growing, the average grayscale of the coronary artery is calculated, and the grayscale distribution of the coronary artery and the coronary artery are calculated . extracting the coronary artery tree along the growth direction of said coronary arteries by comparing the average grayscale of
S04: binarize the image of the coronary artery tree , draw an isosurface image in a three-dimensional coordinate system based on the binarized data , and finally obtain a three-dimensional grid image of the coronary artery,
S05: Calculate the total flow Q _total of the coronary artery ostia in the state of maximum hyperemia = myocardial volume × resting myocardial blood flow × CFR, where CFR is the coronary flow reserve (CFR);
S06: Calculate the blood flow velocity _V1 in the state of maximum hyperemia,
S07: Calculate the pressure loss ΔP from the coronary entrance to the distal coronary artery stenosis and the mean intracoronary artery distal stenosis pressure P _d =P _a −ΔP, where V ₁ is the inlet flow velocity of the coronary artery stenosis, provided that: P _a is the mean aortic pressure of the heart , giving the coronary flow reserve ratio FFR=P _d /P _a
including steps
The step S06 is
S61: Based on the flow-volume scaling law and the three-dimensional reconstruction of the coronary artery tree, the blood flow in any one vessel in the coronary artery tree Q = Q total × ( _V / V _total ) ^3/4 wherein Q _total is the total flow of the coronary ostia at the maximum hyperemic state, V _total is the sum of vessel volumes of all coronary arteries in the coronary tree of the three-dimensional reconstruction; V is the sum of vessel volumes in any one vessel in the coronary artery tree and its downstream vessels;
S62: Determine the blood flow velocity V ₁ =Q/D in any one vessel in the coronary artery tree based on the flow-volume scaling law and the three-dimensional reconstruction of the coronary artery tree, where D is is the average cross-sectional area of the blood vessel,
including
Specifically, the step S07 is
Solve for the vessel 3D grid and solve the continuity and Navier-Stokes equations using numerical methods:

however,

P, ρ, μ are flow velocity, pressure, blood density, blood viscosity, respectively,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia ,
By three-dimensional computational fluid dynamics, the pressure losses ΔP _{1 ,} ΔP 2 _, ΔP ₃ , . ), calculating the mean pressure in the stenotic distal coronary artery Pd = Pa - ΔP, where Pa is the aortic mean,
A method for calculating a coronary flow reserve ratio based on a myocardial blood flow and a CT image, characterized by: In step S02, after obtaining the image containing the coronary ostia, perform connected region analysis on the image containing the coronary ostia, mark each connected region with a different grayscale label, and measure the coronary artery in the aorta. 2. The method for calculating coronary artery flow reserve ratio based on myocardial blood flow and CT images according to claim 1, wherein the position of the mouth is specified. 2. A cardiac aorta image is obtained in said step S02 by extracting the ascending aorta of the heart and its centerline by utilizing the fact that the cross section of the aorta is circular in the heart image . A method for calculating the coronary coronary flow reserve ratio based on the described myocardial blood flow and CT images. Coronary artery image binarization in step S04 includes:
Traversing the voxels in the coronary artery image V1, if a voxel pixel is equal to 0, the pixel value remains unchanged, if not equal to 0, the pixel value is 1, and these pixel values are a new data. including becoming V2,
2. The method for calculating a coronary artery coronary flow reserve ratio based on a myocardial blood flow and a CT image according to claim 1. Resting-state myocardial blood is measured by myocardial contrast echocardiography (MCE), single photon emission tomography (SPECT), positron emission tomography (PET), cardiac magnetic resonance imaging (MRI), or CT perfusion in step S05. determining flow rate and coronary flow reserve (CFR);
2. The method for calculating a coronary artery coronary flow reserve ratio based on a myocardial blood flow and a CT image according to claim 1. The step S07 is
Based on the CT-reconstructed three-dimensional grid image of the coronary artery , the stenosed vessel is straightened, a two-dimensional axisymmetric model is constructed, the two-dimensional grid is divided, and continuity and navigator analysis is performed using numerical methods. Find the solution of the Stokes equation:

where ρ represents the density of blood, u _z and u _r represent the flow velocities in the z and r directions, respectively, μ represents the kinematic viscosity of blood, p represents the strength of blood pressure,
The inlet boundary condition is the inlet flow velocity _V1 of the coronary artery stenosis at maximum hyperemia,
Pressure loss ΔP ₁ , ΔP ₂ , _{ΔP 3} , . . . ), calculating the stenotic distal coronary mean pressure P _d =P _a −ΔP, where P _a is the aortic mean pressure. Method for calculating coronary artery flow reserve ratio based on blood flow and CT images. The step S07 further includes
calculating the pressure drop from the coronary ostium to the distal coronary artery stenosis on a three-dimensional grid image of the coronary artery for different types of curvature of the coronary artery and in light of the two-dimensional axisymmetric model; constructing a database for storing correction factors for the two-dimensional axisymmetric results of each type of curvature;
After obtaining the pressure drop , obtain the corrected pressure drop from the coronary ostium to the distal coronary artery stenosis in light of the correction factors in the database, and then calculate the FFR.
7. The method for calculating the coronary coronary flow reserve ratio based on myocardial blood flow and CT images according to claim 6 , characterized by comprising:

Images