CN112419279B - Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium - Google Patents
Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium Download PDFInfo
- Publication number
- CN112419279B CN112419279B CN202011342206.4A CN202011342206A CN112419279B CN 112419279 B CN112419279 B CN 112419279B CN 202011342206 A CN202011342206 A CN 202011342206A CN 112419279 B CN112419279 B CN 112419279B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- blood vessel
- vessel
- image
- line
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 210000004204 blood vessel Anatomy 0.000 title claims description 120
- 230000002194 synthesizing effect Effects 0.000 title claims description 10
- 230000002792 vascular Effects 0.000 claims abstract description 41
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 14
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 14
- 238000002586 coronary angiography Methods 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 11
- 238000010586 diagram Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000002285 radioactive effect Effects 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 3
- 239000008280 blood Substances 0.000 claims description 2
- 210000004369 blood Anatomy 0.000 claims description 2
- 238000001308 synthesis method Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 7
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 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
- 210000000056 organ Anatomy 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 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
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000000605 extraction Methods 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
- 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
- 238000002601 radiography Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- 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)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Quality & Reliability (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The application provides a two-dimensional image selection and three-dimensional vascular synthesis method and a storage medium, wherein the method comprises the following steps: loading a plurality of groups of coronary artery two-dimensional contrast image sequence groups; selecting a set of said two-dimensional contrast image series of coronary arteries; displaying a series of thumbnails in each group of the coronary artery two-dimensional contrast image sequence groups; selecting said thumbnail of interest from within at least two of said sets of coronary two-dimensional contrast image sequences, said thumbnail showing a magnified image, i.e. obtaining a coronary two-dimensional contrast image for three-dimensional vascular synthesis. The method and the storage medium realize two-dimensional image selection and three-dimensional vascular synthesis, can effectively reduce the difficulty and error of vascular information acquisition, improve the accuracy of vascular information acquisition, and have better practicability.
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a two-dimensional image selection and three-dimensional vascular synthesis method 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.
In the prior art, an image processing algorithm is generally applied to a two-dimensional contrast sequence image, a full-automatic, semi-automatic or manual method is adopted to calculate the path of a blood vessel, and an edge detection method is adopted to obtain the outline of the blood vessel, so that the diameter of the blood vessel is further calculated. However, this method has problems in that various organs around the blood vessel overlap and cross, which all cause uncertainty to the algorithm.
Disclosure of Invention
The invention provides a two-dimensional image selection and three-dimensional vascular synthesis method and a storage medium, which are used for avoiding uncertain factors caused by overlapping and crossing of various organs around blood vessels in two-dimensional contrast sequence images.
In order to achieve the above-mentioned object,
in a first aspect, the present application provides a method of synthesizing a three-dimensional blood vessel, comprising:
loading a plurality of groups of coronary artery two-dimensional contrast image sequence groups;
selecting a set of said two-dimensional contrast image series of coronary arteries;
displaying a series of thumbnails in each group of the coronary artery two-dimensional contrast image sequence groups;
selecting the thumbnail of interest from at least two of the two-dimensional coronary angiography image sequence sets, wherein the thumbnail displays a magnified image, namely, a two-dimensional coronary angiography image for three-dimensional vascular synthesis is obtained;
when each group of coronary artery two-dimensional contrast image sequence groups is displayed on a page, selecting an intermediate frame as a thumbnail by default for display;
acquiring image information of at least two coronary two-dimensional contrast images with different shooting angles and interested coronary two-dimensional contrast images, wherein the image information comprises shooting angles and shooting distances;
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;
synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius;
the acquiring the three-dimensional blood vessel radius according to the image information of the two-dimensional coronary angiography image comprises the following steps:
extracting a two-dimensional vessel centerline from each of the coronary two-dimensional contrast images of interest;
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 of the two-dimensional coronary angiography images of interest according to the two-dimensional vessel contour lines;
the three-dimensional blood vessel radius is obtained according to the two-dimensional blood vessel radius, and the specific formula is as follows:
wherein R represents three-dimensional vessel radius, R 1 、r 2 、r n Representing the two-dimensional vessel radii of the first, second and nth two-dimensional contrast images of interest, respectively;
the acquiring the two-dimensional blood vessel contour line according to the two-dimensional blood vessel 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;
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.
Preferably, the method for acquiring the three-dimensional blood vessel center line according to the image information of the two-dimensional coronary angiography image comprises the following steps:
said extracting a two-dimensional vessel centerline from each of said coronary two-dimensional contrast images of interest;
projecting a radioactive source into a three-dimensional space to form a radioactive point;
the two-dimensional blood vessel center line is projected into a three-dimensional space;
all points in the three-dimensional space are connected with the radiation points, so that a series of crossing points are generated;
and connecting the intersecting points in turn to obtain the three-dimensional blood vessel center line.
Preferably, the method for extracting a two-dimensional vessel centerline from each of the two-dimensional angiographic images of coronary arteries of interest 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.
Preferably, the method for extracting at least one blood vessel local path line 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.
Preferably, the method for meshing the rough vessel map and extracting at least one vessel local path line along the direction from the starting point to the ending 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.
Preferably, the method for selecting one of the vascular path lines 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.
Preferably, 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 computer storage medium containing a computer program for implementing a method of three-dimensional vascular synthesis as claimed in any one of claims 1 to 7 when the computer program is executed by a processor.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the three-dimensional blood vessel synthesizing method has the advantages that in the process of calculating the blood vessel information through the two-dimensional image of the coronary artery, the influence of other parts on the blood vessel is reduced, the radius of the blood vessel can be obtained more easily, the accuracy of the radius of the blood vessel is improved, and the robustness of an algorithm is improved.
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 one embodiment of a method of coronary two-dimensional image selection for three-dimensional vascular synthesis of the present application;
FIG. 2 is a flow chart of another embodiment of a method of coronary two-dimensional image selection for three-dimensional vascular synthesis of the present application;
FIG. 3 is a flow chart of S600 of the present application;
FIG. 4 is a flow chart of S620 of the present application;
FIG. 5 is a flow chart of S630 of the present application;
fig. 6 is a flowchart of S700 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.
In the prior art, the blood vessel contour line is often required to be extracted by calculating the blood vessel evaluation parameters through a blood vessel three-dimensional model, and the blood vessel contour is particularly difficult to extract due to the problems of curling and unclear edges of the blood vessel, and the calculation data is huge and tedious, so how to rapidly extract the blood vessel contour line and the accuracy of extraction are always problems which need to be solved by technicians.
Example 1:
as shown in fig. 1, to solve the above problem, the present application provides a method for selecting a two-dimensional image of a coronary artery for three-dimensional vascular synthesis, including:
s100, loading a plurality of coronary artery two-dimensional contrast image sequence groups;
s200, selecting a group of two-dimensional coronary angiography image sequence groups;
s300, displaying a series of thumbnails in each coronary artery two-dimensional contrast image sequence group;
s400, selecting the thumbnail image of interest from at least two coronary artery two-dimensional contrast image sequence groups, wherein the thumbnail image displays a magnified image, namely, a coronary artery two-dimensional contrast image for three-dimensional vascular synthesis is obtained.
In one embodiment of the present application, when each of the two-dimensional coronary angiography image sequence sets is displayed on a page, an intermediate frame is selected as a thumbnail by default.
Example 2:
as shown in fig. 2, the present application provides a method for synthesizing a three-dimensional blood vessel, including:
s500, according to the method of the embodiment 1, obtaining image information of at least two coronary two-dimensional contrast images of interest with different shooting angles, wherein the image information comprises shooting angles and shooting distances;
s600, as shown in fig. 3, acquiring a three-dimensional blood vessel center line and a three-dimensional blood vessel radius according to the image information of the two-dimensional coronary angiography image, including:
s610, extracting a two-dimensional blood vessel center line from each of the two-dimensional coronary angiography images of interest, including:
s611, reading a coronary artery two-dimensional contrast image;
s612, acquiring a blood vessel segment of interest;
s613, picking up a starting point, a seed point and an ending point of the blood vessel segment of interest;
s614, 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;
s615, extracting at least one local vascular path line from each local vascular zone map, including:
performing image enhancement processing on the local vascular region map to obtain a rough vascular map with strong contrast, wherein the method 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.
Meshing the rough vessel map, extracting at least one vessel local path line along the direction from the starting point to the ending point, and comprising the following steps:
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.
S616, selecting one of the vessel path lines as the two-dimensional vessel center line, including:
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.
S620, as shown in FIG. 4, according to the shooting angle of each coronary artery two-dimensional radiography image, each two-dimensional blood vessel center line is projected into a three-dimensional space, and the three-dimensional blood vessel center line is synthesized, including:
s621, projecting a radioactive source into the three-dimensional space to form a radioactive point;
s622, projecting the two-dimensional blood vessel center line into a three-dimensional space;
s623, connecting all points in the three-dimensional space with the radiation points to generate a series of crossing points;
s624, connecting the intersecting points in sequence to obtain the three-dimensional blood vessel center line;
s630, as shown in FIG. 5, acquiring a two-dimensional blood vessel contour line according to the two-dimensional blood vessel center line comprises:
s631, extracting a two-dimensional blood vessel center line according to the two-dimensional coronary angiography image;
s632, obtaining a straightened blood vessel image according to the two-dimensional blood vessel center line, comprising:
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.
S633, setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
S634, according to said D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
s635, gradually closing the preset contour line of the blood vessel to the center straight line of the blood vessel, and obtaining the contour line of the straightened blood vessel, wherein the step comprises the following steps:
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.
And S636, projecting the contour line of the straightened blood vessel back to the 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.
S640, acquiring a two-dimensional vessel radius in each of the two-dimensional coronary angiography images of interest according to the two-dimensional vessel contour lines;
s650, acquiring the 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.
S700, as shown in fig. 6, synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel centerline and the three-dimensional blood vessel radius, including:
s710, 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;
s720, connecting the points on two adjacent polygons in sequence according to the form of a right triangle, and obtaining the three-dimensional blood vessel.
The present application provides a computer storage medium, which when executed by a processor implements the method for coronary artery two-dimensional image selection for three-dimensional vascular synthesis described above.
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 (8)
1. A method of three-dimensional vascular synthesis, comprising:
loading a plurality of groups of coronary artery two-dimensional contrast image sequence groups;
selecting a set of said two-dimensional contrast image series of coronary arteries;
displaying a series of thumbnails in each group of the coronary artery two-dimensional contrast image sequence groups;
selecting the thumbnail of interest from at least two of the two-dimensional coronary angiography image sequence sets, wherein the thumbnail displays a magnified image, namely, a two-dimensional coronary angiography image for three-dimensional vascular synthesis is obtained;
when each group of coronary artery two-dimensional contrast image sequence groups is displayed on a page, selecting an intermediate frame as a thumbnail by default for display;
acquiring image information of at least two coronary two-dimensional contrast images with different shooting angles and interested coronary two-dimensional contrast images, wherein the image information comprises shooting angles and shooting distances;
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;
synthesizing a three-dimensional blood vessel according to the three-dimensional blood vessel center line and the three-dimensional blood vessel radius;
the acquiring the three-dimensional blood vessel radius according to the image information of the two-dimensional coronary angiography image comprises the following steps:
extracting a two-dimensional vessel centerline from each of the coronary two-dimensional contrast images of interest;
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 of the two-dimensional coronary angiography images of interest according to the two-dimensional vessel contour lines;
the three-dimensional blood vessel radius is obtained according to the two-dimensional blood vessel radius, and the specific formula is as follows:
wherein R represents three-dimensional vessel radius, R 1 、r 2 、r n Representing the two-dimensional vessel radii of the first, second and nth two-dimensional contrast images of interest, respectively;
the acquiring the two-dimensional blood vessel contour line according to the two-dimensional blood vessel 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;
setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
According to said D Threshold value Straight at the center of the blood vesselGenerating a blood vessel preset contour line at two sides of the 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 claim 1, wherein the method for acquiring the three-dimensional vessel centerline from the image information of the two-dimensional coronary angiography image comprises:
extracting a two-dimensional vessel centerline from each of the coronary two-dimensional contrast images of interest;
projecting a radioactive source into a three-dimensional space to form a radioactive point;
the two-dimensional blood vessel center line is projected into a three-dimensional space;
all points in the three-dimensional space are connected with the radiation points, so that a series of crossing points are generated;
and connecting the intersecting points in turn to obtain the three-dimensional blood vessel center line.
3. A method of three-dimensional vascular synthesis according to claim 2, wherein the method of extracting a two-dimensional vascular centerline from each of the two-dimensional coronary angiography images of interest 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.
4. A method of three-dimensional vascular synthesis according to claim 3, wherein the method of extracting at least one local vascular path 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.
5. The method of claim 4, wherein the step of meshing the rough vessel map to extract at least one local vessel 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.
6. The method of claim 5, 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.
7. The method of claim 6, wherein the step of synthesizing a three-dimensional vessel based on the three-dimensional vessel centerline and the three-dimensional 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.
8. A computer storage medium, comprising a computer program which, when executed by a processor, is adapted to carry out a method of three-dimensional vascular synthesis as claimed in any one of claims 1 to 7.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011342206.4A CN112419279B (en) | 2020-11-25 | 2020-11-25 | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium |
| PCT/CN2020/131700 WO2022109904A1 (en) | 2020-11-25 | 2020-11-26 | Method for two-dimensional image selection and three-dimensional blood vessel synthesis and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011342206.4A CN112419279B (en) | 2020-11-25 | 2020-11-25 | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112419279A CN112419279A (en) | 2021-02-26 |
| CN112419279B true CN112419279B (en) | 2024-02-23 |
Family
ID=74842431
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011342206.4A Active CN112419279B (en) | 2020-11-25 | 2020-11-25 | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112419279B (en) |
| WO (1) | WO2022109904A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117058328B (en) * | 2023-10-11 | 2024-01-09 | 杭州脉流科技有限公司 | Coronary vessel tree classification method, apparatus, storage medium and program product |
Citations (4)
| 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 |
| CN108242075A (en) * | 2018-01-05 | 2018-07-03 | 苏州润迈德医疗科技有限公司 | A kind of multi-angle reconstructing blood vessel method based on X ray coronary angiography image |
| 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 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9934566B2 (en) * | 2015-07-14 | 2018-04-03 | Siemens Healthcare Gmbh | 3-D vessel tree surface reconstruction method |
| CN108133478B (en) * | 2018-01-11 | 2021-11-23 | 苏州润迈德医疗科技有限公司 | Method for extracting coronary artery blood vessel central line |
| CN109461138B (en) * | 2018-09-29 | 2020-10-27 | 深圳睿心智能医疗科技有限公司 | Parameter calculation method, system, readable storage medium and computer equipment |
| 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 |
| CN111754506B (en) * | 2020-07-01 | 2024-02-06 | 杭州脉流科技有限公司 | Coronary artery stenosis rate calculation method, device and system based on intra-cavity image and computer storage medium |
-
2020
- 2020-11-25 CN CN202011342206.4A patent/CN112419279B/en active Active
- 2020-11-26 WO PCT/CN2020/131700 patent/WO2022109904A1/en active Application Filing
Patent Citations (4)
| 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 |
| CN108242075A (en) * | 2018-01-05 | 2018-07-03 | 苏州润迈德医疗科技有限公司 | A kind of multi-angle reconstructing blood vessel method based on X ray coronary angiography image |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022109904A1 (en) | 2022-06-02 |
| CN112419279A (en) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112116615B (en) | Method and device for acquiring blood vessel contour line according to blood vessel center line | |
| 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 | |
| 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 | |
| JP5129480B2 (en) | System for performing three-dimensional reconstruction of tubular organ and method for operating blood vessel imaging device | |
| US11076770B2 (en) | System for vascular assessment | |
| JP2004243117A (en) | Method for obtaining physical parameters of physiological structure | |
| WO2022109906A1 (en) | Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium | |
| WO2022109905A1 (en) | Method for adjusting blood vessel contours and centreline, and storage medium | |
| Van Hamersvelt et al. | Diagnostic performance of on-site coronary CT angiography–derived fractional flow reserve based on patient-specific lumped parameter models | |
| WO2021097821A1 (en) | Method and device for extracting vascular centerline from two-dimensional coronary angiographic image | |
| CN112419279B (en) | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium | |
| CN112419277B (en) | Three-dimensional vascular centerline synthesis method, system and storage medium | |
| CN111815584B (en) | Method and system for acquiring heart gravity center based on CT sequence image | |
| EP3886702B1 (en) | Most relevant x-ray image selection for hemodynamic simulation | |
| US10977792B2 (en) | Quantitative evaluation of time-varying data | |
| JP7408381B2 (en) | Image processing device, program and method |
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 |