Method, device, equipment and medium for determining coronary artery blood flow reserve fraction

Abstract

The present invention relates to the field of image processing technology, and discloses a method, device, equipment and medium for determining the coronary blood flow reserve fraction, the method comprising: obtaining a sequence of coronary angiography images to be processed, coronary blood flow parameters and coronary artery ostium pressure parameters; extracting a reference image and a floating image from the sequence of coronary angiography images; performing three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional coronary artery image; extracting a coronary artery stenosis area from the three-dimensional coronary artery image; and determining the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary blood flow parameters and the coronary artery ostium pressure parameters. The method for determining the coronary blood flow reserve fraction provided in this embodiment can objectively and accurately determine the coronary blood flow reserve fraction of the coronary artery stenosis area by determining the coronary blood flow parameters and the coronary artery ostium pressure parameters.

Description

Technical Field

[0001] The present invention relates to the technical field of image processing, and in particular to a method, device, equipment and medium for determining a coronary artery blood flow reserve fraction. Background Art

[0002] Coronary angiography and intravascular ultrasound are both considered the "gold standard" for diagnosing coronary heart disease, but they can only perform imaging evaluation of the degree of stenosis of the lesion, and it is unknown to what extent the stenosis affects distal blood flow; Fractional Flow Reserve (FFR) has now become a recognized indicator for the functional evaluation of coronary artery stenosis, and its most important function is to accurately evaluate the functional consequences of a coronary artery stenosis of unknown effect.

[0003] At present, FFR is often determined by measuring the pressure at the distal end of coronary artery stenosis using a pressure guidewire in angiography and resting states. This method requires the tester to manually measure the coronary artery ostium pressure using a pressure guidewire, but manual measurement may result in errors in the blood flow reserve fraction. Summary of the invention

[0004] In view of this, the present invention provides a method, device, equipment and medium for determining the coronary blood flow reserve fraction to solve the problem that testers manually use a pressure guide wire to measure the coronary artery orifice pressure, but the manual measurement method may be negligent, thus causing errors in the blood flow reserve fraction.

[0005] In a first aspect, the present invention provides a method for determining a coronary blood flow reserve fraction, the method comprising: acquiring a coronary angiography image sequence to be processed, coronary blood flow parameters, and coronary artery orifice pressure parameters; extracting a reference image and a floating image from the coronary angiography image sequence; performing three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional coronary artery image; extracting a coronary artery stenosis area from the three-dimensional coronary artery image; and determining the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameters and the coronary artery orifice pressure parameters.

[0006] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, there may be a location without a lesion in the coronary artery stenosis area, that is, the coronary blood flow reserve fraction is a normal value. According to the coronary blood flow reserve fraction of the coronary artery stenosis area, the location with a lesion and the location without a lesion in the coronary artery stenosis area can be accurately distinguished. Compared with the situation of relying on artificial measurement of the coronary artery ostium pressure using a pressure guide wire, this method can determine the coronary blood flow reserve fraction of the coronary artery stenosis area through coronary blood flow parameters and coronary artery ostium pressure parameters.

[0007] In an optional embodiment, the method of determining coronary blood flow parameters includes: obtaining the volume of the coronary arteries, the length of the coronary arteries, the flow rate of the contrast agent, the start time and the end time of the terminal in the three-dimensional image of the coronary arteries; and determining the coronary blood flow reserve fraction based on the volume of the coronary arteries, the start time, the end time and the flow rate of the contrast agent.

[0008] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, the start time and the end time of the terminal are factors that affect the coronary blood flow reserve fraction, that is, different coronary artery volumes, coronary artery lengths, contrast agent flows, start times and end times will result in different coronary blood flow reserve fractions. Therefore, the coronary blood flow reserve fraction can be accurately determined by the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, the start time and the end time of the terminal.

[0009] In an optional embodiment, the blood flow reserve fraction of the coronary artery stenosis area is determined based on the coronary blood flow parameters and the coronary artery orifice pressure parameters, including: obtaining the density of any vascular point in the coronary artery stenosis area; determining the pressure drop parameters of each vascular point in the coronary artery stenosis area based on the coronary blood flow parameters and the density; determining the coronary blood flow reserve fraction of each vascular point in the coronary artery stenosis area based on the pressure drop parameters and the coronary artery orifice pressure parameters.

[0010] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, since the coronary blood flow parameters corresponding to each vascular point are different, the pressure drop parameters are different, so that the coronary blood flow reserve fraction determined by the pressure drop parameters of each vascular point and the coronary artery orifice pressure parameters are different. This method can accurately determine the pressure drop parameters of each vascular point through the density of any vascular point and the coronary blood flow parameters, thereby accurately determining the coronary blood flow reserve fraction corresponding to each vascular point.

[0011] In an optional embodiment, the method further includes: comparing the coronary blood flow reserve fraction of each vascular point with a preset value, and determining a target vascular point from each vascular point whose coronary blood flow reserve fraction does not meet the preset value based on the comparison result; and determining the coronary artery lesion based on the target vascular point.

[0012] The method for determining the coronary blood flow reserve fraction provided in this embodiment compares the coronary blood flow reserve fraction of each vascular point with a preset value to determine the target vascular point where the coronary blood flow reserve fraction does not meet the preset value. The position corresponding to the vascular point where the coronary blood flow reserve fraction does not meet the preset value is the coronary artery lesion, thereby achieving reasonable determination of the coronary artery lesion.

[0013] In an optional embodiment, a reference image and a floating image are extracted from a coronary angiography image sequence, including: extracting an image of a blood vessel that meets a preset length parameter from the coronary angiography image sequence, and using the image of the blood vessel as a reference image; based on the reference image, extracting images of other blood vessels that meet preset angle parameters with the blood vessel from the coronary angiography image sequence, and using the images of the other blood vessels as floating images.

[0014] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, the lower the angle formed between the reference image and the floating image, the less accurate the result obtained when the reference image and the floating image are subjected to three-dimensional reconstruction processing. By determining the reference image and at least one floating image that meets a preset angle with the reference image to ensure that the reconstruction of the three-dimensional image of the coronary artery is more accurate, the reconstruction quality of the three-dimensional image of the coronary artery is improved.

[0015] In an optional embodiment, the method of obtaining the coronary artery ostium pressure parameter includes: obtaining the radial artery pressure parameter; performing coefficient conversion processing on the radial artery pressure parameter to determine the coronary artery ostium pressure parameter.

[0016] The method for determining the coronary blood flow reserve fraction provided in this embodiment can accurately and conveniently determine the coronary artery ostium pressure parameter by determining the radial artery pressure parameter and performing coefficient conversion processing on the radial artery pressure parameter, because the coronary artery ostium pressure parameter is not convenient to measure directly.

[0017] In an optional embodiment, a three-dimensional reconstruction process is performed on the reference image and the floating image to generate a three-dimensional image of the coronary arteries, including: determining a reference centerline of a first blood vessel in the reference image and a floating centerline of a second blood vessel in the floating image; determining a first radius corresponding to each first blood vessel point based on each first blood vessel point in the reference centerline; determining a second radius corresponding to each second blood vessel point based on each second blood vessel point in the floating centerline; and performing three-dimensional reconstruction process on the reference image and the floating image based on the first radius and the second radius to generate a three-dimensional image of the coronary arteries.

[0018] The method for determining the coronary blood flow reserve fraction provided in this embodiment determines a reference center line and a floating center line, then determines a first radius of a first blood vessel in a reference image according to the reference center line, determines a first radius of a second blood vessel in a floating image according to a floating center point, and then performs three-dimensional reconstruction processing on the reference image and the floating image according to the first radius and the second radius, so as to accurately obtain a three-dimensional image of the coronary arteries.

[0019] In a second aspect, the present invention provides a device for determining a coronary blood flow reserve fraction, the device comprising: an acquisition module for acquiring a coronary angiography image sequence to be processed, coronary blood flow parameters and coronary artery orifice pressure parameters; an extraction module for extracting a reference image and a floating image from the coronary angiography image sequence; a three-dimensional reconstruction module for performing three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary artery; a coronary artery stenosis area determination module for determining a coronary artery stenosis area from the three-dimensional image of the coronary artery; and a blood flow reserve fraction determination module for determining the blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameters and the coronary artery orifice pressure parameters.

[0020] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for determining the coronary blood flow reserve fraction of the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0021] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for determining the coronary blood flow reserve fraction of the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0023] Figure 1 is a flow chart of a method for determining coronary artery blood flow reserve fraction according to an embodiment of the present invention;

[0024] Figure 2 is a schematic diagram of a coronary angiography image sequence according to an embodiment of the present invention;

[0025] Figure 3 is a flow chart of another method for determining coronary artery blood flow reserve fraction according to an embodiment of the present invention;

[0026] Figure 4 is a flow chart of another method for determining coronary artery blood flow reserve fraction according to an embodiment of the present invention;

[0027] Figure 5 is a schematic diagram of a method for determining coronary blood flow reserve fraction according to an embodiment of the present invention;

[0028] Figure 6 is a structural block diagram of a device for determining coronary blood flow reserve fraction according to an embodiment of the present invention;

[0029] Figure 7 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0030] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0031] In the related art, the coronary artery blood flow reserve fraction (FFR) is often determined by measuring the pressure at the distal end of the coronary artery stenosis using a pressure guide wire in the angiography state and the resting state. This method requires the tester to manually measure the coronary artery ostium pressure using a pressure guide wire, but the manual measurement method may be negligent, resulting in errors in the blood flow reserve fraction.

[0032] Based on this, the technical solution of the present invention combines coronary blood flow parameters and coronary artery ostium pressure parameters to determine the coronary blood flow reserve fraction in the coronary artery stenosis area, thereby ensuring accurate determination of the blood flow reserve fraction.

[0033] According to an embodiment of the present invention, an embodiment of a method for determining a coronary blood flow reserve fraction is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0034] In this embodiment, a method for determining the coronary blood flow reserve fraction is provided, which can be used in computer equipment, such as computers, CT scanning equipment, etc. Figure 1 is a flow chart of a method for determining coronary blood flow reserve fraction according to an embodiment of the present invention, such as Figure 1 As shown, the process includes the following steps:

[0035] Step S101, obtaining a sequence of coronary angiography images to be processed, coronary blood flow parameters and coronary ostium pressure parameters.

[0036] Figure 2 A schematic diagram of a coronary angiography image sequence is shown, combined with Figure 2 As shown, the coronary angiography image sequence may be an image formed by multiple blood vessels of the coronary artery, and the coronary angiography image sequence may be obtained by coronary angiography. In the coronary angiography, a catheter is inserted through the femoral artery or other peripheral artery in the thigh, sent to the ascending aorta, and then inserted into the left or right coronary artery orifice, and an iodine contrast agent is injected to visualize the coronary artery, thereby obtaining the coronary angiography image sequence.

[0037] The coronary blood flow parameter may be the flow rate of the coronary blood, which may be calculated from the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, etc. The specific calculation method of the coronary blood flow parameter is described in detail below.

[0038] The coronary artery pressure parameter may be the average pressure of the coronary artery. It is not easy to measure directly by a detection device, but can be calculated by a radial artery pressure parameter. The specific calculation method of the coronary artery pressure parameter is described in detail below.

[0039] Step S102: extracting a reference image and a floating image from the coronary angiography image sequence.

[0040] After the computer device acquires the coronary angiography image sequence, the computer device can extract the reference image and the floating image by means of feature recognition. For example, the coronary angiography image sequence can be input into an image recognition model, and the feature values ​​of blood vessels in the reference image and the floating image can be identified by the image recognition model, and then the reference image and the floating image can be extracted according to the feature values ​​of the blood vessels. The reference image and the floating image can also be extracted by image processing software (such as Photoshop software, PicsArt photo editing software, etc.), which is not specifically limited here.

[0041] Combination Figure 2 As shown, the reference image may be an image of a blood vessel in a coronary angiography image sequence. Figure 3The floating image may be an image of a blood vessel in a coronary angiography image sequence. Figure 3 The blood vessel B may also be the blood vessel C, and no specific limitation is made here.

[0042] Step S103: Perform three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary arteries.

[0043] 3D reconstruction refers to the establishment of a mathematical model suitable for computer representation and processing of 3D objects. It is the basis for processing, operating and analyzing the properties of 3D objects in a computer environment, and is also a key technology for establishing virtual reality in a computer to express the objective world. Computer equipment generates a 3D image of the coronary artery by performing 3D reconstruction on 2D images (i.e., reference images and floating images). Specifically, the specific steps of 3D reconstruction are described in detail below.

[0044] Step S104: determining a coronary artery stenosis region from the coronary artery three-dimensional image.

[0045] The coronary artery stenosis area may be an area where a lesion occurs in the coronary artery. If the coronary artery stenosis area appears in the coronary artery three-dimensional image, it indicates that blood cannot normally enter the heart circulation through the coronary artery. The computer device can determine the pulmonary artery stenosis area from the pulmonary artery three-dimensional image. Specifically, the specific steps for determining the pulmonary artery stenosis area are described in detail below.

[0046] Step S105, determining the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter.

[0047] Coronary artery blood flow reserve fraction (FFR) is the ratio of the coronary artery blood flow parameter provided by the stenotic coronary artery to the myocardium in the innervated area to the coronary artery blood flow parameter provided by the same coronary artery when it is normal. It can be simplified as the ratio of the pressure in the stenotic area (Pd) under the state of maximum myocardial hyperemia to the pressure parameter (Pa) at the coronary artery orifice, that is, FFR = Pd / Pa.

[0048] The computer device can determine the pressure in the stenotic area through the coronary blood flow parameters, and then determine the coronary blood flow reserve fraction through the ratio of the pressure in the stenotic area to the coronary artery ostium pressure parameter.

[0049] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, there may be a location without a lesion in the coronary artery stenosis area, that is, the coronary blood flow reserve fraction is a normal value. According to the coronary blood flow reserve fraction of the coronary artery stenosis area, the location with a lesion and the location without a lesion in the coronary artery stenosis area can be accurately distinguished. Compared with the situation of relying on artificial measurement of the coronary artery ostium pressure using a pressure guide wire, this method can determine the coronary blood flow reserve fraction of the coronary artery stenosis area through coronary blood flow parameters and coronary artery ostium pressure parameters.

[0050] In this embodiment, a method for determining the coronary blood flow reserve fraction is provided, which can be used in computer equipment, such as computers, CT angiography equipment, etc. Figure 3 FIG. 4 is a flow chart of a method for determining coronary blood flow reserve fraction according to an embodiment of the present invention. Figure 3 As shown, the process includes the following steps:

[0051] Step S201, obtaining a sequence of coronary angiography images to be processed, coronary blood flow parameters, and coronary artery ostium pressure parameters. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0052] Step S202: extracting a reference image and a floating image from the coronary angiography image sequence.

[0053] Specifically, the above step S202 includes:

[0054] Step S2021 , extracting an image of a blood vessel that meets a preset length parameter from the coronary angiography image sequence, and using the image of the blood vessel as a reference image.

[0055] The computer device can determine a blood vessel that meets the preset length parameter from the coronary angiography image sequence, then extract the image corresponding to the blood vessel, and use the image corresponding to the blood vessel as the reference image. Specifically, the preset length parameter can be the longest blood vessel among all blood vessels in the coronary angiography image sequence. For example, if the length of blood vessel A is 3 cm, the length of blood vessel B is 5 cm, and the length of blood vessel C is 7 cm, then the computer device can use the image corresponding to blood vessel C as the reference image.

[0056] Step S2022: based on the reference image, extracting images of other blood vessels that meet preset angle parameters with the blood vessel from the coronary angiography image sequence, and using the images of other blood vessels as floating images.

[0057] The number of floating images may be one or more, and is not specifically limited here. After determining the reference image, the computer device may determine images of other blood vessels from the coronary angiography image sequence whose angles with the blood vessels in the reference image are not less than a preset angle parameter. Specifically, the preset angle parameter may be 30°, that is, as long as the angle between the blood vessels in the reference image and the other blood vessels is not less than 30°, they may be used as blood vessels in the floating image.

[0058] Step S203, perform 3D reconstruction on the reference image and the floating image to generate a 3D image of the coronary arteries. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0059] Step S204, determine the coronary artery stenosis area from the coronary artery three-dimensional image. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0060] Step S205, determining the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter.

[0061] Specifically, the above step S205 includes:

[0062] Step S2051, obtaining the density of any blood vessel point in the coronary artery stenosis area.

[0063] The density of the blood vessel point may be the density of the blood at the position corresponding to the blood vessel point. The density of the blood may be detected by a density detection instrument. Specifically, the density detection instrument may be a densitometer, etc., which is not specifically limited here.

[0064] Step S2052: Determine the pressure drop parameter of each blood vessel point in the coronary artery stenosis area based on the coronary artery blood flow parameter and density.

[0065] Formula 1: Formula 2: Where u is the fluid velocity vector, is the pressure drop parameter, is the change of blood flow velocity of coronary blood flow parameters at each vascular point, Find the vector differential for velocity, v is the kinematic viscosity of the fluid, and v is a fixed parameter.

[0066] According to the above formula, if is equal to 0; then It is also equal to 0; is also equal to 0, then formula 3 can be derived from the above formula: According to Formula 3, the pressure drop parameter at each vascular point is related to the change in blood flow velocity of the coronary blood flow parameter at each vascular point. Therefore, the pressure drop parameter at each vascular point can be determined by the change in blood flow velocity of the coronary blood flow parameter at each vascular point.

[0067] Step S2053, determining the coronary blood flow reserve fraction of each vascular point in the coronary artery stenosis area based on the pressure drop parameter and the coronary artery ostium pressure parameter.

[0068] From the above content, it can be seen that FFR = Pd / Pa; among them, Pd is the pressure drop parameter of each vascular point, and Pa is the coronary artery orifice pressure parameter. Therefore, the coronary blood flow reserve fraction of each vascular point can be determined by the ratio of the pressure drop parameter and the coronary artery orifice pressure parameter.

[0069] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, since the coronary blood flow parameters corresponding to each vascular point are different, the pressure drop parameters are different, so that the coronary blood flow reserve fraction determined by the pressure drop parameters of each vascular point and the coronary artery orifice pressure parameters are different. This method can accurately determine the pressure drop parameters of each vascular point through the density of any vascular point and the coronary blood flow parameters, thereby accurately determining the coronary blood flow reserve fraction corresponding to each vascular point.

[0070] In this embodiment, a method for determining the coronary blood flow reserve fraction is provided, which can be used in computer equipment, such as computers, CT angiography equipment, etc. Figure 4 is a flow chart of a method for determining coronary blood flow reserve fraction according to an embodiment of the present invention, such as Figure 4 As shown, the process includes the following steps:

[0071] Step S301, obtaining a sequence of coronary angiography images to be processed, coronary blood flow parameters and coronary ostium pressure parameters.

[0072] Specifically, the above step S301 includes:

[0073] Step S3011, obtaining the volume of the coronary arteries, the length of the coronary arteries, the flow rate of the contrast medium, the start time and the end time of the end in the three-dimensional image of the coronary arteries.

[0074] The flow rate of contrast agent can be used to characterize the amount injected during coronary angiography; the volume of the coronary artery and the length of the coronary artery can be used to characterize the volume and length of the coronary artery in a three-dimensional image of the coronary artery. The volume and length of the coronary artery and the flow rate of contrast agent can be obtained through coronary angiography.

[0075] The start time can be the time when the blood flows to the coronary artery orifice, and the end time can be the time when the blood flows to the coronary artery terminal. For example, if the blood flows to the coronary artery orifice at 10 o'clock, then 10 o'clock is the start time, and if the blood flows to the coronary artery terminal at 10:05, then 10:05 is the end time.

[0076] Step S3012, determining coronary artery blood flow parameters based on the volume of the coronary artery, the start time, the end time, and the flow rate of the contrast agent.

[0077] in, Coronary blood flow parameters are changes in blood flow velocity at each vascular point, t1 is the start time, t2 is the end time, V is the volume of the pulmonary artery, and s is the flow rate of the contrast agent.

[0078] Step S3013, obtaining radial artery pressure parameters.

[0079] The radial artery pressure parameters can be measured by an invasive direct arterial pressure measurement method, wherein the invasive direct arterial pressure measurement method is to place a catheter in a peripheral artery, connect a pressure transducer, and continuously display the arterial pressure waveform and radial artery pressure parameters on a display screen. Other pressure measurement methods can also be used for measurement, which are not specifically limited here.

[0080] Step S3014, performing coefficient conversion processing on the radial artery pressure parameter to determine the coronary artery ostium pressure parameter.

[0081] Pr = A × Pa; where A is a coefficient, Pr is a radial artery pressure parameter, and Pa is a coronary artery pressure parameter. After obtaining the radial artery pressure parameter, the computer device can determine the coronary artery pressure parameter by the ratio of the radial artery pressure parameter and the coefficient.

[0082] Step S3015, obtaining a sequence of coronary angiography images to be processed. Figure 3 Step S201 of the illustrated embodiment will not be described in detail here.

[0083] Step S302: extracting a reference image and a floating image from the coronary angiography image sequence. Figure 3 Step S202 of the illustrated embodiment will not be described in detail here.

[0084] Step S303: Perform three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary arteries.

[0085] Specifically, the above step S303 includes:

[0086] Step S3031, determining a reference center line of a first blood vessel in a reference image and a floating center line of a second blood vessel in a floating image.

[0087] The computer device can divide the blood vessel from the starting segment of the first blood vessel to the end of the first blood vessel into multiple blood vessel segments. The length of each blood vessel segment can be 1-8 mm, which is not specifically limited here. Then, the pixel points in the multiple blood vessel segments are extracted, and the pixel points are fitted to obtain the spherical surface corresponding to each blood vessel segment, and the center point of each spherical surface is determined based on the spherical surface, and the center points are connected to form a reference center point. The method for determining the floating center point is the same as the method for determining the reference center point, which will not be repeated here.

[0088] Step S3032: determining a first radius corresponding to each first blood vessel point based on each first blood vessel point in the reference center line.

[0089] The computer device can obtain the intersection points of the spherical surface and the vascular wall of the first blood vessel, and determine the first radius through the connection lines between each first blood vessel point and the intersection points. The vascular wall of the first blood vessel is detected by a detection device. Specifically, the detection instrument can be a vascular imager, an intravascular ultrasonic diagnostic instrument, etc., which is not specifically limited here.

[0090] Step S3033: determining a second radius corresponding to each second blood vessel point based on each second blood vessel point in the floating center line.

[0091] The method for determining the second radius is the same as the method for determining the first radius, and will not be described in detail here.

[0092] Step S3034: Perform three-dimensional reconstruction processing on the reference image and the floating image based on the first radius and the second radius to generate a three-dimensional image of the coronary arteries.

[0093] Specifically, the two-dimensional planes of the reference image and the floating image are converted from the image plane coordinate system to the projection plane coordinate system, and the three-dimensional coordinate points are obtained using the plane point pairs and distances, rotation matrix, and translation matrix in the projection plane coordinate system, and the three-dimensional image of the coronary artery is obtained based on the three-dimensional coordinate points.

[0094] Step S304, determine the coronary artery stenosis area from the coronary artery three-dimensional image. Figure 3 Step S204 of the illustrated embodiment will not be described in detail here.

[0095] Step S305: Determine the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter. Figure 3 Step S205 of the illustrated embodiment will not be described in detail here.

[0096] The method for determining the coronary blood flow reserve fraction provided in this embodiment can accurately and conveniently determine the coronary artery ostium pressure parameter by determining the radial artery pressure parameter and performing coefficient conversion processing on the radial artery pressure parameter, because the coronary artery ostium pressure parameter is not convenient to measure directly.

[0097] In an optional embodiment, the above method further includes:

[0098] Step a1, comparing the coronary artery blood flow reserve fractions of various vascular points, and determining the target vascular point with the smallest blood flow reserve fraction from among the vascular points based on the comparison result.

[0099] Step a2: determining the lesion of the coronary artery based on the target vessel point.

[0100] After obtaining the coronary blood flow reserve fraction of each vascular point, the computer device can detect the value of the coronary blood flow reserve fraction corresponding to each vascular point in the coronary artery stenosis area, and compare the value of the coronary blood flow reserve fraction with a preset value to determine the target vascular point whose coronary blood flow reserve fraction value is less than the preset value. The position corresponding to the target vascular point is the coronary artery lesion.

[0101] The method for determining the coronary blood flow reserve fraction provided in this embodiment compares the coronary blood flow reserve fraction of each vascular point with a preset value to determine the target vascular point where the coronary blood flow reserve fraction does not meet the preset value. The position corresponding to the vascular point where the coronary blood flow reserve fraction does not meet the preset value is the coronary artery lesion, thereby achieving reasonable determination of the coronary artery lesion.

[0102] Combination Figure 5 As shown, in an optional embodiment, an embodiment of the present invention provides a schematic diagram of a method for determining a coronary blood flow reserve fraction, the method comprising:

[0103] Step b1, measure the radial artery pressure (radial_artery_P), and calculate the radial artery average pressure radial_artery_average_P (i.e. the above-mentioned coronary artery pressure parameter) based on the radial artery pressure. Figure 4 Step S3014 of the illustrated embodiment will not be described in detail here.

[0104] Step b2, obtaining a reference image (img_A) in the coronary angiography (i.e., the above coronary angiography image sequence) and a reference image having an angle greater than a preset angle parameter (i.e., Figure 5 30℃) of floating images (img_B1, img_B2, etc.). For details, see Figure 4 Step S302 of the illustrated embodiment will not be described in detail here.

[0105] Step b3, respectively extract the coronary artery centerline (line_A) of the reference image and the coronary artery centerline (line_B1, line_B2, etc.) of the floating image, and calculate the coronary artery radius (rad_A) based on the coronary artery centerline of the reference image, and calculate the coronary artery radius (rad_B1, rad_B2, etc.) based on the coronary artery centerline of the floating image. For details, please refer to Figure 4 Steps S3031 to S3034 of the illustrated embodiment are not described in detail here.

[0106] Step b4, perform three-dimensional reconstruction (mesh_3D) based on the coronary artery radius of the reference image and the coronary artery radius of the floating image, generate a three-dimensional reconstruction result, and determine the stenosis information narrow_A (i.e., the coronary artery stenosis area) based on the three-dimensional reconstruction result. For details, please refer to Figure 4 Steps S3031 to S3034 of the illustrated embodiment are not described in detail here.

[0107] Step b5, determine velocity velocity_A (i.e. the above coronary artery blood flow parameter). Figure 4 Step S3012 of the illustrated embodiment will not be described in detail here.

[0108] Step b6, determine the coronary quantitative blood flow fraction (i.e. the above-mentioned coronary blood flow reserve fraction) based on the flow velocity, stenosis information and radial artery mean pressure. Figure 4 Step S305 of the illustrated embodiment will not be described in detail here.

[0109] In the method for determining the coronary blood flow reserve fraction provided in this embodiment, there may be a location without a lesion in the coronary artery stenosis area, that is, the coronary blood flow reserve fraction is a normal value. According to the coronary blood flow reserve fraction of the coronary artery stenosis area, the location with a lesion and the location without a lesion in the coronary artery stenosis area can be accurately distinguished. Compared with the situation of relying on artificial measurement of the coronary artery ostium pressure using a pressure guide wire, this method can determine the coronary blood flow reserve fraction of the coronary artery stenosis area through coronary blood flow parameters and coronary artery ostium pressure parameters.

[0110] In this embodiment, a device for determining the coronary blood flow reserve fraction is also provided, and the device is used to implement the above-mentioned embodiments and preferred implementation modes, and the descriptions that have been made are not repeated here. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0111] This embodiment provides a device for determining the coronary blood flow reserve fraction. Figure 6 As shown, including:

[0112] An acquisition module 601 is used to acquire a coronary angiography image sequence to be processed, coronary blood flow parameters and coronary ostium pressure parameters;

[0113] An extraction module 602, for extracting a reference image and a floating image from a coronary angiography image sequence;

[0114] A three-dimensional reconstruction module 603 is used to perform three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary artery;

[0115] A coronary artery stenosis region determination module 604 is used to determine the coronary artery stenosis region from the coronary artery three-dimensional image;

[0116] The blood flow reserve fraction determination module 605 is used to determine the blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter.

[0117] In an optional implementation, the acquisition module 601 includes:

[0118] The first acquisition unit is used to acquire the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, the start time and the end time of the end in the three-dimensional image of the coronary artery.

[0119] The first determination unit is used to determine the coronary artery blood flow parameter based on the volume of the coronary artery, the start time, the end time and the flow rate of the contrast agent.

[0120] In an optional embodiment, the blood flow reserve fraction determination module 605 includes:

[0121] A second acquisition unit is used to acquire the density of any blood vessel point in the coronary artery stenosis area;

[0122] A second determination unit is used to determine the pressure drop parameter of each blood vessel point in the coronary artery stenosis area based on the coronary artery blood flow parameter and density;

[0123] The third determining unit is used to determine the coronary blood flow reserve fraction of each blood vessel point in the coronary artery stenosis area based on the pressure drop parameter and the coronary artery ostium pressure parameter.

[0124] In an optional embodiment, the device further comprises: a comparison module, configured to compare the coronary blood flow reserve fraction of each vascular point with a preset value, and determine, based on the comparison result, a target vascular point from each vascular point whose coronary blood flow reserve fraction does not meet the preset value;

[0125] The lesion determination module is used to determine the lesion of the coronary artery based on the target vessel point.

[0126] In an optional implementation, the extraction module 602 includes:

[0127] A first extraction unit is used to extract an image of a blood vessel that meets a preset length parameter from the coronary angiography image sequence, and use the image of the blood vessel as a reference image;

[0128] The second extraction unit is used to extract images of other blood vessels that meet preset angle parameters with the blood vessel from the coronary angiography image sequence based on the reference image, and use the images of other blood vessels as floating images.

[0129] In an optional implementation, the acquisition module 601 includes:

[0130] The third acquisition unit is used to acquire radial artery pressure parameters.

[0131] The conversion processing unit is used to perform coefficient conversion processing on the radial artery pressure parameter to determine the coronary artery orifice pressure parameter.

[0132] In an optional implementation, the 3D reconstruction module 603 includes:

[0133] a fourth determining unit, configured to determine a reference center line of the first blood vessel in the reference image and a floating center line of the second blood vessel in the floating image;

[0134] a fifth determining unit, configured to determine, based on each first blood vessel point in the reference center line, a first radius corresponding to each first blood vessel point;

[0135] a sixth determining unit, configured to determine, based on each second blood vessel point in the floating center line, a second radius corresponding to each second blood vessel point;

[0136] The reconstruction processing unit is used to perform three-dimensional reconstruction processing on the reference image and the floating image based on the first radius and the second radius to generate a three-dimensional image of the coronary artery.

[0137] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0138] The device for determining the coronary blood flow reserve fraction in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and a memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0139] In the device for determining the coronary blood flow reserve fraction provided in this embodiment, there may be a location without a lesion in the coronary artery stenosis area, that is, the coronary blood flow reserve fraction is a normal value. According to the coronary blood flow reserve fraction of the coronary artery stenosis area, the location where the lesion exists and the location where the lesion does not exist in the coronary artery stenosis area can be accurately distinguished. Compared with the situation of relying on artificial measurement of the coronary artery ostium pressure using a pressure guide wire, this method can determine the coronary blood flow reserve fraction of the coronary artery stenosis area through the coronary artery blood flow parameters and the coronary artery ostium pressure parameters.

[0140] The embodiment of the present invention also provides a computer device having the above Figure 6 An apparatus for determining the coronary blood flow reserve fraction is shown.

[0141] See also Figure 7 , Figure 7 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Figure 7 As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 7 A processor 10 is taken as an example.

[0142] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0143] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0144] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0145] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid state drive; the memory 20 may also include a combination of the above types of memory.

[0146] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0147] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0148] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for determining coronary artery blood flow reserve fraction, characterized in that: The method comprises: Acquiring a sequence of coronary angiography images to be processed, coronary blood flow parameters, and coronary artery ostium pressure parameters; Extracting a reference image and a floating image from the coronary angiography image sequence; wherein the number of floating images is multiple, and the angle between the blood vessel in the reference image and other blood vessels is not less than 30°, and all of them are used as blood vessels in the floating image; Performing three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary arteries; determining a coronary artery stenosis area from the coronary artery three-dimensional image; Determining the coronary blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter; The method of determining the coronary artery blood flow parameter comprises: Acquiring the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, the start time and the end time of the end in the three-dimensional image of the coronary artery; Determining the coronary artery blood flow parameter based on the volume of the coronary artery, the start time, the end time, and the flow rate of the contrast agent; The method of obtaining the coronary artery pressure parameter includes: Obtain radial artery pressure parameters; The radial artery pressure parameter is subjected to coefficient conversion processing to determine the coronary artery ostium pressure parameter.

2. The method for determining the coronary blood flow reserve fraction according to claim 1, characterized in that: Determining the blood flow reserve fraction of the coronary artery stenosis area based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter includes: Acquiring the density of any blood vessel point in the coronary artery stenosis area; Determining the pressure drop parameter of each blood vessel point in the coronary artery stenosis region based on the coronary artery blood flow parameter and the density; The coronary blood flow reserve fraction of each of the vascular points in the coronary artery stenosis region is determined based on the pressure drop parameter and the coronary artery ostium pressure parameter.

3. The method for determining the coronary blood flow reserve fraction according to claim 2, characterized in that: Also includes: comparing the coronary artery blood flow reserve fraction of each of the vascular points with a preset value, and determining, based on the comparison result, a target vascular point from each of the vascular points whose coronary artery blood flow reserve fraction does not meet the preset value; Based on the target vessel point, a lesion of the coronary artery is determined.

4. The method for determining the coronary blood flow reserve fraction according to claim 1, characterized in that: The step of extracting a reference image and a floating image from the coronary angiography image sequence comprises: Extracting an image of a blood vessel that meets a preset length parameter from the coronary angiography image sequence, and using the image of the blood vessel as the reference image; Based on the reference image, images of other blood vessels that meet preset angle parameters with respect to the blood vessel are extracted from the coronary angiography image sequence, and the images of the other blood vessels are used as floating images.

5. The method for determining the coronary blood flow reserve fraction according to claim 1, characterized in that: The three-dimensional reconstruction of the reference image and the floating image to generate a three-dimensional image of the coronary arteries includes: determining a reference centerline of a first blood vessel in the reference image and a floating centerline of a second blood vessel in the floating image; Determine, based on each first blood vessel point in the reference center line, a first radius corresponding to each first blood vessel point; Based on each second blood vessel point in the floating center line, determining a second radius corresponding to each second blood vessel point; The reference image and the floating image are subjected to three-dimensional reconstruction processing based on the first radius and the second radius to generate the coronary artery three-dimensional image.

6. A device for determining coronary artery blood flow reserve fraction, characterized in that: The device comprises: An acquisition module, used for acquiring a sequence of coronary angiography images to be processed, coronary blood flow parameters and coronary artery ostium pressure parameters; An extraction module is used to extract a reference image and a floating image from the coronary angiography image sequence; wherein the number of floating images is multiple, and the angle between the blood vessel in the reference image and other blood vessels is not less than 30°, and all of them are used as blood vessels in the floating image; A three-dimensional reconstruction module, used for performing three-dimensional reconstruction processing on the reference image and the floating image to generate a three-dimensional image of the coronary arteries; A coronary artery stenosis region determination module, used to determine the coronary artery stenosis region from the coronary artery three-dimensional image; A blood flow reserve fraction determination module, configured to determine the blood flow reserve fraction of the coronary artery stenosis region based on the coronary artery blood flow parameter and the coronary artery ostium pressure parameter; The acquisition module includes: a first acquisition unit, which is used to acquire the volume of the coronary artery, the length of the coronary artery, the flow rate of the contrast agent, the start time and the end time of the end in the three-dimensional image of the coronary artery; a first determination unit, which is used to determine the coronary blood flow parameters based on the volume of the coronary artery, the start time, the end time and the flow rate of the contrast agent; a third acquisition unit, which is used to acquire the radial artery pressure parameters; and a conversion processing unit, which is used to perform coefficient conversion processing on the radial artery pressure parameters to determine the coronary artery orifice pressure parameters.

7. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for determining the coronary blood flow reserve fraction according to any one of claims 1 to 5 by executing the computer instructions.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for determining the coronary blood flow reserve fraction according to any one of claims 1 to 5.

Images