CN112419484B - Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium - Google Patents
Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium Download PDFInfo
- Publication number
- CN112419484B CN112419484B CN202011337609.XA CN202011337609A CN112419484B CN 112419484 B CN112419484 B CN 112419484B CN 202011337609 A CN202011337609 A CN 202011337609A CN 112419484 B CN112419484 B CN 112419484B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- blood vessel
- line
- vascular
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000002792 vascular Effects 0.000 title claims abstract description 114
- 210000004351 coronary vessel Anatomy 0.000 title claims abstract description 39
- 238000001308 synthesis method Methods 0.000 title claims abstract description 19
- 238000004458 analytical method Methods 0.000 title claims abstract description 9
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 274
- 238000000034 method Methods 0.000 claims abstract description 71
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 31
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 26
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 26
- 238000000605 extraction Methods 0.000 claims description 24
- 238000002586 coronary angiography Methods 0.000 claims description 20
- 238000010586 diagram Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 9
- 230000002285 radioactive effect Effects 0.000 claims description 9
- 238000002601 radiography Methods 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 abstract description 8
- 230000017531 blood circulation Effects 0.000 abstract description 6
- 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
- 230000006870 function Effects 0.000 description 7
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 230000007812 deficiency Effects 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
- 230000000544 hyperemic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 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
- 210000004165 myocardium Anatomy 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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)
- Public Health (AREA)
- Primary Health Care (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application provides a three-dimensional vascular synthesis method, a three-dimensional vascular synthesis 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 contrast image; and synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius. The method and the device can effectively simulate the vascular state in the real scene, comprise the shape, trend and diameter information of the blood vessel, solve the influence of external conditions on images in the operation process, including the displacement of equipment, the beating of heart and the respiration of a patient, solve the problems of high risk and high cost of guide wire measurement, and provide a basis for calculating blood flow reserve score FFR and other vascular assessment parameters.
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a three-dimensional vascular 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 plaque on the vessel wall, which in turn leads to stenosis of the vessel; especially, the stenosis of blood vessels around the coronary artery will lead to myocardial blood supply deficiency, induce coronary heart disease, angina pectoris and other diseases, and cause serious threat to human health. According to statistics, the number of patients with the existing coronary heart disease in China is about 1100 ten thousand, and the number of patients with cardiovascular interventional operation treatment is increased by more than 10% each year.
Although the conventional medical detection means such as Coronary Angiography (CAG) and Computed Tomography (CT) can show the severity of coronary stenosis of heart, the ischemia of the coronary artery cannot be accurately evaluated. In order to improve the accuracy of coronary blood vessel function evaluation, pijls in 1993 proposed a new index of calculating coronary blood vessel function by pressure measurement, namely fractional flow reserve (Fractional Flow Reserve, FFR), and FFR has become a gold standard for coronary stenosis function evaluation through long-term basic and clinical studies.
Fractional Flow Reserve (FFR) is generally referred to as fractional myocardial flow reserve, defined as the ratio of the maximum blood flow that a diseased coronary can provide to the myocardium to the maximum blood flow at which the coronary is completely normal, and studies have shown that the ratio of blood flow can be replaced with a pressure value at the maximum hyperemic state of the coronary. That is, the FFR value can be measured and then calculated by measuring the pressure at the distal end stenosis of the coronary artery and the proximal end pressure of the coronary artery by the pressure sensor in the maximum congestion state of the coronary artery.
The problems are that: the blood vessel assessment parameters such as fractional flow reserve FFR are obtained through the pressure sensor, a doctor is required to pull the sensor from the artery to a lesion through the guide wire, and a plurality of difficulties in the process need to be solved, such as the invasiveness of operation, the complexity of measurement and the high cost of the pressure guide wire, can be barriers to FFR popularization.
Disclosure of Invention
The invention provides a three-dimensional vascular synthesis method, a three-dimensional vascular synthesis 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 contrast image;
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 above three-dimensional blood vessel synthesizing method, the method for acquiring image information of at least two coronary two-dimensional contrast images with different shooting angles includes:
Acquiring at least two coronary artery two-dimensional contrast image groups with different shooting angles;
reading image information of each coronary artery two-dimensional contrast image group, wherein the image information comprises shooting angles and detection distances;
and selecting a two-dimensional contrast image of interest from each group of the coronary two-dimensional contrast images according to the detection distance.
Optionally, the method for synthesizing a three-dimensional blood vessel includes:
extracting a two-dimensional vessel centerline from each of the two-dimensional contrast images of interest;
and projecting each two-dimensional vessel center line into a three-dimensional space according to the shooting angle and the detection distance of each coronary artery two-dimensional radiography image, and synthesizing the three-dimensional vessel center lines.
Optionally, in the above three-dimensional vessel synthesizing method, the method for projecting each two-dimensional vessel center line into a three-dimensional space according to a shooting angle and a detection distance of each two-dimensional coronary angiography image, the method for synthesizing the three-dimensional vessel center line includes:
using a heart as an origin of coordinates, and establishing a three-dimensional coordinate system;
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 two-dimensional contrast image of interest, wherein the coordinates of each point on the central 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 radial point R, obtaining points on the three-dimensional blood vessel center line from PR connecting lines, and sequentially connecting the points on the three-dimensional blood vessel center line to obtain the three-dimensional blood vessel center line.
Optionally, in the above three-dimensional vascular synthesis method, the method of projecting each point (x, y) into a three-dimensional space to obtain a series of three-dimensional coordinate points P, where the coordinates are (x ", y", z "), includes:
the points (x, y) are rotated around the y axis to obtain (x ', y ', z ') series points, and the specific formula is as follows:
rotating the (x ', y', z ') series points around the x axis to obtain a series of three-dimensional coordinate points P, wherein the coordinates are (x', y ', z');
optionally, the method for synthesizing the three-dimensional blood vessel includes the steps of:
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 coordinates of R in the three-dimensional space are (a, b, c).
Optionally, in the above method for synthesizing a three-dimensional blood vessel, the connecting each three-dimensional coordinate point P with the radiation point R, obtaining a point on the three-dimensional blood vessel centerline from a PR connecting line, and sequentially connecting the points on the three-dimensional blood vessel centerline to obtain the three-dimensional blood vessel centerline, which includes:
correspondingly connecting a series of three-dimensional coordinate points P and radiation points R obtained from the same two-dimensional radiography image of interest to obtain a plurality of PR straight lines;
acquiring the minimum distance point between two PR lines at the same position of the blood vessel, wherein the minimum distance point is a point A and a point B respectively;
connecting the point A with the point B, and acquiring the midpoint of the AB line segment as a point on the central line of the three-dimensional blood vessel;
sequentially connecting the obtained points on the central line of the three-dimensional blood vessel to obtain the central line of the three-dimensional blood vessel.
Optionally, the method for synthesizing a three-dimensional blood vessel includes:
Acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line;
acquiring a two-dimensional vessel radius in each two-dimensional contrast image of interest according to the two-dimensional vessel contour line;
and acquiring the three-dimensional blood vessel radius according to the two-dimensional blood vessel radius.
Optionally, the method for synthesizing a three-dimensional blood vessel includes:
wherein R represents three-dimensional vessel radius, R 1 、r 2 、r n The two-dimensional vessel radii of the first, second, and nth two-dimensional contrast images of interest are represented, respectively.
Optionally, in the above three-dimensional vessel synthesis method, the method for extracting a two-dimensional vessel centerline from each of the two-dimensional coronary angiography images includes:
reading a two-dimensional coronary angiography image;
acquiring a vessel segment of interest;
picking up a start point, a seed point and an end point of the vessel segment of interest;
dividing two-dimensional contrast images between two adjacent points of a starting point, a seed point and an ending point respectively to obtain at least two local vessel region diagrams;
extracting at least one local vascular path line from each local vascular zone map;
Connecting the corresponding blood vessel local route lines on each local blood vessel region graph to obtain at least one blood vessel route line;
and selecting one vascular path line as the two-dimensional vascular center line.
Optionally, the method for synthesizing a three-dimensional blood vessel, wherein the method for extracting at least one local path line of the blood vessel from each local blood vessel region map comprises the following steps:
performing image enhancement processing on the local vascular region map to obtain a rough vascular map with strong contrast;
and meshing the rough blood vessel graph, and extracting at least one local path line of the blood vessel along the direction from the starting point to the ending point.
Optionally, in the above three-dimensional vascular synthesis method, the method for performing image enhancement processing on the local vascular region map to obtain a rough vascular map with strong contrast includes:
in each local vessel region graph, the vessel segment of interest is taken as a foreground, other regions are taken as a background, the foreground is strengthened, the background is weakened, and the rough vessel graph with strong contrast is obtained.
Optionally, in the above three-dimensional vascular synthesis method, the method for meshing the rough vascular map, and extracting at least one local vascular path along the direction from the start point to the end point includes:
Grid dividing the rough blood vessel map;
searching a shortest time path of the intersection points on the starting point and the n grids on the periphery along the extending direction of the blood vessel from the starting point to the ending point as a second point, searching the shortest time path of the intersection points on the second point and the n grids on the periphery as a third point, and repeating the steps until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and according to the searching sequence, connecting lines from the starting point to the ending point in the extending direction of the blood vessel, and obtaining at least one local path line of the blood vessel.
Optionally, in the above three-dimensional vascular synthesis method, the method for selecting one vascular path line as the two-dimensional vascular centerline includes:
summing the time taken from the start point to the end point for each vessel path line if the vessel path line is two or more;
the vessel path line at least when taken is taken as the two-dimensional vessel centerline.
Optionally, the method for obtaining a two-dimensional vascular contour line according to the vascular centerline in the three-dimensional vascular synthesis method is characterized by comprising the following steps:
Extracting a two-dimensional blood vessel center line according to the two-dimensional coronary angiography image;
obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line;
setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
According to said D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
gradually converging the preset contour line of the blood vessel towards the center 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 an image for extracting the central line of the two-dimensional blood vessel, so as to obtain the contour line of the two-dimensional blood vessel.
Optionally, in the above three-dimensional blood vessel synthesizing method, the method for obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line includes:
straightening the two-dimensional blood vessel center line to obtain a blood vessel center line;
dividing the local vessel region map into x units along the vessel extending direction from the starting point to the ending point, wherein x is a positive integer;
correspondingly arranging the two-dimensional blood vessel center line of each unit along the blood vessel center line;
and the image after corresponding setting is the straightened blood vessel image.
Optionally, in the above three-dimensional vascular synthesis method, the method for gradually converging the preset contour line of the blood vessel toward the center line of the blood vessel to obtain the contour line of the straightened blood vessel includes:
Dividing the vascular preset contour line into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on a preset contour line of each blood vessel;
respectively converging z points towards the central straight line of the blood vessel in a grading way along the direction perpendicular to the central straight line of the blood vessel to generate z converging points, wherein z is a positive integer;
setting RGB difference threshold to delta RGB Threshold value Comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line along the direction perpendicular to the blood vessel center straight line every time the blood vessel center straight line is close, wherein the difference value is less than or equal to delta RGB Threshold value When the blood vessel is closed, the closing point stops closing to the center line of the blood vessel;
acquiring the close points as contour points;
and connecting the contour points in sequence 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 center line and the three-dimensional blood vessel radius comprises the following steps:
drawing a picture in the three-dimensional space along the radius of the corresponding three-dimensional blood vessel at the point on the central line of each three-dimensional blood vessel to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon similar to a circle;
And sequentially connecting the points on two adjacent polygons according to the form of a right triangle to obtain the three-dimensional blood vessel.
In a second aspect, the present application provides a three-dimensional vascular 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 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 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 and the three-dimensional vessel center line transmitted by the three-dimensional vessel center line acquisition device, and acquiring the three-dimensional vessel radius according to the image information and the three-dimensional vessel center line;
The three-dimensional vascular synthesis device is used for receiving the three-dimensional vascular center line transmitted by the three-dimensional vascular center line acquisition device and receiving the three-dimensional vascular radius transmitted by the three-dimensional vascular radius acquisition device, and synthesizing a three-dimensional vascular according to the three-dimensional vascular center line and the three-dimensional vascular radius.
Optionally, in the three-dimensional vascular synthesis system, the three-dimensional vascular centerline acquisition device includes: the two-dimensional blood vessel central line extraction structure is connected with the image reading device and the three-dimensional blood vessel central line acquisition structure;
the two-dimensional blood vessel center line extraction structure is used for receiving the coronary artery two-dimensional contrast images sent by the image reading device and extracting a two-dimensional blood vessel center line from each two-dimensional contrast image of interest;
the three-dimensional blood vessel center line acquisition structure is used for receiving the two-dimensional blood vessel center lines sent by the two-dimensional blood vessel center line extraction structure, receiving the two-dimensional blood vessel center lines sent by the image reading device, and projecting each two-dimensional blood vessel center line into a three-dimensional space according to the shooting angle of each coronary artery two-dimensional radiography image to synthesize the three-dimensional blood vessel center lines.
Optionally, in the above three-dimensional blood vessel synthesizing method, the two-dimensional blood vessel centerline extracting structure includes: the central line extraction 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 extraction unit is connected with the image reading device and is used for extracting a blood vessel central line according to the two-dimensional coronary angiography 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 extraction unit;
the first blood vessel contour line unit is used for setting a blood vessel diameter threshold D on the straightened blood vessel image sent by the straightening unit Threshold value The method comprises the steps of carrying out a first treatment on the surface of the According to said D Threshold value Generating a blood vessel preset contour line on two sides of the blood vessel center straight line; gradually converging the preset contour line of the blood vessel towards the center straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
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 center line to obtain a blood vessel contour line.
In a third aspect, the present application provides a coronary artery analysis system comprising: the three-dimensional vascular synthesis system described above.
In a fourth aspect, the present application provides a computer storage medium, which when executed by a processor implements the above-described three-dimensional vascular synthesis method.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the three-dimensional vascular synthesis method can effectively simulate the vascular state in a real scene, comprises the shape, trend and diameter information of the blood vessel, solves the problems of influence of external conditions on images, including displacement of equipment, beating of heart and respiration of a patient, high risk of guide wire measurement and high cost in the surgical process, and provides a basis for calculating blood flow reserve score FFR and other vascular assessment parameters.
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 do not constitute a limitation on the invention. In the drawings:
reference numerals are described below:
FIG. 1 is a flow chart of a method of synthesizing a three-dimensional blood vessel of the present application;
FIG. 2 is a flow chart of S100 of the present application;
FIG. 3 is a flow chart of S200 of the present application;
FIG. 4 is a flowchart of S210 of the present application;
FIG. 5 is a flow chart 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 flow chart 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 flow chart of S235 of the present application;
FIG. 13 is a flowchart of S300 of the present application;
FIG. 14 is a block diagram of a three-dimensional vascular synthesis system of the present application;
FIG. 15 is another block diagram of the three-dimensional vascular synthesis system of the present application;
fig. 16 is a block diagram of a two-dimensional vessel centerline extraction structure 210 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 specific embodiments of the present invention and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Various embodiments of the invention are disclosed in the following drawings, in which details of the practice are set forth in the following description for the purpose of clarity. However, it should be understood that these practical details are not to be taken as limiting the invention. That is, in some embodiments of the invention, these practical details are unnecessary. Moreover, for the purpose of simplifying the drawings, some conventional structures and components are shown in the drawings in a simplified schematic manner.
Problems of the prior art: the blood vessel assessment parameters such as fractional flow reserve FFR are obtained through the pressure sensor, a doctor is required to pull the sensor from the artery to a lesion through the guide wire, and a plurality of difficulties in the process need to be solved, such as the invasiveness of operation, the complexity of measurement and the high cost of the pressure guide wire, can be barriers to FFR popularization.
Example 1:
as shown in fig. 1, in order to solve the above-mentioned problems, the present application provides a three-dimensional vascular synthesis method, which includes:
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 two-dimensional coronary angiography image;
S300, synthesizing the three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius.
The three-dimensional vascular synthesis method can effectively simulate the vascular state in a real scene, comprises the shape, trend and diameter information of the blood vessel, solves the problems of influence of external conditions on images, including displacement of equipment, beating of heart and respiration of a patient, high risk of guide wire measurement and high cost in the surgical process, and provides a basis for calculating blood flow reserve score FFR and other vascular assessment parameters.
Example 2:
further optimizing the present embodiment on the basis of embodiment 1;
as shown in fig. 1, there is provided a three-dimensional vascular synthesis method comprising:
s100, as shown in FIG. 2, acquiring image information of at least two coronary two-dimensional contrast images with different shooting angles, wherein the image information comprises:
s110, acquiring at least two coronary artery two-dimensional contrast image groups with different shooting angles;
s120, reading image information of each coronary two-dimensional contrast image group, including shooting angles and detection distances;
s130, selecting a two-dimensional contrast image of interest from each group of coronary two-dimensional contrast images according to the detection distance.
S200, as shown in FIG. 3, acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to image information of a two-dimensional coronary angiography image comprises:
s210, as shown in FIG. 4, extracting a two-dimensional vessel centerline from each two-dimensional contrast image of interest, comprising:
s211, reading a coronary artery two-dimensional contrast image;
s212, acquiring a blood vessel segment of interest;
s213, picking up a starting point, a seed point and an ending point of the blood vessel segment of interest;
s214, respectively dividing two-dimensional contrast images between two adjacent points of a starting point, a seed point and an ending point to obtain at least two local vessel region diagrams;
s215, as shown in fig. 5, extracting at least one local vascular path line from each local vascular zone map, including:
s2151, carrying out image enhancement processing on the local blood vessel region graph to obtain a rough blood vessel graph with strong contrast, comprising: in each local vessel region diagram, the vessel segment of interest is taken as a foreground, other regions are taken as a background, the foreground is strengthened, the background is weakened, and a rough vessel diagram with strong contrast is obtained.
S2152, as shown in fig. 6, grids the rough vessel map, extracting at least one vessel local path line along the direction from the start point to the end point, including:
S21521, performing grid division on the rough blood vessel map;
s21522, searching the shortest time path of the intersection points on the starting point and the n grids at the periphery as a second point along the extending direction of the blood vessel from the starting point to the ending point, searching the shortest time path of the intersection points on the second point and the n grids at the periphery as a third point, and repeating the steps until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
s21523, obtaining at least one local path line of the blood vessel by connecting the blood vessel extending directions from the start point to the end point in the search order.
S216, connecting the corresponding blood vessel local route lines on each local blood vessel region graph to obtain at least one blood vessel route line;
s217, as shown in FIG. 7, selecting a blood vessel path line as a two-dimensional blood vessel center line comprises:
s2171, if there are two or more blood vessel path lines, summing the time taken from the start point to the end point for each blood vessel path line;
s2172, the smallest blood vessel route line is taken as the two-dimensional blood vessel center line.
S220, as shown in FIG. 8, according to the shooting angle and the detection distance of each coronary two-dimensional radiography image, projecting each two-dimensional vessel center line into a three-dimensional space, and synthesizing the three-dimensional vessel center line, comprising:
S221, establishing a three-dimensional coordinate system by taking a heart as an origin of coordinates;
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 radiography image, wherein coordinates of each point on a two-dimensional blood vessel central line are (x, y);
s223, 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 '), and the method comprises the following steps:
the points (x, y) are rotated around the y axis to obtain (x ', y ', z ') series points, and the specific formula is as follows:
rotating the (x ', y', z ') series points around the x axis to obtain a series of three-dimensional coordinate points P, wherein the coordinates are (x', y ', z');
s224, projecting the radioactive source into a three-dimensional space to form a radioactive point R, comprising:
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 coordinates of R in the three-dimensional space are (a, b, c).
S225, as shown in fig. 9, connecting each three-dimensional coordinate point P with a radial point R, obtaining a point on a three-dimensional vessel centerline from a PR connection, and sequentially connecting the points on the three-dimensional vessel centerline to obtain a three-dimensional vessel centerline, including:
s2251, a series of three-dimensional coordinate points P and radiation points R obtained from the same two-dimensional radiography image of interest are correspondingly connected to obtain a plurality of PR straight lines;
S2252, obtaining the minimum distance point between two PR lines at the same position of the blood vessel, namely an A point and a B point;
s2253, connecting the point A with the point B, and obtaining the midpoint of the AB line segment as a point on the central line of the three-dimensional blood vessel;
s2254, sequentially connecting points on the central line of the obtained series of three-dimensional blood vessels to obtain the central line of the three-dimensional blood vessels.
In S230, as shown in fig. 10, in the prior art, the blood vessel contour line is often required to be extracted when the blood vessel evaluation parameter is calculated by the blood vessel three-dimensional model, because the blood vessel has the problems of curling and unclear edges, the blood vessel contour extraction is particularly difficult, and the calculation data is huge and tedious, so how to quickly extract the blood vessel contour line, and the accuracy of the extraction are always problems that the technician needs to solve, in order to solve the above problems, the present application further implements S230, including:
s231, extracting a two-dimensional blood vessel center line according to the two-dimensional coronary angiography image;
s232, as shown in FIG. 11, obtaining a straightened vessel image from a two-dimensional vessel centerline, comprising:
s2321, straightening a two-dimensional blood vessel center line to obtain a blood vessel center line;
S2322, dividing the local vascular zone map into x units along the vascular extension direction from the starting point to the ending point, wherein x is a positive integer;
s2323, correspondingly arranging the two-dimensional blood vessel center line of each unit along the blood vessel center line;
s2324, the image after corresponding setting is a straightened blood vessel image.
S233, setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
S234, according to D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
s235, as shown in FIG. 12, the preset contour line of the blood vessel is gradually closed to the center line of the blood vessel, and the contour line of the straightened blood vessel is obtained, which comprises:
s2351, dividing the blood vessel preset contour line 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, respectively converging z points towards the center line of the blood vessel in a grading manner along the direction perpendicular to the center line of the blood vessel to generate z converging points, wherein z is a positive integer;
s2354, setting RGB difference value threshold to ΔRGB Threshold value Each time the blood vessel is closed along the direction perpendicular to the straight line of the center of the blood vessel, the RGB value of the closed point is compared with the RGB value of the point on the straight line of the center of the blood vessel, and when the difference value is less than or equal to delta RGB Threshold value When the blood vessel is closed, the closing point stops closing to the center of the blood vessel linearly;
s2355, acquiring a close point as a contour point;
s2356, the smooth curve formed by connecting the contour points in sequence is the contour line of the straightened blood vessel.
S236, projecting the contour line of the straightened blood vessel back to the image for extracting the center line of the two-dimensional blood vessel, and obtaining the contour line of the two-dimensional blood vessel.
S240, acquiring a two-dimensional vessel radius in each interested two-dimensional contrast image according to the two-dimensional vessel contour line;
s250, acquiring a three-dimensional blood vessel radius according to the two-dimensional blood vessel radius, wherein the specific formula is as follows:
wherein R represents three-dimensional vessel radius, R 1 、r 2 、r n The two-dimensional vessel radii of the first, second, and nth two-dimensional contrast images of interest are represented, respectively.
In the operation process, two-dimensional contrast images with the angle difference of 30 degrees are usually selected to synthesize a three-dimensional blood vessel, so that the three-dimensional blood vessel radius formula adopted here is
S300, as shown in fig. 13, synthesizing a three-dimensional blood vessel according to a three-dimensional blood vessel center line and a three-dimensional blood vessel radius, comprising:
s310, drawing points on the central line of each three-dimensional blood vessel in a three-dimensional space along the radius of the corresponding three-dimensional blood vessel to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon similar to a circle;
S320, connecting points on two adjacent polygons in sequence according to the form of a right triangle, and obtaining the three-dimensional blood vessel.
According to the blood vessel center line, a straightened blood vessel image is obtained; setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value The method comprises the steps of carrying out a first treatment on the surface of the According to said D Threshold value Generating a blood vessel preset contour line on two sides of the blood vessel center straight line; gradually converging the preset contour line of the blood vessel towards the center 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 central line of the blood vessel to obtain the contour line of the blood vessel; the extraction of the blood vessel contour line is quick and accurate.
As shown in fig. 14, the present application provides a three-dimensional vascular synthesis system comprising: the three-dimensional blood vessel radius acquisition device 300 is connected with the image reading device 100, and the three-dimensional blood vessel center line acquisition device 200, the three-dimensional blood vessel radius acquisition device 300 and the three-dimensional blood vessel synthesis device 400 are connected in sequence; the image reading apparatus 100 is configured to acquire image information of at least two coronary two-dimensional contrast images with different imaging angles; the three-dimensional vessel centerline obtaining device 200 is configured to receive the image information of the two-dimensional coronary angiography image transmitted by the image reading device, and obtain a three-dimensional vessel centerline according to the image information; the three-dimensional vessel radius acquisition device 300 is used for receiving the image information of the coronary artery two-dimensional contrast image transmitted by the image reading device and receiving the three-dimensional vessel center line transmitted by the three-dimensional vessel center line acquisition device, and acquiring the three-dimensional vessel radius according to the image information and the three-dimensional vessel center line; the three-dimensional vascular synthesis device 400 is configured to receive the three-dimensional vascular centerline transmitted by the three-dimensional vascular centerline acquisition device and the three-dimensional vascular radius transmitted by the three-dimensional vascular radius acquisition device, and synthesize a three-dimensional vascular according to the three-dimensional vascular centerline and the three-dimensional vascular radius.
As shown in fig. 15, in one embodiment of the present application, the three-dimensional blood 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 the two-dimensional coronary angiography images sent by the image reading device, and extract a two-dimensional vessel centerline from each two-dimensional angiography image of interest; the three-dimensional vessel centerline obtaining structure 220 is configured to receive the two-dimensional vessel centerlines sent by the two-dimensional vessel centerline extracting structure, and receive the two-dimensional vessel centerlines sent by the image reading device, and project each two-dimensional vessel centerline into the three-dimensional space according to the shooting angle of each two-dimensional coronary angiography image, so as to synthesize the three-dimensional vessel centerlines.
As shown in fig. 16, in one embodiment of the present application, the two-dimensional vessel centerline extraction structure 210 includes: a center line extraction 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 sequentially connected; 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 two-dimensional coronary angiography image; the straightening unit 212 is configured to obtain a straightened blood vessel image according to the blood vessel center line extracted by the center line extraction unit 211; the first blood vessel contour line unit 213 is used for setting a blood vessel diameter threshold D on the straightened blood vessel image sent by the straightening unit 212 Threshold value The method comprises the steps of carrying out a first treatment on the surface of the According to D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line; gradually converging a preset contour line of a blood vessel towards a straight line at the center of the blood vessel to obtain a contour line of the straightened blood vessel; the second blood vessel contour line unit 214 is configured to project the contour line of the straightened blood vessel sent by the first blood vessel contour line unit 213 back onto the image of the blood vessel center line, so as to obtain a blood vessel contour line.
The present application provides a coronary artery analysis system comprising: the three-dimensional vascular synthesis system described above.
The application provides a computer storage medium, and a computer program realizes the three-dimensional vascular synthesis method when being executed by a processor.
Those skilled in the art will appreciate that the various aspects of the present invention may be implemented as a system, method, or computer program product. Accordingly, aspects of the invention may be embodied in the following forms, namely: an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining hardware and software aspects may all generally be referred to herein as a "circuit," module "or" system. Furthermore, in some embodiments, aspects of the invention may also be implemented 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 methods and/or systems 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 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 method and/or system as herein, such as a computing platform for executing a plurality of instructions, are performed by a data processor. 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 are optionally also provided.
Any combination of one or more computer readable may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any 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 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.
The computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, either in baseband or as part of a carrier wave. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. 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 remote computers, the remote computer may be connected to the user computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (e.g., connected 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 (article of manufacture).
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 device or other devices provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
The foregoing embodiments of the present invention have been described in some detail by way of illustration of the principles of the invention, and it is to be understood that the invention is not limited to the specific embodiments of the invention but is intended to cover modifications, equivalents, alternatives and modifications within the spirit and principles of the invention.
Claims (21)
1. A three-dimensional vascular synthesis method is characterized by comprising the following steps of
Acquiring at least two coronary artery two-dimensional contrast image groups with different shooting angles;
reading image information of each coronary artery two-dimensional contrast image group, wherein the image information comprises shooting angles and detection distances;
selecting a two-dimensional contrast image of interest from each group of the coronary two-dimensional contrast images according to the detection distance;
acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to the coronary artery two-dimensional radiography image;
synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius;
the method for acquiring the three-dimensional blood vessel center line and the three-dimensional blood vessel radius according to the coronary artery two-dimensional contrast image comprises the following steps:
extracting a two-dimensional vessel centerline from each of the two-dimensional contrast images of interest;
projecting each two-dimensional blood vessel center line into a three-dimensional space according to the shooting angle and the detection distance of each coronary artery two-dimensional contrast image, and synthesizing the three-dimensional blood vessel center lines;
acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line;
acquiring a two-dimensional vessel radius in each two-dimensional contrast image of interest according to the two-dimensional vessel contour line;
Acquiring the three-dimensional blood vessel radius according to the two-dimensional blood vessel radius;
the method for acquiring the two-dimensional vascular contour line according to the two-dimensional vascular center line comprises the following steps:
extracting a two-dimensional blood vessel center line according to the two-dimensional coronary angiography image;
obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line;
on the straightened vessel image, settingVascular diameter threshold D Threshold value ；
According to said D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
gradually converging the preset contour line of the blood vessel towards the center 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 an image for extracting the central line of the two-dimensional blood vessel, so as to obtain the contour line of the two-dimensional blood vessel.
2. The method of synthesizing a three-dimensional blood vessel according to claim 1, wherein the method of projecting each of the two-dimensional blood vessel centerlines into a three-dimensional space based on a photographing angle and a detection distance of each of the two-dimensional coronary angiography images, the method comprising:
using a heart as an origin of coordinates, and establishing a three-dimensional coordinate system;
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 two-dimensional contrast image of interest, wherein the coordinates of each point on the central 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 radial point R, obtaining points on the three-dimensional blood vessel center line from PR connecting lines, and sequentially connecting the points on the three-dimensional blood vessel center line to obtain the three-dimensional blood vessel center line.
3. A method of three-dimensional vascular synthesis according to claim 2, wherein the projection of points (x, y) into a three-dimensional space yields a series of three-dimensional coordinate points P, with coordinates (x ', y ', z '), comprising:
the points (x, y) are rotated around the y axis to obtain (x ', y ', z ') series points, and the specific formula is as follows:
rotating the (x ', y', z ') series points around the x axis to obtain a series of three-dimensional coordinate points P, wherein the coordinates are (x', y ', z'), and the specific formula is as follows:
4. a method of three-dimensional vascular synthesis according to claim 3, wherein the method of projecting a radiation 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 coordinates of R in the three-dimensional space are (a, b, c).
5. The method of synthesizing a three-dimensional blood vessel according to claim 4, wherein the connecting each of the three-dimensional coordinate points P with the radiation point R, obtaining a point on the three-dimensional blood vessel centerline from a PR connection, sequentially connecting the points on the three-dimensional blood vessel centerline, and obtaining the three-dimensional blood vessel centerline comprises:
correspondingly connecting a series of three-dimensional coordinate points P and radiation points R obtained from the same two-dimensional radiography image of interest to obtain a plurality of PR straight lines;
acquiring the minimum distance point between two PR lines at the same position of the blood vessel, wherein the minimum distance point is a point A and a point B respectively;
connecting the point A with the point B, and acquiring the midpoint of an AB line segment as a point on the central line of the three-dimensional blood vessel;
sequentially connecting the obtained points on the central line of the three-dimensional blood vessel to obtain the central line of the three-dimensional blood vessel.
6. The method of three-dimensional vascular synthesis according to claim 1, wherein the method of acquiring a vessel centerline and a three-dimensional vessel radius from the coronary two-dimensional contrast image comprises:
Acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line;
acquiring a two-dimensional vessel radius in each two-dimensional contrast image of interest according to the two-dimensional vessel contour line;
and acquiring the three-dimensional blood vessel radius according to the two-dimensional blood vessel radius.
7. The method of three-dimensional vascular synthesis according to claim 6, wherein the method of obtaining the three-dimensional vascular radius from the two-dimensional vascular radius comprises:
wherein R represents three-dimensional vessel radius, R 1 、r 2 、r n The two-dimensional vessel radii of the first, second, and nth two-dimensional contrast images of interest are represented, respectively.
8. The method of three-dimensional vascular synthesis according to claim 1, wherein the method of extracting a two-dimensional vascular centerline from each of the two-dimensional coronary angiography images comprises:
reading a two-dimensional coronary angiography image;
acquiring a vessel segment of interest;
picking up a start point, a seed point and an end point of the vessel segment of interest;
dividing two-dimensional contrast images between two adjacent points of a starting point, a seed point and an ending point respectively to obtain at least two local vessel region diagrams;
extracting at least one local vascular path line from each local vascular zone map;
Connecting the corresponding blood vessel local route lines on each local blood vessel region graph to obtain at least one blood vessel route line;
and selecting one vascular path line as the two-dimensional vascular center line.
9. The method of three-dimensional vascular synthesis according to claim 8, wherein the extracting at least one local vascular path line from each of the local vascular zone maps comprises:
performing image enhancement processing on the local vascular region map to obtain a rough vascular map with strong contrast;
and meshing the rough blood vessel graph, and extracting at least one local path line of the blood vessel along the direction from the starting point to the ending point.
10. The method for synthesizing a three-dimensional blood vessel according to claim 9, wherein the method for obtaining a rough blood vessel map with strong contrast by performing image enhancement processing on the local blood vessel region map comprises the following steps:
in each local vessel region graph, the vessel segment of interest is taken as a foreground, other regions are taken as a background, the foreground is strengthened, the background is weakened, and the rough vessel graph with strong contrast is obtained.
11. The method of three-dimensional vascular synthesis according to claim 10, wherein the step of meshing the rough vascular map to extract at least one local vascular path line along the direction from the start point to the end point comprises:
Grid dividing the rough blood vessel map;
searching a shortest time path of the intersection points on the starting point and the n grids on the periphery along the extending direction of the blood vessel from the starting point to the ending point as a second point, searching the shortest time path of the intersection points on the second point and the n grids on the periphery as a third point, and repeating the steps until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and according to the searching sequence, connecting lines from the starting point to the ending point in the extending direction of the blood vessel, and obtaining at least one local path line of the blood vessel.
12. The method of claim 11, wherein selecting one of the vessel path lines as the two-dimensional vessel centerline comprises:
summing the time taken from the start point to the end point for each vessel path line if the vessel path line is two or more;
the vessel path line at least when taken is taken as the two-dimensional vessel centerline.
13. The method of three-dimensional vascular synthesis according to claim 6, wherein the method of acquiring a two-dimensional vascular profile from the vascular centerline comprises:
Extracting a two-dimensional blood vessel center line according to the two-dimensional coronary angiography image;
obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line;
setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
According to said D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
gradually converging the preset contour line of the blood vessel towards the center 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 an image for extracting the central line of the two-dimensional blood vessel, so as to obtain the contour line of the two-dimensional blood vessel.
14. The method of claim 13, wherein the method of obtaining a straightened vessel image from the two-dimensional vessel centerline comprises:
straightening the two-dimensional blood vessel center line to obtain a blood vessel center line;
dividing a local vascular region map of the vascular center line into x units along the vascular extending direction from the starting point to the ending point of the vascular center line, wherein x is a positive integer;
correspondingly arranging the two-dimensional blood vessel center line of each unit along the blood vessel center line;
and the image after corresponding setting is the straightened blood vessel image.
15. The method of claim 14, wherein the step-wise approaching the preset contour line of the blood vessel to the straight line of the center of the blood vessel, and the step-wise approaching the contour line of the straightened blood vessel comprises the steps of:
dividing the vascular preset contour line into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on a preset contour line of each blood vessel;
respectively converging z points towards the central straight line of the blood vessel in a grading way along the direction perpendicular to the central straight line of the blood vessel to generate z converging points, wherein z is a positive integer;
setting RGB difference threshold to delta RGB Threshold value Comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line along the direction perpendicular to the blood vessel center straight line every time the blood vessel center straight line is close, wherein the difference value is less than or equal to delta RGB Threshold value When the blood vessel is closed, the closing point stops closing to the center line of the blood vessel;
acquiring the close points as contour points;
and connecting the contour points in sequence to form a smooth curve which is the contour line of the straightened blood vessel.
16. The method of synthesizing a three-dimensional blood vessel according to claim 1, wherein the method of 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 radius of the corresponding three-dimensional blood vessel at the point on the central line of each three-dimensional blood vessel to obtain a plurality of edge points, and sequentially connecting the edge points to obtain a polygon similar to a circle;
and sequentially connecting the points on two adjacent polygons according to the form of a right triangle to obtain the three-dimensional blood vessel.
17. A three-dimensional vascular synthesis system for use in the method of three-dimensional vascular synthesis of any of claims 1 to 16, 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 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 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 and the three-dimensional vessel center line transmitted by the three-dimensional vessel center line acquisition device, and acquiring the three-dimensional vessel radius according to the image information and the three-dimensional vessel center line;
The three-dimensional vascular synthesis device is used for receiving the three-dimensional vascular center line transmitted by the three-dimensional vascular center line acquisition device and receiving the three-dimensional vascular radius transmitted by the three-dimensional vascular radius acquisition device, and synthesizing a three-dimensional vascular according to the three-dimensional vascular center line and the three-dimensional vascular radius.
18. The three-dimensional vascular synthesis system of claim 17, wherein the three-dimensional vascular centerline acquisition device comprises: the two-dimensional blood vessel central line extraction structure is connected with the image reading device and the three-dimensional blood vessel central line acquisition structure;
the two-dimensional blood vessel center line extraction structure is used for receiving the coronary artery two-dimensional contrast images sent by the image reading device and extracting a two-dimensional blood vessel center line from each two-dimensional contrast image of interest;
the three-dimensional blood vessel center line acquisition structure is used for receiving the two-dimensional blood vessel center lines sent by the two-dimensional blood vessel center line extraction structure, receiving the two-dimensional blood vessel center lines sent by the image reading device, and projecting each two-dimensional blood vessel center line into a three-dimensional space according to the shooting angle of each coronary artery two-dimensional radiography image to synthesize the three-dimensional blood vessel center lines.
19. The three-dimensional vascular synthesis system of claim 18, wherein the two-dimensional vascular centerline extraction structure comprises: the central line extraction 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 extraction unit is connected with the image reading device and is used for extracting a blood vessel central line according to the two-dimensional coronary angiography 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 extraction unit;
the first blood vessel contour line unit is used for setting a blood vessel diameter threshold D on the straightened blood vessel image sent by the straightening unit Threshold value The method comprises the steps of carrying out a first treatment on the surface of the According to said D Threshold value Generating a blood vessel preset contour line on two sides of the blood vessel center straight line; gradually converging the preset contour line of the blood vessel towards the center straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
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 center line to obtain a blood vessel contour line.
20. A coronary artery analysis system, comprising: the three-dimensional vascular synthesis system of any of claims 17-19.
21. A computer storage medium, comprising a computer program which, when executed by a processor, implements the method of synthesizing a three-dimensional blood vessel according to any one of claims 1 to 16.
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 CN112419484A (en) | 2021-02-26 |
| CN112419484B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114391792B (en) * | 2021-09-13 | 2023-02-24 | 南京诺源医疗器械有限公司 | Tumor prediction method and device based on narrow-band imaging and imaging endoscope |
| CN116012522B (en) * | 2022-11-03 | 2023-11-14 | 成都寻创信息技术有限公司 | Three-dimensional imaging system for head, neck, jaw and face soft tissues, bones and blood vessels |
| CN116712157B (en) * | 2023-01-16 | 2024-03-08 | 天津市鹰泰利安康医疗科技有限责任公司 | Irreversible perforation system for heart vessel |
| CN117036640B (en) * | 2023-10-10 | 2023-12-19 | 杭州脉流科技有限公司 | Coronary artery blood vessel model reconstruction method, device, equipment and storage medium |
| CN117058328B (en) * | 2023-10-11 | 2024-01-09 | 杭州脉流科技有限公司 | Coronary vessel tree classification method, apparatus, storage medium and program product |
Citations (5)
| 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 |
| 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 (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1665130A4 (en) * | 2003-09-25 | 2009-11-18 | Paieon Inc | System and method for three-dimensional reconstruction of a tubular organ |
| 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 (5)
| 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 |
| 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 |
|---|
| 基于造影图像的冠状动脉血管三维表面重建方法;陈晓冬 等;工程图学学报(第03期);111-116 * |
| 基于造影图像的能量场反投影合成冠脉三维重建;丛伟建 等;中国科学:信息科学;第43卷(第08期);947-963 * |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022109903A1 (en) | 2022-06-02 |
| CN112419484A (en) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
| US11896416B2 (en) | Method for calculating coronary artery fractional flow reserve on basis of myocardial blood flow and CT images | |
| CN112164020B (en) | Method, device, analysis system and storage medium for accurately extracting blood vessel center line | |
| CN112116615B (en) | Method and device for acquiring blood vessel contour line according to blood vessel center line | |
| JP5129480B2 (en) | System for performing three-dimensional reconstruction of tubular organ and method for operating blood vessel imaging device | |
| JP2022169579A (en) | Diagnostically useful results in real time | |
| US9375191B2 (en) | Method and apparatus for determining three-dimensional reconstruction of an object | |
| CN112132949B (en) | Three-dimensional modeling method, device and system for blood vessel with narrow lesion region | |
| Vukicevic et al. | Three-dimensional reconstruction and NURBS-based structured meshing of coronary arteries from the conventional X-ray angiography projection images | |
| US11017531B2 (en) | Shell-constrained localization of vasculature | |
| CN112419462B (en) | Rendering synthesis method, system and storage medium of three-dimensional blood vessel | |
| WO2017148332A1 (en) | Perfusion analysis method and device | |
| CN111476791B (en) | Image processing method, image processing apparatus, and non-transitory computer readable medium | |
| WO2022109905A1 (en) | Method for adjusting blood vessel contours and centreline, and storage medium | |
| CN114596311B (en) | Blood vessel function evaluation method and blood vessel function evaluation device based on blood vessel image | |
| CN112446866A (en) | Blood flow parameter calculation method, device, equipment and storage medium | |
| JP2024515635A (en) | System and method for reconstructing 3D images from ultrasound and camera images - Patents.com | |
| CN112132882A (en) | Method and device for extracting blood vessel central line from coronary artery two-dimensional contrast image | |
| CN112419277B (en) | Three-dimensional vascular centerline synthesis method, system and storage medium | |
| CN112669449A (en) | CAG and IVUS accurate linkage analysis method and system based on 3D reconstruction technology | |
| CN112419279B (en) | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium | |
| CN111815584B (en) | Method and system for acquiring heart gravity center based on CT sequence image | |
| US10977792B2 (en) | Quantitative evaluation of time-varying data | |
| US20240070855A1 (en) | Method and System of Calculating the Cross-section Area and Included Angle of Three-dimensional Blood Vessel Branch | |
| CN112487342B (en) | Method, system and storage medium for accurately acquiring narrow lesion interval |
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 |