CN110338830B - Method for automatically extracting head and neck blood vessel central path in CTA image - Google Patents
Method for automatically extracting head and neck blood vessel central path in CTA image Download PDFInfo
- Publication number
- CN110338830B CN110338830B CN201910693725.6A CN201910693725A CN110338830B CN 110338830 B CN110338830 B CN 110338830B CN 201910693725 A CN201910693725 A CN 201910693725A CN 110338830 B CN110338830 B CN 110338830B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- region
- candidate
- artery
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 171
- 238000000034 method Methods 0.000 title claims abstract description 29
- 210000001715 carotid artery Anatomy 0.000 claims abstract description 64
- 210000002385 vertebral artery Anatomy 0.000 claims abstract description 61
- 238000001914 filtration Methods 0.000 claims abstract description 45
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 210000002376 aorta thoracic Anatomy 0.000 claims abstract description 10
- 238000010606 normalization Methods 0.000 claims abstract description 6
- 230000005484 gravity Effects 0.000 claims description 31
- 210000001519 tissue Anatomy 0.000 claims description 23
- 210000001841 basilar artery Anatomy 0.000 claims description 19
- 210000005013 brain tissue Anatomy 0.000 claims description 12
- 210000000988 bone and bone Anatomy 0.000 claims description 10
- 210000004004 carotid artery internal Anatomy 0.000 claims description 10
- 238000012216 screening Methods 0.000 claims description 10
- 230000002792 vascular Effects 0.000 claims description 10
- 210000001367 artery Anatomy 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 6
- 210000003625 skull Anatomy 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 238000007917 intracranial administration Methods 0.000 claims description 2
- 210000004072 lung Anatomy 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims 3
- 230000004807 localization Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 4
- 210000003270 subclavian artery Anatomy 0.000 description 4
- 210000000711 cavernous sinus Anatomy 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 1
- 210000000269 carotid artery external Anatomy 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 210000001951 dura mater Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 210000000103 occipital bone Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 231100000216 vascular lesion Toxicity 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/501—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Vascular Medicine (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention provides a method for automatically extracting a central path of a head-neck blood vessel in a CTA image, which is suitable for CTA image data scanned from an aortic arch to a cranial vertex, and comprises the following steps: s1: automatically positioning the starting point and the end point of the carotid artery and the vertebral artery respectively according to the shape, the gray scale and the position characteristics of the blood vessel in the cross section of the CTA image; s2: performing multi-scale gradient normalization blood vessel filtering based on Raycasting on the region according with the gray features of the blood vessel between the start point and the end point of the blood vessel positioned in the step S1 to construct a blood vessel filtering weight value graph; s3: in the vessel filtering weight value map constructed in the step S2, the extraction of the central paths of the carotid artery and the vertebral artery is finished by adopting a Dijkstra optimal path extraction algorithm. The invention has the advantages that: the automatic positioning of the blood vessel is realized without manual participation of a user; the extracted central line is accurate, and the later adjustment of a user is not needed; the method has high running speed and meets the real-time requirement.
Description
Technical Field
The invention relates to a method for automatically extracting a head and neck blood vessel central path in a CTA image. The invention belongs to the technical field of medical image processing.
Background
Cerebrovascular diseases are one of the main diseases threatening human health at present, and clinically, the three-dimensional structure information of blood vessels is obtained by analyzing images after CTA (non-invasive vascular imaging examination) images and extracting the central path of the blood vessels to assist doctors in judging the degree of vascular lesions. However, since the distribution of blood vessels in a human body is very complex, the blood vessels often pass through bones, and both the bones and the blood vessels present high-brightness gray features in an image generated after a CTA image, which brings great difficulty to the extraction of a blood vessel center path and lesion analysis.
Currently, when a blood vessel analysis is performed on an image generated after a CTA image, a subtraction technique is usually used to remove bones in the CTA image to obtain a blood vessel region, and then a starting point and a terminating point of a blood vessel are manually defined to obtain a central path of the blood vessel. The disadvantages of this method are: because the subtraction operation requires two scans (CT and CTA scans) of the patient, the patient is likely to shift during the interval between the two scans, which may cause the condition of missing or breaking of the blood vessel structure after subtraction, and further affect the subsequent lesion analysis of the blood vessel.
Although the method for acquiring the central path of the blood vessel by using the blood vessel enhancement method based on Hessian can directly extract the central path of the blood vessel in a CTA image, the method is not ideal in general effect on the blood vessel with complex shape and surrounding tissues, such as the encephalic segment of a vertebra and an internal carotid artery, and also needs a doctor to manually define the starting point and the ending point of the blood vessel, so that certain burden is brought to the analysis work of the doctor.
Disclosure of Invention
In view of the above, the present invention provides a method for automatically extracting central paths of the head and neck blood vessels, i.e., the vertebral artery blood vessels and the carotid artery blood vessels, in the CTA image.
In order to achieve the purpose, the invention adopts the following technical scheme: a method for automatically extracting a central path of a head-neck blood vessel in a CTA image is suitable for CTA image data scanned from an aortic arch to a cranial vertex, and is characterized in that: it comprises the following steps:
s1: automatically positioning the starting point and the end point of the carotid artery and the vertebral artery respectively according to the shape, the gray scale and the position characteristics of the blood vessel in the cross section of the CTA image;
s2: performing Raycasting-based multi-scale gradient normalized blood vessel filtering on the region according with the gray level characteristics of the blood vessel between the starting point and the end point of the blood vessel positioned in the step S1 to construct a blood vessel filtering weight value graph;
s3: in the vessel filtering weight value map constructed in the step S2, the extraction of the central paths of the carotid artery and the vertebral artery is finished by adopting a Dijkstra optimal path extraction algorithm.
The invention has the advantages that: the automatic positioning of the blood vessel is realized without manual participation of a user; the extracted central line is accurate, and the later adjustment of a user is not needed; the method has high running speed and meets the real-time requirement.
Drawings
Fig. 1 is a flowchart of a method for automatically extracting a central path of a head and neck blood vessel according to the present invention;
FIG. 2 is a graph of multi-scale gradient normalized vascular filtering based on Raycasting (ray casting);
FIG. 3A is the origin (green marked point) of the carotid and vertebral arteries located after step S1.1.2;
FIG. 3B is the vertebral artery end point (green marker point) located after step S1.2.1.2;
FIG. 3C is the carotid endpoint (green marker point) located after step S1.2.2.2;
fig. 3D is a graph of the right carotid artery central path extracted after step S3.2 and its corresponding CPR (curved reconstruction) results;
fig. 3E is a graph of the left carotid artery central path extracted after step S3.2 and its corresponding CPR (curved reconstruction) results;
FIG. 3F is a graph of the right vertebral artery central path extracted after step S3.2 and its corresponding CPR (curved reconstruction) results;
fig. 3G is a graph of the left vertebral artery central path extracted after step S3.2 and its corresponding CPR (curved reconstruction) results.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and examples. It should be noted that various modifications can be made to the embodiments disclosed herein, and therefore, the embodiments disclosed in the specification should not be construed as limiting the present invention, but merely as exemplifications of embodiments thereof, which are intended to make the features of the present invention obvious.
As shown in fig. 1, the method for automatically extracting the central path of the head and neck blood vessels (vertebral artery blood vessels and carotid artery blood vessels) in the CTA image provided by the present invention is suitable for CTA image data (hereinafter, described as CTA image data meeting the requirement) scanned from the aortic arch to the cranial vertex. The method provided by the invention mainly comprises three key steps:
s1: according to the shape, gray scale and position characteristics of the blood vessel in the cross section of the CTA image, automatically positioning the starting point and the end point of the carotid artery and the vertebral artery respectively;
s2: performing Raycasting-based multi-scale gradient normalized blood vessel filtering on the region according with the gray level characteristics of the blood vessel between the starting point and the end point of the blood vessel positioned in the step S1, and constructing a blood vessel filtering weight value graph according to the calculated filtering value;
s3: in the vessel filtering weight value map constructed in the step S2, the extraction of the central paths of the carotid artery and the vertebral artery is finished by adopting a Dijkstra optimal path extraction algorithm.
The human vertebral artery starts from the subclavian artery, passes through the transverse apopore of the vertebra, bends towards the inner side at the back of the lateral mass of the atlas, passes through the macropore of the occipital bone and enters the cranial cavity, and is combined with the contralateral vertebral artery to form the basilar artery at the lower edge of the brain bridge, so the invention needs to position the starting point of the left and right vertebral arteries and simultaneously needs to position the basilar artery point as the termination point of the left and right vertebral arteries.
The carotid artery originates from the aortic arch and subclavian artery and branches upwardly into the internal and external carotid arteries. Wherein the internal carotid artery extends upwards to enter the cranial cavity through the rupture hole and enters the cavernous sinus through the subkaryoid sinus, an S-shaped bend is formed in the cavernous sinus and extends forwards to extend through the dura mater, and the branch has a posterior traffic branch, a cerebral anterior and a middle artery to participate in forming a cerebral artery loop. The invention needs to separately locate the starting point of the left and right carotid artery and the termination point of the internal carotid artery.
Step S1 of the invention: the method comprises the following steps of automatically positioning the starting point and the end point of a carotid artery and a vertebral artery respectively according to the shape and the gray scale of a blood vessel in the cross section of a CTA image and the position characteristics, and comprises the following specific steps:
s1.1 automatic positioning of starting points of vertebral artery blood vessels and carotid artery blood vessels
S1.1.1 locating candidate levels of carotid and vertebral artery origins
Firstly, performing binarization preprocessing on CTA image data meeting requirements to distinguish human tissues and background areas in an image; then identifying a hole area with a gray value less than-450 HU in the human tissue area, scanning layer by layer from bottom to top (from the autonomous artery arch to the cranial vertex), and when the area of a single hole is more than 60mm 2 And less than 1200mm 2 When the current level is considered as the upper edge of the lung parenchyma, the level is taken as a candidate level of the starting points of the vertebral artery and the carotid artery, and the scanning is stopped;
s1.1.2 locating the starting point of vertebral artery and carotid artery
And (3) performing label extraction on all regions with the gray values larger than 150HU in the candidate layer to obtain candidate regions of vertebral artery and carotid artery blood vessels, and screening the candidate regions by using the following rules:
1) Calculating the gravity center of the marked tissue, and removing the area of the gravity center which is not marked;
2) On the basis of the step 1), calculating the adjacent edge ratio of the minimum circumscribed rectangle of the remaining candidate region, and removing the region with the ratio larger than 3 or smaller than 1/3;
3) According to the size characteristics of the normal vertebral artery and the carotid artery, on the basis of the steps 1) and 2), further removing the candidate region with the area larger than 60mm 2 And less than 3mm 2 The area of (a);
4) Calculating the gray level mean value and the variance of the remaining candidate regions one by one, taking the gravity center of the current candidate region as a seed point, taking the 2-time variance as a growth step length, and carrying out three-dimensional region growth in the upper and lower 3-layer image ranges by taking the layer of the current candidate region as a reference; judging the region grown on each layer again according to the adjacent edge ratio in the step 2), and removing the candidate regions which do not meet the conditions;
5) The cross section of the blood vessel presents a circular-like characteristic in the CTA image, and the distribution in the human body is approximately centrosymmetric, so that the circularity, central Symmetry, bounding box aspect Ratio and Area are respectively calculated for the region screened in the step 4), wherein the central Symmetry represents that whether a candidate region exists in the range of 5mm above, below, left and right of the opposite side is searched by taking the center (x direction) of the human tissue after binarization as a reference; if the distance can be found, calculating the distance between the gravity centers of all the candidate blood vessel regions in the search range and the gravity center of the current region in the y direction, and taking the corresponding minimum distance as the central symmetry degree of the candidate blood vessel region; if not, setting the central symmetry of the current blood vessel candidate region as the side length of the whole image;
6) Sorting the 4 metric parameters calculated in step 5): the circularity and the area of the region are sorted from large to small, the horizontal-longitudinal ratio of the bounding boxes is sorted from small to large according to the difference value between the circularity and the area of the region and the area of the bounding boxes, the central symmetry is sorted from small to large, and the sorted order is used as the fraction of the value of the region metric: the circularity fraction is denoted S circlarity And the symmetry fraction is represented as S symmetry The aspect ratio fraction is represented as S ratio The area fraction is represented as S area ；
7) Defining a vessel cross-section similarity score: score = w1 × S circlarity +w2*S symmetry +w3*S ratio +w4*S area W1, w2, w3 and w4 respectively represent weights corresponding to the circularity fraction, the symmetry fraction, the aspect ratio fraction and the area fraction of the candidate blood vessel region; the smaller the vessel similarity score is, the more likely the region is to be a vessel region;
8) Sorting the blood vessel similarity Score according to a sequence from small to large to obtain a candidate region ranked in the first four; firstly, defining two regions with larger areas as the initial regions of the left and right carotid arteries, and defining the other two regions as the initial regions of the left and right vertebral arteries; then distinguishing left and right carotid arteries and left and right vertebral arteries according to the obtained x coordinates of the centers of gravity of the carotid artery and vertebral artery initial regions; finally, the barycenter of the blood vessel region is used as the starting point of the left and right carotid arteries (such as points A and B in FIG. 3A) and the starting point of the left and right vertebral arteries (such as points C and D in FIG. 3A).
9) If the screened areas are smaller than 4, or the screened areas are positioned at one side of the center (x direction) of the human tissue, the scanning upwards is continued by taking the current candidate layer as the reference, and simultaneously the starting points of the carotid artery and the vertebral artery are continuously searched by returning to S1.1.1.
Since the vertebral artery originates from the subclavian artery, the carotid artery originates from the aortic arch and the subclavian artery, and the starting points of the vertebral artery and carotid artery are located in the lower half of the CTA image data that meets the requirements, the positioning of the starting points of the vertebral artery and carotid artery is performed until the lower half of the CTA data is scanned layer by layer.
S1.2 automatic definition of the end points of vertebral artery and carotid artery blood vessels
S1.2.1 automatic definition of vertebral artery vascular endpoint (i.e. basilar artery point)
S1.2.1.1 locating basilar artery point candidate slice
Scanning layer by layer (from the vertex to the aortic arch direction) from top to bottom, extracting real brain tissue by taking a 120HU gray value as a threshold value, and simultaneously calculating the real brain tissue and the body tissue bounding box after binarization processing in S1.1.1; when the bounding box area size of the brain tissue is 1/3 of the volume tissue bounding box, the scanning is stopped, and the layer is taken as a candidate layer of the basilar artery.
S1.2.1.2 locates the basilar artery point
Performing label extraction on the region with the gray value of 150-750HU in the candidate layer to obtain a candidate region of the basilar artery, and screening the candidate region by using the following rules:
1) Calculating the gravity center of the marked tissue, and removing the area of the gravity center which is not marked;
2) Adopting a blood vessel enhancement function based on a Hessian matrix to perform enhancement processing on the candidate region filtered in the step 1);
3) Removing the candidate area corresponding to the threshold value by taking 0.6 as the threshold value of the blood vessel enhancement;
4) Taking the gravity center of the candidate region filtered in the step 3) as a seed point and 100HU as a step length, and performing two-dimensional region growth;
5) According to the size characteristics of normal basilar artery, the removal area is less than 1mm 2 Greater than 35mm 2 The blood vessel candidate region of (a);
6) Calculating a bounding box of each candidate region, and removing regions with the bounding box aspect ratio larger than 4;
7) Extracting the boundary of the blood vessel candidate region, calculating the maximum and minimum distances from the gravity point to the boundary point, and removing the blood vessel candidate region of which the ratio of the maximum distance to the minimum distance is more than 8;
8) Taking the gravity center of the blood vessel candidate region as a seed point and 100HU as a step length, carrying out three-dimensional region growth on the upper layer image and the lower layer image, extracting a blood vessel region adjacent to the current blood vessel candidate region, and calculating a bounding box of the adjacent blood vessel region; removing the blood vessel candidate area with the aspect ratio of the bounding box being more than 4 or the longest side being more than 10 mm;
9) Calculating the centering degree, the height and the area of the blood vessel candidate region screened in the step 8), wherein the centering degree represents the distance between the gravity center of the blood vessel candidate region and the center of the brain tissue; the height represents the coordinate of the gravity center y direction of the blood vessel candidate region;
10 ) the 3 metric parameters calculated in step 9) are sorted: wherein the intermediate levels are ordered from small to large, height andthe areas are sorted from big to small, and the sorted order is used as the fraction of the metric value of the region: the median score is denoted S center Height fraction is represented as S y And the area fraction is represented as S area ；
11 Construction of vascular similarity Score = w1 × S) center +w2*S y +w3*S area W1, w1 and w3 respectively represent the score of the median of the candidate blood vessel region, the height score and the weight corresponding to the area; taking the gravity center of the candidate blood vessel region with the smallest blood vessel similarity score as a basilar artery point, such as point E in fig. 3B;
12 If no blood vessel candidate area exists after the screening in the step 8), continuing to scan downwards by taking the current layer as a reference, and returning to S1.2.1.1 to continue searching for the base artery point.
For satisfactory CTA image data, the basilar artery points are in the upper half of the data, so the basilar artery points are positioned until the upper half of the CTA data is scanned layer-by-layer.
S1.2.2 locates the carotid endpoint, i.e., the internal carotid endpoint
S1.2.2.1 locating carotid artery endpoint candidate level
Scanning layer by layer from top to bottom (from the cranial vertex to the aortic arch direction), calculating a skull bounding box, constructing an intracranial tissue bounding box by using an area which is contracted inwards by 1/2 of the skull bounding box, and counting the number of pixel points of which the internal gray value is more than 550 HU; if all pixel points in the range of the bounding box are more than 550HU, judging that the current layer is a cranial vertex region, not judging at the moment, and continuing to scan downwards; when the gray values of all pixels in the bounding box are less than 550HU, the scanning is considered to enter the brain tissue area; and after the brain tissue area appears, judging the proportion of the number of the bone pixel points in the bounding box to the number of all the pixel points in the tissue bounding box, if the proportion is more than 0.005, stopping searching, and taking the layer as a candidate layer of the internal carotid artery endpoint.
S1.2.2.2 locating carotid endpoint
1) Calculating the mean value CT of the gray level of the blood vessel region where the carotid artery starting point positioned in S1.1 is located mean Sum standard deviation CT std Using carotid artery starting point as seed point and CT mean +CT std 、CT mean -CT std Performing unidirectional upward three-dimensional region growth as an upper threshold and a lower threshold until a candidate layer of the internal carotid artery endpoint is reached, and marking a candidate blood vessel region in the candidate layer;
2) In the layer where the internal carotid artery endpoint is located, searching pixel points with the gray value higher than 750HU, taking the pixel points as bone tissue seed points, extracting all tissues communicated with the seed points with the threshold value within the range of 120-3071HU in a two-dimensional region growing mode, and removing the candidate blood vessel region adhered to the bone in the step 1);
3) Calculating the central Symmetry, the bounding box aspect Ratio and the Area of the blood vessel candidate region obtained after screening in the step 2), wherein the definition of the central Symmetry is the same as that of the step 5) in the S1.1.2;
4) Sorting the 3 metric parameters calculated in step 3): the central symmetry degrees are sorted from small to large, the horizontal-vertical ratio of the bounding boxes is sorted from small to large according to the difference value of the horizontal-vertical ratio of the bounding boxes and 1, and the area of the area is sorted from large to small. And taking the sorted order as the score of the region metric value: the central symmetry fraction is denoted S symmetry And the horizontal-vertical ratio fraction of the bounding box is expressed as S ratio And the area fraction of the region is represented as S area ；
5) Defining a vessel cross-section similarity score: score = w1 × S symmetry +w2*S ratio +w3*S area . Wherein, w1, w2 and w3 respectively represent the central symmetry fraction of the candidate blood vessel region, the horizontal-vertical ratio fraction of the bounding box and the weight corresponding to the region area fraction. The smaller the vessel similarity score is, the more likely the region is a vessel region;
6) Sorting the blood vessel section similarity scores Score in a descending order, and defining two blood vessel candidate regions with the highest scores as the corresponding gravity centers as the end points of the left and right carotid arteries according to the x coordinate direction, such as points F and G in FIG. 3C;
7) If the candidate blood vessel regions screened in the step 2) are less than 2, continuing to scan downwards by taking the current layer as a reference, and simultaneously returning to the step 2) to continue searching for the stiff artery end point.
For satisfactory CTA image data, the carotid endpoint is located at the top half of the data, so the carotid endpoint is located until the top half of the CTA data is scanned layer-by-layer. .
Step S2 of the invention: and between the starting point and the end point of the blood vessel positioned in the step S1, carrying out Raycasting-based multi-scale gradient normalized blood vessel filtering on the region conforming to the gray level characteristics of the blood vessel, and constructing a blood vessel filtering weight value graph.
S2.1 vascular Pre-extraction
Respectively calculating the mean value CT of the gray level of the corresponding region according to the carotid artery and vertebral artery blood vessel regions positioned in S1.1 mean And standard deviation CT std Using carotid and vertebral artery starting points located in S1.1 as seed points, CT mean +2*CT std 、CT mean -2*CT std Performing unidirectional upward three-dimensional region growth as an upper threshold and a lower threshold; the carotid artery grows to the level of the carotid artery terminal point, and the vertebral artery grows to the level of the vertebral artery terminal point;
s2.2 Raycasting-based multi-scale gradient normalization vascular filtering
In the CTA image, the cross section of the blood vessel has a feature of a shape similar to a circle, and the gray scale change has a feature of a gaussian-like distribution (the gray scale value at the center of the blood vessel is high, and the gray scale value gradually decreases with the increase of the radius with the center as the center of the circle). The multi-scale gradient normalization blood vessel filtering based on Raycasting just utilizes the two characteristics of blood vessels to calculate the blood vessel filtering value of each pixel point in the pre-extracted blood vessel region in S2.1, and the specific steps are as follows:
1) Since the cross-sectional dimension of the blood vessel varies in its course, it is necessary to calculate the filtered value of the blood vessel at different radius scales (as shown in fig. 2), R e (R ∈) min ,R max ](in the subsequent step with R min And R max Represent the minimum and maximum radius dimensions, no longer distinguish between carotid and vertebral);
2) Any blood vessel pre-extracted in S2.1Taking the region point as the center of a circle, respectively calculating the gradient response V of the projected radius scale R and the center of the circle point R Minimum V of gradient response on the 0 to R radius scale R,min ，R min To R max Maximum value V at radius scale R,max ；
3) Defining the normalized gradient response at the radius scale R as E R ＝(V R -V R,min )/V R,max If the current point is the center point of the blood vessel, E R About equal to 1; if it is a non-vascular point, E R About 0;
4) Repeating steps 2) -3), calculating a plurality of projection radius scales R min To R max E of R And will be at different scales E R The maximum value of (2) is used as a final blood vessel filtering value of the point;
s2.3 construction of a vessel filtering weight map
For each blood vessel point extracted in S2.1, a blood vessel filtering value can be calculated through S2.2, and a blood vessel filtering weight value of the blood vessel point is defined as:
wherein the content of the first and second substances,
and
respectively, the reciprocal of the blood vessel filtering value calculated in S2.2 for two adjacent points, dist is the physical distance between the two points.
And constructing a blood vessel filtering value weight map according to the calculated blood vessel filtering weights of all the blood vessel points.
Step S3 of the invention: in the vessel filtering weight value map constructed in the step S2, the extraction of the central paths of the carotid artery and the vertebral artery is finished by adopting a Dijkstra optimal path extraction algorithm.
Obtaining the shortest central path between the starting point and the end point of the blood vessel by adopting an optimal path extraction algorithm of Dijkstra, and firstly, defining an energy function between the starting point and the end point of the blood vessel positioned in the step S1:
E＝∫(w(s)+ε)ds
wherein w (S) is the vessel filtering enhancement weight calculated in S2.3, epsilon is a regular term, and S is a path between the starting point and the end point.
And then, selecting a path with the minimum blood vessel filtering weight sum in all paths between the starting point and the end point of the blood vessel as a blood vessel center path.
Fig. 3D is a diagram showing the extracted right carotid artery central path and its corresponding CPR (curved surface reconstruction) result, fig. 3E is a diagram showing the extracted left carotid artery central path and its corresponding CPR (curved surface reconstruction) result, fig. 3F is a diagram showing the extracted right vertebral artery central path and its corresponding CPR (curved surface reconstruction) result, and fig. 3G is a diagram showing the extracted left vertebral artery central path and its corresponding CPR (curved surface reconstruction) result.
Compared with the prior art, the invention has the advantages that:
1. the automatic positioning of the blood vessel is realized without manual participation of a user;
2. the extracted central line is accurate, and the later adjustment of a user is not needed;
3. the method has high running speed and meets the real-time requirement.
Finally, it should be noted that: the above-mentioned embodiments are only used for illustrating the technical solution of the present invention, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (6)
1. A method for automatically extracting the central path of a head-neck blood vessel in a CTA image is suitable for CTA image data scanned from an aortic arch to a cranial vertex, and is characterized in that: it comprises the following steps:
s1: automatically positioning the starting point and the end point of the carotid artery and the vertebral artery respectively according to the shape, the gray scale and the position characteristics of the blood vessel in the cross section of the CTA image; the method specifically comprises the following steps:
s1.1, automatically positioning the starting points of carotid and vertebral artery blood vessels:
s1.1.1, positioning candidate layers of starting points of a carotid artery and a vertebral artery;
carrying out binarization preprocessing on CTA image data, and distinguishing human body tissues and background areas in the image; then identifying the hole area with the gray value less than-450 HU in the human tissue area, and scanning layer by layer from bottom to top, when the area of a single hole is more than 60mm 2 And less than 1200mm 2 When the current level is considered as the upper edge of the lung parenchyma, the level is taken as a candidate level of the starting points of the vertebral artery and the carotid artery, and the scanning is stopped;
s1.1.2 locates the starting point of vertebral artery and carotid artery;
marking and extracting all regions with the gray values larger than 150HU in the candidate layer to obtain candidate regions of the vertebral artery and the carotid artery blood vessel, screening the candidate regions and determining the starting points of the vertebral artery blood vessel and the carotid artery blood vessel;
s1.2, automatically defining the end points of vertebral artery blood vessels and carotid artery blood vessels;
s1.2.1, automatically positioning vertebral artery terminal:
s1.2.1.1 locating a basilar artery point candidate slice;
scanning layer by layer from top to bottom, namely from the cranial vertex to the aortic arch, extracting real brain tissue by taking a 120HU gray value as a threshold value, and simultaneously calculating the real brain tissue and the body tissue bounding box subjected to binarization processing in S1.1.1; when the bounding box area size of the brain tissue is 1/3 of the volume tissue bounding box, stopping scanning, and taking the layer as a candidate layer of the basilar artery;
s1.2.1.2 locating the basilar artery point, the terminal point of the vertebral artery;
marking and extracting a region with the gray value of 150-750HU in a candidate layer to obtain a candidate region of the basilar artery, and screening the candidate region to determine a basilar artery point;
s1.2.2 automatic carotid endpoint localization:
s1.2.2.1 locating carotid artery endpoint candidate level
Scanning layer by layer from top to bottom, namely from the vertex of the skull to the aortic arch, calculating a skull bounding box, constructing an intracranial tissue bounding box by using a region which is shrunk inwards by 1/2 of the skull bounding box, and counting the number of pixel points of which the internal gray value is greater than 550 HU; if all the pixel points in the bounding box range are more than 550HU, judging that the current layer is a cranial vertex area, and continuing to scan downwards without judging; when the gray values of all pixels in the bounding box are less than 550HU, the scanning is considered to enter the brain tissue area; after a brain tissue area appears, judging the proportion of the number of bone pixel points in the bounding box to all the pixel points in the tissue bounding box in which the bone pixel points are located layer by layer, if the proportion is more than 0.005, stopping searching, and taking the layer as a candidate layer of the internal carotid artery endpoint;
s1.2.2.2 position carotid endpoint
1) Calculating the mean value CT of the gray level of the blood vessel region where the carotid artery starting point positioned in S1.1 is located mean Sum standard deviation CT std Taking the origin of carotid artery as a seed point and CT mean +CT std 、CT mean -CT std Performing unidirectional upward three-dimensional region growth as an upper threshold and a lower threshold until reaching a candidate layer of the internal carotid artery terminal point, and marking a candidate blood vessel region in the candidate layer;
2) In the layer where the internal carotid artery endpoint is located, searching pixel points with the gray value higher than 750HU, taking the pixel points as bone tissue seed points, extracting all tissues communicated with the seed points with the threshold value within the range of 120-3071HU in a two-dimensional region growing mode, and removing the candidate blood vessel region adhered to the bone in the step 1);
3) Calculating the central Symmetry, the bounding box aspect Ratio and the Area of the blood vessel candidate region obtained after screening in the step 2), wherein the Symmetry is the same as the definition of the step 5) in the step S1.1.2;
4) Sorting the 3 metric parameters calculated in step 3): wherein, the central symmetry degrees are sorted from small to large, and the horizontal-vertical ratio of the bounding box is from small to small according to the difference value with 1Sorting the areas according to the size, and taking the sorted order as the fraction of the measurement value of the area: the central symmetry fraction is denoted S symmetry And the horizontal-vertical ratio fraction of the bounding box is expressed as S ratio And the area fraction of the region is represented as S area ；
5) Defining a vessel cross-section similarity score: score = w1 × S symmetry +w2*S ratio +w3*S area (ii) a W1, w2 and w3 respectively represent the central symmetry fraction of the candidate blood vessel region, the horizontal-vertical ratio fraction of the bounding box and the weight corresponding to the region area fraction; the smaller the vessel similarity score is, the more likely the region is a vessel region;
6) Sorting the vessel section similarity scores Score in a descending order, and defining two vessel candidate regions with the highest scores as corresponding gravity centers as end points of the left and right carotid arteries according to the x coordinate direction;
7) If the number of the candidate blood vessel regions screened in the step 2) is less than 2, continuing to scan downwards by taking the current layer as a reference, and simultaneously returning to the step 2) to continue searching for a stiff artery end point;
s2: performing Raycasting-based multi-scale gradient normalized blood vessel filtering on the region according with the gray level characteristics of the blood vessel between the starting point and the end point of the blood vessel positioned in the step S1 to construct a blood vessel filtering weight value graph;
s3: in the vessel filtering weight value map constructed in the step S2, the extraction of the central paths of the carotid artery and the vertebral artery is finished by adopting a Dijkstra optimal path extraction algorithm.
2. The method for automatically extracting the central path of the head and neck blood vessels from the CTA image as recited in claim 1, wherein: the method for screening the candidate regions in the step S1.1.2 is as follows:
1) Calculating the gravity center of the marked tissue, and removing the area of the gravity center which is not marked;
2) On the basis of the step 1), calculating the adjacent edge ratio of the minimum circumscribed rectangle of the remaining candidate region, and removing the region with the ratio larger than 3 or smaller than 1/3;
3) According to the normal vertebraThe size characteristics of arteries and carotid arteries are further removed, and the area of the candidate region is larger than 60mm on the basis of the steps 1) and 2) 2 And less than 3mm 2 The area of (a);
4) Calculating the gray level mean value and the variance of the remaining candidate regions one by one, taking the gravity center of the current candidate region as a seed point, taking the 2-time variance as a growth step length, and carrying out three-dimensional region growth in the upper and lower 3-layer image ranges by taking the layer of the current candidate region as a reference; judging the region grown on each layer again according to the adjacent edge ratio in the step 2), and removing the candidate regions which do not meet the conditions;
5) The cross section of the blood vessel presents a similar circular characteristic in a CTA image, and the distribution in the human body is approximately centrosymmetric, so that the circularity, central Symmetry, bounding box aspect Ratio and Area are respectively calculated for the region screened in the step 4), wherein the central Symmetry represents whether a candidate region exists in the range of 5mm above, below, left and right of the opposite side by taking the center of the human tissue after binarization as a reference; if the distance can be found, calculating the distance between the gravity centers of all the candidate blood vessel regions in the search range and the gravity center of the current region in the y direction, and taking the corresponding minimum distance as the central symmetry degree of the candidate blood vessel region; if not, setting the central symmetry of the current blood vessel candidate region as the side length of the whole image;
6) Sorting the 4 metric parameters calculated in step 5): the circularity and the area are sorted from large to small, the transverse-longitudinal ratio is sorted from small to large according to the difference value of 1, the symmetry is sorted from small to large, and the sorted order is used as the fraction of the measurement value of the region: the circularity fraction is denoted S circlarity And the symmetry fraction is represented as S symmetry The aspect ratio fraction is represented as S ratio The area fraction is represented as S area ；
7) Defining a vessel cross-section similarity score: score = w1 × S circlarity +w2*S symmetry +w3*S ratio +w4*S area Wherein w1, w2, w3 and w4 respectively represent the circularity fraction, symmetry fraction, aspect ratio fraction and surface of the candidate blood vessel regionThe weight value corresponding to the integral number; the smaller the vessel similarity score is, the more likely the region is a vessel region;
8) Sorting the blood vessel similarity Score according to a sequence from small to large to obtain a candidate region ranked in the first four; firstly, defining two regions with larger areas as the initial regions of the left and right carotid arteries, and defining the other two regions as the initial regions of the left and right vertebral arteries; then distinguishing the left and right carotid arteries and the left and right vertebral arteries according to the obtained x coordinates of the centers of gravity of the carotid artery and vertebral artery initial regions; finally, the gravity center of the blood vessel region is used as the starting points of the left and right carotid arteries and the left and right vertebral arteries;
9) If the screened areas are smaller than 4 or the screened areas are positioned at one side of the center of the human tissue, the scanning upwards is continued by taking the current candidate layer as the reference, and simultaneously the scanning returns to S1.1.1 to continue searching for the starting points of the carotid artery and the vertebral artery.
3. The method of claim 2 for automatically extracting the central path of the head and neck blood vessels in the CTA image, wherein: the method for screening the candidate regions to determine the basilar artery points in the step S1.2.1.2 is as follows:
1) Calculating the gravity center of the marked tissue, and removing the area of the gravity center which is not marked;
2) Adopting a blood vessel enhancement function based on a Hessian matrix to perform enhancement processing on the candidate region filtered in the step 1);
3) Taking 0.6 as a threshold value of blood vessel enhancement, and removing candidate regions corresponding to the threshold value;
4) Taking the gravity center of the candidate region filtered in the step 3) as a seed point and 100HU as a step length, and performing two-dimensional region growth;
5) According to the size characteristics of normal basilar artery, the removal area is less than 1mm 2 Greater than 35mm 2 The blood vessel candidate region of (a);
6) Calculating a bounding box of each candidate region, and removing regions with the bounding box aspect ratio larger than 4;
7) Extracting the boundary of the blood vessel candidate region, calculating the maximum and minimum distances from the gravity point to the boundary point, and removing the blood vessel candidate region of which the ratio of the maximum distance to the minimum distance is more than 8;
8) Taking the gravity center of the blood vessel candidate region as a seed point and 100HU as a step length, carrying out three-dimensional region growth on the upper layer image and the lower layer image, extracting a blood vessel region adjacent to the current blood vessel candidate region, and calculating a bounding box of the adjacent blood vessel region; removing the blood vessel candidate area with the aspect ratio of the bounding box being more than 4 or the longest side being more than 10 mm;
9) Calculating the centering degree, the height and the area of the blood vessel candidate region screened in the step 8), wherein the centering degree represents the distance between the gravity center of the blood vessel candidate region and the center of the brain tissue; the height represents the coordinate of the gravity center y direction of the blood vessel candidate region;
10 ) the 3 metric parameters calculated in step 9) are sorted: the centering degree is sorted from small to large, the height and the area are sorted from large to small, and the sorted order is used as the fraction of the measurement value of the region: the median score is denoted S center Height fraction Sy and area fraction S area ；
11 ) construct a vascular similarity Score = w 1S center +w2*S y +w3*S area W1, w1 and w3 respectively represent the score of the median of the candidate blood vessel region, the height score and the weight corresponding to the area; taking the gravity center of the blood vessel candidate region with the minimum blood vessel similarity score as a basilar artery point;
12 If no blood vessel candidate area exists after the screening in the step 8), continuing to scan downwards by taking the current layer as a reference, and returning to S1.2.1.1 to continue searching for the base artery point.
4. The method of claim 3 for automatically extracting the central path of the head and neck blood vessels in the CTA image, wherein:
in the step S2, a multi-scale gradient normalization blood vessel filtering based on Raycasting is performed on the region conforming to the gray level feature of the blood vessel between the starting point and the end point of the blood vessel, and a method for constructing a blood vessel filtering weight map is as follows:
s2.1 vascular Pre-extraction
Carotid and vertebral arteries according to localization in S1.1Blood vessel region, calculating the mean value CT of the gray levels of the corresponding regions mean Sum standard deviation CT std Taking carotid and vertebral artery starting points located in S1.1 as seed points and CT mean +2*CT std 、CT mean -2*CT std Performing unidirectional upward three-dimensional region growth as an upper threshold and a lower threshold; the carotid artery grows to the level of the carotid artery terminal point, and the vertebral artery grows to the level of the vertebral artery terminal point;
s2.2, based on Raycasting multi-scale gradient normalization blood vessel filtering, calculating a blood vessel filtering value of each pixel point in the pre-extracted blood vessel region in S2.1;
s2.3 construction of a vessel filtering weight map
For each blood vessel point extracted in S2.1, a blood vessel filtering value can be calculated through S2.2, and a blood vessel filtering weight value of the blood vessel point is defined as:
wherein the content of the first and second substances,
and
respectively calculating the reciprocal of the blood vessel filtering value calculated in S2.2 for two adjacent points, and Dist is the physical distance between the two points;
and constructing a blood vessel filtering value weight map according to the calculated blood vessel filtering weights of all the blood vessel points.
5. The method of claim 4 for automatically extracting the central path of the head and neck blood vessels in the CTA image, wherein: s2.2, based on Raycasting multi-scale gradient normalization blood vessel filtering, calculating a blood vessel filtering value of each pixel point in a pre-extracted blood vessel region in S2.1, and specifically comprising the following steps:
1) Since the cross-sectional dimension of the blood vessel varies during its courseCalculating the filtered value of the blood vessel under different radius scales as required, wherein R belongs to (R) min ,R max )；
2) Taking any one blood vessel region point pre-extracted in S2.1 as a circle center, respectively calculating the projection radius scale as R and the gradient response V of the circle center point R Minimum V of gradient response on the 0 to R radius scale R,min ，R min To R max Maximum value V at radius scale R,max ；
3) Defining the normalized gradient response at the radius scale R as E R ＝(V R -V R,min )/V R,max If the current point is the center point of the blood vessel, E R Approximately equal to 1; if it is a non-vascular point, E R About 0;
4) Repeating steps 2) -3), calculating a plurality of projection radius scales R min To R max E of R And will be at different scales E R The maximum value of (a) is used as the final blood vessel filtering value at the point.
6. The method of claim 5 for automatically extracting the central path of the head and neck blood vessels in the CTA image, wherein: the step S3: in the vessel filtering weight map constructed in step S2, the method for extracting the central paths of the carotid artery and the vertebral artery by using Dijkstra optimal path extraction algorithm is as follows:
s3.1: defining an energy function between the start and end points of the blood vessels located in step S1:
E＝∫(w(s)+ε)ds
wherein w (S) is the blood vessel filtering enhancement weight calculated in S2.3, epsilon is a regular term, and S is a path between a starting point and an end point;
and S3.2, selecting the path with the minimum blood vessel filtering weight sum in all paths between the starting point and the end point of the blood vessel as the blood vessel central path.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910693725.6A CN110338830B (en) | 2019-07-30 | 2019-07-30 | Method for automatically extracting head and neck blood vessel central path in CTA image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910693725.6A CN110338830B (en) | 2019-07-30 | 2019-07-30 | Method for automatically extracting head and neck blood vessel central path in CTA image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110338830A CN110338830A (en) | 2019-10-18 |
| CN110338830B true CN110338830B (en) | 2022-12-13 |
Family
ID=68179062
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910693725.6A Active CN110338830B (en) | 2019-07-30 | 2019-07-30 | Method for automatically extracting head and neck blood vessel central path in CTA image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110338830B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113554742B (en) * | 2020-04-26 | 2024-02-02 | 上海联影医疗科技股份有限公司 | Three-dimensional image reconstruction method, device, equipment and storage medium |
| CN111815587A (en) * | 2020-06-29 | 2020-10-23 | 苏州润心医疗器械有限公司 | Method and system for picking up points on aorta centerline based on CT sequence image |
| CN112308846B (en) * | 2020-11-04 | 2021-07-13 | 赛诺威盛科技(北京)股份有限公司 | Blood vessel segmentation method and device and electronic equipment |
| CN112837288B (en) * | 2021-02-01 | 2021-11-23 | 数坤(北京)网络科技股份有限公司 | Blood vessel centerline extraction method and device and readable storage medium |
| CN113689388B (en) * | 2021-08-03 | 2024-02-06 | 慧影医疗科技(北京)股份有限公司 | Three-dimensional center line starting point positioning method and device for aortic curved surface reconstruction |
| CN114972859A (en) * | 2022-05-19 | 2022-08-30 | 推想医疗科技股份有限公司 | Pixel classification method, model training method, device, equipment and medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017116074A1 (en) * | 2015-12-29 | 2017-07-06 | 서울대학교산학협력단 | Method and apparatus for measuring ct image atherosclerotic plaque volume using adaptive threshold value image segmentation |
| CN108133478A (en) * | 2018-01-11 | 2018-06-08 | 苏州润心医疗器械有限公司 | A kind of method for extracting central line of coronary artery vessel |
| CN109035194A (en) * | 2018-02-22 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | A kind of vessel extraction method and device |
| CN109448004A (en) * | 2018-10-26 | 2019-03-08 | 强联智创（北京）科技有限公司 | A kind of intercept method and system of the intracranial vessel image based on center line |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050256400A1 (en) * | 2003-12-03 | 2005-11-17 | Bhargav Raman | Method to identify arterial and venous vessels |
| US9727968B2 (en) * | 2012-08-13 | 2017-08-08 | Koninklijke Philips N.V. | Tubular structure tracking |
-
2019
- 2019-07-30 CN CN201910693725.6A patent/CN110338830B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017116074A1 (en) * | 2015-12-29 | 2017-07-06 | 서울대학교산학협력단 | Method and apparatus for measuring ct image atherosclerotic plaque volume using adaptive threshold value image segmentation |
| CN108133478A (en) * | 2018-01-11 | 2018-06-08 | 苏州润心医疗器械有限公司 | A kind of method for extracting central line of coronary artery vessel |
| CN109035194A (en) * | 2018-02-22 | 2018-12-18 | 青岛海信医疗设备股份有限公司 | A kind of vessel extraction method and device |
| CN109448004A (en) * | 2018-10-26 | 2019-03-08 | 强联智创（北京）科技有限公司 | A kind of intercept method and system of the intracranial vessel image based on center line |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110338830A (en) | 2019-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110338830B (en) | Method for automatically extracting head and neck blood vessel central path in CTA image | |
| Grefkes et al. | Human somatosensory area 2: observer-independent cytoarchitectonic mapping, interindividual variability, and population map | |
| CN108335304B (en) | Aortic aneurysm segmentation method of abdominal CT scanning sequence image | |
| CN108961261B (en) | Optic disk region OCT image hierarchy segmentation method based on space continuity constraint | |
| US8139837B2 (en) | Bone number determination apparatus and recording medium having stored therein program | |
| CN106846346B (en) | Method for rapidly extracting pelvis outline of sequence CT image based on key frame mark | |
| Septiarini et al. | AUTOMATIC GLAUCOMA DETECTION BASED ON THE TYPE OF FEATURES USED: A REVIEW. | |
| CN111374712B (en) | Ultrasonic imaging method and ultrasonic imaging equipment | |
| CN114757960B (en) | Tooth segmentation and reconstruction method based on CBCT image and storage medium | |
| CN116982077A (en) | System and method for automatically segmenting patient-specific anatomy for pathology-specific measurements | |
| CN106780491B (en) | Initial contour generation method adopted in segmentation of CT pelvic image by GVF method | |
| CN108986114A (en) | A kind of abdominal CT sequence image liver automatic division method based on level set and shape description symbols | |
| CN112132837A (en) | Chest bone automatic extraction method and system, electronic equipment and storage medium | |
| CN106780492B (en) | Method for extracting key frame of CT pelvic image | |
| WO2005048844A1 (en) | Automatic identification of the anterior and posterior commissure landmarks | |
| US20150206300A1 (en) | Method and apparatus for extraction and quantification of hematoma from a brain scan such as computed tomography data | |
| Tan et al. | Automatic feature extraction of optical coherence tomography for lamina cribrosa detection | |
| Zhang et al. | Automatic rib segmentation in chest CT volume data | |
| Memiş et al. | Computerized 2D detection of the multiform femoral heads in magnetic resonance imaging (MRI) sections with the integro-differential operator | |
| Zhao et al. | Auto-threshold bone segmentation based on CT image and its application on CTA bone-subtraction | |
| CN110378868B (en) | Pericardium segmentation method based on CT sequence and method for roughly positioning pericardium area from pericardium sequence center layer slice | |
| US11481905B2 (en) | Atlas for automatic segmentation of retina layers from OCT images | |
| Bozkurt et al. | Comparison of seeded region growing and random walk methods for vessel and bone segmentation in CTA images | |
| JP2004174219A (en) | Apparatus and method for processing image and recording medium for storing program used for causing computer to execute the method | |
| JP2004174217A (en) | Image processing apparatus, image processing method, and recording medium which stores program according to which computer perform pertinent image processing |
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 | ||
| CB02 | Change of applicant information |
Address after: 100176 1st floor, building 8, 11 Kangding street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant after: Sinovision Technology (Beijing) Co.,Ltd. Address before: 100176 1st floor, building 8, 11 Kangding street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Applicant before: SAINUO WEISHENG TECHNOLOGY (BEIJING) Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |