CN108564574B - Method, computer device and computer-readable storage medium for determining fractional flow reserve - Google Patents
Method, computer device and computer-readable storage medium for determining fractional flow reserve Download PDFInfo
- Publication number
- CN108564574B CN108564574B CN201810320618.4A CN201810320618A CN108564574B CN 108564574 B CN108564574 B CN 108564574B CN 201810320618 A CN201810320618 A CN 201810320618A CN 108564574 B CN108564574 B CN 108564574B
- Authority
- CN
- China
- Prior art keywords
- surface mesh
- coronary artery
- blood vessel
- mesh
- coronary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Abstract
The embodiment of the invention provides a method, computer equipment and a computer readable storage medium for determining fractional flow reserve, wherein the method comprises the following steps: acquiring a CTA image of a subject; constructing a vessel surface mesh and a coronary artery centerline from the CTA image; determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line; cutting the blood vessel surface mesh at the cutting position to generate an intermediate surface mesh; generating a body mesh model of the blood vessel according to the intermediate surface mesh; and obtaining a fractional flow reserve result of the coronary artery according to the body mesh model and the boundary condition of the blood vessel. The technical scheme provided by the embodiment of the invention can automatically determine the cutting position, avoid the problem of long time consumption and simultaneously improve the accuracy of the calculation of the blood flow reserve fraction of the blood vessel.
Description
[ technical field ] A method for producing a semiconductor device
The present invention relates to the field of image processing technologies, and in particular, to a method, a computer device, and a computer-readable storage medium for determining fractional flow reserve.
[ background of the invention ]
Fractional Flow Reserve (FFR) is a functional indicator for assessing the degree of myocardial ischemia resulting from coronary artery stenosis. With the development of computer technology and medical imaging technology, a Computed Tomography Angiography (CTA) image is combined with computational fluid dynamics, and FFR results can be obtained noninvasively by constructing a blood vessel model and fluid mechanics simulation.
In order to solve the problem of constructing a blood vessel model, a scheme adopted in the prior art is to directly generate a blood vessel volume mesh model, wherein the blood vessel comprises an aorta and a coronary artery, and specifically, after the aorta and the coronary artery are segmented from a CTA image, the volume mesh model is directly generated without any treatment on the coronary artery.
However, in the prior art, the blood vessel volume mesh model is directly generated, errors of the generated volume mesh caused by uneven thickness of each position of the coronary artery are not considered, and the obtained FFR result has low accuracy. In addition, in the prior art, the blood vessel mesh model is directly generated, and the cutting position is manually determined in the generation process, so that the cutting position cannot be automatically determined, and thus, the detection time is too long, and the operation is complex. In view of this, there is a need for an improvement of existing FFR detection methods.
[ summary of the invention ]
In view of this, embodiments of the present invention provide a method, a computer device, and a computer-readable storage medium for determining a fractional flow reserve, so as to solve the problem in the prior art that an error exists in a generated volume mesh due to non-uniform thickness of each position of a coronary artery, and an obtained FFR result is low in accuracy when a blood vessel volume mesh model is directly generated.
In one aspect, an embodiment of the present invention provides a method for determining a fractional flow reserve, where the method includes:
acquiring a CTA image of a subject;
constructing a vessel surface mesh and coronary artery centerlines from the CTA image;
determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line;
cutting the vessel surface mesh at the cutting location, generating an intermediate surface mesh;
generating a body mesh model of the blood vessel according to the intermediate surface mesh;
and obtaining a fractional flow reserve result of the coronary artery according to the volume mesh model and the boundary condition of the blood vessel.
The above aspect and any possible implementation further provide an implementation, before obtaining the fractional flow reserve result of the coronary artery according to the volume mesh model and the boundary condition of the blood vessel, the method further includes:
acquiring the boundary condition, wherein the boundary condition comprises at least one of an aorta inlet speed boundary condition, an aorta outlet pressure boundary condition and a coronary artery outlet flow resistance boundary condition.
The above aspect and any possible implementation further provide an implementation manner, wherein obtaining the coronary artery outlet flow resistance boundary condition comprises:
obtaining at least one of cardiac output and myocardial mass;
obtaining the total coronary flow based on at least one of the cardiac output and myocardial mass;
and acquiring the boundary condition of the coronary artery outlet flow resistance according to the total coronary artery flow and the parameters of the coronary artery.
The above aspect and any possible implementation further provide an implementation in which the parameter of the coronary artery includes at least one of a coronary bifurcation cross-sectional area, a coronary artery exit cross-sectional area, and a coronary artery intraluminal density attenuation gradient.
The above aspect and any possible implementation further provide an implementation in which, when the blood vessel is an aorta or a coronary artery, an area of a unit cell in the first surface mesh is larger than an area of a unit cell in the second surface mesh;
the first surface mesh is an outer surface mesh of the aorta;
the second surface mesh is an outer surface mesh of the coronary artery.
The above-described aspect and any possible implementation further provide an implementation in which a cutting position of the vessel surface mesh is determined from the vessel surface mesh and a coronary artery centerline, the method further comprising:
and determining the cutting position of the blood vessel surface mesh according to the roundness of an intersection line ring formed by the cross section determined by the point on the centerline of the coronary artery and the blood vessel surface mesh.
The above aspect and any possible implementation further provide an implementation, before determining a cutting position of the vessel surface mesh according to the vessel surface mesh and a coronary artery centerline, the method further comprising:
the coronary centerline is examined and corrected.
In another aspect, an embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor is configured to implement the following operations when executing the computer program:
acquiring a CTA image of a subject;
constructing a vessel surface mesh and coronary artery centerlines from the CTA image;
determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line;
cutting the vessel surface mesh at the cutting location, generating an intermediate surface mesh;
generating a body mesh model of the blood vessel according to the intermediate surface mesh;
and obtaining a fractional flow reserve result of the coronary artery according to the volume mesh model and the boundary condition of the blood vessel.
The above-described aspects and any possible implementation further provide an implementation, and the processor, when executing the computer program, further implements the following operations:
in response to the fractional flow reserve result for the coronary artery being less than a set threshold, a stenosis location is determined at the coronary artery.
In another aspect, an embodiment of the present invention provides a computer-readable storage medium, including: computer-executable instructions which, when executed, perform any of the above-described methods of determining fractional flow reserve.
In the embodiment of the invention, the surface grids of the blood vessels and the coronary artery central lines in the computed tomography angiography CTA image are constructed to automatically cut and position, so that the problem of overlong time consumption caused by manually positioning the cutting position is effectively avoided, and the whole detection operation is simplified; the method comprises the steps of taking a point on a centerline of a coronary artery as a reference point, enabling a cross section passing through the point to be intersected with a surface mesh of a blood vessel to form a ring, determining whether a position corresponding to the cross section is a required cutting position according to the maximum distance between the point on the ring and the reference point, and avoiding the influence of the position of excessive stenosis of the coronary artery on a calculation result. Therefore, compared with the scheme of directly generating the blood vessel body mesh model in the prior art, the embodiment of the invention can solve the problems that the blood vessel body mesh model is directly generated in the prior art, errors exist in the generated body mesh caused by non-uniform thickness of each position of a coronary artery, and the obtained FFR result is low in accuracy.
[ description of the drawings ]
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for determining fractional flow reserve according to an embodiment of the present invention;
FIG. 2 is a schematic illustration of a blood vessel surface mesh constructed in accordance with an embodiment of the present invention;
FIG. 3 is a schematic diagram of an intermediate surface mesh provided by an embodiment of the present invention;
FIG. 4 is a schematic diagram of a volume grid provided by an embodiment of the present invention;
FIG. 5 is a schematic flow chart for obtaining boundary conditions of coronary artery outlet flow resistance according to an embodiment of the present invention;
FIG. 6 is a graphical representation of fractional flow reserve results provided by an embodiment of the present invention;
fig. 7 is a schematic diagram of a computer device according to an embodiment of the present invention.
[ detailed description ] embodiments
For better understanding of the technical solutions of the present invention, the following detailed descriptions of the embodiments of the present invention are provided with reference to the accompanying drawings.
It should be understood that the described embodiments are only some embodiments of the invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The terminology used in the embodiments of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the examples of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It should be understood that the term "and/or" as used herein is merely one type of association that describes an associated object, meaning that three relationships may exist, e.g., a and/or B may mean: a exists alone, A and B exist simultaneously, and B exists alone. In addition, the character "/" herein generally indicates that the former and latter related objects are in an "or" relationship.
The word "if" as used herein may be interpreted as "at … …" or "when … …" or "in response to a determination" or "in response to a detection", depending on the context. Similarly, the phrases "if determined" or "if detected (a stated condition or event)" may be interpreted as "when determined" or "in response to a determination" or "when detected (a stated condition or event)" or "in response to a detection (a stated condition or event)", depending on the context.
Example one
An embodiment of the present invention provides a method for determining fractional flow reserve, as shown in fig. 1, which is a schematic flow chart of the method for determining fractional flow reserve provided in the embodiment of the present invention, and the method may specifically include the following steps:
s101, a CTA image of the subject is acquired.
Alternatively, a computed tomography angiography CTA image may be obtained non-invasively, which can effectively avoid the risk of invasive detection to the subject.
S102, from the CTA image, a vessel surface mesh and a coronary artery centerline are constructed.
In one particular implementation, the surface mesh of the vessel can be constructed simultaneously with the coronary centerline, allowing for increased processing speed of the CTA image. In another embodiment, optimization processing may also be performed on the extracted coronary artery centerlines: detecting whether the length of the coronary artery central line exceeds a set threshold value, and if the length of the coronary artery central line does not exceed the set threshold value, reserving the central line of the part; otherwise, the part is deleted. Through the above operation, the influence of the vessel bifurcation on the subsequent operation can be eliminated.
S103, determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line.
It should be noted that, in the prior art, when a blood vessel mesh model is directly generated, the cutting position is usually determined manually and cannot be determined automatically, which results in a long time consuming FFR detection process. In the embodiment of the invention, the cutting position of the blood vessel surface mesh can be automatically determined according to the blood vessel surface mesh and the coronary artery central line, so that the time for manually determining the cutting position is effectively reduced, and the detection efficiency is improved; meanwhile, the information of the blood vessel surface grid and the coronary artery central line is jointly utilized, so that human errors can be effectively avoided, and the determined cutting position is more accurate.
S104, cutting the blood vessel surface mesh at the cutting position to generate an intermediate surface mesh.
And S105, generating a body mesh model of the blood vessel according to the intermediate surface mesh.
And S106, obtaining a fractional flow reserve result of the coronary artery according to the body mesh model and the boundary condition of the blood vessel.
For ease of understanding, the following describes implementation of the above steps.
In executing step S102, first, a three-dimensional image of a blood vessel needs to be extracted from the CTA image. Specifically, by acquiring computer-scanned (CT) data of a blood vessel in a CTA image, the CT data is matched with blood vessel theoretical model data to obtain an initial point and a traveling direction of a minimum cost path. In the process of advancing the minimum cost path, the vessel theoretical model is corrected in a segmented mode, and the advancing direction is optimized. And when the minimum cost path is finished, searching the minimum point from the end point in a reverse mode, generating a center line of the blood vessel, and checking and correcting the center line to ensure the correctness of the result. Wherein, the blood vessel may include, but is not limited to, at least one of a main vein, a coronary vein, an aorta, and a coronary artery.
Due to the difference of heart individuals and the error of registration, the blood vessel theoretical model is corrected in a segmented manner in different heart CT data of blood vessels to be extracted so as to improve the matching degree of the blood vessel theoretical model and the actual blood vessels, and the blood vessel three-dimensional image can be acquired fully automatically without human-computer interaction operation in the process.
Secondly, constructing surface grids of the blood vessels in the blood vessel three-dimensional image, wherein the area of unit grids in each blood vessel surface grid can be the same or different. Please refer to fig. 2, which is a schematic diagram of a blood vessel surface mesh constructed according to an embodiment of the present invention.
As shown in fig. 2, when the blood vessels are an aorta and a coronary artery, a first surface mesh is constructed on an outer surface of the aorta, and a second surface mesh is constructed on an outer surface of the coronary artery, and an area of a unit mesh in the first surface mesh is larger than an area of a unit mesh in the second surface mesh. The coronary artery is fine, unit grids with small areas are needed, the aorta is large, unit grids with relatively large areas are used, the accuracy of the calculation result can be guaranteed, and meanwhile, the calculation amount is reduced.
In executing the step S103, the cutting position of the blood vessel surface mesh is determined based on the roundness of the line of intersection formed by the cross section determined by the point on the centerline of the coronary artery and the blood vessel surface mesh. Specifically, points on the centerline of the coronary artery are taken as reference points, a cross section passing through the points intersects with a surface mesh of a blood vessel to form a ring, so that the distance between each point on the ring and the reference points is calculated to obtain the maximum distance, namely the roundness of the ring, the roundness of the ring is compared with a preset roundness threshold, and if the roundness of the ring is the preset roundness threshold, the position corresponding to the cross section is determined to be a required cutting position, thereby avoiding the influence of the position of the excessive stenosis of the coronary artery on the calculation result. Wherein the preset roundness threshold value is any value within the range of 0-1.
Specifically, the preset roundness threshold value is 0.6, the point 1 and the point 2 on the centerline of the coronary artery near the outlet are taken as reference points, the cross section passing through the point 1 intersects with the surface mesh of the blood vessel to form a ring 1, the cross section passing through the point 2 intersects with the surface mesh of the blood vessel to form a ring 2, at this time, the maximum distance between the point on the ring 1 and the reference point 1 is 0.3 and is not equal to the preset roundness threshold value, and the maximum distance between the point on the ring 2 and the reference point 2 is 0.6 and is equal to the preset roundness threshold value, so that the position corresponding to the cross section passing through the point 2 is determined as the required cutting position.
Based on this, the surface mesh of the blood vessel is cut at the position corresponding to the cross section passing through the point 2 to generate the middle surface mesh, and the entrance and exit of the middle surface mesh obtained by cutting are marked, so that the cutting of the model end face can be automatically realized without manual operation in the process. For example, the position corresponding to the cross section passing through point 2 is labeled as exit a. As shown in fig. 3, which is a schematic diagram of the middle surface mesh provided by the embodiment of the present invention, it can be seen that the outlet a of the coronary artery has been cut flat.
On the basis of the intermediate surface mesh, as shown in fig. 4, a mesh is also generated inside the blood vessel, that is, a volume mesh model of the blood vessel can be obtained, and can be used for simulating the blood flow conditions inside the blood vessel and on the surface, so that the calculation result has high accuracy.
In one embodiment, generating a volumetric mesh model of the blood vessel from the intermediate surface mesh may comprise the steps of: performing boundary protection on the middle surface grid by using a relaxation method, and performing volume sampling on the three-dimensional uniform grid to obtain a volume sampling point set; and extracting a body mesh model of the generated blood vessel from the obtained body sampling point set by using a three-dimensional Dirony triangulation method. In this embodiment, boundary protecting the intermediate surface mesh using a relaxation method and performing volume sampling on the three-dimensional uniform mesh to obtain a volume sampling point set may include: relaxing the surface sampling point set; and carrying out volume sampling on the three-dimensional uniform grid by using the relaxed surface point set to obtain a volume sampling point set.
In a specific application scenario, a fractional flow reserve FFR result of a coronary artery can be obtained according to a volume mesh model of a blood vessel and boundary conditions, and the result is used for representing the degree of myocardial ischemia.
First, boundary conditions need to be acquired. Wherein the boundary conditions may include, but are not limited to, one or more of aortic inlet velocity boundary conditions, aortic outlet pressure boundary conditions, and coronary outlet flow resistance boundary conditions in combination.
The aortic inlet velocity boundary is calculated from cardiac output and the aortic outlet pressure boundary is calculated from blood pressure. The cardiac output is obtained directly by means of ultrasound measurement, or by multiplying the heart rate by the difference between the volume of the left ventricle at end systole and end diastole.
Referring to fig. 5, a schematic flow chart of obtaining a boundary condition of a coronary artery outlet flow resistance according to an embodiment of the present invention includes the following steps:
s501, at least one of cardiac output and myocardial mass is obtained.
Wherein the myocardial mass is the mass of the left ventricular myocardium wall.
S502, acquiring total coronary flow according to at least one of cardiac output and myocardial mass.
The total coronary flow is obtained according to the relationship that the total coronary flow is in direct proportion to cardiac output or the relationship that the total coronary flow is in direct proportion to myocardial mass.
And S503, acquiring a boundary condition of the flow resistance of the coronary artery outlet according to the total flow of the coronary artery and the parameters of the coronary artery.
The parameters of the coronary artery may include, but are not limited to, one or more of a coronary bifurcation cross-sectional area, a coronary outlet cross-sectional area, and a coronary artery luminal density Attenuation Gradient (TAG). For example: the parameters of the coronary artery select the coronary artery exit cross-sectional area, which is generally proportional to the FFR detection results, i.e.: the larger the cross-sectional area, the greater the blood flow rate therethrough, whereas the smaller the cross-sectional area, the smaller the blood flow rate therethrough. When TAG is chosen as the coronary parameter, it is generally inversely proportional to the FFR detection result, i.e.: the larger the cross-sectional area, the smaller the blood flow rate therethrough, whereas the smaller the cross-sectional area, the larger the blood flow rate therethrough.
After the boundary condition determination is completed, the hyperemic state is simulated in two ways:
first, the flow at the aorta entrance end is constant, and the flow at each coronary artery exit end is proportionally increased, and the velocity at each coronary artery exit end is proportionally increased.
Secondly, the equivalent flow resistance at the outlet end of each coronary artery is reduced in an equal proportion, and the flow at the outlet end of each coronary artery cannot be increased in an equal proportion, so that the method is closer to the actual situation, and the accuracy of the calculation result is high.
Based on this, according to the blood vessel volume mesh model and the boundary conditions, the blood vessel congestion state is simulated, and further a fractional flow reserve FFR result is obtained, please refer to fig. 6, which is a schematic diagram of the fractional flow reserve result provided by the embodiment of the present invention, as shown in fig. 6, which is a gray scale of the fractional flow reserve FFR result, where the variation range of all pixels in the whole gray scale fluctuates between 0.75 to 1, where: the higher the gray value is, the higher the blood flow reserve fraction of the corresponding part is, and the higher the blood flow of the part is; conversely, a lower gray scale value indicates a lower fractional flow reserve at the corresponding region, and the blood flow at that region is lower, and a position with a fractional flow reserve FFR result of 0.75 or less is generally considered medically ischemic. As shown in fig. 6, there is a relatively severe ischemic event at the end of the left coronary artery, which corresponds to a clinically constricted region of a blood vessel.
On the basis of the above embodiments, the present embodiment provides a computer-readable storage medium.
In particular, the computer-readable storage medium comprises computer-executable instructions, which when executed, are capable of performing the method of determining fractional flow reserve presented in this embodiment.
In an embodiment of the invention, the cut location is performed by constructing the surface mesh of the vessel and the coronary artery centerline in a computed tomography angiography CTA image, based on which, taking a point on the centerline of the coronary artery as a reference point, a cross section passing through the point can intersect with the surface mesh of the blood vessel to form a ring, whether the position corresponding to the cross section is the required cutting position can be determined according to the maximum distance between the point on the ring and the reference point, the influence of the position of the coronary artery over-stenosis on the calculation result is avoided, and therefore, cutting the surface mesh of the blood vessel at the cutting position to generate an intermediate surface mesh, and further generating a mesh inside the blood vessel based on the intermediate surface mesh, the body mesh model of the blood vessel can be obtained, and can be used for simulating the blood flow condition inside and on the surface of the blood vessel, so that the calculation result has high accuracy. Therefore, compared with the scheme of directly generating the blood vessel body mesh model in the prior art, the embodiment of the invention can solve the problems that the blood vessel body mesh model is directly generated in the prior art, errors exist in the generated body mesh caused by non-uniform thickness of each position of a coronary artery, and the obtained FFR result is low in accuracy; the cutting position is automatically determined, and the problem of long time consumption is avoided.
Example two
Based on the method for determining fractional flow reserve provided in the first embodiment, an embodiment of the present invention further provides a computer device for implementing the method. Specifically, please refer to fig. 7, which is a schematic diagram illustrating a computer device according to an embodiment of the present invention. As shown in fig. 7, the computer device 71 comprises a memory 701, a processor 702 and a computer program 703 (not shown in fig. 7) stored on the memory 701 and executable on the processor 702.
A memory 701 for storing a computer program 703 executable on a processor 702;
a processor 702 configured to perform the following operations when executing the computer program 703:
acquiring a CTA image of a subject;
constructing a vessel surface mesh and a coronary artery centerline from the CTA image;
determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line;
cutting the blood vessel surface mesh at the cutting position to generate an intermediate surface mesh;
generating a body mesh model of the blood vessel according to the intermediate surface mesh;
and obtaining a fractional flow reserve result of the coronary artery according to the body mesh model and the boundary condition of the blood vessel.
In a specific implementation, the processor 702, when configured to execute the computer program 703, further performs the following operations:
in response to the fractional flow reserve result for a coronary artery being less than a set threshold, a stenosis location is determined for the coronary artery.
In a particular implementation, the processor 702 is further configured to construct a first surface mesh on an outer surface of the aorta when the vessels are the aorta and the coronary arteries;
constructing a second surface mesh on the outer surface of the coronary artery;
the area of the unit cells in the first surface mesh is larger than the area of the unit cells in the second surface mesh.
In one particular implementation, the processor 702 is further configured to determine a cutting location of the blood vessel surface mesh based on a roundness of an intersection line loop formed by the cross-section and the blood vessel surface mesh determined from a point on the centerline of the coronary artery.
In a specific implementation, the processor 702 is further configured to check and correct the coronary artery centerline before determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery centerline to ensure the correctness of the result.
In one particular implementation, the processor 702 is further configured to obtain boundary conditions, wherein the boundary conditions may include, but are not limited to, one or more combinations of aortic inlet velocity boundary conditions, aortic outlet pressure boundary conditions, and coronary outlet flow resistance boundary conditions.
In one particular implementation, the processor 702 is further configured to obtain at least one of cardiac output and myocardial mass;
obtaining total coronary flow based on at least one of cardiac output and myocardial mass;
and acquiring a boundary condition of the flow resistance of the coronary artery outlet according to the total flow of the coronary artery and the parameters of the coronary artery. The parameters of the coronary artery may include, but are not limited to, one or more of a coronary bifurcation cross-sectional area, a coronary outlet cross-sectional area, and a coronary intraluminal Attenuation Gradient (TAG).
In one embodiment, the coronary parameters select TAG, which is the linear regression coefficient between the radioactive decay value in the lumen of the coronary artery and the axial distance from the coronary ostium to that point, i.e., the amount of change in CT value (HU) per 10mm unit length interval from the coronary ostium. Generally, the larger the TAG, the more likely a stenosis in the vessel is present; smaller TAG indicates better vascular uniformity. In this embodiment, the method of measuring TAG may comprise: firstly, reconstructing a cross-sectional image of the coronary artery perpendicular to the centerline of the blood vessel; then, the cross-sectional area from the coronary ostia to the vessel endings is larger than the set area (e.g. 2 mm)2) Measuring variables such as the cross-sectional area of the lumen, the average diameter of the lumen and the radioactive attenuation value of the lumen at certain intervals (for example, 5 mm); the contour of the target region and the centerline of the vessel may be corrected manually if necessary. TAG may represent the amount of change in CT attenuation (HU) per 10mm of intra-coronary artery. In this embodiment, a boundary condition is determined from the TAG and applied to the calculation of the FFR, taking into account that the TAG is at oneThe degree of the stenosis of the blood vessel is represented in a fixed degree, and the calculation precision of the FFR can be improved.
Since each unit in this embodiment can execute the method shown in the first embodiment, reference may be made to the related description of the first embodiment for a part of this embodiment that is not described in detail.
In an embodiment of the invention, the cut location is performed by constructing a surface mesh of the vessel and the coronary artery centerline in a computed tomography angiography CTA image, based on which points on the coronary artery centerline are taken. The cross section passing through the point is a reference point and can be intersected with the surface mesh of the blood vessel to form a ring, whether the position corresponding to the cross section is a required cutting position or not can be determined according to the maximum distance between the point on the ring and the reference point, and the influence of the position of the over-narrow coronary artery on a calculation result is avoided. Therefore, compared with the scheme of directly generating the blood vessel body mesh model in the prior art, the embodiment of the invention can solve the problems that the blood vessel body mesh model is directly generated in the prior art, errors exist in the generated body mesh caused by non-uniform thickness of each position of a coronary artery, and the obtained FFR result is low in accuracy.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the embodiments provided in the present invention, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions in actual implementation, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, or in a form of hardware plus a software functional unit.
The integrated unit implemented in the form of a software functional unit may be stored in a computer readable storage medium. The software functional unit is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device) or a Processor (Processor) to execute some steps of the methods according to the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (9)
1. A method of determining fractional flow reserve, the method comprising:
acquiring a CTA image of a subject;
constructing a vessel surface mesh and coronary artery centerlines from the CTA image;
determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line;
cutting the vessel surface mesh at the cutting location, generating an intermediate surface mesh;
generating a body mesh model of the blood vessel according to the intermediate surface mesh;
obtaining a fractional flow reserve result of the coronary artery according to the volume mesh model and the boundary condition of the blood vessel;
the blood vessel is an aorta and a coronary artery, the blood vessel surface mesh comprises a first surface mesh and a second surface mesh, and the area of the unit mesh in the first surface mesh is larger than that of the unit mesh in the second surface mesh;
the first surface mesh is an outer surface mesh of the aorta;
the second surface mesh is an outer surface mesh of the coronary artery.
2. The method of claim 1, wherein prior to obtaining the fractional flow reserve result for the coronary artery based on the volumetric mesh model and the boundary conditions of the vessel, the method further comprises:
acquiring the boundary condition, wherein the boundary condition comprises at least one of an aorta inlet speed boundary condition, an aorta outlet pressure boundary condition and a coronary artery outlet flow resistance boundary condition.
3. The method of claim 2, wherein obtaining the coronary outlet flow resistance boundary condition comprises:
obtaining at least one of cardiac output and myocardial mass;
obtaining the total coronary flow based on at least one of the cardiac output and myocardial mass;
and acquiring the boundary condition of the coronary artery outlet flow resistance according to the total coronary artery flow and the parameters of the coronary artery.
4. The method of claim 3, wherein the coronary parameter comprises at least one of a coronary bifurcation cross-sectional area, a coronary exit cross-sectional area, and a coronary intraluminal density decay gradient.
5. The method of claim 1, wherein determining the cutting location of the vessel surface mesh from the vessel surface mesh and a coronary artery centerline comprises:
and determining the cutting position of the blood vessel surface mesh according to the roundness of an intersection line ring formed by the cross section determined by the point on the centerline of the coronary artery and the blood vessel surface mesh.
6. The method of claim 1, wherein prior to determining the cutting location of the vessel surface mesh from the vessel surface mesh and a coronary artery centerline, the method further comprises:
the coronary centerline is examined and corrected.
7. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor is configured to perform the following operations when the computer program is executed:
acquiring a CTA image of a subject;
constructing a vessel surface mesh and coronary artery centerlines from the CTA image;
determining the cutting position of the blood vessel surface mesh according to the blood vessel surface mesh and the coronary artery central line;
cutting the vessel surface mesh at the cutting location, generating an intermediate surface mesh;
generating a body mesh model of the blood vessel according to the intermediate surface mesh;
obtaining a fractional flow reserve result of the coronary artery according to the volume mesh model and the boundary condition of the blood vessel;
the blood vessel is an aorta and a coronary artery, and the area of the unit grid in the first surface grid is larger than that of the unit grid in the second surface grid;
the first surface mesh is an outer surface mesh of the aorta;
the second surface mesh is an outer surface mesh of the coronary artery.
8. The computer device of claim 7, wherein the processor, when executing the computer program, is further configured to:
in response to the fractional flow reserve result for the coronary artery being less than a set threshold, a stenosis location is determined at the coronary artery.
9. A computer-readable storage medium, comprising: computer-executable instructions which, when executed, perform the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810320618.4A CN108564574B (en) | 2018-04-11 | 2018-04-11 | Method, computer device and computer-readable storage medium for determining fractional flow reserve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810320618.4A CN108564574B (en) | 2018-04-11 | 2018-04-11 | Method, computer device and computer-readable storage medium for determining fractional flow reserve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108564574A CN108564574A (en) | 2018-09-21 |
| CN108564574B true CN108564574B (en) | 2021-04-20 |
Family
ID=63534473
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810320618.4A Active CN108564574B (en) | 2018-04-11 | 2018-04-11 | Method, computer device and computer-readable storage medium for determining fractional flow reserve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108564574B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110706770B (en) * | 2019-09-30 | 2020-08-04 | 上海杏脉信息科技有限公司 | Cardiac data processing apparatus, cardiac data processing method, and computer-readable storage medium |
| CN112116615B (en) * | 2019-11-19 | 2023-12-05 | 苏州润迈德医疗科技有限公司 | Method and device for acquiring blood vessel contour line according to blood vessel center line |
| CN111738982A (en) * | 2020-05-28 | 2020-10-02 | 数坤（北京）网络科技有限公司 | Vascular cavity concentration gradient extraction method and device and readable storage medium |
| CN111754522B (en) * | 2020-06-19 | 2021-08-03 | 上海杏脉信息科技有限公司 | Method and device for acquiring coronary artery hemodynamic data |
| CN112001893B (en) * | 2020-07-31 | 2024-02-20 | 上海联影医疗科技股份有限公司 | Calculation method, device and equipment of vascular parameters and storage medium |
| CN112690814B (en) * | 2020-11-06 | 2022-10-14 | 杭州阿特瑞科技有限公司 | Low-error coronary artery blood flow reserve fraction measuring method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101152104A (en) * | 2006-04-21 | 2008-04-02 | 美国西门子医疗解决公司 | System and method for semi-automatic aortic aneurysm analysis |
| CN106023202A (en) * | 2016-05-20 | 2016-10-12 | 苏州润心医疗科技有限公司 | Coronary artery fractional flow reserve calculation method based on heart CT image |
| CN106327487A (en) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | Coronary artery blood flow reserve fraction calculation method based on X ray coronary artery angiographic image |
| CN106659400A (en) * | 2014-06-30 | 2017-05-10 | 皇家飞利浦有限公司 | Apparatus for determining a fractional flow reserve value |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8315812B2 (en) * | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
-
2018
- 2018-04-11 CN CN201810320618.4A patent/CN108564574B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101152104A (en) * | 2006-04-21 | 2008-04-02 | 美国西门子医疗解决公司 | System and method for semi-automatic aortic aneurysm analysis |
| CN106659400A (en) * | 2014-06-30 | 2017-05-10 | 皇家飞利浦有限公司 | Apparatus for determining a fractional flow reserve value |
| CN106023202A (en) * | 2016-05-20 | 2016-10-12 | 苏州润心医疗科技有限公司 | Coronary artery fractional flow reserve calculation method based on heart CT image |
| CN106327487A (en) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | Coronary artery blood flow reserve fraction calculation method based on X ray coronary artery angiographic image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108564574A (en) | 2018-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108564574B (en) | Method, computer device and computer-readable storage medium for determining fractional flow reserve | |
| US11896416B2 (en) | Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images | |
| US20230360803A1 (en) | Methods and systems for predicting sensitivity of blood flow calculations to changes in anatomical geometry | |
| EP2966615B1 (en) | Method and system for improved hemodynamic computation in coronary arteries | |
| US10463336B2 (en) | Method and system for purely geometric machine learning based fractional flow reserve | |
| CN111754506B (en) | Coronary artery stenosis rate calculation method, device and system based on intra-cavity image and computer storage medium | |
| US10134129B2 (en) | Method and system for hemodynamic computation in coronary arteries | |
| US10595806B2 (en) | Fractional flow reserve (FFR) index with adaptive boundary condition parameters | |
| EP3394829B1 (en) | Heart model guided coronary artery segmentation | |
| EP3140757B1 (en) | Method and system for non-invasive functional assessment of coronary artery stenosis using flow computations in models based on diseased patients and hypothetical normal anatomical models | |
| JP7236769B2 (en) | Calculation method of microcirculatory resistance index based on myocardial blood flow and CT images | |
| WO2019169779A1 (en) | Method for measuring fractional flow reserve without vasodilator | |
| CN110638438A (en) | Method and system for machine learning-based assessment of fractional flow reserve | |
| EP3244790B1 (en) | Instantaneous wave-free ratio (ifr) computer tomography (ct) | |
| CN113040795B (en) | Detection method for non-guide wire FFR, non-guide wire IMR and non-guide wire CFR | |
| EP3820357B1 (en) | Patient-specific modeling of hemodynamic parameters in coronary arteries | |
| CN116090364A (en) | Method for obtaining coronary blood flow reserve fraction based on CTA image and readable storage medium | |
| US10354349B2 (en) | Systems and methods for using geometry sensitivity information for guiding workflow | |
| CN111652881A (en) | Coronary artery reconstruction and fractional flow reserve calculation method, device and equipment based on deep learning and readable storage medium | |
| CN113180614A (en) | Detection method for non-guide wire FFR, non-guide wire IMR and non-guide wire CFR | |
| KR102000615B1 (en) | A method for automatically extracting a starting point of coronary arteries, and an apparatus thereof | |
| JP2015217113A (en) | Blood vessel analysis device, medical image diagnostic device, blood vessel analysis method, and blood vessel analysis program | |
| CN117036530B (en) | Cross-modal data-based coronary artery fractional flow reserve prediction method and device | |
| US20240070855A1 (en) | Method and System of Calculating the Cross-section Area and Included Angle of Three-dimensional Blood Vessel Branch | |
| CN116313101A (en) | Method, system, equipment and medium for determining fractional flow reserve of coronary artery |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |