CN112116711A - Synthetic method and device of circular truncated cone blood vessel mathematical model for hydrodynamics analysis - Google Patents
Synthetic method and device of circular truncated cone blood vessel mathematical model for hydrodynamics analysis Download PDFInfo
- Publication number
- CN112116711A CN112116711A CN201911237196.5A CN201911237196A CN112116711A CN 112116711 A CN112116711 A CN 112116711A CN 201911237196 A CN201911237196 A CN 201911237196A CN 112116711 A CN112116711 A CN 112116711A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- model
- vessel
- dimensional
- point
- 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
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 249
- 238000013178 mathematical model Methods 0.000 title claims abstract description 63
- 238000004458 analytical method Methods 0.000 title claims abstract description 61
- 238000010189 synthetic method Methods 0.000 title claims description 3
- 238000000034 method Methods 0.000 claims abstract description 84
- 239000002356 single layer Substances 0.000 claims abstract description 50
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 47
- 239000012530 fluid Substances 0.000 claims abstract description 22
- 239000010410 layer Substances 0.000 claims abstract description 21
- 210000004351 coronary vessel Anatomy 0.000 claims description 55
- 230000002792 vascular Effects 0.000 claims description 47
- 239000008280 blood Substances 0.000 claims description 16
- 210000004369 blood Anatomy 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 239000000843 powder Substances 0.000 claims description 6
- 238000005192 partition Methods 0.000 claims description 3
- 238000002601 radiography Methods 0.000 claims description 2
- 208000031481 Pathologic Constriction Diseases 0.000 abstract description 5
- 208000037804 stenosis Diseases 0.000 abstract description 5
- 230000036262 stenosis Effects 0.000 abstract description 5
- 230000017531 blood circulation Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000002586 coronary angiography Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000000605 extraction Methods 0.000 description 5
- 238000002583 angiography Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000003709 image segmentation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000013517 stratification Methods 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000001308 synthesis method Methods 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 206010011089 Coronary artery stenosis Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 150000002632 lipids Chemical group 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
-
- 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/11—Region-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/66—Analysis of geometric attributes of image moments or centre of gravity
-
- 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/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The application provides a method and a device for synthesizing a circular truncated cone blood vessel mathematical model for hydrodynamics analysis, wherein the method comprises the following steps: according to the real-time diameter D of the vesseltCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model; carrying out N-edge type mesh division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer mesh model; and carrying out surface layering treatment on the single-layer mesh model to form a double-layer mesh model, namely a blood vessel mathematical model. The method solves the problem that no blood vessel three-dimensional grid model for fluid mechanics analysis exists in the prior art, and fills the blank of the industry; since the wall of a blood vessel has a certain thickness and a stenosis problem occurs mainly in the inner wall of the blood vessel,therefore, the double-layer mesh model is built by the blood vessel mathematical model, the double-layer mesh model has a certain thickness, and the inner-layer mesh model has elasticity, can correct the blood flow speed and is closer to the real state of the blood vessel.
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a method, a device and a system for synthesizing a circular truncated cone blood vessel mathematical model for hydrodynamics analysis.
Background
The deposition of lipids and carbohydrates in human blood on the vessel wall will form plaques on the vessel wall, which in turn leads to vessel stenosis; especially, the blood vessel stenosis near the coronary artery of the heart can cause insufficient blood supply of cardiac muscle, induce diseases such as coronary heart disease, angina pectoris and the like, and cause serious threat to the health of human beings. According to statistics, about 1100 million patients with coronary heart disease in China currently have the number of patients treated by cardiovascular interventional surgery increased by more than 10% every year.
Although conventional medical detection means such as coronary angiography CAG and computed tomography CT can display the severity of coronary stenosis of the heart, the ischemia of the coronary cannot be accurately evaluated. In order to improve the accuracy of coronary artery function evaluation, Pijls in 1993 proposes a new index for estimating coronary artery function through pressure measurement, namely Fractional Flow Reserve (FFR), and the FFR becomes the gold standard for coronary artery stenosis function evaluation through long-term basic and clinical research.
FFR is one of the coronary artery blood vessel evaluation parameters, and microcirculation resistance index IMR and the like belong to the coronary artery blood vessel evaluation parameters.
In coronary angiography images, the coronary vessel assessment parameters need to be calculated in combination with fluid mechanics analysis (CFD), and no vessel mathematical model for fluid mechanics analysis exists in the prior art.
Disclosure of Invention
The invention provides a method, a device and a system for synthesizing a circular truncated cone mathematical model for fluid mechanics analysis, and aims to solve the problem that no mathematical model for fluid mechanics analysis exists in the prior art.
In order to achieve the above object, in a first aspect, the present application provides a method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, including:
according to the real-time diameter D of the vesseltCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model;
carrying out N-edge type grid division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer grid model, wherein N is more than or equal to 6;
and carrying out surface layering treatment on the single-layer grid model to form a double-layer grid model, namely a circular truncated cone blood vessel mathematical model.
Optionally, in the above method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the N-edge mesh partition is performed along the circumferential surface of the three-dimensional vascular model to form a single-layer mesh model, where N is greater than or equal to 6, and the method includes:
carrying out mesh division by taking a triangle as a minimum unit along the circumferential surface of the three-dimensional blood vessel model;
according to the sequence, every N triangles are combined and converted into 1N-sided polygon, and an N-sided polygon initial grid is formed;
and deleting the connecting lines inside each N-polygon in the N-polygon initial grids to form a single-layer N-polygon grid model, wherein N is more than or equal to 6.
Optionally, in the method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the method for meshing with a triangle as a minimum unit along the circumferential surface of the three-dimensional vascular model includes:
segmenting the three-dimensional vessel model into K segments,
and performing mesh division on the circumferential surface of each section of the three-dimensional blood vessel model by taking a triangle as a minimum unit.
Optionally, in the above method for synthesizing a circular truncated cone mathematical model for hydrodynamics analysis, the triangle as the minimum unit is an isosceles triangle.
Optionally, in the method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamic analysis, the method for performing surface stratification on the single-layer mesh model to form a double-layer mesh model, that is, a mathematical vascular model, includes:
obtaining the wall thickness h of the blood vessel;
according to the wall thickness h and the initial diameter D of the blood vesselGet upEnd of vessel diameter DPowderAnd in blood vesselsPerforming three-dimensional modeling on the length L of the core line, and forming a circular truncated cone three-dimensional model on the inner surface or the outer surface of the single-layer grid model;
according to the acquisition method of the single-layer grid model, carrying out N-edge type grid division along the circumferential surface of the circular truncated cone three-dimensional model to form another single-layer grid model;
and forming the double-layer mesh model, namely the blood vessel mathematical model, by using the two layers of the single-layer mesh model and the blood vessel wall thickness h.
Optionally, in the method for synthesizing the circular truncated cone mathematical vascular model for fluid mechanics analysis, the real-time diameter D of the blood vessel is used as the basistThe method for forming the three-dimensional blood vessel model by three-dimensional modeling of the length L of the central line of the blood vessel comprises the following steps:
acquiring coronary artery two-dimensional contrast images of at least two body positions;
obtaining the real-time diameter D of the blood vessel according to the coronary artery two-dimensional contrast imagetAnd the length L of the straightened blood vessel center line;
according to said DtAnd L three-dimensional modeling to form a three-dimensional model of the circular truncated cone.
Optionally, in the method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the real-time diameter D of the blood vessel is obtained according to the coronary artery two-dimensional angiography imagetAnd the method for the length L of the straightened blood vessel central line comprises the following steps:
extracting a blood vessel central line from the coronary artery two-dimensional contrast image of each posture along the direction from a coronary artery entrance to the coronary artery tail end;
acquiring a straightened vessel image according to the coronary artery two-dimensional radiography image and the vessel central line;
acquiring a straightened blood vessel contour line according to the straightened blood vessel center line and the straightened blood vessel image;
acquiring geometric information of the straightened vessel, comprising: real-time diameter D of blood vesseltAnd the length of the straightened blood vessel central line is the length L of the central straight line.
Alternatively, the above-mentioned for hydrodynamic divisionA method for synthesizing a circular truncated cone mathematical model of a blood vessel according to said DtAnd L three-dimensional modeling, wherein the method for forming the circular truncated cone three-dimensional model comprises the following steps:
performing three-dimensional modeling according to the geometric information, the central line and the contour line to obtain a three-dimensional blood vessel model;
real-time diameter D from the vesseltObtaining the starting diameter D of the blood vessel internallyGet upAnd end diameter D of vesselPowder；
According to said DGet up、DPowderAnd L, performing three-dimensional modeling to form the three-dimensional model of the circular truncated cone.
Optionally, in the above method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, after the acquiring of the two-dimensional coronary angiography images of at least two body positions, the real-time diameter D of the blood vessel is obtained according to the two-dimensional coronary angiography imagestObtaining the starting diameter D of the blood vessel internallyGet upAnd end diameter D of vesselPowderAnd the length L after the blood vessel central line is straightened further comprises:
acquiring a vessel segment of interest from the coronary artery two-dimensional angiogram;
picking up a starting point and an end point of the vessel segment of interest;
and segmenting local vessel region maps corresponding to the starting points and the end points from the coronary artery two-dimensional contrast image.
Optionally, in the method for synthesizing the circular truncated cone mathematical vascular model for hydrodynamic analysis, the method for segmenting the local vascular region map corresponding to the starting point and the ending point from the two-dimensional coronary angiography image further includes:
picking up at least one seed point of the vessel segment of interest;
and respectively segmenting the two-dimensional contrast images between the adjacent two points of the starting point, the seed point and the ending point to obtain at least two local blood vessel region images.
Optionally, in the method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the method for extracting a vessel centerline from the two-dimensional coronary angiography image of each posture along the direction from the coronary artery entrance to the coronary artery end comprises:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
performing mesh division on the rough blood vessel map, and extracting at least one blood vessel path line along the direction from the starting point to the end point;
and selecting one blood vessel path line as the blood vessel central line.
Optionally, in the above method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the method for meshing the rough blood vessel map and extracting at least one blood vessel path line along the direction from the starting point to the ending point includes:
gridding the rough vessel map;
searching the shortest time path between the starting point and the intersection points on the peripheral n grids along the extending direction of the blood vessels from the starting point to the ending point to serve as a second point, searching the shortest time path between the second point and the intersection points on the peripheral n grids to serve as a third point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and connecting the extending directions of the blood vessels from the starting point to the ending point according to the searching sequence to obtain at least one blood vessel path line.
Optionally, in the method for synthesizing the circular truncated cone mathematical vascular model for hydrodynamic analysis, the selecting one of the vascular path lines as the vascular centerline includes:
if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
the vessel path line that is the least in time is taken as the vessel centerline.
Optionally, in the method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, the method for extracting a vessel centerline from a coronary artery two-dimensional contrast image of each posture along a direction from a coronary artery entrance to a coronary artery end includes:
performing image processing on the local blood vessel region image to obtain a rough trend line of the blood vessel between the starting point and the ending point;
obtaining rough blood vessel edge lines, wherein images between the rough blood vessel edge lines including the rough blood vessel trend lines are blood vessel skeletons;
extracting the vessel centerline from the vascular skeleton.
Optionally, in the method for synthesizing a circular truncated vessel mathematical model for hydrodynamics analysis, the method for extracting the vessel centerline from the vascular skeleton includes:
performing grid division on the processed region image;
searching the blood vessel skeleton according to RGB values along the direction from the starting point to the ending point, searching a point where the minimum value of RGB difference values of the starting point and intersection points on m grids around the starting point is located as a second point, searching a point where the minimum value of RGB difference values of the second point and intersection points on m grids around the second point is located as a third point, and repeating the steps until the third point reaches the ending point, wherein m is a positive integer greater than or equal to 1;
obtaining at least one connecting line from the starting point to the end point according to the searching sequence;
and if the connecting line is two or more, selecting one connecting line as the center line of the blood vessel.
Optionally, in the method for synthesizing the circular truncated cone mathematical vascular model for hydrodynamics analysis, the method for obtaining the blood vessel contour line after straightening according to the blood vessel centerline after straightening and the blood vessel image after straightening includes:
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating a vascular preset on both sides of the central line of said vesselContour lines;
and gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the straightened contour line of the blood vessel.
In a second aspect, the present application provides an apparatus for synthesizing a mathematical model of a blood vessel, comprising: the three-dimensional model structure of the blood vessel, the structure of the single-layer grid model and the structure of the mathematical model of the blood vessel are connected in sequence, and the structure of the mathematical model of the blood vessel is connected with the structure of the three-dimensional model;
the three-dimensional blood vessel model structure is used for real-time diameter D according to blood vesselstCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model;
the single-layer grid model structure is used for carrying out N-edge type grid division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer grid model, wherein N is more than or equal to 6;
the vessel mathematical model structure is used for carrying out surface layering treatment on the single-layer grid model to form a double-layer grid model, namely the vessel mathematical model.
In a third aspect, the present application provides a coronary artery analysis system comprising: the device for synthesizing the blood vessel mathematical model is described above.
In a fourth aspect, the present application provides a computer storage medium, and a computer program, when executed by a processor, implements the method for synthesizing a circular truncated vessel mathematical model for hydrodynamics analysis described above.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the application provides a synthesis method of a circular truncated cone mathematical vascular model for hydrodynamics analysis, solves the problem that no vascular mathematical model for hydrodynamics analysis exists in the prior art, and fills the blank of the industry. Because the vascular wall has certain thickness, and the stenosis problem can appear mainly at the vascular inner wall, consequently this application will be through building into double-deck grid model with vascular mathematical model, has certain thickness, and inlayer grid model has elasticity, can play the corrective action to blood flow speed, and outer grid model can play the fixed effect of shape to inlayer grid model, combines mechanical analysis can effectually alleviate the deformation of vascular inner wall, more approaches the narrow condition of true blood vessel. Further, the minimum unit of the single-layer grid model is set to be a polygon with the number of sides being more than or equal to 6, because the deformation capacity of the triangle is poor, when one side is impacted by external force, the other side can also deform, the deformation of the triangle is larger, and the hexagon is impacted by external force, only two sides deform, and the other 4 sides can not deform, therefore, the deformation of the hexagon is smaller, and the double-layer grid model can form a hexagonal prism, the hexagonal prism is more stable relative to the triangular prism, and the hexagon has the advantages of less sampling points and higher sampling efficiency relative to the triangle, on the basis of keeping the original blood vessel shape, the calculation efficiency is improved when the fluid mechanics analysis CFD is calculated, and the calculation time is greatly shortened.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and not to limit the invention. In the drawings:
FIG. 1 is a schematic structural diagram of a circular truncated cone mathematical vascular model for hydrodynamics analysis according to the present application;
FIG. 2 is a flow chart of a method of synthesizing a frustoconical vascular mathematical model for hydrodynamic analysis according to the present application;
FIG. 3 is a schematic structural diagram of a single-layer mesh model of the present application;
FIG. 4 is a flowchart of S02 of the present application;
fig. 5 is a flowchart of S03 of the present application;
fig. 6 is a flowchart of S01 of the present application;
fig. 7 is a flowchart of S400 of the present application;
fig. 8 is a flowchart of S500 of the present application;
FIG. 9 is a flow chart of a first method of S510 of the present application;
fig. 10 is a flowchart of S520 of the present application;
FIG. 11 is a flow chart of a second method of S510 of the present application;
FIG. 12 is a flow chart of S530' of the present application;
fig. 13 is a flowchart of S600 of the present application;
fig. 14 is a flowchart of S700 of the present application;
fig. 15 is a flowchart of S730 of the present application;
fig. 16 is a flowchart of S900 of the present application;
FIG. 17 is a three-dimensional vessel model;
FIG. 18 is a block diagram of an apparatus for synthesizing a circular truncated vessel mathematical model;
FIG. 19 is a block diagram of a single-layer mesh model structure;
FIG. 20 is a block diagram of one embodiment of a three-dimensional vessel model structure 1 of the present application;
fig. 21 is another structural block diagram of an embodiment of the three-dimensional blood vessel model structure 1 of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the specific embodiments of the present invention and the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. 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.
In the following description, for purposes of explanation, numerous implementation details are set forth in order to provide a thorough understanding of the various embodiments of the present invention. It should be understood, however, that these implementation details are not to be interpreted as limiting the invention. That is, in some embodiments of the invention, such implementation details are not necessary. In addition, some conventional structures and components are shown in simplified schematic form in the drawings.
In coronary angiography images, coronary vessel assessment parameters need to be calculated in combination with fluid mechanics analysis, while no vessel mathematical model for fluid mechanics analysis exists in the prior art.
Example 1:
as shown in fig. 2, the present application provides a method for synthesizing a circular truncated cone mathematical vascular model for hydrodynamics analysis, comprising:
s01, according to the real-time diameter D of the blood vesseltCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model;
s02, carrying out N-edge type mesh division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer mesh model shown in figure 3, wherein N is more than or equal to 6;
s03, performing surface stratification on the single-layer mesh model to form a double-layer mesh model, i.e. the mathematical model of the blood vessel shown in fig. 1.
The application provides a synthesis method of a circular truncated cone blood vessel mathematical model for hydrodynamics analysis, solves the problem that no blood vessel three-dimensional grid model for hydrodynamics analysis exists in the prior art, and fills the blank of the industry. Because the vascular wall has certain thickness, and the stenosis problem can appear mainly at the vascular inner wall, consequently this application will be through building into double-deck grid model with vascular mathematical model, more is close to vascular true condition, and outer grid model can play the fixed effect of shape to inner grid model, combines mechanical analysis can effectually alleviate the deformation of vascular inner wall, more is close to the narrow condition of true blood vessel. Further, the minimum unit of the single-layer grid model is set to be a polygon with the number of sides being more than or equal to 6, because the deformation capacity of the triangle is poor, when one side is impacted by external force, the other side can also deform, the deformation of the triangle is larger, and the hexagon is impacted by external force, only two sides deform, and the other 4 sides can not deform, therefore, the deformation of the hexagon is smaller, and the double-layer grid model can form a hexagonal prism, the hexagonal prism is more stable relative to the triangular prism, and the hexagon has the advantages of less sampling points and higher sampling efficiency relative to the triangle, on the basis of keeping the original blood vessel shape, the calculation efficiency is improved when the fluid mechanics analysis CFD is calculated, and the calculation time is greatly shortened.
Example 2:
the application provides a method for synthesizing a circular truncated cone blood vessel mathematical model for fluid mechanics analysis, as shown in fig. 2, the method comprises the following steps:
s01, according to the real-time diameter D of the blood vesseltAnd the length L of the central line of the blood vessel is subjected to three-dimensional modeling to form a three-dimensional blood vessel model, as shown in figure 6, the three-dimensional model comprises the following steps:
s100, acquiring coronary artery two-dimensional contrast images of at least two body positions; preferably, the included angle between the two body positions is more than or equal to 30 degrees;
s200, obtaining a blood vessel section of interest from the coronary artery two-dimensional contrast image;
s300, picking up a starting point and an end point of the interested blood vessel section;
s400, segmenting a local blood vessel region map corresponding to the starting point and the ending point from the coronary artery two-dimensional contrast image, as shown in fig. 7, including:
s410, picking up at least one seed point of the blood vessel section of interest;
s420, segmenting the two-dimensional contrast image between two adjacent points of the starting point, the seed point and the ending point respectively to obtain at least two local blood vessel region images;
s500, obtaining the real-time diameter D of the blood vessel according to the coronary artery two-dimensional contrast imagetAnd the length L of the straightened blood vessel center line comprises:
extracting a vessel central line from the coronary artery two-dimensional contrast image of each posture along the direction from the coronary artery entrance to the coronary artery end, comprising two methods:
as shown in fig. 8, the first method is:
s510, respectively extracting at least one blood vessel path line from the local blood vessel region map of each posture, as shown in fig. 9, including:
s511, in each local blood vessel region image, the blood vessel section of interest is used as a foreground, other regions are used as backgrounds, the foreground is strengthened, the backgrounds are weakened, and a rough blood vessel image with strong contrast is obtained;
s512, performing mesh division on the rough blood vessel map;
s513, along the extending direction of the blood vessel from the starting point to the ending point, searching the shortest time path between the starting point and the intersection points on the peripheral n grids as a second point, searching the shortest time path between the second point and the intersection points on the peripheral n grids as a third point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
s514, connecting the extending directions of the blood vessels from the starting point to the end point according to the searching sequence to obtain at least one blood vessel route line;
s520, selecting a blood vessel path line as a blood vessel centerline, as shown in fig. 10, including:
s521, if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
s522, the least blood vessel route line is taken as the blood vessel central line.
As shown in fig. 11, the second method is:
s510', image processing is carried out on the local blood vessel region image, and a rough trend line of the blood vessel between the starting point and the ending point is obtained;
s520', obtaining rough blood vessel edge lines, wherein the image between the rough blood vessel edge lines including the rough blood vessel trend lines is the blood vessel skeleton;
s530', extracting a vessel centerline from the vascular skeleton, as shown in fig. 12, including:
s531', performing mesh division on the processed region image;
s532', searching the blood vessel skeleton according to the RGB values along the direction from the starting point to the ending point, searching a point where the minimum value of the RGB difference values of the starting point and the intersection points on the m grids around the starting point is located as a second point, searching a point where the minimum value of the RGB difference values of the second point and the intersection points on the m grids around the second point is located as a third point, and repeating the steps for the third point until the ending point is reached, wherein m is a positive integer greater than or equal to 1;
s533', obtaining at least one connecting line from the starting point to the end point according to the searching sequence;
s534', if the connecting line is two or more, selecting one connecting line as the center line of the blood vessel.
S600, acquiring a straightened vessel image according to the coronary artery two-dimensional angiography image and the vessel centerline, as shown in fig. 13, including:
s610, straightening the center line of the blood vessel to obtain a center straight line of the blood vessel;
s620, dividing the local blood vessel region map into x units along the extending direction of the blood vessel from the starting point to the ending point, wherein x is a positive integer;
s630, correspondingly arranging the center line of the blood vessel of each unit along the center straight line of the blood vessel;
s640, the correspondingly set image is a straightened blood vessel image;
s700, obtaining a blood vessel contour line after straightening according to the blood vessel centerline after straightening and the blood vessel image after straightening, as shown in fig. 14, including:
s710, setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
S720, according to DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel;
s730, gradually drawing the preset contour line of the blood vessel toward the central straight line of the blood vessel to obtain the straightened contour line of the blood vessel, as shown in fig. 15, including:
s731, dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
s732, acquiring z points of each unit, which are positioned on a preset contour line of each blood vessel;
s733, along the direction vertical to the blood vessel center straight line, z points are respectively closed to the blood vessel center straight line in a grading mode to generate z closed points, wherein z is a positive integer;
s734, setting the RGB difference threshold value as delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, the RGB value of the close point is matched with the blood vessel center straight line in each closingWhen the difference value is less than or equal to Δ RGBThreshold(s)When the blood vessel is closed, the closing point stops closing towards the center line of the blood vessel;
s735, acquiring a close point as a contour point;
s736, sequentially connecting the contour points to form a smooth curve which is the contour line of the blood vessel;
s800, acquiring geometric information of the straightened blood vessel, wherein the geometric information comprises the following steps: real-time diameter D of blood vesseltAnd the length of the straightened blood vessel central line, namely the length L of the central straight line, is as follows:
(1) real-time diameter D of blood vesseltAnd (2) the length L of the central straight line of the blood vessel;
(1) real-time diameter D of blood vesseltThe obtaining method comprises the following steps:
obtaining the distance between all the oppositely arranged contour points along the direction perpendicular to the center line of the blood vessel, namely the real-time diameter D of the blood vesselt。
S900, according to the DtAnd L three-dimensional modeling to form a three-dimensional model of the circular truncated cone, as shown in fig. 16, including:
s910, from the vessel real-time diameter DtObtaining the starting diameter D of the blood vessel internallyGet upAnd end diameter D of vesselPowderAnd a vessel center straight length L;
s920, according to DGet upAnd DPowderAnd L three-dimensional modeling, forming a three-dimensional model of the circular truncated cone as shown in FIG. 17;
s02, performing N-edge mesh partition along the circumferential surface of the three-dimensional blood vessel model to form a single-layer mesh model, where N is greater than or equal to 6, as shown in fig. 4, including:
s021 performing mesh division with a triangle as a minimum unit along a circumferential surface of the three-dimensional blood vessel model, including:
segmenting the three-dimensional vessel model into K segments,
performing mesh division on the circumferential surface of each segmented three-dimensional blood vessel model by taking a triangle as a minimum unit, preferably, the triangle as the minimum unit is an isosceles triangle;
s022, converting every N triangle combinations into 1N-sided shape according to the sequence to form an N-sided shape initial grid;
s023, deleting connecting lines inside each N-polygon in the N-polygon initial grids to form a single-layer N-polygon grid model, wherein N is more than or equal to 6;
s03, performing surface stratification process on the single-layer mesh model to form a double-layer mesh model, i.e. a mathematical model of blood vessels, as shown in fig. 5, including:
s031, obtaining the wall thickness h of the blood vessel; preferably, h is 0.2mm to 2 mm;
s032, according to the wall thickness h of the blood vessel, the starting diameter D of the blood vessel and the ending diameter D of the blood vesselPowderPerforming three-dimensional modeling on the length L of the central line of the blood vessel, and forming a circular truncated cone three-dimensional model on the inner surface or the outer surface of the single-layer grid model;
s033, according to the method for obtaining the single-layer mesh model, performing N-edge mesh division along the circumferential surface of the circular truncated cone three-dimensional model to form another single-layer mesh model;
s034, forming the double-layer mesh model, namely the blood vessel mathematical model, by two layers of the single-layer mesh model and the blood vessel wall thickness h.
Example 3:
as shown in fig. 18, the present application provides an apparatus for synthesizing a mathematical model of a blood vessel, comprising: the three-dimensional model structure comprises a three-dimensional blood vessel model structure 1, a single-layer grid model structure 2 and a circular truncated cone mathematical model structure 3 which are sequentially connected, wherein the blood vessel mathematical model structure 3 is connected with the three-dimensional model structure 1; the three-dimensional blood vessel model structure 1 is used for measuring the real-time diameter D of the blood vesseltCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model; the single-layer grid model structure 2 is used for carrying out N-edge type grid division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer grid model, wherein N is more than or equal to 6; the vessel mathematical model structure 3 is used for performing surface layering processing on the single-layer mesh model to form a double-layer mesh model, namely a vessel mathematical model.
As shown in fig. 20, the three-dimensional blood vessel model structure 1 further includes: the device comprises a central line extracting unit 100, a straightening unit 200, a contour line unit 300, a geometric information unit 400 and a three-dimensional modeling unit 500 which are connected in sequence; the straightening unit 200 is connected with the geometric information unit 400, and the three-dimensional modeling unit 500 is connected with the straightening unit 200 and the contour line unit 300; the centerline extraction unit 100 is configured to extract a blood vessel centerline from the two-dimensional coronary angiography images of at least two body positions along a direction from the coronary artery entrance to the coronary artery end; the straightening unit 200 is configured to receive the blood vessel center line sent by the center line extracting unit 100, and obtain a straightened blood vessel image according to the coronary artery two-dimensional angiography image and the blood vessel center line; the contour line unit 300 is configured to receive the straightened blood vessel image sent by the straightening unit 200, and obtain a straightened blood vessel contour line according to the straightened blood vessel center line and the straightened blood vessel image; the geometric information unit 400 is configured to receive the straightened blood vessel image sent by the straightening unit 200 and the blood vessel contour line sent by the contour line unit 300, and obtain geometric information of the straightened blood vessel; the three-dimensional modeling unit 500 is configured to receive the straightened blood vessel image sent by the straightening unit 200, the blood vessel contour line sent by the contour line unit 300, and the geometric information of the blood vessel sent by the geometric information unit 400, perform three-dimensional modeling according to the geometric information, the center line, and the contour line, and obtain a three-dimensional blood vessel model.
As shown in fig. 21, in an embodiment of the present application, the apparatus further includes: an image segmentation unit 600 connected to the center line extraction unit 100; the image segmentation unit 600 is configured to segment a local blood vessel region map corresponding to a start point and an end point from a two-dimensional coronary angiography image, or segment a two-dimensional angiography image between two adjacent points of the start point, the seed point, and the end point to obtain at least two local blood vessel region maps.
As shown in fig. 21, in an embodiment of the present application, the centerline extraction unit 100 further includes: the blood vessel path module 110 and the blood vessel center line extraction module 120 are connected in sequence, and the blood vessel path module 110 is connected with the image segmentation unit 600; a vessel path module 110, configured to extract at least one vessel path line from the local vessel region map of each body position; a blood vessel centerline extraction module 120, configured to select one of the blood vessel path lines sent by the blood vessel path module 110 as a blood vessel centerline.
As shown in fig. 19, in an embodiment of the present application, the single-layer mesh model structure 2 further includes: the device comprises a triangular mesh dividing unit 21, an N-edge mesh dividing unit 22 and a single-layer mesh model unit 23 which are connected in sequence, wherein the single-layer mesh model unit 23 is connected with a blood vessel mathematical model structure 3; the triangular mesh division unit 21 is connected with the three-dimensional modeling unit 500, and is used for carrying out mesh division by taking a triangle as a minimum unit along the circumferential surface of the three-dimensional blood vessel model; the N-polygon mesh dividing unit 22 is configured to convert every N triangle combinations into 1N polygons in sequence to form N-polygon initial meshes; the single-layer grid model unit 23 is configured to delete a connection line inside each N-polygon in the N-polygon initial grid to form a single-layer N-polygon grid model, where N is greater than or equal to 6.
The present application provides a coronary artery analysis system comprising: the device for synthesizing the blood vessel mathematical model is described above.
The present application provides a computer storage medium, and a computer program, when executed by a processor, implements the method for synthesizing a circular truncated cone mathematical model for fluid mechanics analysis.
As will be appreciated by one skilled in the art, aspects of the present invention may be embodied as a system, method or computer program product. Thus, various aspects of the invention may be embodied in the form of: an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining hardware and software aspects that may all generally be referred to herein as a "circuit," module "or" system. Furthermore, in some embodiments, aspects of the invention may also be embodied in the form of a computer program product in one or more computer-readable media having computer-readable program code embodied therein. Implementation of the method and/or system of embodiments of the present invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
For example, hardware for performing selected tasks according to embodiments of the invention could be implemented as a chip or a circuit. As software, selected tasks according to embodiments of the invention could be implemented as a plurality of software instructions being executed by a computer using any suitable operating system. In an exemplary embodiment of the invention, one or more tasks according to exemplary embodiments of the methods and/or systems as described herein are performed by a data processor, such as a computing platform for executing a plurality of instructions. Optionally, the data processor comprises volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, e.g. a magnetic hard disk and/or a removable medium. Optionally, a network connection is also provided. A display and/or a user input device, such as a keyboard or mouse, is optionally also provided.
Any combination of one or more computer readable media may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following:
an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
For example, computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer-readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.
The computer program instructions may also be loaded onto a computer (e.g., a coronary artery analysis system) or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer, other programmable apparatus or other devices provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
The above embodiments of the present invention have been described in further detail for the purpose of illustrating the invention, and it should be understood that the above embodiments are only illustrative of the present invention and are not to be construed as limiting the present 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 (19)
1. The synthetic method of the circular truncated cone blood vessel mathematical model for the fluid mechanics analysis is characterized by comprising the following steps:
according to the real-time diameter D of the vesseltCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model;
carrying out N-edge type grid division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer grid model, wherein N is more than or equal to 6;
and carrying out surface layering treatment on the single-layer grid model to form a double-layer grid model, namely a circular truncated cone blood vessel mathematical model.
2. The method for synthesizing a circular truncated vessel mathematical model for hydrodynamics analysis as claimed in claim 1, wherein said N-edge mesh partition along the circumferential surface of said three-dimensional vessel model forms a single-layer mesh model, wherein the method of N ≧ 6 includes:
carrying out mesh division by taking a triangle as a minimum unit along the circumferential surface of the three-dimensional blood vessel model;
according to the sequence, every N triangles are combined and converted into 1N-sided polygon, and an N-sided polygon initial grid is formed;
and deleting the connecting lines inside each N-polygon in the N-polygon initial grids to form a single-layer N-polygon grid model, wherein N is more than or equal to 6.
3. The method of synthesizing a frustoconical vascular mathematical model for hydrodynamic analysis of claim 2, wherein the method of meshing with a triangle as a minimum unit along a circumferential surface of the three-dimensional vascular model comprises:
segmenting the three-dimensional vessel model into K segments,
and performing mesh division on the circumferential surface of each section of the three-dimensional blood vessel model by taking a triangle as a minimum unit.
4. The method of synthesizing a frustoconical vascular mathematical model for hydrodynamic analysis of claim 2, wherein the triangle as the smallest unit is an isosceles triangle.
5. The method of synthesizing a frustoconical mathematical vascular model for fluid mechanics analysis of claim 1, wherein the method of surface-layering the single-layer mesh model to form a two-layer mesh model, i.e., a mathematical vascular model, comprises:
obtaining the wall thickness h of the blood vessel;
according to the wall thickness h and the initial diameter D of the blood vesselGet upEnd of vessel diameter DPowderPerforming three-dimensional modeling on the length L of the central line of the blood vessel, and forming a circular truncated cone three-dimensional model on the inner surface or the outer surface of the single-layer grid model;
according to the acquisition method of the single-layer grid model, carrying out N-edge type grid division along the circumferential surface of the circular truncated cone three-dimensional model to form another single-layer grid model;
and forming the double-layer grid model, namely the circular truncated cone blood vessel mathematical model, by the two layers of the single-layer grid model and the blood vessel wall thickness h.
6. The method for synthesizing a frustoconical vessel mathematical model for fluid mechanics analysis of claim 1, wherein the real-time vessel diameter D is measured in terms of vessel diametertThe method for forming the three-dimensional blood vessel model by three-dimensional modeling of the length L of the central line of the blood vessel comprises the following steps:
acquiring coronary artery two-dimensional contrast images of at least two body positions;
obtaining the real-time diameter D of the blood vessel according to the coronary artery two-dimensional contrast imagetAnd the length L of the straightened blood vessel center line;
according to said DtAnd L three-dimensional modeling to form another round table three-dimensional model, namely a three-dimensional blood vessel model.
7. The method for synthesizing the circular truncated vessel mathematical model for fluid mechanics analysis of claim 6, wherein the obtaining of the real-time vessel diameter D from the coronary artery two-dimensional contrast imagetAnd the method for the length L of the straightened blood vessel central line comprises the following steps:
extracting a blood vessel central line from the coronary artery two-dimensional contrast image of each posture along the direction from a coronary artery entrance to the coronary artery tail end;
acquiring a straightened vessel image according to the coronary artery two-dimensional radiography image and the vessel central line;
acquiring a straightened blood vessel contour line according to the straightened blood vessel center line and the straightened blood vessel image;
acquiring geometric information of the straightened vessel, comprising: real-time diameter D of blood vesseltAnd the length of the straightened blood vessel central line is the length L of the central straight line.
8. The method of synthesizing a frustoconical vascular mathematical model for hydrodynamic analysis of claim 6, wherein the method is based on the DtAnd L three-dimensional modeling, wherein the method for forming the circular truncated cone three-dimensional model comprises the following steps:
performing three-dimensional modeling according to the geometric information, the central line and the contour line to obtain a three-dimensional blood vessel model;
real-time diameter D from the vesseltObtaining the starting diameter D of the blood vessel internallyGet upAnd end diameter D of vesselPowder；
According to said DGet up、DPowderAnd L, performing three-dimensional modeling to form the three-dimensional model of the circular truncated cone.
9. The method of claim 6, wherein after said acquiring the coronary artery two-dimensional angiographic images of at least two body positions, obtaining a real-time vessel diameter D from said coronary artery two-dimensional angiographic imagestObtaining the starting diameter D of the blood vessel internallyGet upAnd end diameter D of vesselPowderAnd the length L after the blood vessel central line is straightened further comprises:
acquiring a vessel segment of interest from the coronary artery two-dimensional angiogram;
picking up a starting point and an end point of the vessel segment of interest;
and segmenting local vessel region maps corresponding to the starting points and the end points from the coronary artery two-dimensional contrast image.
10. The method for synthesizing a circular truncated vessel mathematical model for fluid mechanics analysis of claim 9, wherein the method for segmenting the local vessel region map corresponding to the starting point and the ending point from the coronary artery two-dimensional contrast image further comprises:
picking up at least one seed point of the vessel segment of interest;
and respectively segmenting the two-dimensional contrast images between the adjacent two points of the starting point, the seed point and the ending point to obtain at least two local blood vessel region images.
11. The method for synthesizing a circular truncated vessel mathematical model for fluid mechanics analysis of claim 7 wherein said method for extracting a vessel centerline from said coronary artery two-dimensional angiographic image of each body position along the direction from coronary artery entrance to coronary artery end comprises:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
performing mesh division on the rough blood vessel map, and extracting at least one blood vessel path line along the direction from the starting point to the end point;
and selecting one blood vessel path line as the blood vessel central line.
12. The method of synthesizing a frustoconical mathematical vascular model for fluid dynamics analysis as claimed in claim 11, wherein the step of meshing the rough vessel map and extracting at least one vessel path line along the direction from the starting point to the ending point comprises:
gridding the rough vessel map;
searching the shortest time path between the starting point and the intersection points on the peripheral n grids along the extending direction of the blood vessels from the starting point to the ending point to serve as a second point, searching the shortest time path between the second point and the intersection points on the peripheral n grids to serve as a third point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and connecting the extending directions of the blood vessels from the starting point to the ending point according to the searching sequence to obtain at least one blood vessel path line.
13. The method of synthesizing a frustoconical vessel mathematical model for fluid mechanics analysis of claim 12, wherein said method of selecting one of said vessel path lines as said vessel centerline comprises:
if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
the vessel path line that is the least in time is taken as the vessel centerline.
14. The method for synthesizing a circular truncated vessel mathematical model for fluid mechanics analysis of claim 7 wherein said method for extracting a vessel centerline from a coronary artery two-dimensional angiographic image of each body position along the direction from coronary artery entrance to coronary artery end comprises:
performing image processing on the local blood vessel region image to obtain a rough trend line of the blood vessel between the starting point and the ending point;
obtaining rough blood vessel edge lines, wherein images between the rough blood vessel edge lines including the rough blood vessel trend lines are blood vessel skeletons;
extracting the vessel centerline from the vascular skeleton.
15. The method of synthesizing a frustoconical vessel mathematical model for hydrodynamic analysis of claim 14, wherein the method of extracting the vessel centerline from the vascular skeleton comprises:
performing grid division on the processed region image;
searching the blood vessel skeleton according to RGB values along the direction from the starting point to the ending point, searching a point where the minimum value of RGB difference values of the starting point and intersection points on m grids around the starting point is located as a second point, searching a point where the minimum value of RGB difference values of the second point and intersection points on m grids around the second point is located as a third point, and repeating the steps until the third point reaches the ending point, wherein m is a positive integer greater than or equal to 1;
obtaining at least one connecting line from the starting point to the end point according to the searching sequence;
and if the connecting line is two or more, selecting one connecting line as the center line of the blood vessel.
16. The method of synthesizing a frustoconical vessel mathematical model for fluid dynamics analysis of claim 15, wherein the method of obtaining a straightened vessel contour from the straightened vessel centerline and the straightened vessel image comprises:
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line;
and gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the straightened contour line of the blood vessel.
17. The device for synthesizing the circular truncated cone-shaped blood vessel mathematical model is used for the method for synthesizing the circular truncated cone-shaped blood vessel mathematical model for hydrodynamic analysis according to any one of claims 1 to 16, and is characterized by comprising the following steps: the three-dimensional model structure of the blood vessel, the structure of the single-layer grid model and the structure of the mathematical model of the blood vessel are connected in sequence, and the structure of the mathematical model of the blood vessel is connected with the structure of the three-dimensional model;
the three-dimensional blood vessel model structure is used for real-time diameter D according to blood vesselstCarrying out three-dimensional modeling on the length L of the central line of the blood vessel to form a three-dimensional blood vessel model;
the single-layer grid model structure is used for carrying out N-edge type grid division along the circumferential surface of the three-dimensional blood vessel model to form a single-layer grid model, wherein N is more than or equal to 6;
the vessel mathematical model structure is used for carrying out surface layering treatment on the single-layer grid model to form a double-layer grid model, namely the vessel mathematical model.
18. A coronary artery analysis system, comprising: the apparatus for synthesizing a frustoconical vascular mathematical model of claim 17.
19. A computer storage medium, wherein a computer program when executed by a processor implements the method of synthesizing a circular truncated vessel mathematical model for hydrodynamic analysis of any one of claims 1 to 16.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911237196.5A CN112116711B (en) | 2019-12-05 | 2019-12-05 | Synthesizing method and device for truncated cone blood vessel mathematical model for hydrodynamic analysis |
| PCT/CN2019/123615 WO2021109116A1 (en) | 2019-12-05 | 2019-12-06 | Synthesis method of circular truncated blood vessel mathematical model for hydromechanics analysis and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911237196.5A CN112116711B (en) | 2019-12-05 | 2019-12-05 | Synthesizing method and device for truncated cone blood vessel mathematical model for hydrodynamic analysis |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112116711A true CN112116711A (en) | 2020-12-22 |
| CN112116711B CN112116711B (en) | 2024-01-23 |
Family
ID=73795370
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911237196.5A Active CN112116711B (en) | 2019-12-05 | 2019-12-05 | Synthesizing method and device for truncated cone blood vessel mathematical model for hydrodynamic analysis |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112116711B (en) |
| WO (1) | WO2021109116A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115422856A (en) * | 2022-09-05 | 2022-12-02 | 青岛埃米博创医疗科技有限公司 | CFD teaching-oriented teaching blood vessel model generation method |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102010009701A1 (en) * | 2010-03-01 | 2011-09-01 | Siemens Aktiengesellschaft | Method for automatically identifying narrow points of coronary blood vessel in heart of patient, involves identifying narrow points in hollow vessel based on automated evaluation of profile course curve of hollow vessel |
| US20120203530A1 (en) * | 2011-02-07 | 2012-08-09 | Siemens Corporation | Method and System for Patient-Specific Computational Modeling and Simulation for Coupled Hemodynamic Analysis of Cerebral Vessels |
| CN105096388A (en) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | Computational Fluid Dynamics (CFD) based coronary artery blood flow simulating system and method |
| CN106327487A (en) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | Coronary artery blood flow reserve fraction calculation method based on X ray coronary artery angiographic image |
| CN108550180A (en) * | 2018-03-09 | 2018-09-18 | 南京信息工程大学 | Blood vessel modeling method based on interior point set region constraint and Gaussian process parameter optimization |
| CN109461138A (en) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | Calculation method of parameters, system, readable storage medium storing program for executing and computer equipment |
| CN110367965A (en) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | The method, apparatus and system of convenient measurement coronary artery assessment parameters |
| CN110384494A (en) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | The method for measuring microcirculation drag index |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8170304B2 (en) * | 2007-04-03 | 2012-05-01 | Siemens Aktiengesellschaft | Modeling cerebral aneurysms in medical images |
| CN105095615B (en) * | 2014-04-23 | 2019-05-21 | 北京冠生云医疗技术有限公司 | The method and system that blood flow data in blood vessel is handled |
-
2019
- 2019-12-05 CN CN201911237196.5A patent/CN112116711B/en active Active
- 2019-12-06 WO PCT/CN2019/123615 patent/WO2021109116A1/en active Application Filing
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102010009701A1 (en) * | 2010-03-01 | 2011-09-01 | Siemens Aktiengesellschaft | Method for automatically identifying narrow points of coronary blood vessel in heart of patient, involves identifying narrow points in hollow vessel based on automated evaluation of profile course curve of hollow vessel |
| US20120203530A1 (en) * | 2011-02-07 | 2012-08-09 | Siemens Corporation | Method and System for Patient-Specific Computational Modeling and Simulation for Coupled Hemodynamic Analysis of Cerebral Vessels |
| CN105096388A (en) * | 2014-04-23 | 2015-11-25 | 北京冠生云医疗技术有限公司 | Computational Fluid Dynamics (CFD) based coronary artery blood flow simulating system and method |
| CN106327487A (en) * | 2016-08-18 | 2017-01-11 | 苏州润心医疗科技有限公司 | Coronary artery blood flow reserve fraction calculation method based on X ray coronary artery angiographic image |
| CN108550180A (en) * | 2018-03-09 | 2018-09-18 | 南京信息工程大学 | Blood vessel modeling method based on interior point set region constraint and Gaussian process parameter optimization |
| CN110367965A (en) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | The method, apparatus and system of convenient measurement coronary artery assessment parameters |
| CN110384494A (en) * | 2018-09-19 | 2019-10-29 | 苏州润迈德医疗科技有限公司 | The method for measuring microcirculation drag index |
| CN109461138A (en) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | Calculation method of parameters, system, readable storage medium storing program for executing and computer equipment |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115422856A (en) * | 2022-09-05 | 2022-12-02 | 青岛埃米博创医疗科技有限公司 | CFD teaching-oriented teaching blood vessel model generation method |
| CN115422856B (en) * | 2022-09-05 | 2023-08-08 | 青岛埃米博创医疗科技有限公司 | CFD teaching-oriented teaching blood vessel model generation method |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2021109116A1 (en) | 2021-06-10 |
| CN112116711B (en) | 2024-01-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210272030A1 (en) | Systems and methods for using geometry sensitivity information for guiding workflow | |
| CN112132949B (en) | Three-dimensional modeling method, device and system for blood vessel with narrow lesion region | |
| WO2019210553A1 (en) | Microcirculation resistance index calculation method based on angiogram image and hydrodynamics model | |
| CN111612743B (en) | CT image-based coronary artery central line extraction method | |
| CN112116615B (en) | Method and device for acquiring blood vessel contour line according to blood vessel center line | |
| CN110889896A (en) | Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis | |
| US20220319116A1 (en) | Method and device for synthesizing mathematical model of blood vessel having stenotic lesion | |
| CN112164020A (en) | Method, device, analysis system and storage medium for accurately extracting blood vessel center line | |
| CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
| CN110786842B (en) | Method, device, system and storage medium for measuring diastolic blood flow velocity | |
| WO2022109906A1 (en) | Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium | |
| CN116090364A (en) | Method for obtaining coronary blood flow reserve fraction based on CTA image and readable storage medium | |
| CN110786840B (en) | Method, device and storage medium for obtaining blood vessel evaluation parameter based on physiological parameter | |
| CN112132882A (en) | Method and device for extracting blood vessel central line from coronary artery two-dimensional contrast image | |
| CN112419276A (en) | Method for regulating blood vessel contour and central line and storage medium | |
| US10354349B2 (en) | Systems and methods for using geometry sensitivity information for guiding workflow | |
| CN112116711B (en) | Synthesizing method and device for truncated cone blood vessel mathematical model for hydrodynamic analysis | |
| KR101514795B1 (en) | Quantification method of vessel | |
| CN112419280B (en) | Method for precisely acquiring narrow lesion interval and storage medium | |
| CN112419277A (en) | Three-dimensional blood vessel center line synthesis method, system and storage medium | |
| CN112419279A (en) | Method and storage medium for two-dimensional image selection and three-dimensional blood vessel synthesis | |
| CN115222665B (en) | Plaque detection method and device, electronic equipment and readable storage medium | |
| CN116342608A (en) | Medical image-based stent adherence measurement method, device, equipment and medium | |
| Shen et al. | Evaluation and analysis of cardiovascular function in intensive care unit patients by ultrasound image segmentation based on deep learning | |
| CN115035064A (en) | Myocardial bridge depth determination system, myocardial bridge depth determination device and storage medium |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |