CN112419276A - Method for regulating blood vessel contour and central line and storage medium - Google Patents
Method for regulating blood vessel contour and central line and storage medium Download PDFInfo
- Publication number
- CN112419276A CN112419276A CN202011337588.1A CN202011337588A CN112419276A CN 112419276 A CN112419276 A CN 112419276A CN 202011337588 A CN202011337588 A CN 202011337588A CN 112419276 A CN112419276 A CN 112419276A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- contour
- point
- line
- adjusting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 264
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000001105 regulatory effect Effects 0.000 title description 4
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 22
- 238000002601 radiography Methods 0.000 claims abstract description 3
- 238000012937 correction Methods 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 11
- 239000013598 vector Substances 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000003993 interaction Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 201000000057 Coronary Stenosis Diseases 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 206010011089 Coronary artery stenosis Diseases 0.000 description 1
- 206010020565 Hyperaemia Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000002586 coronary angiography Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 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
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Images
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 method for adjusting the contour and the central line of a blood vessel and a storage medium, comprising the following steps: obtaining a blood vessel central line and a blood vessel contour line of coronary artery two-dimensional radiography; selecting an adjustment area of the interested blood vessel contour line; and adjusting the contour line of the blood vessel in the adjusting area to be close to or far from the center line of the blood vessel to obtain a corrected contour line of the blood vessel. According to the method for adjusting the blood vessel contour and the center line and the storage medium, manual interaction is introduced, the influence of the whole blood vessel information on a single adjusting point is combined, the automatic calculation result of the blood vessel contour and the center line can be corrected, and the accuracy of calculation of the blood vessel contour and the center line is improved.
Description
Technical Field
The present invention relates to the technical field of coronary artery medicine, and particularly to a method for adjusting a blood vessel contour and a centerline and a storage medium.
Background
The deposition of lipids and carbohydrates in human blood on the vessel wall will form plaques on the vessel wall, which in turn leads to vessel stenosis; especially, the blood vessel stenosis near the coronary artery of the heart can cause insufficient blood supply of cardiac muscle, induce diseases such as coronary heart disease, angina pectoris and the like, and cause serious threat to the health of human beings. According to statistics, about 1100 million patients with coronary heart disease in China currently have the number of patients treated by cardiovascular interventional surgery increased by more than 10% every year.
Although conventional medical detection means such as coronary angiography CAG and computed tomography CT can display the severity of coronary stenosis of the heart, the ischemia of the coronary cannot be accurately evaluated. In order to improve the accuracy of coronary artery function evaluation, Pijls in 1993 proposes a new index for estimating coronary artery function through pressure measurement, namely Fractional Flow Reserve (FFR), and the FFR becomes the gold standard for coronary artery stenosis function evaluation through long-term basic and clinical research.
The Fractional Flow Reserve (FFR) generally refers to the fractional flow reserve of myocardium, and is defined as the ratio of the maximum blood flow provided by a diseased coronary artery to the maximum blood flow when the coronary artery is completely normal. Namely, the FFR value can be measured and calculated by measuring the pressure at the position of the coronary stenosis and the pressure at the position of the coronary stenosis under the maximal hyperemia state of the coronary artery through a pressure sensor.
In the prior art, an accurate blood vessel contour and a central line are obtained by judging image pixels through an algorithm, but in the actual process, certain errors exist in the acquisition of the image pixels by an imaging device, 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 precision of details cannot be determined, and if the image pixels are simply relied on, the details of the blood vessel contour and the central line are greatly different from the actual condition.
Disclosure of Invention
The invention 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 the accurate blood vessel contour and central line cannot be obtained by judging image pixels through an algorithm.
To achieve the above object, in a first aspect, the present application provides a method for adjusting a blood vessel contour, comprising:
obtaining a blood vessel central line and a blood vessel contour line of coronary artery two-dimensional radiography;
selecting an adjustment area of the interested blood vessel contour line;
and adjusting the contour line of the blood vessel in the adjusting area to be close to or far from the center line of the blood vessel to obtain a corrected contour line of the blood vessel.
Optionally, in the method for adjusting a blood vessel contour, the method for adjusting the blood vessel contour line in the adjustment region to approach or separate from the blood vessel centerline to obtain a corrected blood vessel contour line includes:
acquiring a single adjustment step length m;
selecting any point in 1/3-2/3 areas in the adjusting area as an original contour point;
acquiring points corresponding to the original contour points and located on the center line of the blood vessel as corresponding center points;
acquiring unit vectors and Euclidean distances between the original contour points and the corresponding central points;
moving the original contour point along the unit vector by a plurality of steps to obtain a contour adjusting point;
and in the adjusting area, picking up the starting point, the contour adjusting point, the plurality of seed points and the tail point of the blood vessel contour line, and fitting by adopting a cubic B-spline fitting curve algorithm to obtain a blood vessel correction contour line.
Optionally, the method for adjusting a blood vessel contour described above further includes: and correcting and adjusting the blood vessel correction contour line.
Optionally, in the method for adjusting a blood vessel contour, the method for correcting and adjusting the blood vessel modified contour line includes:
acquiring an included angle alpha of an original contour point, a corresponding central point and a contour adjusting point on the blood vessel contour line;
acquiring an Euclidean distance L between the corresponding central point and the contour adjusting point;
if said alpha is > 90 DEG or L < DminAcquiring a contour replacement point according to the unit vector of the original contour point and the corresponding central point, the coordinates of the original contour point, the minimum vessel diameter of the blood vessel and the coordinates of the corresponding central point;
if alpha is less than or equal to 90 degrees and L is more than or equal to DminTaking the contour adjusting point as a wheel buckle replacing point;
and in the adjusting area, picking up the starting point, the contour replacing point, the plurality of seed points and the tail point of the blood vessel contour line, and fitting by adopting a cubic B-spline fitting curve algorithm to obtain a new blood vessel correction contour line.
Optionally, the method for adjusting the vessel contour described above, the method for acquiring the vessel centerline and the vessel contour line of the coronary artery two-dimensional contrast includes:
extracting a blood vessel central line according to the coronary artery two-dimensional contrast image;
obtaining a straightened blood vessel image according to the blood vessel central line;
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line;
gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and projecting the contour line of the straightened blood vessel back to the image of the center line of the blood vessel to obtain the contour line of the blood vessel.
Optionally, in the method for adjusting a blood vessel contour, the method for obtaining a straightened blood vessel image according to the blood vessel centerline includes:
straightening the central line of the blood vessel to obtain a central straight line of the blood vessel;
dividing the local blood vessel region map into x units along the extending direction of the blood vessel from the starting point to the ending point, wherein x is a positive integer;
correspondingly arranging the center line of the blood vessel of each unit along the center line of the blood vessel;
and the correspondingly set image is the straightened blood vessel image.
Optionally, in the method for adjusting a blood vessel contour, the method for gradually drawing the preset contour line of the blood vessel toward the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel includes:
dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on each preset blood vessel contour line;
respectively closing the z points to the blood vessel center straight line in a grading way along the direction vertical to the blood vessel center straight line to generate z closing points, wherein z is a positive integer;
setting RGB difference threshold to delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)When the blood vessel is in the closed state, the closing point stops closing towards the center line of the blood vessel;
acquiring the approach point as a contour point;
and sequentially connecting the contour points to form a smooth curve which is the contour line of the straightened blood vessel.
Optionally, the method for adjusting a blood vessel contour as described above, the method for extracting a blood vessel centerline from a coronary artery two-dimensional contrast image, includes:
reading a coronary artery two-dimensional contrast image;
obtaining a vessel segment of interest;
picking up a starting point, a seed point and an end point of the vessel segment of interest;
respectively segmenting two-dimensional contrast images between two adjacent points of a starting point, a seed point and an end point to obtain at least two local blood vessel region images;
extracting at least one blood vessel local path line from each local blood vessel region map;
connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
and selecting one blood vessel path line as the blood vessel central line.
Optionally, the method for adjusting a blood vessel contour described above, the method for extracting at least one blood vessel local path line from the local blood vessel region map respectively includes:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
and performing grid division on the rough blood vessel map, and extracting at least one blood vessel local path line along the direction from the starting point to the end point.
Optionally, in the method for adjusting a blood vessel contour, the method for performing image enhancement processing on the local blood vessel region map to obtain a rough blood vessel map with a strong contrast includes:
in each local blood vessel region image, the blood vessel section of interest is used as a foreground, other regions are used as backgrounds, the foreground is strengthened, the backgrounds are weakened, and the rough blood vessel image with strong contrast is obtained.
Optionally, in the method for adjusting a blood vessel contour, the method for meshing the rough blood vessel map and extracting at least one blood vessel local path line along a direction from the starting point to the ending point includes:
gridding the rough vessel map;
searching the shortest time path between the starting point and the intersection points on the peripheral n grids along the extending direction of the blood vessels from the starting point to the ending point to serve as a second point, searching the shortest time path between the second point and the intersection points on the peripheral n grids to serve as a third point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and connecting the extending directions of the blood vessels from the starting point to the ending point according to the searching sequence to obtain at least one blood vessel local path line.
Optionally, in the method for adjusting a blood vessel contour, the selecting one of the blood vessel route lines as the blood vessel center line includes:
if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
the vessel path line that is the least in time is taken as the vessel centerline.
In a second aspect, the present application provides a method of adjusting a centerline of a vessel, comprising:
the method for regulating the blood vessel contour;
and acquiring central curves of the two blood vessel correction contour lines, namely the blood vessel correction central lines.
In a third aspect, the present application provides a computer storage medium, a computer program being executed by a processor for implementing the above-mentioned method of adjusting a blood vessel contour.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
according to the method and the device, the automatic calculation results of the blood vessel contour and the center line are improved, the blood vessel contour and the center line are manually adjusted and driven, the actual condition of the blood vessel is combined, dynamic correction is performed according to the experience of a user, and the calculation accuracy of the blood vessel contour and the center line 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 not to limit the invention. In the drawings:
the following reference numerals are used for the description:
FIG. 1 is a flow chart of one embodiment of a method of adjusting a blood vessel contour of the present application;
fig. 2 is a flowchart of S100 of the present application;
fig. 3 is a flowchart of S110 of the present application;
fig. 4 is a flowchart of S120 of the present application;
FIG. 5 is a flow chart 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 adjusting a blood vessel profile of the present application;
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 invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the specific embodiments of the present invention and the accompanying drawings. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the following description, for purposes of explanation, numerous implementation details are set forth in order to provide a thorough understanding of the various embodiments of the present invention. It should be understood, however, that these implementation details are not to be interpreted as limiting the invention. That is, in some embodiments of the invention, such implementation details are not necessary. In addition, some conventional structures and components are shown in simplified schematic form in the drawings.
As shown in fig. 1, to solve the above problem, the present application provides a method of adjusting a blood vessel contour, including:
s100, as shown in fig. 2, acquiring a vessel centerline and a vessel contour line of a coronary artery two-dimensional angiography, including:
s110, as shown in fig. 3, extracting a vessel centerline from the coronary artery two-dimensional contrast image includes:
s111, reading a coronary artery two-dimensional contrast image;
s112, obtaining the interested blood vessel section;
s113, picking up a starting point, a seed point and an end point of the interested blood vessel section;
s114, segmenting the two-dimensional contrast image between two adjacent points of the starting point, the seed point and the ending point respectively to obtain at least two local blood vessel region images;
s115, extracting at least one blood vessel local path line from each local blood vessel region map, including:
in each local blood vessel region image, an interested blood vessel section is used as a foreground, other regions are used as backgrounds, the foreground is strengthened, the backgrounds are weakened, and a rough blood vessel image with strong contrast is obtained;
meshing the rough blood vessel map;
searching the shortest time path between the starting point and the intersection points on the peripheral n grids as a second point and searching the shortest time path between the second point and the intersection points on the peripheral n grids as a third point along the extending direction of the blood vessel from the starting point to the ending point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
according to the search sequence, connecting the extending directions of the blood vessels from the starting point to the end point to obtain at least one blood vessel local path line;
s116, connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
s117, selecting a blood vessel path line as a blood vessel central line, including:
if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
the least vascular path line when used 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 center line of the blood vessel to obtain a center straight line of the blood vessel;
s122, dividing the local blood vessel region map into x units along the extending direction of the blood vessel from the starting point to the ending point, wherein x is a positive integer;
s123, correspondingly arranging the center line of the blood vessel of each unit along the center straight line of the blood vessel;
and S124, the correspondingly set image is a straightened blood vessel image.
S130, setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
S140, according to DThreshold(s)Generating preset contour lines of the blood vessels on two sides of the central straight line of the blood vessel;
s150, as shown in fig. 5, gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel, including:
s151, dividing the preset contour line of the blood vessel 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, along the direction perpendicular to the blood vessel center straight line, closing z points to the blood vessel center straight line in a grading mode respectively to generate z closing points, wherein z is a positive integer;
s154, setting the RGB difference value threshold value to delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)When the blood vessel is closed, the closing point stops closing towards the center line of the blood vessel;
s155, acquiring a close point as a contour point;
and S156, sequentially connecting the contour points to form a smooth curve, namely 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 centerline of the blood vessel to obtain the contour line of the blood vessel.
S200, selecting an adjustment area of the interested blood vessel contour line;
s300, as shown in fig. 6, adjusting the contour line of the blood vessel in the adjustment region to approach or depart from the centerline of the blood vessel to obtain a corrected contour line of the blood vessel, including:
s310, obtaining a single adjustment step length m;
s320, selecting any point in 1/3-2/3 areas in the adjusting area as an original contour point;
s330, acquiring points on the center line of the blood vessel corresponding to the original contour points as corresponding center points;
s340, acquiring unit vectors and Euclidean distances between the original contour points and corresponding central points;
s350, moving the original contour point along the unit vector by a plurality of steps to obtain a contour adjusting point;
and S360, in the adjusting area, picking up the starting point, the contour adjusting point, the plurality of seed points and the tail point of the blood vessel contour line, 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, the embodiment 1 further includes: s400, correcting and adjusting the blood vessel correction contour, as shown in fig. 8, includes:
s410, acquiring an included angle alpha of an original contour point, a corresponding central point and a contour adjusting point on a blood vessel contour line;
s420, acquiring the Euclidean distance L between the corresponding center point and the contour adjusting point;
s430, if alpha is more than 90 degrees or L is less than DminAcquiring a contour replacement point according to the unit vector of the original contour point and the corresponding central point, the coordinates of the original contour point, the minimum vessel diameter of the blood vessel and the coordinates of the corresponding central point;
s440, if alpha is less than or equal to 90 degrees and L is more than or equal to DminTaking the contour adjusting point as a wheel buckle replacing point;
s450, in the adjusting area, the starting point, the contour replacing point, the plurality of seed points and the tail point of the blood vessel contour line are picked up, and a new blood vessel correction contour line is obtained by adopting cubic B-spline fitting curve algorithm fitting.
The present application provides a method of adjusting a vessel centerline, comprising:
the method for regulating the blood vessel contour;
and acquiring central curves of the two blood vessel correction contour lines, namely the blood vessel correction central lines.
The present application provides a computer storage medium, a computer program being executed by a processor for implementing the above-mentioned method of adjusting a blood vessel contour.
As will be appreciated by one skilled in the art, aspects of the present invention may be embodied as a system, method or computer program product. Thus, various aspects of the invention may be embodied in the form of: an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining hardware and software aspects that may all generally be referred to herein as a "circuit," module "or" system. Furthermore, in some embodiments, aspects of the invention may also be embodied in the form of a computer program product in one or more computer-readable media having computer-readable program code embodied therein. Implementation of the method and/or system of embodiments of the present invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
For example, hardware for performing selected tasks according to embodiments of the invention could be implemented as a chip or a circuit. As software, selected tasks according to embodiments of the invention could be implemented as a plurality of software instructions being executed by a computer using any suitable operating system. In an exemplary embodiment of the invention, one or more tasks according to exemplary embodiments of the methods and/or systems as described herein are performed by a data processor, such as a computing platform for executing a plurality of instructions. Optionally, the data processor comprises volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, e.g. a magnetic hard disk and/or a removable medium. Optionally, a network connection is also provided. A display and/or a user input device, such as a keyboard or mouse, is optionally also provided.
Any combination of one or more computer readable media may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following:
an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
A computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, for example, in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
For example, computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + + or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer-readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks.
The computer program instructions may also be loaded onto a computer (e.g., a coronary artery analysis system) or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer, other programmable apparatus or other devices provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
The above embodiments of the present invention have been described in further detail for the purpose of illustrating the invention, and it should be understood that the above embodiments are only illustrative of the present invention and are not to be construed as limiting the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (14)
1. A method of adjusting a blood vessel contour, comprising:
obtaining a blood vessel central line and a blood vessel contour line of coronary artery two-dimensional radiography;
selecting an adjustment area of the interested blood vessel contour line;
and adjusting the contour line of the blood vessel in the adjusting area to be close to or far from the center line of the blood vessel to obtain a corrected contour line of the blood vessel.
2. The method of claim 1, wherein the adjusting the contour line of the blood vessel to be within the adjustment region toward or away from the center line of the blood vessel to obtain the modified contour line of the blood vessel comprises:
acquiring a single adjustment step length m;
selecting any point in 1/3-2/3 areas in the adjusting area as an original contour point;
acquiring points corresponding to the original contour points and located on the center line of the blood vessel as corresponding center points;
acquiring unit vectors and Euclidean distances between the original contour points and the corresponding central points;
moving the original contour point along the unit vector by a plurality of steps to obtain a contour adjusting point;
and in the adjusting area, picking up the starting point, the contour adjusting point, the plurality of seed points and the tail point of the blood vessel contour line, and fitting by adopting a cubic B-spline fitting curve algorithm to obtain a blood vessel correction contour line.
3. The method of adjusting a vessel contour as defined in claim 2, further comprising: and correcting and adjusting the blood vessel correction contour line.
4. The method for adjusting the contour of a blood vessel according to claim 3, wherein the method for correcting and adjusting the blood vessel correction contour line comprises:
acquiring an included angle alpha of an original contour point, a corresponding central point and a contour adjusting point on the blood vessel contour line;
acquiring an Euclidean distance L between the corresponding central point and the contour adjusting point;
if said alpha is > 90 DEG or L < DminAcquiring a contour replacement point according to the unit vector of the original contour point and the corresponding central point, the coordinates of the original contour point, the minimum vessel diameter of the blood vessel and the coordinates of the corresponding central point;
if alpha is less than or equal to 90 degrees and L is more than or equal to DminTaking the contour adjusting point as a wheel buckle replacing point;
and in the adjusting area, picking up the starting point, the contour replacing point, the plurality of seed points and the tail point of the blood vessel contour line, and fitting by adopting a cubic B-spline fitting curve algorithm to obtain a new blood vessel correction contour line.
5. The method for adjusting the contour of a blood vessel as claimed in claim 1, wherein the method for acquiring the centerline and contour of the blood vessel of the coronary artery two-dimensional contrast comprises:
extracting a blood vessel central line according to the coronary artery two-dimensional contrast image;
obtaining a straightened blood vessel image according to the blood vessel central line;
setting a blood vessel diameter threshold value D on the straightened blood vessel imageThreshold(s)；
According to said DThreshold(s)Generating preset blood vessel contour lines on two sides of the blood vessel central straight line;
gradually drawing the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel;
and projecting the contour line of the straightened blood vessel back to the image of the center line of the blood vessel to obtain the contour line of the blood vessel.
6. The method for adjusting the contour of a blood vessel as claimed in claim 5, wherein the method for obtaining the straightened blood vessel image according to the center line of the blood vessel comprises:
straightening the central line of the blood vessel to obtain a central straight line of the blood vessel;
dividing the local blood vessel region map into x units along the extending direction of the blood vessel from the starting point to the ending point, wherein x is a positive integer;
correspondingly arranging the center line of the blood vessel of each unit along the center line of the blood vessel;
and the correspondingly set image is the straightened blood vessel image.
7. The method for adjusting the contour of a blood vessel according to claim 6, wherein the step-by-step approach of the preset contour line of the blood vessel to the central straight line of the blood vessel to obtain the contour line of the straightened blood vessel comprises:
dividing the preset contour line of the blood vessel into y units, wherein y is a positive integer;
acquiring z points of each unit, which are positioned on each preset blood vessel contour line;
respectively closing the z points to the blood vessel center straight line in a grading way along the direction vertical to the blood vessel center straight line to generate z closing points, wherein z is a positive integer;
setting RGB difference threshold to delta RGBThreshold(s)Along the direction perpendicular to the blood vessel center straight line, comparing the RGB value of the close point with the RGB value of the point on the blood vessel center straight line every time of closing, and when the difference value is less than or equal to delta RGBThreshold(s)When the blood vessel is in the closed state, the closing point stops closing towards the center line of the blood vessel;
acquiring the approach point as a contour point;
and sequentially connecting the contour points to form a smooth curve which is the contour line of the straightened blood vessel.
8. The method for adjusting the contour of a blood vessel as claimed in claim 7, wherein the method for extracting the centerline of the blood vessel from the coronary artery two-dimensional contrast image comprises:
reading a coronary artery two-dimensional contrast image;
obtaining a vessel segment of interest;
picking up a starting point, a seed point and an end point of the vessel segment of interest;
respectively segmenting two-dimensional contrast images between two adjacent points of a starting point, a seed point and an end point to obtain at least two local blood vessel region images;
extracting at least one blood vessel local path line from each local blood vessel region map;
connecting corresponding blood vessel local path lines on each local blood vessel region map to obtain at least one blood vessel path line;
and selecting one blood vessel path line as the blood vessel central line.
9. The method for adjusting blood vessel contour according to claim 8, wherein the method for extracting at least one blood vessel local path line from the local blood vessel region map respectively comprises:
performing image enhancement processing on the local blood vessel region image to obtain a rough blood vessel image with strong contrast;
and performing grid division on the rough blood vessel map, and extracting at least one blood vessel local path line along the direction from the starting point to the end point.
10. The method for adjusting the contour of a blood vessel as claimed in claim 9, wherein said image enhancement processing on said local blood vessel region map to obtain a coarse blood vessel map with strong contrast comprises:
in each local blood vessel region image, the blood vessel section of interest is used as a foreground, other regions are used as backgrounds, the foreground is strengthened, the backgrounds are weakened, and the rough blood vessel image with strong contrast is obtained.
11. The method of claim 10, wherein the step of gridding the rough vessel map to extract at least one local path line of the vessel along the direction from the starting point to the ending point comprises:
gridding the rough vessel map;
searching the shortest time path between the starting point and the intersection points on the peripheral n grids along the extending direction of the blood vessels from the starting point to the ending point to serve as a second point, searching the shortest time path between the second point and the intersection points on the peripheral n grids to serve as a third point, and repeating the steps at the third point until the shortest time path reaches the ending point, wherein n is a positive integer greater than or equal to 1;
and connecting the extending directions of the blood vessels from the starting point to the ending point according to the searching sequence to obtain at least one blood vessel local path line.
12. The method of adjusting a blood vessel contour according to claim 11, wherein the method of selecting one of the blood vessel route lines as the blood vessel center line comprises:
if the number of the blood vessel path lines is two or more, summing the time from the starting point to the end point of each blood vessel path line;
the vessel path line that is the least in time is taken as the vessel centerline.
13. A method of adjusting a centerline of a blood vessel, comprising:
a method of modulating blood vessel contours as defined in any one of claims 1 to 12;
and acquiring central curves of the two blood vessel correction contour lines, namely the blood vessel correction central lines.
14. A computer storage medium, characterized in that a computer program is executed by a processor for implementing the method of adjusting a blood vessel contour according to any one of claims 1 to 12.
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 true CN112419276A (en) | 2021-02-26 |
| CN112419276B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113256707A (en) * | 2021-06-25 | 2021-08-13 | 西南交通大学 | 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 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115005846B (en) * | 2022-08-08 | 2022-10-11 | 江苏盛恩祥生物技术有限公司 | Data processing method and device suitable for blood vessel development |
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 (5)
| 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 |
| US20130184570A1 (en) * | 2011-07-19 | 2013-07-18 | Toshiba Medical Systems Corporation | Moving object contour extraction apparatus, left ventricle image segmentation apparatus, moving object contour extraction method and left ventricle image segmentation method |
| 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 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113256707A (en) * | 2021-06-25 | 2021-08-13 | 西南交通大学 | Rail center line extraction method based on rail top surface seed point automatic tracking |
| 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112419276B (en) | 2023-12-05 |
| WO2022109905A1 (en) | 2022-06-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210272030A1 (en) | Systems and methods for using geometry sensitivity information for guiding workflow | |
| EP3394829B1 (en) | Heart model guided coronary artery segmentation | |
| JP2022169579A (en) | Diagnostically useful results in real time | |
| US9424637B2 (en) | Vessel segmentation method and apparatus using multiple thresholds values | |
| 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 | |
| CN112419276A (en) | Method for regulating blood vessel contour and central line and storage medium | |
| CN112116615B (en) | Method and device for acquiring blood vessel contour line according to blood vessel center line | |
| CN107580470B (en) | Instantaneous flow reserve-computed tomography | |
| CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
| CN112132949A (en) | Blood vessel three-dimensional modeling method, device and system with narrow lesion interval | |
| CN117115150B (en) | Method, computing device and medium for determining branch vessels | |
| WO2022109906A1 (en) | Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium | |
| CN117036530B (en) | Cross-modal data-based coronary artery fractional flow reserve prediction method and device | |
| US11010899B2 (en) | Image processing apparatus, control method of image processing apparatus, and storage medium | |
| CN112419277B (en) | Three-dimensional vascular centerline synthesis method, system and storage medium | |
| CN112472112A (en) | Method, system and storage medium for regulating vascular stenosis region | |
| US10600182B2 (en) | Image processing apparatus and image processing method, that determine a conformable image | |
| CN112419279B (en) | Method for selecting two-dimensional image and synthesizing three-dimensional blood vessel and storage medium | |
| 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 |