CN112419276B - Method for adjusting contour and center line of blood vessel and storage medium - Google Patents
Method for adjusting contour and center line of blood vessel and storage medium Download PDFInfo
- Publication number
- CN112419276B CN112419276B CN202011337588.1A CN202011337588A CN112419276B CN 112419276 B CN112419276 B CN 112419276B CN 202011337588 A CN202011337588 A CN 202011337588A CN 112419276 B CN112419276 B CN 112419276B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- point
- line
- vessel
- contour
- 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
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 190
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000003860 storage Methods 0.000 title claims abstract description 18
- 238000012937 correction Methods 0.000 claims abstract description 20
- 238000002586 coronary angiography Methods 0.000 claims abstract description 13
- 230000002792 vascular Effects 0.000 claims description 21
- 238000010586 diagram Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 210000004351 coronary vessel Anatomy 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 7
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 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
- 238000011156 evaluation Methods 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
- 238000004458 analytical method 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
- 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
- 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
- 238000002601 radiography Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000003313 weakening effect Effects 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
Abstract
The application provides a method for adjusting the contour and the central line of a blood vessel and a storage medium, wherein the method comprises the following steps: acquiring a blood vessel central line and a blood vessel contour line of a two-dimensional coronary angiography; selecting an adjustment region of the vessel contour line of interest; and adjusting the blood vessel contour line in the adjusting area to be close to or far from the blood vessel center line to obtain a blood vessel correction contour line. According to the method and the storage medium for adjusting the blood vessel contour and the central line, provided by the application, by introducing manual interaction and combining the influence of the blood vessel integral information on a single adjusting point, the automatic calculation result of the blood vessel contour and the central line can be corrected, and the accuracy of calculating the blood vessel contour and the central line is improved.
Description
Technical Field
The application relates to the technical field of coronary artery medicine, in particular to a method for adjusting the contour and the central line of a blood vessel 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 accurate blood vessel contour and a central line are obtained by attempting to judge image pixels through an algorithm, but in the actual process, the imaging equipment has a certain error in obtaining the image pixels, the distance between the single image pixels is also a fixed parameter, the image pixels can only reflect the whole condition of the blood vessel contour and the central line, the accuracy on details cannot be determined, and if the image pixels are simply relied on, the blood vessel contour and the central line can be greatly different from the actual condition on details.
Disclosure of Invention
The application provides a method for adjusting a blood vessel contour and a central line and a storage medium, which are used for solving the problem that accurate blood vessel contour and central line can not be obtained through judgment of image pixels by an algorithm.
To achieve the above object, in a first aspect, the present application provides a method for adjusting a contour of a blood vessel, comprising:
acquiring a blood vessel central line and a blood vessel contour line of a two-dimensional coronary angiography;
selecting an adjustment region of the vessel contour line of interest;
and adjusting the blood vessel contour line in the adjusting area to be close to or far from the blood vessel center line to obtain a blood vessel correction contour line.
Optionally, in the method for adjusting a blood vessel contour, the method for adjusting the blood vessel contour line in the adjusting area to be close to or far from the blood vessel center line to obtain the corrected blood vessel contour line includes:
obtaining a single adjustment step length m;
selecting any point in a 1/3-2/3 area in the adjustment area as an original contour point;
acquiring a point which corresponds to the original contour point and is positioned on the blood vessel central line as a corresponding central point;
obtaining a unit vector and Euclidean distance between the original contour point and the corresponding center point;
the original contour point moves a plurality of steps along the unit vector to obtain a contour adjusting point;
and in the adjustment area, picking up a starting point, a contour adjustment point, a plurality of seed points and an end point of the blood vessel contour line, and adopting a cubic B spline fitting curve algorithm to fit to obtain a blood vessel correction contour line.
Optionally, the method for adjusting the profile of a blood vessel further comprises: and correcting and adjusting the blood vessel correction contour line.
Optionally, the method for adjusting the profile of the blood vessel, the method for correcting and adjusting the profile line of the blood vessel comprises the following steps:
acquiring an included angle alpha of an original contour point, a corresponding center point and a contour adjusting point on the contour line of the blood vessel;
acquiring Euclidean distance L between the corresponding center point and the profile adjusting point;
if said alpha > 90 DEG or L < D min Acquiring a contour replacement point according to the unit vector of the original contour point and the corresponding center point, the original contour point coordinate, the minimum vessel diameter and the corresponding center point coordinate;
if said alpha is not more than 90 DEG and L is not less than D min Taking the profile adjusting point as a wheel buckle replacing point;
and in the adjustment area, picking up a starting point, a profile replacement point, a plurality of seed points and an end point of the blood vessel profile line, and adopting a cubic B spline fitting curve algorithm to fit to obtain a new blood vessel correction profile line.
Optionally, the method for adjusting the blood vessel profile, the method for acquiring the blood vessel centerline and the blood vessel profile of the coronary two-dimensional radiography comprises the following steps:
extracting a blood vessel center line according to the two-dimensional coronary angiography image;
obtaining a straightened blood vessel image according to the 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 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;
and 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.
Optionally, the method for adjusting the contour of a blood vessel, wherein the method for obtaining a straightened blood vessel image according to the blood vessel center line comprises the following steps:
straightening the central line of the blood vessel to obtain a central straight line of the blood vessel;
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 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 method for adjusting a contour of a blood vessel, the method for gradually converging the preset contour line of the blood vessel toward the center line of the blood vessel to obtain a 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 adjusting the blood vessel profile 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 blood vessel path line as the blood vessel central line.
Optionally, the method for adjusting the blood vessel profile includes:
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 method for adjusting a blood vessel profile, the method for performing image enhancement processing on the local blood vessel region map to obtain a rough blood vessel 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 method for adjusting a blood vessel profile, the method for meshing the rough blood vessel map and extracting at least one local path line of the blood vessel 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 method for adjusting a blood vessel profile, the method for selecting one of the blood vessel path lines as the blood vessel center line 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 which the minimum is taken as the vessel centerline.
In a second aspect, the present application provides a method of adjusting a vessel centerline comprising:
the above-described method of modulating a vascular profile;
and obtaining the central curves of the two blood vessel correction contour lines, namely the blood vessel correction central line.
In a third aspect, the present application provides a computer storage medium, which when executed by a processor implements the above-described method of adjusting a vessel profile.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
the application improves the automatic calculation result of the blood vessel contour and the central line, and combines the actual condition of the blood vessel by manually adjusting and driving the blood vessel contour and the central line, and dynamically corrects the blood vessel contour and the central line according to the experience of a user, thereby improving the calculation accuracy of the blood vessel contour and the central line.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
reference numerals are described below:
FIG. 1 is a flow chart of one embodiment of a method of the present application for adjusting a vessel profile;
FIG. 2 is a flow chart of S100 of the present application;
FIG. 3 is a flow chart of S110 of the present application;
FIG. 4 is a flowchart of S120 of the present application;
FIG. 5 is a flowchart of S150 of the present application;
FIG. 6 is a flowchart of S300 of the present application;
FIG. 7 is a flow chart of another embodiment of a method of the present application for adjusting a vessel profile;
fig. 8 is a flowchart of S400 of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments of the present application and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
Various embodiments of the application 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 application. That is, in some embodiments of the application, 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.
As shown in fig. 1, in order to solve the above problems, the present application provides a method for adjusting a profile of a blood vessel, comprising:
s100, as shown in FIG. 2, acquiring a vessel centerline and a vessel contour line of a two-dimensional coronary angiography comprises:
s110, as shown in fig. 3, extracting a vessel centerline from the two-dimensional coronary angiography image, including:
s111, reading a two-dimensional coronary angiography image;
s112, acquiring a blood vessel segment of interest;
s113, picking up a starting point, a seed point and an ending point of the blood vessel segment of interest;
s114, 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;
s115, extracting at least one blood vessel local path line from each local blood vessel region graph, wherein the step comprises the following steps:
in each local vessel region diagram, taking a vessel segment of interest as a foreground, taking other regions as a background, strengthening the foreground, weakening the background, and obtaining a rough vessel diagram with strong contrast;
grid division is carried out on the rough blood vessel map;
searching a shortest time path of the intersection points on the starting point and the n grids at 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 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;
according to the searching sequence, connecting lines from the starting point to the ending point in the extending direction of the blood vessel to obtain at least one local path line of the blood vessel;
s116, connecting the corresponding blood vessel local route lines on each local blood vessel region graph to obtain at least one blood vessel route line;
s117, selecting a blood vessel path line as a blood vessel center line, comprising:
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 smallest vessel path line is taken as the vessel centerline.
S120, as shown in FIG. 4, obtaining a straightened vessel image according to the vessel centerline, including:
s121, straightening the central line of the blood vessel to obtain a central straight line of the blood vessel;
s122, 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;
s123, correspondingly arranging the blood vessel center line of each unit along the blood vessel center line;
s124, the image after corresponding setting is a straightened blood vessel image.
S130, setting a blood vessel diameter threshold D on the straightened blood vessel image Threshold value ；
S140, according to D Threshold value Generating a blood vessel preset contour line on two sides of a blood vessel center straight line;
s150, as shown in FIG. 5, 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:
s151, dividing a blood vessel preset contour line into y units, wherein y is a positive integer;
s152, acquiring z points of each unit, which are positioned on a preset contour line of each blood vessel;
s153, 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;
s154, setting RGB difference value threshold as delta 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;
s155, acquiring a close point as a contour point;
and S156, connecting the contour points in sequence to form a smooth curve which is the contour line of the straightened blood vessel.
And S160, 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.
S200, selecting an adjustment area of a blood vessel contour line of interest;
s300, as shown in FIG. 6, the blood vessel contour line in the adjusting area is adjusted to be close to or far from the blood vessel center line, so as to obtain a blood vessel correction contour line, which comprises the following steps:
s310, obtaining a single adjustment step length m;
s320, selecting any point in a 1/3-2/3 area in the adjustment area as an original contour point;
s330, acquiring a point which corresponds to the original contour point and is positioned on the blood vessel central line as a corresponding central point;
s340, obtaining a unit vector and Euclidean distance between an original contour point and a corresponding center point;
s350, the original contour point moves a plurality of steps along the unit vector to obtain a contour adjusting point;
s360, picking up a starting point, a contour adjusting point, a plurality of seed points and an end point of a blood vessel contour line in the adjusting area, and fitting by adopting a cubic B spline fitting curve algorithm to obtain a blood vessel correction contour line.
Example 2:
as shown in fig. 7, on the basis of embodiment 1, further comprising: s400, correcting and adjusting the blood vessel correction contour line, as shown in FIG. 8, comprises:
s410, acquiring an included angle alpha of an original contour point, a corresponding center point and a contour adjusting point on a blood vessel contour line;
s420, acquiring Euclidean distance L between the corresponding center point and the profile adjusting point;
s430, if alpha > 90 DEG or L < D min Acquiring a contour replacement point according to the unit vector of the original contour point and the corresponding center point, the original contour point coordinate, the minimum vessel diameter and the corresponding center point coordinate;
s440, if alpha is not more than 90 DEG and L is not less than D min Taking the profile adjusting point as a wheel buckle replacing point;
s450, picking up a starting point, a profile replacement point, a plurality of seed points and an end point of a blood vessel profile line in the adjustment area, and adopting a cubic B spline fitting curve algorithm to fit to obtain a new blood vessel correction profile line.
The application provides a method for adjusting the center line of a blood vessel, which comprises the following steps:
the above-described method of modulating a vascular profile;
and obtaining the central curves of the two blood vessel correction contour lines, namely the blood vessel correction central line.
The present application provides a computer storage medium which when executed by a processor implements the method of adjusting a vessel profile described above.
Those skilled in the art will appreciate that the various aspects of the present application may be implemented as a system, method, or computer program product. Accordingly, aspects of the application 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 application 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 application 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 application could be implemented as a chip or circuit. As software, selected tasks according to embodiments of the application 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 application, 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 application 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 application have been described in some detail by way of illustration of the principles of the application, and it is to be understood that the application is not limited to the specific embodiments of the application but is intended to cover modifications, equivalents, alternatives and modifications within the spirit and principles of the application.
Claims (13)
1. A method of modulating a contour of a blood vessel, comprising:
acquiring a blood vessel central line and a blood vessel contour line of a two-dimensional coronary angiography;
selecting an adjustment region of the vessel contour line of interest;
the blood vessel contour line in the adjusting area is adjusted to be close to or far from the blood vessel center line, so that a blood vessel correction contour line is obtained;
the method for adjusting the blood vessel contour line in the adjusting area to be close to or far from the blood vessel center line to obtain the corrected blood vessel contour line comprises the following steps:
obtaining a single adjustment step length m;
selecting any point in a 1/3-2/3 area in the adjustment area as an original contour point;
acquiring a point which corresponds to the original contour point and is positioned on the blood vessel central line as a corresponding central point;
obtaining a unit vector and Euclidean distance between the original contour point and the corresponding center point;
the original contour point moves a plurality of steps along the unit vector to obtain a contour adjusting point;
and in the adjustment area, picking up a starting point, a contour adjustment point, a plurality of seed points and an end point of the blood vessel contour line, and adopting a cubic B spline fitting curve algorithm to fit to obtain a blood vessel correction contour line.
2. The method of modulating a vessel profile of claim 1, further comprising: and correcting and adjusting the blood vessel correction contour line.
3. A method of adjusting a vessel profile as claimed in claim 2, wherein the method of correcting and adjusting the vessel correction profile comprises:
acquiring an included angle alpha of an original contour point, a corresponding center point and a contour adjusting point on the contour line of the blood vessel;
acquiring Euclidean distance L between the corresponding center point and the profile adjusting point;
if said alpha > 90 DEG or L < D min Acquiring a contour replacement point according to the unit vector of the original contour point and the corresponding center point, the original contour point coordinate, the minimum vessel diameter and the corresponding center point coordinate;
if said alpha is not more than 90 DEG and L is not less than D min Taking the profile adjusting point as a wheel buckle replacing point;
and in the adjustment area, picking up a starting point, a profile replacement point, a plurality of seed points and an end point of the blood vessel profile line, and adopting a cubic B spline fitting curve algorithm to fit to obtain a new blood vessel correction profile line.
4. The method of claim 1, wherein the method of acquiring a vessel centerline and a vessel contour line of a two-dimensional contrast of a coronary artery comprises:
extracting a blood vessel center line according to the two-dimensional coronary angiography image;
obtaining a straightened blood vessel image according to the 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 the blood vessel center line;
gradually converging the preset contour line of the blood vessel towards the central line of the blood vessel to obtain the contour line of the straightened blood vessel;
and projecting the contour line of the straightened blood vessel back to the image of the central line of the blood vessel to obtain the contour line of the blood vessel.
5. The method of adjusting a vessel profile of claim 4, wherein the method of obtaining a straightened vessel image from the vessel centerline comprises:
straightening the central line of the blood vessel to obtain a central straight line of the blood vessel;
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;
correspondingly arranging the 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.
6. The method of claim 5, wherein the step-wise approaching the preset vascular contour line to the central line of the blood vessel, and the step-wise approaching the preset vascular contour line to the central line of the 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.
7. The method of claim 6, wherein the method of extracting a vessel centerline from a two-dimensional coronary angiography image 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 blood vessel path line as the blood vessel central line.
8. The method of claim 7, wherein the method of extracting at least one local vascular path line from the local vascular zone map, respectively, 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.
9. The method of claim 8, wherein the image enhancement processing is performed on the local vessel region map to obtain a rough vessel map with strong contrast, and the method comprises:
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.
10. The method of claim 9, wherein the meshing the rough vessel map to extract at least one vessel local 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.
11. The method of claim 10, wherein selecting one of the vessel path lines as the 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 which the minimum is taken as the vessel centerline.
12. A method of adjusting a vessel centerline, comprising:
a method of modulating a vascular profile as claimed in any one of claims 1 to 11;
and obtaining the central curves of the two blood vessel correction contour lines, namely the blood vessel correction central line.
13. A computer storage medium, characterized in that a computer program, when executed by a processor, implements the method of adjusting a vessel profile according to any one of claims 1-11.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337588.1A CN112419276B (en) | 2020-11-25 | 2020-11-25 | Method for adjusting contour and center line of blood vessel and storage medium |
| PCT/CN2020/131701 WO2022109905A1 (en) | 2020-11-25 | 2020-11-26 | Method for adjusting blood vessel contours and centreline, and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011337588.1A CN112419276B (en) | 2020-11-25 | 2020-11-25 | Method for adjusting contour and center line of blood vessel and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112419276A CN112419276A (en) | 2021-02-26 |
| CN112419276B true CN112419276B (en) | 2023-12-05 |
Family
ID=74843794
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011337588.1A Active CN112419276B (en) | 2020-11-25 | 2020-11-25 | Method for adjusting contour and center line of blood vessel and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112419276B (en) |
| WO (1) | WO2022109905A1 (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113256707B (en) * | 2021-06-25 | 2021-09-28 | 西南交通大学 | Rail center line extraction method based on rail top surface seed point automatic tracking |
| CN114898036A (en) * | 2022-04-26 | 2022-08-12 | 青岛埃米博创医疗科技有限公司 | Liver blood vessel model generation method based on exploration operator |
| CN115005846B (en) * | 2022-08-08 | 2022-10-11 | 江苏盛恩祥生物技术有限公司 | Data processing method and device suitable for blood vessel development |
| CN116758050A (en) * | 2023-07-12 | 2023-09-15 | 强联智创（北京）科技有限公司 | Method and product for blind completion of central line of intracranial Wills ring blood vessel |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102890823A (en) * | 2011-07-19 | 2013-01-23 | 株式会社东芝 | Motion object outline extraction and left ventricular image segmentation method and device |
| CN104318557A (en) * | 2014-10-17 | 2015-01-28 | 重庆大学 | Vascular skeleton line reconstruction and precise vessel diameter calculation method |
| CN109461138A (en) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | Calculation method of parameters, system, readable storage medium storing program for executing and computer equipment |
| CN110889896A (en) * | 2019-11-11 | 2020-03-17 | 苏州润迈德医疗科技有限公司 | Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111325761A (en) * | 2020-02-26 | 2020-06-23 | 医慧(嘉兴)科技有限公司 | Intravascular ultrasound image segmentation and measurement method based on cubic B-spline interpolation |
-
2020
- 2020-11-25 CN CN202011337588.1A patent/CN112419276B/en active Active
- 2020-11-26 WO PCT/CN2020/131701 patent/WO2022109905A1/en active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102890823A (en) * | 2011-07-19 | 2013-01-23 | 株式会社东芝 | Motion object outline extraction and left ventricular image segmentation method and device |
| CN104318557A (en) * | 2014-10-17 | 2015-01-28 | 重庆大学 | Vascular skeleton line reconstruction and precise vessel diameter calculation method |
| CN109461138A (en) * | 2018-09-29 | 2019-03-12 | 深圳睿心智能医疗科技有限公司 | Calculation method of parameters, system, readable storage medium storing program for executing and computer equipment |
| 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 |
|---|---|
| CN112419276A (en) | 2021-02-26 |
| WO2022109905A1 (en) | 2022-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112419276B (en) | Method for adjusting contour and center line of blood vessel and storage medium | |
| US10134129B2 (en) | Method and system for hemodynamic computation in coronary arteries | |
| CN112164020B (en) | Method, device, analysis system and storage medium for accurately extracting blood vessel center line | |
| JP2022169579A (en) | Diagnostically useful results in real time | |
| JP5129480B2 (en) | System for performing three-dimensional reconstruction of tubular organ and method for operating blood vessel imaging device | |
| US9424637B2 (en) | Vessel segmentation method and apparatus using multiple thresholds values | |
| US9842401B2 (en) | Segmentation apparatus for interactively segmenting blood vessels in angiographic image data | |
| US20220277447A1 (en) | Method and apparatus for acquiring contour line of blood vessel according to centerline of blood vessel | |
| JP7313284B2 (en) | Customize, calibrate and/or train fractional coronary flow reserve simulation parameters | |
| CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
| CN107580470B (en) | Instantaneous flow reserve-computed tomography | |
| WO2017187269A1 (en) | Enhanced personalized evaluation of coronary artery disease using an integration of multiple medical imaging techniques | |
| CN112132949A (en) | Blood vessel three-dimensional modeling method, device and system with narrow lesion interval | |
| CN112446866A (en) | Blood flow parameter calculation method, device, equipment and storage medium | |
| CN112494020B (en) | Method for acquiring aortic pressure curve of interest and storage medium | |
| WO2022109906A1 (en) | Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium | |
| JP2022544052A (en) | Method and apparatus for correcting blood flow velocity based on contrast-enhanced image interval time | |
| CN112132882A (en) | Method and device for extracting blood vessel central line from coronary artery two-dimensional contrast image | |
| CN112472112B (en) | Method, system and storage medium for regulating vascular stenosis | |
| CN112419277B (en) | Three-dimensional vascular centerline synthesis method, system and storage medium | |
| US11481957B2 (en) | Medical image processing apparatus and storage medium | |
| Cunha et al. | A method for segmentation of dental implants and crestal bone | |
| CN112419279B (en) | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium | |
| CN110929604B (en) | Method, device, system and storage medium for screening flow velocity based on contrast image | |
| CN110494893B (en) | FFR-based interactive monitoring of non-invasive imaging |
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 |