CN112419484A - Three-dimensional blood vessel synthesis method and system, coronary artery analysis system and storage medium - Google Patents
Three-dimensional blood vessel synthesis method and system, coronary artery analysis system and storage medium Download PDFInfo
- Publication number
- CN112419484A CN112419484A CN202011337609.XA CN202011337609A CN112419484A CN 112419484 A CN112419484 A CN 112419484A CN 202011337609 A CN202011337609 A CN 202011337609A CN 112419484 A CN112419484 A CN 112419484A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- dimensional
- line
- image
- 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 414
- 210000004351 coronary vessel Anatomy 0.000 title claims abstract description 66
- 238000004458 analytical method Methods 0.000 title claims abstract description 9
- 238000001308 synthesis method Methods 0.000 title abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 85
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 62
- 238000002601 radiography Methods 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 11
- 238000003786 synthesis reaction Methods 0.000 claims description 11
- 230000002285 radioactive effect Effects 0.000 claims description 10
- 241000282414 Homo sapiens Species 0.000 claims description 8
- 238000002583 angiography Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000005855 radiation Effects 0.000 claims description 8
- 238000002586 coronary angiography Methods 0.000 claims description 7
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 2
- 230000002792 vascular Effects 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 8
- 230000008569 process Effects 0.000 abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 6
- 238000010009 beating Methods 0.000 abstract description 3
- 238000006073 displacement reaction Methods 0.000 abstract description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 201000000057 Coronary Stenosis Diseases 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 206010011089 Coronary artery stenosis Diseases 0.000 description 1
- 206010020565 Hyperaemia Diseases 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
- 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
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 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
-
- 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/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
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Epidemiology (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application provides a three-dimensional blood vessel synthesis method and system, a coronary artery analysis system and a storage medium. Acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles; acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to the image information of the coronary artery two-dimensional radiography image; and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel central line and the three-dimensional blood vessel radius. The method can effectively simulate the state of the blood vessel in a real scene, including the shape, the trend and the diameter information of the blood vessel, solves the influence of external conditions on the image in the operation process, including the displacement of equipment, the beating of the heart and the respiration of a patient, solves the problems of high risk and high cost of guide wire measurement, and provides a basis for calculating blood vessel evaluation parameters such as the Fractional Flow Reserve (FFR).
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a three-dimensional blood vessel synthesis method and system, a coronary artery analysis system and a storage medium.
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.
The Fractional Flow Reserve (FFR) generally refers to the fractional flow reserve of myocardium, and is defined as the ratio of the maximum blood flow provided by a diseased coronary artery to the maximum blood flow when the coronary artery is completely normal. Namely, the FFR value can be measured and calculated by measuring the pressure at the position of the coronary stenosis and the pressure at the position of the coronary stenosis under the maximal hyperemia state of the coronary artery through a pressure sensor.
The problems that exist are that: the mode of obtaining blood vessel evaluation parameters such as Fractional Flow Reserve (FFR) and the like through the pressure sensor needs a doctor to draw the sensor to a lesion part through a guide wire from an artery, and in the process, a plurality of difficulties need to be solved, such as invasiveness of operation, complexity of measurement and high cost of the pressure guide wire, can become obstacles for popularization of the FFR.
Disclosure of Invention
The invention provides a three-dimensional blood vessel synthesis method and system, a coronary artery analysis system and a storage medium, which are used for synthesizing a three-dimensional blood vessel through image simulation without adopting pressure guide wire measurement, and solve the use problem of the pressure guide wire in the prior art.
To achieve the above object, in a first aspect, a method for synthesizing a three-dimensional blood vessel includes:
acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles;
acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to the image information of the coronary artery two-dimensional radiography image;
and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel central line and the three-dimensional blood vessel radius.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for acquiring image information of at least two-dimensional coronary angiography images with different capturing angles includes:
acquiring at least two groups of coronary artery two-dimensional contrast image groups with different shooting angles;
reading image information of each group of coronary artery two-dimensional contrast image group, wherein the image information comprises a shooting angle and a detection distance;
and respectively selecting an interested two-dimensional contrast image from each group of the coronary artery two-dimensional contrast images according to the detection distance.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for acquiring a three-dimensional blood vessel centerline from the coronary artery two-dimensional contrast image includes:
extracting a two-dimensional vessel centerline from each of the two-dimensional angiographic images of interest;
and projecting each two-dimensional blood vessel central line into a three-dimensional space according to the shooting angle and the detection distance of each two-dimensional coronary artery angiography image, and synthesizing the three-dimensional blood vessel central lines.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for projecting each two-dimensional blood vessel centerline into a three-dimensional space according to the capturing angle and the detection distance of each two-dimensional coronary angiography image includes:
establishing a three-dimensional coordinate system by taking the heart as a coordinate origin;
acquiring a left-right angle alpha, a front-back angle beta and a distance S between a human body and a flat panel detector of each interested two-dimensional radiography image, wherein the coordinates of each point on the center line of the two-dimensional blood vessel are (x, y);
projecting each point (x, y) into a three-dimensional space to obtain a series of three-dimensional coordinate points P, wherein the coordinates are (x ', y ', z ');
projecting a radioactive source into the three-dimensional space to form a radioactive point R;
and connecting each three-dimensional coordinate point P with the radiation point R, obtaining points on the central line of the three-dimensional blood vessel from the PR connection line, and sequentially connecting the points on the central line of the three-dimensional blood vessel to obtain the central line of the three-dimensional blood vessel.
Optionally, in the above method for synthesizing a three-dimensional blood vessel, the method for projecting each point (x, y) into a three-dimensional space to obtain a series of three-dimensional coordinate points P with coordinates (x ", y", z ") includes:
rotating each point (x, y) around the y axis to obtain a series of (x ', y ', z ') points, wherein the specific formula is as follows:
rotating the series of points (x ', y ', z ') around an x axis to obtain a series of three-dimensional coordinate points P with coordinates (x ", y", z ");
optionally, the method for synthesizing a three-dimensional blood vessel as described above, the method for projecting a radiation source into the three-dimensional space to form a radiation point R, includes:
obtaining a distance S' between a human body and a radioactive source in each two-dimensional contrast image of interest;
according to the formula
Wherein, the coordinate of R in the three-dimensional space is (a, b, c).
Optionally, in the method for synthesizing a three-dimensional blood vessel, the connecting line between each three-dimensional coordinate point P and the radiation point R is used to obtain a point on a centerline of the three-dimensional blood vessel from a PR connecting line, and the connecting lines are sequentially connected to obtain the centerline of the three-dimensional blood vessel, and the method includes:
correspondingly connecting a series of three-dimensional coordinate points P and radial points R obtained from the same two-dimensional contrast image of interest to obtain a plurality of PR straight lines;
acquiring points of the minimum distance between two PR straight lines at the same position of the blood vessel, wherein the points are respectively a point A and a point B;
connecting the point A with the point B, and acquiring the midpoint of the line segment AB as a point on the central line of the three-dimensional blood vessel;
and sequentially connecting the obtained points on the central line of the series of three-dimensional blood vessels to obtain the central line of the three-dimensional blood vessel.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for acquiring a blood vessel centerline and a three-dimensional blood vessel radius from the coronary artery two-dimensional contrast image includes:
acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line;
acquiring the two-dimensional vessel radius in each two-dimensional contrast image of interest according to the two-dimensional vessel contour line;
and obtaining the three-dimensional vessel radius according to the two-dimensional vessel radius.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for obtaining the three-dimensional blood vessel radius according to the two-dimensional blood vessel radius includes:
wherein R represents the three-dimensional vessel radius, R1、r2、rnThe two-dimensional vessel radii of the first, second and nth two-dimensional contrast images of interest are represented, respectively.
Optionally, in the method for synthesizing three-dimensional blood vessels, the method for extracting a two-dimensional blood vessel centerline from each of the two-dimensional coronary angiography images includes:
reading a coronary artery two-dimensional contrast image;
obtaining a vessel segment of interest;
picking up a starting point, a seed point and an end point of the vessel segment of interest;
respectively segmenting two-dimensional contrast images between two adjacent points of a starting point, a seed point and an end point to obtain at least two local blood vessel region images;
extracting at least one blood vessel local path line from each local blood vessel region map;
connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
and selecting one blood vessel path line as the two-dimensional blood vessel central line.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for extracting at least one blood vessel local path line from each local blood vessel region map includes:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
and performing grid division on the rough blood vessel map, and extracting at least one blood vessel local path line along the direction from the starting point to the end point.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for performing image enhancement processing on the local blood vessel region map to obtain a rough blood vessel map with a strong contrast includes:
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 the rough blood vessel image with strong contrast is obtained.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for meshing the rough blood vessel map and extracting at least one blood vessel local 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 local path line.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the selecting one of the blood vessel path lines as the two-dimensional blood vessel center line 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 two-dimensional vessel centerline.
Optionally, the method for synthesizing a three-dimensional blood vessel described above, wherein the method for obtaining a two-dimensional blood vessel contour line from the blood vessel centerline includes:
extracting a two-dimensional blood vessel central line according to the coronary artery two-dimensional radiography image;
obtaining a straightened blood vessel image according to the two-dimensional blood vessel central line;
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel;
gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and projecting the contour line of the straightened blood vessel back to the image for extracting the center line of the two-dimensional blood vessel to obtain the contour line of the two-dimensional blood vessel.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line includes:
straightening the center line of the two-dimensional blood vessel to obtain a blood vessel center straight line;
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;
correspondingly arranging the two-dimensional blood vessel central line of each unit along the blood vessel central straight line;
and the correspondingly set image is the straightened blood vessel image.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the method for gradually drawing the preset contour line of the blood vessel toward the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel includes:
dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on each preset blood vessel contour line;
respectively closing the z points to the blood vessel center straight line in a grading way along the direction vertical to the blood vessel center straight line to generate z closing points, wherein z is a positive integer;
setting RGB difference threshold to delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)Then the closing point stops moving to the center of the blood vesselThe straight lines are closed;
acquiring the approach point as a contour point;
and sequentially connecting the contour points to form a smooth curve which is the contour line of the straightened blood vessel.
Optionally, the method for synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel centerline and the three-dimensional blood vessel radius includes:
drawing a picture in the three-dimensional space along the corresponding three-dimensional blood vessel radius to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon approximate to a circle;
and sequentially connecting points on two adjacent polygons according to a right-angle triangle form to obtain the three-dimensional blood vessel.
In a second aspect, the present application provides a three-dimensional vessel synthesis system comprising: the three-dimensional blood vessel radius acquisition device is connected with the image reading device and the three-dimensional blood vessel center line acquisition device;
the image reading device is used for acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles;
the three-dimensional blood vessel center line acquisition device is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device and acquiring a three-dimensional blood vessel center line according to the image information;
the three-dimensional blood vessel radius acquisition device is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device, receiving the three-dimensional blood vessel center line transmitted by the three-dimensional blood vessel center line acquisition device, and acquiring a three-dimensional blood vessel radius according to the image information and the three-dimensional blood vessel center line;
the three-dimensional blood vessel synthesizing device is used for receiving the three-dimensional blood vessel center line transmitted by the three-dimensional blood vessel center line acquiring device and receiving the three-dimensional blood vessel radius transmitted by the three-dimensional blood vessel radius acquiring device, and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius.
Optionally, in the three-dimensional vessel synthesis system, the three-dimensional vessel centerline obtaining device includes: the two-dimensional blood vessel center line extracting structure is connected with the image reading device, and the three-dimensional blood vessel center line acquiring structure is connected with the two-dimensional blood vessel center line extracting structure;
the two-dimensional vessel centerline extraction structure is used for receiving the coronary artery two-dimensional contrast images sent by the image reading device and extracting a two-dimensional vessel centerline from each two-dimensional contrast image of interest;
the three-dimensional blood vessel center line obtaining structure is used for receiving the two-dimensional blood vessel center line sent by the two-dimensional blood vessel center line extracting structure, receiving the two-dimensional blood vessel center line sent by the image reading device, projecting each two-dimensional blood vessel center line into a three-dimensional space according to the shooting angle of each two-dimensional coronary artery angiography image, and synthesizing the three-dimensional blood vessel center lines.
Optionally, in the method for synthesizing a three-dimensional blood vessel, the two-dimensional blood vessel centerline extraction structure includes: the center line extracting unit, the straightening unit, the first blood vessel contour line unit and the second blood vessel contour line unit are sequentially connected;
the central line extracting unit is connected with the image reading device and used for extracting a blood vessel central line according to a coronary artery two-dimensional contrast image;
the straightening unit is used for obtaining a straightened blood vessel image according to the blood vessel central line extracted by the central line extracting unit;
the first blood vessel contour line unit is used for setting a blood vessel diameter threshold value D on the straightened blood vessel image sent by the straightening unitThreshold(s)(ii) a According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line; gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and the second blood vessel contour line unit is used for projecting the contour line of the straightened blood vessel sent by the first blood vessel contour line unit back to the image of the blood vessel central line to obtain the blood vessel contour line.
In a third aspect, the present application provides a coronary artery analysis system comprising: the three-dimensional vessel synthesis system.
In a fourth aspect, the present application provides a computer storage medium, and a computer program, when executed by a processor, implements the above-described three-dimensional blood vessel synthesis method.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the application provides a three-dimensional blood vessel synthesis method, which can effectively simulate the state of blood vessels in a real scene, including the shape, the trend and the diameter information of the blood vessels, solves the influence of external conditions on images in the operation process, including the displacement of equipment, the beating of the heart and the respiration of a patient, solves the problems of high risk and high cost of guide wire measurement, and provides a basis for calculating blood vessel evaluation parameters such as blood flow reserve fraction (FFR).
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:
the following reference numerals are used for the description:
FIG. 1 is a flow chart of a method of synthesizing a three-dimensional blood vessel according to the present application;
fig. 2 is a flowchart of S100 of the present application;
fig. 3 is a flowchart of S200 of the present application;
fig. 4 is a flowchart of S210 of the present application;
FIG. 5 is a flowchart of S215 of the present application;
fig. 6 is a flowchart of S2152 of the present application;
FIG. 7 is a flowchart of S217 of the present application;
FIG. 8 is a flowchart of S220 of the present application;
FIG. 9 is a flowchart of S225 of the present application;
FIG. 10 is a flowchart of S230 of the present application;
fig. 11 is a flowchart of S232 of the present application;
fig. 12 is a flowchart of S235 of the present application;
fig. 13 is a flowchart of S300 of the present application;
FIG. 14 is a block diagram of the structure of the three-dimensional vessel synthesis system of the present application;
FIG. 15 is another block diagram of the three-dimensional vessel synthesis system of the present application;
fig. 16 is a block diagram of the two-dimensional blood vessel centerline extraction structure 210 according to 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.
The problems of the prior art are as follows: the mode of obtaining blood vessel evaluation parameters such as Fractional Flow Reserve (FFR) and the like through the pressure sensor needs a doctor to draw the sensor to a lesion part through a guide wire from an artery, and in the process, a plurality of difficulties need to be solved, such as invasiveness of operation, complexity of measurement and high cost of the pressure guide wire, can become obstacles for popularization of the FFR.
Example 1:
as shown in fig. 1, in order to solve the above problem, the present application provides a method for synthesizing a three-dimensional blood vessel, including:
s100, acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles;
s200, acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to image information of a coronary artery two-dimensional contrast image;
and S300, synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel central line and the three-dimensional blood vessel radius.
The application provides a three-dimensional blood vessel synthesis method, which can effectively simulate the state of blood vessels in a real scene, including the shape, the trend and the diameter information of the blood vessels, solves the influence of external conditions on images in the operation process, including the displacement of equipment, the beating of the heart and the respiration of a patient, solves the problems of high risk and high cost of guide wire measurement, and provides a basis for calculating blood vessel evaluation parameters such as blood flow reserve fraction (FFR).
Example 2:
on the basis of the embodiment 1, the embodiment is further optimized;
as shown in fig. 1, there is provided a method for synthesizing a three-dimensional blood vessel, comprising:
s100, as shown in fig. 2, acquiring image information of at least two coronary artery two-dimensional contrast images with different imaging angles, including:
s110, acquiring at least two coronary artery two-dimensional contrast image groups with different shooting angles;
s120, reading image information of each group of coronary artery two-dimensional contrast image group, wherein the image information comprises a shooting angle and a detection distance;
and S130, respectively selecting an interested two-dimensional contrast image from each group of coronary artery two-dimensional contrast images according to the detection distance.
S200, as shown in fig. 3, acquiring a three-dimensional blood vessel centerline and a three-dimensional blood vessel radius according to image information of a two-dimensional coronary angiography image, including:
s210, as shown in fig. 4, extracting a two-dimensional vessel centerline from each two-dimensional contrast image of interest, including:
s211, reading a coronary artery two-dimensional contrast image;
s212, obtaining a blood vessel section of interest;
s213, picking up the starting point, the seed point and the end point of the interested blood vessel section;
s214, 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;
s215, as shown in fig. 5, extracting at least one blood vessel local path line from each local blood vessel region map, including:
s2151, performing image enhancement on the local blood vessel region map to obtain a coarse blood vessel map with strong contrast, including: 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.
S2152, as shown in fig. 6, performing mesh division on the rough blood vessel map, and extracting at least one blood vessel local path line along a direction from the starting point to the ending point, including:
s21521, grid division is conducted on the rough blood vessel map;
s21522, searching a shortest time path between the starting point and intersection points on n peripheral grids along the extending direction of the blood vessel from the starting point to the ending point to serve as a second point, searching a shortest time path between the second point and the intersection points on the n peripheral 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;
s21523, 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 local path line.
S216, connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
s217, as shown in fig. 7, selecting a blood vessel path line as a two-dimensional blood vessel center line, including:
s2171, 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;
s2172, the least vessel path line is taken as the two-dimensional vessel centerline.
S220, as shown in fig. 8, projecting each two-dimensional blood vessel centerline into a three-dimensional space according to the shooting angle and the detection distance of each two-dimensional coronary artery angiography image, and synthesizing the three-dimensional blood vessel centerlines, including:
s221, establishing a three-dimensional coordinate system by taking the heart as a coordinate origin;
s222, acquiring a left-right angle alpha, a front-back angle beta and a distance S between a human body and a flat panel detector of each interested two-dimensional contrast image, wherein the coordinates of each point on a two-dimensional blood vessel central line are (x, y);
s223, projecting each point (x, y) into the three-dimensional space to obtain a series of three-dimensional coordinate points P, the coordinates being (x ", y", z "), including:
rotating each point (x, y) around the y axis to obtain a series of (x ', y ', z ') points, wherein the specific formula is as follows:
rotating the series of points (x ', y ', z ') around an x axis to obtain a series of three-dimensional coordinate points P with coordinates (x ", y", z ");
s224, projecting the radioactive source into a three-dimensional space to form a radioactive point R, wherein the method comprises the following steps:
obtaining the distance S' between the human body and the radioactive source in each interested two-dimensional radiography image;
according to the formula
Wherein, the coordinate of R in the three-dimensional space is (a, b, c).
S225, as shown in fig. 9, connecting each three-dimensional coordinate point P with the radiation point R, obtaining a point on the three-dimensional blood vessel centerline from the PR connection line, and sequentially connecting the points on the three-dimensional blood vessel centerline to obtain the three-dimensional blood vessel centerline, including:
s2251, correspondingly connecting a series of three-dimensional coordinate points P and radial points R obtained from the same two-dimensional contrast image of interest to obtain a plurality of PR straight lines;
s2252, acquiring points with the minimum distance between two PR straight lines at the same position of the blood vessel, namely points A and B;
s2253, connecting the point A with the point B, and acquiring the middle point of the line segment AB as a point on the center line of the three-dimensional blood vessel;
s2254, the points on the center lines of the obtained series of three-dimensional blood vessels are sequentially connected to obtain a three-dimensional blood vessel center line.
S230, obtaining a two-dimensional blood vessel contour line according to the two-dimensional blood vessel centerline in S210, as shown in fig. 10, in the prior art, calculating blood vessel evaluation parameters through a blood vessel three-dimensional model often requires extracting a blood vessel contour line, and since a blood vessel has a problem of curl and an unclear edge, the blood vessel contour line is particularly difficult to extract, and the operation data is huge and tedious, so how to quickly extract a blood vessel contour line, and the accuracy of extraction is always a problem that a technician needs to solve, in order to solve the above problem, the present application further implements S230, including:
s231, extracting a two-dimensional blood vessel center line according to the coronary artery two-dimensional contrast image;
s232, as shown in fig. 11, obtaining a straightened vessel image according to the two-dimensional vessel centerline, including:
s2321, straightening the center line of the two-dimensional blood vessel to obtain a blood vessel center straight line;
s2322, 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;
s2323, correspondingly arranging the two-dimensional blood vessel center lines of each unit along a blood vessel center straight line;
s2324, the correspondingly set image is a straightened blood vessel image.
S233, setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
S234, according to DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel;
s235, as shown in fig. 12, gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel, including:
s2351, dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
s2352, acquiring z points of each unit, which are positioned on a preset contour line of each blood vessel;
s2353, along the direction vertical to the central straight line of the blood vessel, closing the z points to the central straight line of the blood vessel in a grading way to generate z closing points, wherein z is a positive integer;
s2354, setting RGB difference threshold value as delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)When the blood vessel is closed, the closing point stops closing towards the center line of the blood vessel;
s2355, acquiring a close point as a contour point;
s2356, the smooth curves formed by connecting the contour points in sequence are the contour lines of the straightened blood vessel.
And S236, projecting the contour line of the straightened blood vessel back to the image for extracting the two-dimensional blood vessel center line to obtain the two-dimensional blood vessel contour line.
S240, acquiring the two-dimensional vessel radius in each interested two-dimensional contrast image according to the two-dimensional vessel contour line;
s250, obtaining the radius of the three-dimensional blood vessel according to the radius of the two-dimensional blood vessel, wherein the specific formula is as follows:
wherein R represents the three-dimensional vessel radius, R1、r2、rnThe two-dimensional vessel radii of the first, second and nth two-dimensional contrast images of interest are represented, respectively.
Two-dimensional contrast images with an angle difference of 30 degrees are usually selected to synthesize a three-dimensional blood vessel in the operation process, so that a three-dimensional blood vessel radius formula commonly adopted here is
S300, as shown in fig. 13, synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel centerline and the three-dimensional blood vessel radius, including:
s310, drawing a picture in a three-dimensional space along the corresponding three-dimensional blood vessel radius at each point on the center line of the three-dimensional blood vessel to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon approximate to a circle;
and S320, sequentially connecting points on two adjacent polygons according to a right-angle triangle form to obtain the three-dimensional blood vessel.
Obtaining a straightened blood vessel image according to the blood vessel central line; setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)(ii) a According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line; gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel; projecting the contour line of the straightened blood vessel back to the image of the center line of the blood vessel to obtain the contour line of the blood vessel; the extraction of the blood vessel contour line is rapid and accurate.
As shown in fig. 14, the present application provides a three-dimensional vessel synthesis system comprising: the three-dimensional blood vessel radius measuring device comprises an image reading device 100, a three-dimensional blood vessel center line obtaining device 200, a three-dimensional blood vessel radius obtaining device 300 and a three-dimensional blood vessel synthesizing device 400 which are connected in sequence, wherein the three-dimensional blood vessel radius obtaining device 300 is connected with the image reading device 100; the image reading apparatus 100 is configured to acquire image information of at least two coronary artery two-dimensional contrast images with different imaging angles; the three-dimensional blood vessel center line obtaining device 200 is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device and obtaining a three-dimensional blood vessel center line according to the image information; the three-dimensional vessel radius acquiring device 300 is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device, receiving the three-dimensional vessel center line transmitted by the three-dimensional vessel center line acquiring device, and acquiring the three-dimensional vessel radius according to the image information and the three-dimensional vessel center line; the three-dimensional blood vessel synthesizing device 400 is used for receiving the three-dimensional blood vessel center line transmitted by the three-dimensional blood vessel center line acquiring device and receiving the three-dimensional blood vessel radius transmitted by the three-dimensional blood vessel radius acquiring device, and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius.
As shown in fig. 15, in one embodiment of the present application, a three-dimensional vessel centerline acquisition device 200 includes: a two-dimensional blood vessel centerline extraction structure 210 and a three-dimensional blood vessel centerline acquisition structure 220 connected to the image reading apparatus 100, the two-dimensional blood vessel centerline extraction structure 210 being connected to the three-dimensional blood vessel centerline acquisition structure 220; the two-dimensional vessel centerline extraction structure 210 is configured to receive a two-dimensional angiography image of a coronary artery sent by the image reading apparatus, and extract a two-dimensional vessel centerline from each two-dimensional angiography image of interest; the three-dimensional blood vessel centerline acquisition structure 220 is configured to receive the two-dimensional blood vessel centerline sent by the two-dimensional blood vessel centerline extraction structure, and receive the two-dimensional blood vessel centerline sent by the image reading device, and project each two-dimensional blood vessel centerline into a three-dimensional space according to the shooting angle of each two-dimensional coronary artery angiography image, so as to synthesize a three-dimensional blood vessel centerline.
As shown in fig. 16, in one embodiment of the present application, a two-dimensional vessel centerline extraction structure 210 comprises: a central line extracting unit 211, a straightening unit 212, a first blood vessel contour line unit 213 and a second blood vessel contour line unit 214 which are connected in sequence; the center line extraction unit 211 is connected to the image reading apparatus 100, and is configured to extract a blood vessel center line from the coronary artery two-dimensional contrast image; the straightening unit 212 is used for obtaining a straightened blood vessel image according to the blood vessel central line extracted by the central line extraction unit 211; first blood vessel outline sheetElement 213 is used to set a blood vessel diameter threshold value D on the straightened blood vessel image sent by straightening unit 212Threshold(s)(ii) a According to DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel; gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel; the second blood vessel contour unit 214 is configured to project the contour of the straightened blood vessel sent by the first blood vessel contour unit 213 back to the image of the blood vessel centerline to obtain the blood vessel contour.
The present application provides a coronary artery analysis system comprising: the three-dimensional vessel synthesis system.
The present application provides a computer storage medium, and a computer program, when executed by a processor, implements the above-described three-dimensional blood vessel synthesis method.
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 (23)
1. A method of synthesizing a three-dimensional blood vessel, comprising:
acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles;
acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to the image information of the coronary artery two-dimensional radiography image;
and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel central line and the three-dimensional blood vessel radius.
2. The method for synthesizing three-dimensional blood vessels according to claim 1, wherein the method for acquiring image information of at least two coronary artery two-dimensional contrast images with different capturing angles comprises:
acquiring at least two groups of coronary artery two-dimensional contrast image groups with different shooting angles;
reading image information of each group of coronary artery two-dimensional contrast image group, wherein the image information comprises a shooting angle and a detection distance;
and respectively selecting an interested two-dimensional contrast image from each group of the coronary artery two-dimensional contrast images according to the detection distance.
3. The method for synthesizing three-dimensional vessel according to claim 2, wherein the method for acquiring the three-dimensional vessel centerline from the coronary artery two-dimensional contrast image comprises:
extracting a two-dimensional vessel centerline from each of the two-dimensional angiographic images of interest;
and projecting each two-dimensional blood vessel central line into a three-dimensional space according to the shooting angle and the detection distance of each two-dimensional coronary artery angiography image, and synthesizing the three-dimensional blood vessel central lines.
4. The method for synthesizing three-dimensional blood vessels according to claim 3, wherein the method for synthesizing the three-dimensional blood vessel center lines projects each of the two-dimensional blood vessel center lines into a three-dimensional space according to the capturing angle and the detection distance of each of the two-dimensional coronary angiography images, and comprises:
establishing a three-dimensional coordinate system by taking the heart as a coordinate origin;
acquiring a left-right angle alpha, a front-back angle beta and a distance S between a human body and a flat panel detector of each interested two-dimensional radiography image, wherein the coordinates of each point on the center line of the two-dimensional blood vessel are (x, y);
projecting each point (x, y) into a three-dimensional space to obtain a series of three-dimensional coordinate points P, wherein the coordinates are (x ', y ', z ');
projecting a radioactive source into the three-dimensional space to form a radioactive point R;
and connecting each three-dimensional coordinate point P with the radiation point R, obtaining points on the central line of the three-dimensional blood vessel from the PR connection line, and sequentially connecting the points on the central line of the three-dimensional blood vessel to obtain the central line of the three-dimensional blood vessel.
5. The method for synthesizing three-dimensional blood vessel according to claim 4, wherein the method for projecting each point (x, y) into three-dimensional space to obtain a series of three-dimensional coordinate points P with coordinates (x ", y", z ") comprises:
rotating each point (x, y) around the y axis to obtain a series of (x ', y ', z ') points, wherein the specific formula is as follows:
rotating the series of points (x ', y ', z ') around an x axis to obtain a series of three-dimensional coordinate points P with coordinates (x ", y", z ");
6. the method for synthesizing three-dimensional blood vessel according to claim 5, wherein the method for projecting a radioactive source into the three-dimensional space to form a radiation point R comprises:
obtaining a distance S' between a human body and a radioactive source in each two-dimensional contrast image of interest;
according to the formula
Wherein, the coordinate of R in the three-dimensional space is (a, b, c).
7. The method for synthesizing a three-dimensional blood vessel according to claim 6, wherein the method for connecting each three-dimensional coordinate point P with the radiation point R, obtaining a point on the centerline of the three-dimensional blood vessel from a PR connection line, and sequentially connecting the points on the centerline of the three-dimensional blood vessel to obtain the centerline of the three-dimensional blood vessel comprises:
correspondingly connecting a series of three-dimensional coordinate points P and radial points R obtained from the same two-dimensional contrast image of interest to obtain a plurality of PR straight lines;
acquiring points of the minimum distance between two PR straight lines at the same position of the blood vessel, wherein the points are respectively a point A and a point B;
connecting the point A with the point B, and acquiring the midpoint of the line segment AB as a point on the central line of the three-dimensional blood vessel;
and sequentially connecting the obtained points on the central line of the series of three-dimensional blood vessels to obtain the central line of the three-dimensional blood vessel.
8. The method for synthesizing three-dimensional blood vessel according to claim 3, wherein the method for obtaining the centerline and the radius of the three-dimensional blood vessel from the coronary artery two-dimensional contrast image comprises:
acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line;
acquiring the two-dimensional vessel radius in each two-dimensional contrast image of interest according to the two-dimensional vessel contour line;
and obtaining the three-dimensional vessel radius according to the two-dimensional vessel radius.
9. The method for synthesizing three-dimensional blood vessel according to claim 8, wherein the method for obtaining the three-dimensional blood vessel radius according to the two-dimensional blood vessel radius comprises:
wherein R represents the three-dimensional vessel radius, R1、r2、rnThe two-dimensional vessel radii of the first, second and nth two-dimensional contrast images of interest are represented, respectively.
10. The method for synthesizing three-dimensional blood vessels according to claim 3, wherein the method for extracting a two-dimensional blood vessel centerline from each of the two-dimensional coronary angiography images comprises:
reading a coronary artery two-dimensional contrast image;
obtaining a vessel segment of interest;
picking up a starting point, a seed point and an end point of the vessel segment of interest;
respectively segmenting two-dimensional contrast images between two adjacent points of a starting point, a seed point and an end point to obtain at least two local blood vessel region images;
extracting at least one blood vessel local path line from each local blood vessel region map;
connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
and selecting one blood vessel path line as the two-dimensional blood vessel central line.
11. The method for synthesizing three-dimensional blood vessels according to claim 10, wherein the method for extracting at least one blood vessel local path line from each local blood vessel region map comprises:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
and performing grid division on the rough blood vessel map, and extracting at least one blood vessel local path line along the direction from the starting point to the end point.
12. The method for synthesizing three-dimensional blood vessels according to claim 11, wherein the method for performing image enhancement processing on the local blood vessel region map to obtain a coarse blood vessel map with strong contrast comprises:
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 the rough blood vessel image with strong contrast is obtained.
13. The method for synthesizing three-dimensional blood vessels according to claim 12, wherein the method for gridding the rough blood vessel map and extracting at least one local path line of blood vessels 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 local path line.
14. The method for synthesizing a three-dimensional blood vessel according to claim 13, wherein the step of selecting one of the blood vessel route lines as the two-dimensional blood vessel center line 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 two-dimensional vessel centerline.
15. The method for synthesizing a three-dimensional blood vessel according to claim 8, wherein the method for obtaining a two-dimensional blood vessel contour line from the blood vessel centerline comprises:
extracting a two-dimensional blood vessel central line according to the coronary artery two-dimensional radiography image;
obtaining a straightened blood vessel image according to the two-dimensional blood vessel central line;
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel;
gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and projecting the contour line of the straightened blood vessel back to the image for extracting the center line of the two-dimensional blood vessel to obtain the contour line of the two-dimensional blood vessel.
16. The method for synthesizing three-dimensional blood vessel according to claim 15, wherein the method for obtaining the straightened blood vessel image according to the two-dimensional blood vessel central line comprises the following steps:
straightening the center line of the two-dimensional blood vessel to obtain a blood vessel center straight line;
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;
correspondingly arranging the two-dimensional blood vessel central line of each unit along the blood vessel central straight line;
and the correspondingly set image is the straightened blood vessel image.
17. The method for synthesizing a three-dimensional blood vessel according to claim 16, wherein the step of gradually drawing the preset contour line of the blood vessel toward the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel comprises the steps of:
dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on each preset blood vessel contour line;
respectively closing the z points to the blood vessel center straight line in a grading way along the direction vertical to the blood vessel center straight line to generate z closing points, wherein z is a positive integer;
setting RGB difference threshold to delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)When the blood vessel is in the closed state, the closing point stops closing towards the center line of the blood vessel;
acquiring the approach point as a contour point;
and sequentially connecting the contour points to form a smooth curve which is the contour line of the straightened blood vessel.
18. The method for synthesizing a three-dimensional blood vessel according to claim 1, wherein the method for synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel centerline and the three-dimensional blood vessel radius comprises:
drawing a picture in the three-dimensional space along the corresponding three-dimensional blood vessel radius to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon approximate to a circle;
and sequentially connecting points on two adjacent polygons according to a right-angle triangle form to obtain the three-dimensional blood vessel.
19. A three-dimensional blood vessel synthesis system for use in a method for synthesizing a three-dimensional blood vessel according to any one of claims 1 to 18, comprising: the three-dimensional blood vessel radius acquisition device is connected with the image reading device and the three-dimensional blood vessel center line acquisition device;
the image reading device is used for acquiring image information of at least two coronary artery two-dimensional contrast images with different shooting angles;
the three-dimensional blood vessel center line acquisition device is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device and acquiring a three-dimensional blood vessel center line according to the image information;
the three-dimensional blood vessel radius acquisition device is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device, receiving the three-dimensional blood vessel center line transmitted by the three-dimensional blood vessel center line acquisition device, and acquiring a three-dimensional blood vessel radius according to the image information and the three-dimensional blood vessel center line;
the three-dimensional blood vessel synthesizing device is used for receiving the three-dimensional blood vessel center line transmitted by the three-dimensional blood vessel center line acquiring device and receiving the three-dimensional blood vessel radius transmitted by the three-dimensional blood vessel radius acquiring device, and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius.
20. The three-dimensional vessel synthesis system according to claim 19, wherein the three-dimensional vessel centerline acquisition device comprises: the two-dimensional blood vessel center line extracting structure is connected with the image reading device, and the three-dimensional blood vessel center line acquiring structure is connected with the two-dimensional blood vessel center line extracting structure;
the two-dimensional vessel centerline extraction structure is used for receiving the coronary artery two-dimensional contrast images sent by the image reading device and extracting a two-dimensional vessel centerline from each two-dimensional contrast image of interest;
the three-dimensional blood vessel center line obtaining structure is used for receiving the two-dimensional blood vessel center line sent by the two-dimensional blood vessel center line extracting structure, receiving the two-dimensional blood vessel center line sent by the image reading device, projecting each two-dimensional blood vessel center line into a three-dimensional space according to the shooting angle of each two-dimensional coronary artery angiography image, and synthesizing the three-dimensional blood vessel center lines.
21. The three-dimensional vessel synthesis system of claim 20, wherein the two-dimensional vessel centerline extraction structure comprises: the center line extracting unit, the straightening unit, the first blood vessel contour line unit and the second blood vessel contour line unit are sequentially connected;
the central line extracting unit is connected with the image reading device and used for extracting a blood vessel central line according to a coronary artery two-dimensional contrast image;
the straightening unit is used for obtaining a straightened blood vessel image according to the blood vessel central line extracted by the central line extracting unit;
the first blood vessel contour line unit is used for setting a blood vessel diameter threshold value D on the straightened blood vessel image sent by the straightening unitThreshold(s)(ii) a According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line; gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and the second blood vessel contour line unit is used for projecting the contour line of the straightened blood vessel sent by the first blood vessel contour line unit back to the image of the blood vessel central line to obtain the blood vessel contour line.
22. A coronary artery analysis system, comprising: the three-dimensional vascular synthesis system of any one of claims 19 to 21.
23. A computer storage medium, wherein a computer program is executed by a processor to implement the method of synthesizing a three-dimensional blood vessel according to any one of claims 1 to 18.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337609.XA CN112419484B (en) | 2020-11-25 | 2020-11-25 | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium |
| PCT/CN2020/131699 WO2022109903A1 (en) | 2020-11-25 | 2020-11-26 | Three-dimensional blood vessel synthesis method and system, coronary artery analysis system and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337609.XA CN112419484B (en) | 2020-11-25 | 2020-11-25 | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112419484A true CN112419484A (en) | 2021-02-26 |
| CN112419484B CN112419484B (en) | 2024-03-22 |
Family
ID=74842083
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011337609.XA Active CN112419484B (en) | 2020-11-25 | 2020-11-25 | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112419484B (en) |
| WO (1) | WO2022109903A1 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114391792A (en) * | 2021-09-13 | 2022-04-26 | 南京诺源医疗器械有限公司 | Tumor prediction method and device based on narrow-band imaging and imaging endoscope |
| CN116712157A (en) * | 2023-01-16 | 2023-09-08 | 天津市鹰泰利安康医疗科技有限责任公司 | Irreversible perforation system for heart vessel |
| CN117036640A (en) * | 2023-10-10 | 2023-11-10 | 杭州脉流科技有限公司 | Coronary artery blood vessel model reconstruction method, device, equipment and storage medium |
| CN117058328A (en) * | 2023-10-11 | 2023-11-14 | 杭州脉流科技有限公司 | Coronary vessel tree classification method, apparatus, storage medium and program product |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116012522B (en) * | 2022-11-03 | 2023-11-14 | 成都寻创信息技术有限公司 | Three-dimensional imaging system for head, neck, jaw and face soft tissues, bones and blood vessels |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
| US20070116342A1 (en) * | 2003-09-25 | 2007-05-24 | Michael Zarkh | System and method for three-dimensional reconstruction of a tubular organ |
| CN109872321A (en) * | 2019-02-26 | 2019-06-11 | 数坤（北京）网络科技有限公司 | A kind of hemadostewnosis detection method and equipment |
| CN110287956A (en) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | Vessel centerline automatic matching method and device |
| CN110367965A (en) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | The method, apparatus and system of convenient measurement coronary artery assessment parameters |
| CN110889896A (en) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060036167A1 (en) * | 2004-07-03 | 2006-02-16 | Shina Systems Ltd. | Vascular image processing |
| CN109727242B (en) * | 2018-12-28 | 2021-02-26 | 上海联影医疗科技股份有限公司 | Blood vessel center line extraction method, device, computer equipment and storage medium |
| CN111161342B (en) * | 2019-12-09 | 2023-08-29 | 杭州脉流科技有限公司 | Method, apparatus, device, system and readable storage medium for obtaining fractional flow reserve based on coronary angiography image |
| CN111127552B (en) * | 2020-03-26 | 2020-07-14 | 南京景三医疗科技有限公司 | Blood vessel three-dimensional reconstruction method, medical device and storage medium |
-
2020
- 2020-11-25 CN CN202011337609.XA patent/CN112419484B/en active Active
- 2020-11-26 WO PCT/CN2020/131699 patent/WO2022109903A1/en active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6047080A (en) * | 1996-06-19 | 2000-04-04 | Arch Development Corporation | Method and apparatus for three-dimensional reconstruction of coronary vessels from angiographic images |
| US20070116342A1 (en) * | 2003-09-25 | 2007-05-24 | Michael Zarkh | System and method for three-dimensional reconstruction of a tubular organ |
| CN110367965A (en) * | 2018-09-19 | 2019-10-25 | 苏州润迈德医疗科技有限公司 | The method, apparatus and system of convenient measurement coronary artery assessment parameters |
| CN109872321A (en) * | 2019-02-26 | 2019-06-11 | 数坤（北京）网络科技有限公司 | A kind of hemadostewnosis detection method and equipment |
| CN110287956A (en) * | 2019-06-13 | 2019-09-27 | 北京理工大学 | Vessel centerline automatic matching method and device |
| CN110889896A (en) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis |
Non-Patent Citations (2)
| Title |
|---|
| 丛伟建 等: "基于造影图像的能量场反投影合成冠脉三维重建", 中国科学:信息科学 * |
| 陈晓冬 等: "基于造影图像的冠状动脉血管三维表面重建方法", 工程图学学报 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114391792A (en) * | 2021-09-13 | 2022-04-26 | 南京诺源医疗器械有限公司 | Tumor prediction method and device based on narrow-band imaging and imaging endoscope |
| CN114391792B (en) * | 2021-09-13 | 2023-02-24 | 南京诺源医疗器械有限公司 | Tumor prediction method and device based on narrow-band imaging and imaging endoscope |
| CN116712157A (en) * | 2023-01-16 | 2023-09-08 | 天津市鹰泰利安康医疗科技有限责任公司 | Irreversible perforation system for heart vessel |
| CN116712157B (en) * | 2023-01-16 | 2024-03-08 | 天津市鹰泰利安康医疗科技有限责任公司 | Irreversible perforation system for heart vessel |
| CN117036640A (en) * | 2023-10-10 | 2023-11-10 | 杭州脉流科技有限公司 | Coronary artery blood vessel model reconstruction method, device, equipment and storage medium |
| CN117036640B (en) * | 2023-10-10 | 2023-12-19 | 杭州脉流科技有限公司 | Coronary artery blood vessel model reconstruction method, device, equipment and storage medium |
| CN117058328A (en) * | 2023-10-11 | 2023-11-14 | 杭州脉流科技有限公司 | Coronary vessel tree classification method, apparatus, storage medium and program product |
| CN117058328B (en) * | 2023-10-11 | 2024-01-09 | 杭州脉流科技有限公司 | Coronary vessel tree classification method, apparatus, storage medium and program product |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022109903A1 (en) | 2022-06-02 |
| CN112419484B (en) | 2024-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11896416B2 (en) | Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images | |
| CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
| CN112164020B (en) | Method, device, analysis system and storage medium for accurately extracting blood vessel center line | |
| US20220277447A1 (en) | Method and apparatus for acquiring contour line of blood vessel according to centerline of blood vessel | |
| CN107451983A (en) | The three-dimensional fusion method and system of CT images | |
| CN112132949A (en) | Blood vessel three-dimensional modeling method, device and system with narrow lesion interval | |
| WO2017148332A1 (en) | Perfusion analysis method and device | |
| CN111476791B (en) | Image processing method, image processing apparatus, and non-transitory computer readable medium | |
| WO2022109906A1 (en) | Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium | |
| CN111260703B (en) | Method, system, medium and storage medium for obtaining spinal straightening image set | |
| US20200118264A1 (en) | Shell-constrained localization of vasculature | |
| CN112419276A (en) | Method for regulating blood vessel contour and central line and storage medium | |
| US20220319116A1 (en) | Method and device for synthesizing mathematical model of blood vessel having stenotic lesion | |
| CN112132882A (en) | Method and device for extracting blood vessel central line from coronary artery two-dimensional contrast image | |
| CN113888566B (en) | Target contour curve determination method and device, electronic equipment and storage medium | |
| WO2022109902A1 (en) | Three-dimensional vascular centerline synthesis method and system, and storage medium | |
| CN114596311A (en) | Blood vessel function evaluation method and blood vessel function evaluation device based on blood vessel image | |
| US11317875B2 (en) | Reconstruction of flow data | |
| WO2022109904A1 (en) | Method for two-dimensional image selection and three-dimensional blood vessel synthesis and storage medium | |
| CN111815584B (en) | Method and system for acquiring heart gravity center based on CT sequence image | |
| CN111798468B (en) | Image processing method and device, storage medium and electronic terminal | |
| JP7446645B2 (en) | System for acquiring aortic images based on deep learning | |
| US20200193590A1 (en) | Quantitative evaluation of time-varying data | |
| WO2018173273A1 (en) | Living body model generation device, living body model generation method, and living body model generation program | |
| CN114049282B (en) | Coronary artery construction method, device, terminal 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 |