CN108564574A - Determine method, computer equipment and the computer readable storage medium of blood flow reserve score - Google Patents
Determine method, computer equipment and the computer readable storage medium of blood flow reserve score Download PDFInfo
- Publication number
- CN108564574A CN108564574A CN201810320618.4A CN201810320618A CN108564574A CN 108564574 A CN108564574 A CN 108564574A CN 201810320618 A CN201810320618 A CN 201810320618A CN 108564574 A CN108564574 A CN 108564574A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- coronary artery
- surface grid
- center line
- blood
- 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.)
- Granted
Links
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
An embodiment of the present invention provides a kind of method, computer equipment and the computer readable storage medium of determining blood flow reserve score, the method includes：Obtain the CTA images of subject；According to CTA images, blood vessel surface grid and coronary artery center line are constructed；According to blood vessel surface grid and coronary artery center line, the cutting position of blood vessel surface grid is determined；In location of cut cutting vessel surface mesh, intermediate surface grid is generated；According to intermediate surface grid, the volume mesh model of blood vessel is generated；According to the volume mesh model and boundary condition of blood vessel, blood flow reserve fractional result coronarius is obtained.Technical solution provided in an embodiment of the present invention can automatically determine cutting position, avoid taking longer problem, while improve the precision that vascular flow deposit score calculates.
Description
【Technical field】
The present invention relates to technical field of image processing more particularly to a kind of method of determining blood flow reserve score, computers
Equipment and computer readable storage medium.
【Background technology】
Coronary flow reserve score (Fractional Flow Reserve, FFR) is that a kind of assessment coronary artery is narrow
The functional attributes of degree of myocardial ischemia caused by narrow.With the development of computer technology and Medical Imaging Technology, will calculate
Machine tomography angiography (Computed Tomographic Angiography, CTA) image and Fluid Mechanics Computation phase
In conjunction with can non-invasively obtain FFR results by building vascular pattern and fluid mechanical emulation.
For this problem of structure vascular pattern, the scheme used in the prior art is to directly generate blood vessel volume mesh mould
Type, the blood vessel include aorta and coronary artery, specifically, after being partitioned into aorta and coronary artery in CTA images, no
Any processing is done to coronary artery, just directly generates volume mesh model.
However, in the prior art, directly generating vascular bodies grid model, do not consider that each position thickness of coronary artery is uneven
Even to cause the volume mesh generated there are error, the FFR result accuracy of acquisition is low.In addition, in the prior art, directly generating blood vessel
Volume mesh model and the manually determined cutting position in generating process can not automatically determine cutting position, thus cause
Detection is time-consuming long, complicated for operation.In consideration of it, it is necessary to be improved to existing FFR detection methods.
【Invention content】
In view of this, an embodiment of the present invention provides a kind of method, computer equipment and the meters of determining blood flow reserve score
Calculation machine readable storage medium storing program for executing does not consider that coronary artery is each to solve to directly generate vascular bodies grid model in the prior art
Position thickness is uneven to cause the volume mesh generated there are error, the low problem of the FFR result accuracy of acquisition.
On the one hand, an embodiment of the present invention provides a kind of methods of determining blood flow reserve score, which is characterized in that the side
Method includes：
Obtain the CTA images of subject；
According to the CTA images, blood vessel surface grid and coronary artery center line are constructed；
According to the blood vessel surface grid and coronary artery center line, the cutting position of the blood vessel surface grid is determined；
The blood vessel surface grid is cut in the location of cut, generates intermediate surface grid；
According to the intermediate surface grid, the volume mesh model of blood vessel is generated；
According to the volume mesh model and boundary condition of the blood vessel, the blood flow reserve score knot coronarius is obtained
Fruit.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, according to the blood
The volume mesh model and boundary condition of pipe, before obtaining the blood flow reserve fractional result coronarius, the method is also wrapped
It includes：
The boundary condition is obtained, the boundary condition includes aorta entrance velocity boundary condition, aorta outlet pressure
At least one of force boundary condition and coronary artery outlet flow resistance boundary condition.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method obtains the hat
Shape artery exports flow resistance boundary condition, including：
Obtain at least one of cardiac output and myocardial mass；
According at least one of the cardiac output and myocardial mass, the coronary artery total flow is obtained；
According to the coronary artery total flow and the parameter coronarius, coronary artery outlet flow resistance side is obtained
Boundary's condition.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, it is described coronal dynamic
The parameter of arteries and veins includes coronary artery bifurcated mouth cross-sectional area, coronary artery exit cross-sectional area and coronary artery intracavitary density decay
At least one of gradient.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, when the blood vessel
When for aorta with coronary artery, the area of unit grids is more than unit grids in second surface grid in first surface grid
Area；
The first surface grid is the appearance surface grids of the aorta；
The second surface grid is the appearance surface grids coronarius.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, according to the blood
Pipe surface grid and coronary artery center line determine that the cutting position of the blood vessel surface grid, the method further include：
According to the friendship of cross section and the blood vessel surface grid protocol determined by the point on the coronary artery center line
The circularity of wire loop determines the cutting position of the blood vessel surface grid.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, according to the blood
Pipe surface grid and coronary artery center line, before the cutting position for determining the blood vessel surface grid, the method further includes：
The coronary artery center line is checked and corrected.
On the other hand, it an embodiment of the present invention provides a kind of computer equipment, including memory, processor and is stored in
On reservoir and the computer program that can run on a processor, realized such as when the processor is for executing the computer program
Lower operation:
Obtain the CTA images of subject；
According to the CTA images, blood vessel surface grid and coronary artery center line are constructed；
According to the blood vessel surface grid and coronary artery center line, the cutting position of the blood vessel surface grid is determined；
The blood vessel surface grid is cut in the location of cut, generates intermediate surface grid；
According to the intermediate surface grid, the volume mesh model of blood vessel is generated；
According to the volume mesh model and boundary condition of the blood vessel, the blood flow reserve score knot coronarius is obtained
Fruit.
The aspect and any possible implementation manners as described above, it is further provided a kind of realization method, the processor
Following operation is also realized when for executing the computer program:
It is less than given threshold in response to the blood flow reserve fractional result coronarius, is determined in the coronary artery narrow
Narrow position.
In another aspect, an embodiment of the present invention provides a kind of computer readable storage mediums, including：Computer is executable to be referred to
It enables, the method that any of the above-described kind of determining blood flow reserve score is executed when the computer executable instructions are run.
In the embodiment of the present invention, by the surface mesh for constructing computer tomography angiography CTA image medium vessels
It was taken caused by effectively avoiding result in manual positioning cutting position with coronary artery center line to carry out automatic cutting positioning
Long problem simplifies entire detection operation；Using the point on coronary artery center line as reference point, by the cross section meeting of the point and blood
The surface mesh of pipe is crossed to form annular, can determine that cross section institute is right according to the maximum distance of point and reference point in annular
Whether the position answered is the cutting position needed, avoids influence of the excessively narrow position of coronary artery to result of calculation, from
And in the surface mesh of the location of cut cutting vessel, intermediate surface grid is generated, in turn, according to intermediate surface grid,
Internal blood vessel also generates grid, you can obtains the volume mesh model of blood vessel, can be used for the blood inside simulated blood vessel with surface
Mobility status so that result of calculation accuracy is high.Therefore, compared to the prior art in directly generate the side of vascular bodies grid model
Case, the embodiment of the present invention can solve to directly generate vascular bodies grid model in the prior art, not consider each position of coronary artery
Set that thickness is uneven to cause the volume mesh generated there are error, the low problem of the FFR result accuracy of acquisition.
【Description of the drawings】
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this field
For those of ordinary skill, without creative efforts, other attached drawings are can also be obtained according to these attached drawings.
Fig. 1 is the flow diagram of the method for the determination blood flow reserve score that the embodiment of the present invention is provided；
Fig. 2 is the schematic diagram for the blood vessel surface grid that the embodiment of the present invention is constructed；
Fig. 3 is the intermediate surface grid schematic diagram that the embodiment of the present invention is provided；
Fig. 4 is the volume mesh schematic diagram that the embodiment of the present invention is provided；
Fig. 5 is the flow diagram for obtaining coronary artery outlet flow resistance boundary condition that the embodiment of the present invention is provided；
Fig. 6 is the schematic diagram for the blood flow reserve fractional result that the embodiment of the present invention is provided；
Fig. 7 is the composition schematic diagram for the computer equipment that the embodiment of the present invention is provided.
【Specific implementation mode】
For a better understanding of the technical solution of the present invention, being retouched in detail to the embodiment of the present invention below in conjunction with the accompanying drawings
It states.
It will be appreciated that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, those of ordinary skill in the art obtained without creative efforts it is all its
Its embodiment, shall fall within the protection scope of the present invention.
The term used in embodiments of the present invention is not intended to be limiting merely for the sake of the purpose of description specific embodiment
The present invention.In the embodiment of the present invention and "an" of singulative used in the attached claims, " described " and "the"
It is also intended to including most forms, unless context clearly shows that other meanings.
It should be appreciated that term "and/or" used herein is only a kind of incidence relation of description affiliated partner, indicate
There may be three kinds of relationships, for example, A and/or B, can indicate：Individualism A, exists simultaneously A and B, individualism B these three
Situation.In addition, character "/" herein, it is a kind of relationship of "or" to typically represent forward-backward correlation object.
Depending on context, word as used in this " if " can be construed to " ... when " or " when ...
When " or " in response to determination " or " in response to detection ".Similarly, depend on context, phrase " if it is determined that " or " if detection
(condition or event of statement) " can be construed to " when determining " or " in response to determination " or " when the detection (condition of statement
Or event) when " or " in response to detection (condition or event of statement) ".
Embodiment one
An embodiment of the present invention provides a kind of methods of determining blood flow reserve score, as shown in Figure 1, it is implemented for the present invention
The flow diagram of the method for the determination blood flow reserve score that example is provided, this method can specifically include following steps：
S101 obtains the CTA images of subject.
Optionally, computer tomography angiography CTA images noninvasive can acquire, this can be effectively
Avoid the invasive risk for detecting and being brought to subject.
S102 constructs blood vessel surface grid and coronary artery center line according to CTA images.
During a concrete implementation, construct the surface mesh of blood vessel and coronary artery center line can simultaneously into
Row, can improve the processing speed to CTA images.In another embodiment, may be used also for the coronary artery center line of extraction
Optimize processing：Whether the length for detecting coronary artery center line is more than given threshold, if the length of coronary artery center line
Degree is less than given threshold, then retains the center line of the part；Otherwise, the part is deleted.By aforesaid operations, blood vessel can be excluded
Influence of the bifurcated for subsequent operation.
S103 determines the cutting position of blood vessel surface grid according to blood vessel surface grid and coronary artery center line.
It should be pointed out that in the prior art, directly generate vascular bodies grid model, usual manually determined cutting position and
Cutting position can not be automatically determined, cause FFR detection process to take long.And the present invention is in embodiment, it can basis
Blood vessel surface grid and coronary artery center line, automatically determine the cutting position of blood vessel surface grid, effectively reduce artificial true
Determine the time of cutting position, improves detection efficiency；Meanwhile combining the letter that blood vessel surface grid and coronary artery center line is utilized
Breath, can effectively avoid human error, determining cutting position is more accurate.
S104 generates intermediate surface grid in location of cut cutting vessel surface mesh.
S105 generates the volume mesh model of blood vessel according to intermediate surface grid.
S106 obtains blood flow reserve fractional result coronarius according to the volume mesh model and boundary condition of blood vessel.
In order to make it easy to understand, hereinafter, the realization method to above steps illustrates.
When executing S102 steps, firstly, it is necessary to according to CTA images, the 3-D view of blood vessel is extracted.Specifically, lead to
Computer scanning (Computed Tomography, the CT) data for obtaining CTA image medium vessels are crossed, by the CT data and blood vessel
Theoretical model data are matched, to obtain the initial point and direction of travel of minimal cost path.In the row of minimal cost path
Into in the process, blood vessel theoretical model is segmented correction, and direction of travel is optimized.After minimal cost path is advanced, from
Terminal reverse search smallest point generates the center line of blood vessel, and the center line is checked and corrected, to ensure result just
True property.Wherein, blood vessel can include but is not limited at least one of cardinal vein, coronary vein, aorta and coronary artery.
Due to the error of the difference and registration of heart individual, in the cardiac CT data of different blood vessels to be extracted, blood vessel
Theoretical model will be segmented correction, and to improve the matching degree of blood vessel theoretical model and actual vessel, this process is not necessarily to man-machine friendship
Interoperability, you can fully automatically obtain blood vessel 3-D view.
Secondly, the surface mesh for constructing blood vessel 3-D view medium vessels, the face of unit grids in each blood vessel surface grid
Product can be identical, can not also be identical.Referring to FIG. 2, the signal of its blood vessel surface grid constructed by the embodiment of the present invention
Figure.
As shown in Fig. 2, when blood vessel is aorta and coronary artery, first surface net is constructed on the outer surface of aorta
Lattice construct second surface grid on outer surface coronarius, and the area of unit grids is more than second in first surface grid
The area of unit grids in surface mesh.Since coronary artery is tiny, need using the unit grids compared with small area, and aorta
It is coarse, then the unit grids of relatively large area are used, can ensure the precision of result of calculation in this way, while reducing operand.
When executing S103 steps, according to cross section determined by the point on coronary artery center line and blood vessel surface grid
The circularity of the intersection ring of formation, determines the cutting position of blood vessel surface grid.Specifically, it is with the point on coronary artery center line
Reference point can be crossed to form annular by the cross section of the point with the surface mesh of blood vessel, to calculate the every bit in annular
At a distance from reference point, maximum distance, the as circularity of the annular are obtained, the circularity of the annular and preset roundness threshold are done
Compare, if the circularity of the annular is preset roundness threshold, it is determined that the position corresponding to the cross section is the cutting needed
Position avoids influence of the excessively narrow position of coronary artery to result of calculation.Wherein, preset roundness threshold is 0~1 model
Enclose any one interior numerical value.
Specifically, it is 0.6 to preset roundness threshold, and with the point 1 close to exit on coronary artery center line, point 2 is reference
The cross section of point, passing point 1 can be crossed to form annular 1 with the surface mesh of blood vessel, and the cross section of passing point 2 can be with the table of blood vessel
Surface grids are crossed to form annular 2, at this point, the maximum distance of the point in annular 1 and reference point 1 is 0.3, not with default roundness threshold
It is equal, and the maximum distance of the point in annular 2 and reference point 2 is 0.6, it is equal with default roundness threshold, therefore, it is determined that pass through
Position corresponding to the cross section of point 2 is the cutting position needed.
Based on this, the surface mesh of cutting vessel at the position corresponding to the cross section of passing point 2 generates intermediate surface
Grid, and make marks to the entrance for cutting obtained intermediate surface grid, this process does not need manual operation, it will be able to
The automatically cutting of implementation model end face.For example, exit a will be labeled as at the position corresponding to the cross section of passing point 2.Such as
Shown in Fig. 3, the intermediate surface grid schematic diagram that is provided by the embodiment of the present invention, it can be seen that exit a coronarius
It has been cut flat with.
On the basis of intermediate surface grid, as shown in figure 4, also generating grid in the blood vessels, you can obtain blood vessel
Volume mesh model can be used for the blood flow situation inside simulated blood vessel with surface so that result of calculation accuracy is high.
In one embodiment, according to intermediate surface grid, the volume mesh model for generating blood vessel may include following steps：Make
Boundary protection is carried out to intermediate surface mesh with relaxation method, and carries out sampler body on three-dimensional uniform grid, obtains sampler body
Point set；The volume mesh model for generating blood vessel is extracted using three-dimensional Di Luoni Triangulation Algorithms to obtained sampler body point set.Herein
In embodiment, boundary protection is carried out to intermediate surface mesh using relaxation method, and sampler body is carried out on three-dimensional uniform grid,
Sampler body point set is obtained to may include：It relaxes to surface sampling point set；Using the surface point after relaxation, in three-dimensional uniform net
Sampler body is carried out on lattice, obtains sampler body point set.
In a specific application scenarios, it can be obtained coronal dynamic according to the volume mesh model and boundary condition of blood vessel
The blood flow reserve score FFR of arteries and veins is as a result, the result is used to characterize degree of myocardial ischemia.
Firstly, it is necessary to obtain boundary condition.Wherein, boundary condition can include but is not limited to aorta entrance velocity boundary
The combination of condition, aorta outlet pressure boundary condition and one or more of coronary artery outlet flow resistance boundary condition.
Aorta entrance velocity boundary condition is calculated to obtain by cardiac output, and aorta outlet pressure boundary condition is by blood pressure
It is calculated.Wherein, cardiac output is directly obtained by the method for ultrasonic measurement, alternatively, according to left ventricle in end-systole
The difference of the volume of volume and diastasis is multiplied to obtain with heart rate.
Coronary artery exports the acquisition methods of flow resistance boundary condition, referring to FIG. 5, what it was provided by the embodiment of the present invention
The flow diagram for obtaining coronary artery outlet flow resistance boundary condition, can specifically include following steps：
S501 obtains at least one of cardiac output and myocardial mass.
Wherein, myocardial mass is the quality of myocardium of left ventricle wall.
S502 obtains coronary artery total flow according at least one of cardiac output and myocardial mass.
Wherein, according to the coronary artery total flow relationship directly proportional to cardiac output, alternatively, according to coronary artery total flow
The relationship directly proportional to myocardial mass obtains coronary artery total flow.
S503 obtains coronary artery and exports flow resistance boundary condition according to coronary artery total flow and parameter coronarius.
Wherein, parameter coronarius can include but is not limited to coronary artery bifurcated mouth cross-sectional area, coronary artery goes out
In mouth cross-sectional area and coronary artery intracavitary density decay gradient (Transluminal Attenuation Gradient, TAG)
One or more.Such as：Parameter selection coronary artery exit cross-sectional area coronarius, it is usual with FFR testing results
It is directly proportional, i.e.,：Cross-sectional area is bigger, by blood flow flow it is bigger, conversely, cross-sectional area is smaller, by blood flow flow get over
It is small.When parameter selection TAG coronarius, usually it is inversely proportional with FFR testing results, i.e.,：Cross-sectional area is bigger, by blood
It is smaller to flow flow, conversely, cross-sectional area is smaller, by blood flow flow it is bigger.
After the completion of boundary condition determines, congestive state is simulated by the following two kinds mode：
The first, the flow of aorta arrival end is constant, the flow of each coronary artery outlet end of raising of equal proportion, this
When, the speed of each coronary artery outlet end also can equal proportion increase.
Second, the effective resistance of each coronary artery outlet end of reduction of equal proportion, since each coronary artery exports
The flow at end will not equal proportion raising, closer in actual conditions, result of calculation accuracy is high.
Based on this, according to vascular bodies grid model and boundary condition, simulated blood vessel congestive state, and then blood flow reserve is obtained
Score FFR is as a result, referring to FIG. 6, the schematic diagram of its blood flow reserve fractional result provided by the embodiment of the present invention, such as Fig. 6
The gray-scale map of blood flow reserve score FFR results is shown, between the variation range of all pixels is 0.75~1 in entire gray-scale map
Fluctuation, wherein：Gray value is higher, indicates that the blood flow reserve score of corresponding position is higher, the blood flow at the position is higher；Instead
It, gray value is lower, indicates that the blood flow reserve score of corresponding position is lower, the blood flow at the position is lower, under normal circumstances
Blood flow reserve score FFR results are medically considered ischemic in 0.75 position below.As shown in fig. 6, left coronal in figure
There are more serious ischemic phenomenon, which is corresponding hemadostewnosis region clinically for the end of artery.
On the basis of the above embodiments, the present embodiment provides a kind of computer readable storage medium.
Specifically, computer readable storage medium includes computer executable instructions, when the computer executable instructions
When being run, the method that is able to carry out the determination blood flow reserve score gone out given in the present embodiment.
In the embodiment of the present invention, by the surface mesh for constructing computer tomography angiography CTA image medium vessels
It is based on this, using the point on coronary artery center line as reference point, by this with coronary artery center line to carry out cutting positioning
The cross section of point can be crossed to form annular with the surface mesh of blood vessel, can be with according to the maximum distance of point and reference point in annular
It determines whether the position corresponding to the cross section is the cutting position needed, avoids the excessively narrow position of coronary artery to meter
The influence for calculating result, to, in the surface mesh of the location of cut cutting vessel, intermediate surface grid is generated, in turn, according to
Intermediate surface grid, also generates grid in the blood vessels, you can obtains the volume mesh model of blood vessel, can be used in simulated blood vessel
The blood flow situation in portion and surface so that result of calculation accuracy is high.Therefore, compared to the prior art in directly generate blood vessel
The scheme of volume mesh model, the embodiment of the present invention can solve to directly generate vascular bodies grid model in the prior art, not consider
The each position thickness of coronary artery is uneven to cause the volume mesh generated there are error, and the FFR result accuracy of acquisition is low to ask
Topic；Cutting position is automatically determined, avoids taking longer problem.
Embodiment two
Based on the method for the determination blood flow reserve score that above-described embodiment one is provided, the embodiment of the present invention further provides
Realize a kind of computer equipment of the above method.Specifically, referring to FIG. 7, its computer for being provided by the embodiment of the present invention
The composition schematic diagram of equipment.As shown in fig. 7, the computer equipment 71 includes memory 701, processor 702 and it is stored in storage
On device 701 and the computer program 703 (being not shown in Fig. 7) that can be run on processor 702.
Memory 701, for storing the computer program 703 that can be run on processor 702；
Processor 702 realizes following operation when for executing computer program 703:
Obtain the CTA images of subject；
According to CTA images, blood vessel surface grid and coronary artery center line are constructed；
According to blood vessel surface grid and coronary artery center line, the cutting position of blood vessel surface grid is determined；
In location of cut cutting vessel surface mesh, intermediate surface grid is generated；
According to intermediate surface grid, the volume mesh model of blood vessel is generated；
According to the volume mesh model and boundary condition of blood vessel, blood flow reserve fractional result coronarius is obtained.
During a concrete implementation, processor 702 is for when executing computer program 703, also realizing following behaviour
Make:
It is less than given threshold in response to blood flow reserve fractional result coronarius, narrow position is determined in the coronary artery
It sets.
During a concrete implementation, processor 702 is additionally operable to when blood vessel is aorta and coronary artery, in master
First surface grid is constructed on the outer surface of artery；
Second surface grid is constructed on outer surface coronarius；
The area of unit grids is more than the area of unit grids in second surface grid in first surface grid.
During a concrete implementation, processor 702 is additionally operable to be determined according to the point on coronary artery center line
Cross section and blood vessel surface grid protocol intersection ring circularity, determine the cutting position of blood vessel surface grid.
During a concrete implementation, processor 702 is additionally operable to according to blood vessel surface grid and coronary artery center
Line before the cutting position for determining blood vessel surface grid, is checked and is corrected to coronary artery center line, to ensure result
Correctness.
During a concrete implementation, processor 702 is additionally operable to obtain boundary condition, wherein boundary condition can be with
Including but not limited to aorta entrance velocity boundary condition, aorta outlet pressure boundary condition and coronary artery outlet flow resistance side
The combination of one or more of boundary's condition.
During a concrete implementation, processor 702 is additionally operable to obtain in cardiac output and myocardial mass at least
One；
According at least one of cardiac output and myocardial mass, coronary artery total flow is obtained；
According to coronary artery total flow and parameter coronarius, obtains coronary artery and export flow resistance boundary condition.Wherein,
Parameter coronarius can include but is not limited to coronary artery bifurcated mouth cross-sectional area, coronary artery exit cross-sectional area and hat
One or more of attenuation gradient (Transluminal Attenuation Gradient, TAG) in shape lumen of artery.
In one embodiment, parameter selection TAG coronarius, be coronary artery tube chamber in radiation decay depreciation with
It is to the linear regression coeffficient between the axial distance since coronary ostium, i.e., long per 10mm units away from coronary ostium
The variable quantity of degree interval CT values (HU).In general TAG is bigger, indicates that there may be narrow for blood vessel；TAG is smaller, indicates blood vessel
Uniformity is better.In this embodiment, the method for measuring TAG may include：First, coronary artery is rebuild perpendicular to vessel centerline
Cross sectional image；Then, setting area (such as 2mm is more than from coronary ostium to blood vessel distal cross-sectional area2) at, every
The variables such as tube chamber cross-sectional area, tube chamber average diameter and tube chamber radiation decay depreciation are measured every certain spacing (such as 5mm)；Mesh
The center line of the shape and blood vessel of marking region when necessary can manual correction.TAG can indicate CT pad values in every 10mm coronary arteries
(HU) variable quantity.In this embodiment, boundary condition is determined according to TAG, and the boundary condition is applied to the calculating of FFR,
In view of TAG characterizes the degree of hemadostewnosis to a certain extent, the precision of FFR calculating can be improved.
Method shown in embodiment one is able to carry out by each unit in this present embodiment, what the present embodiment was not described in detail
Part can refer to the related description to embodiment one.
In the embodiment of the present invention, by the surface mesh for constructing computer tomography angiography CTA image medium vessels
It is based on this, to carry out cutting positioning with the point on coronary artery center line with coronary artery center line.For reference point, by this
The cross section of point can be crossed to form annular with the surface mesh of blood vessel, can be with according to the maximum distance of point and reference point in annular
It determines whether the position corresponding to the cross section is the cutting position needed, avoids the excessively narrow position of coronary artery to meter
The influence for calculating result, to, in the surface mesh of the location of cut cutting vessel, intermediate surface grid is generated, in turn, according to
Intermediate surface grid, also generates grid in the blood vessels, you can obtains the volume mesh model of blood vessel, can be used in simulated blood vessel
The blood flow situation in portion and surface so that result of calculation accuracy is high.Therefore, compared to the prior art in directly generate blood vessel
The scheme of volume mesh model, the embodiment of the present invention can solve to directly generate vascular bodies grid model in the prior art, not consider
The each position thickness of coronary artery is uneven to cause the volume mesh generated there are error, and the FFR result accuracy of acquisition is low to ask
Topic.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided by the present invention, it should be understood that disclosed system, device and method can be with
It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit
It divides, only a kind of division of logic function, formula that in actual implementation, there may be another division manner, for example, multiple units or group
Part can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown
Or the mutual coupling, direct-coupling or communication connection discussed can be by some interfaces, device or unit it is indirect
Coupling or communication connection can be electrical, machinery or other forms.
The unit illustrated as separating component may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, you can be located at a place, or may be distributed over multiple
In network element.Some or all of unit therein can be selected according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, it can also
It is that each unit physically exists alone, it can also be during two or more units be integrated in one unit.Above-mentioned integrated list
The form that hardware had both may be used in member is realized, can also be realized in the form of hardware adds SFU software functional unit.
The above-mentioned integrated unit being realized in the form of SFU software functional unit can be stored in one and computer-readable deposit
In storage media.Above-mentioned SFU software functional unit is stored in a storage medium, including some instructions are used so that a computer
It is each that device (can be personal computer, server or network equipment etc.) or processor (Processor) execute the present invention
The part steps of embodiment the method.And storage medium above-mentioned includes：USB flash disk, mobile hard disk, read-only memory (Read-
Only Memory, ROM), random access memory (Random Access Memory, RAM), magnetic disc or CD etc. it is various
The medium of program code can be stored.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention
With within principle, any modification, equivalent substitution, improvement and etc. done should be included within the scope of protection of the invention god.
Claims (10)
1. a kind of method of determining blood flow reserve score, which is characterized in that the method includes：
Obtain the CTA images of subject；
According to the CTA images, blood vessel surface grid and coronary artery center line are constructed；
According to the blood vessel surface grid and coronary artery center line, the cutting position of the blood vessel surface grid is determined；
The blood vessel surface grid is cut in the location of cut, generates intermediate surface grid；
According to the intermediate surface grid, the volume mesh model of blood vessel is generated；
According to the volume mesh model and boundary condition of the blood vessel, the blood flow reserve fractional result coronarius is obtained.
2. according to the method described in claim 1, it is characterized in that, according to the volume mesh model and boundary condition of the blood vessel,
Before obtaining the blood flow reserve fractional result coronarius, the method further includes：
The boundary condition is obtained, the boundary condition includes aorta entrance velocity boundary condition, aorta outlet pressure side
At least one of boundary's condition and coronary artery outlet flow resistance boundary condition.
3. according to the method described in claim 2, it is characterized in that, obtaining coronary artery outlet flow resistance boundary condition, packet
It includes：
Obtain at least one of cardiac output and myocardial mass；
According at least one of the cardiac output and myocardial mass, the coronary artery total flow is obtained；
According to the coronary artery total flow and the parameter coronarius, coronary artery outlet flow resistance perimeter strip is obtained
Part.
4. according to the method described in claim 3, it is characterized in that, the parameter coronarius includes coronary artery bifurcated mouth
At least one of cross-sectional area, coronary artery exit cross-sectional area and coronary artery intracavitary density decay gradient.
5. according to the method described in claim 1, it is characterized in that, when the blood vessel is aorta and coronary artery, first
The area of unit grids is more than the area of unit grids in second surface grid in surface mesh；
The first surface grid is the appearance surface grids of the aorta；
The second surface grid is the appearance surface grids coronarius.
6. according to the method described in claim 1, it is characterized in that, described according in the blood vessel surface grid and coronary artery
Heart line determines that the cutting position of the blood vessel surface grid includes：
According to the intersection ring of cross section and the blood vessel surface grid protocol determined by the point on the coronary artery center line
Circularity, determine the cutting position of the blood vessel surface grid.
7. according to the method described in claim 1, it is characterized in that, according to the blood vessel surface grid and coronary artery center
Line, before the cutting position for determining the blood vessel surface grid, the method further includes：
The coronary artery center line is checked and corrected.
8. a kind of computer equipment, including memory, processor and storage are on a memory and the meter that can run on a processor
Calculation machine program, which is characterized in that following operation is realized when the processor is for executing the computer program:
Obtain the CTA images of subject；
According to the CTA images, blood vessel surface grid and coronary artery center line are constructed；
According to the blood vessel surface grid and coronary artery center line, the cutting position of the blood vessel surface grid is determined；
The blood vessel surface grid is cut in the location of cut, generates intermediate surface grid；
According to the intermediate surface grid, the volume mesh model of blood vessel is generated；
According to the volume mesh model and boundary condition of the blood vessel, the blood flow reserve fractional result coronarius is obtained.
9. computer equipment according to claim 8, which is characterized in that the processor is for executing the computer journey
Following operation is also realized when sequence:
It is less than given threshold in response to the blood flow reserve fractional result coronarius, narrow position is determined in the coronary artery
It sets.
10. a kind of computer readable storage medium, which is characterized in that including：Computer executable instructions, when the computer can
It executes instruction and executes method as described in any one of claim 1 to 7 when being run.
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 true CN108564574A (en) | 2018-09-21 |
| CN108564574B 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) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110706770A (en) * | 2019-09-30 | 2020-01-17 | 上海杏脉信息科技有限公司 | Cardiac data processing apparatus, cardiac data processing method, and computer-readable storage medium |
| CN111738982A (en) * | 2020-05-28 | 2020-10-02 | 数坤（北京）网络科技有限公司 | Vascular cavity concentration gradient extraction method and device and readable storage medium |
| CN111754522A (en) * | 2020-06-19 | 2020-10-09 | 上海杏脉信息科技有限公司 | Method and device for acquiring coronary artery hemodynamic data |
| CN111815598A (en) * | 2020-06-30 | 2020-10-23 | 上海联影医疗科技有限公司 | Blood vessel parameter calculation method, device, equipment and storage medium |
| CN112001893A (en) * | 2020-07-31 | 2020-11-27 | 上海联影医疗科技股份有限公司 | Blood vessel parameter calculation method, device, equipment and storage medium |
| CN112116615A (en) * | 2019-11-19 | 2020-12-22 | 苏州润迈德医疗科技有限公司 | Method and device for obtaining blood vessel contour line according to blood vessel central line |
| CN112690814A (en) * | 2020-11-06 | 2021-04-23 | 杭州阿特瑞科技有限公司 | Low-error coronary artery blood flow reserve fraction measuring method |
Citations (5)
| 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 |
| US20160232667A1 (en) * | 2010-08-12 | 2016-08-11 | Heartflow, Inc. | Method and system for image processing and patient-specific modeling of blood flow |
| 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 |
-
2018
- 2018-04-11 CN CN201810320618.4A patent/CN108564574B/en active Active
Patent Citations (5)
| 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 |
| US20160232667A1 (en) * | 2010-08-12 | 2016-08-11 | Heartflow, Inc. | Method and system for image processing and patient-specific modeling of blood flow |
| 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 |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110706770A (en) * | 2019-09-30 | 2020-01-17 | 上海杏脉信息科技有限公司 | Cardiac data processing apparatus, cardiac data processing method, and computer-readable storage medium |
| CN110706770B (en) * | 2019-09-30 | 2020-08-04 | 上海杏脉信息科技有限公司 | Cardiac data processing apparatus, cardiac data processing method, and computer-readable storage medium |
| CN112116615A (en) * | 2019-11-19 | 2020-12-22 | 苏州润迈德医疗科技有限公司 | Method and device for obtaining blood vessel contour line according to blood vessel central line |
| 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 |
| CN111754522A (en) * | 2020-06-19 | 2020-10-09 | 上海杏脉信息科技有限公司 | Method and device for acquiring coronary artery hemodynamic data |
| CN111754522B (en) * | 2020-06-19 | 2021-08-03 | 上海杏脉信息科技有限公司 | Method and device for acquiring coronary artery hemodynamic data |
| CN111815598A (en) * | 2020-06-30 | 2020-10-23 | 上海联影医疗科技有限公司 | Blood vessel parameter calculation method, device, equipment and storage medium |
| CN111815598B (en) * | 2020-06-30 | 2024-04-26 | 上海联影医疗科技股份有限公司 | Calculation method, device and equipment of vascular parameters and storage medium |
| CN112001893A (en) * | 2020-07-31 | 2020-11-27 | 上海联影医疗科技股份有限公司 | Blood vessel parameter calculation method, device, equipment and storage medium |
| CN112001893B (en) * | 2020-07-31 | 2024-02-20 | 上海联影医疗科技股份有限公司 | Calculation method, device and equipment of vascular parameters and storage medium |
| CN112690814A (en) * | 2020-11-06 | 2021-04-23 | 杭州阿特瑞科技有限公司 | Low-error coronary artery blood flow reserve fraction measuring method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108564574B (en) | 2021-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108564574A (en) | Determine method, computer equipment and the computer readable storage medium of blood flow reserve score | |
| US11896416B2 (en) | Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images | |
| US10580526B2 (en) | System and method for calculating vessel flow parameters based on angiography | |
| Muraru et al. | New speckle-tracking algorithm for right ventricular volume analysis from three-dimensional echocardiographic data sets: validation with cardiac magnetic resonance and comparison with the previous analysis tool | |
| CN110367965B (en) | Method, device and system for conveniently measuring coronary artery blood vessel evaluation parameters | |
| CN110226923B (en) | Method for measuring fractional flow reserve without vasodilator | |
| US10134129B2 (en) | Method and system for hemodynamic computation in coronary arteries | |
| EP3394829B1 (en) | Heart model guided coronary artery segmentation | |
| US10373700B2 (en) | Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance | |
| US9924869B2 (en) | Method and apparatus for determining blood flow required, method and apparatus for producing blood flow image, and method and apparatus for processing myocardial perfusion image | |
| CN108122616B (en) | Method for generating individual-specific cardiovascular models and use thereof | |
| CN105326486A (en) | Method and system for calculating blood vessel pressure difference and fractional flow reserve | |
| CN108992057A (en) | A kind of method and apparatus of determining coronary flow reserve score FFR | |
| CN108665449B (en) | Deep learning model and device for predicting blood flow characteristics on blood flow vector path | |
| JP2018509971A (en) | Method and apparatus for quantitative flow analysis | |
| CN112384136A (en) | Method, device and system for obtaining blood vessel evaluation parameters based on radiography images | |
| CN109009061B (en) | Calculation method and device for obtaining blood flow characteristic value based on blood pressure correction | |
| CN111091913A (en) | Microcirculation resistance calculation method based on fractional flow reserve and coronary artery CT (computed tomography) contrast images | |
| CN108742570B (en) | Device for acquiring blood vessel pressure difference based on coronary artery advantage type | |
| CN112089433B (en) | Coronary artery blood flow reserve fraction measuring method based on CTA and DSA | |
| CN109907772A (en) | The method and apparatus for obtaining coronary blood flow and blood flow velocity | |
| Schmidt et al. | Anatomic validation of a novel method for left ventricular volume and mass measurements with use of real-time 3-dimensional echocardiography | |
| CN112704505B (en) | Method for measuring coronary artery flow reserve fraction by using CTA (computed tomography angiography) and DSA (digital signal amplification) | |
| EP3886702B1 (en) | Most relevant x-ray image selection for hemodynamic simulation | |
| JP7227400B2 (en) | Method and device for acquiring coronary blood flow and blood flow velocity |
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 | ||
| 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. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |