Next Article in Journal
Significance of Postprocedural Contrast Medium Injection after CT-Guided Abscess Drainage
Previous Article in Journal
Doppler Examination of the Testicular Artery of Beagle-Breed Dogs from Birth to Puberty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Path-Finding Algorithms on the Labeling of the Centerlines of Circle of Willis Arteries

Division of Digital Healthcare, College of Software and Digital Healthcare Convergence, Yonsei University, Wonju 26493, Republic of Korea
*
Author to whom correspondence should be addressed.
Submission received: 19 May 2023 / Revised: 21 July 2023 / Accepted: 21 July 2023 / Published: 24 July 2023

Abstract

:
Quantitative analysis of intracranial vessel segments typically requires the identification of the vessels’ centerlines, and a path-finding algorithm can be used to automatically detect vessel segments’ centerlines. This study compared the performance of path-finding algorithms for vessel labeling. Three-dimensional (3D) time-of-flight magnetic resonance angiography (MRA) images from the publicly available dataset were considered for this study. After manual annotations of the endpoints of each vessel segment, three path-finding methods were compared: (Method 1) depth-first search algorithm, (Method 2) Dijkstra’s algorithm, and (Method 3) A* algorithm. The rate of correctly found paths was quantified and compared among the three methods in each segment of the circle of Willis arteries. In the analysis of 840 vessel segments, Method 2 showed the highest accuracy (97.1%) of correctly found paths, while Method 1 and 3 showed an accuracy of 83.5% and 96.1%, respectively. The AComm artery was highly inaccurately identified in Method 1, with an accuracy of 43.2%. Incorrect paths by Method 2 were noted in the R-ICA, L-ICA, and R-PCA-P1 segments. The Dijkstra and A* algorithms showed similar accuracy in path-finding, and they were comparable in the speed of path-finding in the circle of Willis arterial segments.

1. Introduction

Abnormal morphological characteristics of the cerebral blood vessels are associated with vascular diseases. For example, stenosis of any artery due to intracranial atherosclerosis can cause ischemia in the brain tissue, resulting in stroke and other cognitive brain disorders [1,2]. Development of intracranial aneurysms [3,4] may be associated with high blood pressure on the weakened vessel wall and potentially could result in ruptures and cerebral hemorrhages if left untreated [5]. Three-dimensional (3D) computed tomography angiography (CTA) and 3D time-of-flight MR angiography (TOF-MRA) are commonly used to noninvasively obtain data on the cerebral vessels [6,7]. The geometric information such as vessel tortuosity [8] and distributions of the diameters of the vessels [9] is extracted and evaluated to study correlations between vessel abnormalities and diseases. Quantitative evaluation based on the geometry of the cerebral arteries has been reported in the literature [9,10,11]. For example, a cumulative distribution function of the cross-sectional diameters in a vessel segment was evaluated [9]. Tortuosity descriptors such as the sum of angle metrics, inflection count metric, triangular index, relative length, and product of angle distance were used to find the relationship between diseases and morphological features [12,13]. These quantities are typically calculated on a vessel segment after identifying the vessel’s centerline, and a path-finding algorithm can automate the identification of the centerline.
Labeling of the arteries in the circle of Willis (CoW) has been of interest to researchers since the CoW is essential in maintaining the circulation of blood flow even in the occlusion of an artery and has important bifurcations that anatomically separate the intracranial arteries of interest [14,15,16]. For example, the internal carotid artery (ICA) branches into the anterior cerebral artery (ACA) and middle cerebral artery (MCA). In the posterior circulation, the basilar artery (BA) is connected to the left and right posterior cerebral arteries (PCA). Tracking of the vessel’s centerline can be performed via a path-finding algorithm after manual annotations of two endpoints of a vessel segment. This procedure is relatively simple and straightforward when compared with deep learning-based segmentation of the intracranial vessel segments [17]. The development of deep learning segmentation models requires manual segmentations of individual arterial segments for the generation of training data. The manual segmentation process involves slice-wise manual tracing in a 3D volume consisting of hundreds of slices, and it is thus time-consuming and laborious [18]. Hence, the centerline tracking approaches can improve the efficiency of vessel labeling and its subsequent quantification without the need for deep learning model training data generation. Previous studies have demonstrated methods adopting path-finding algorithms such as Dijkstra’s algorithm in analyzing 3D angiography data [19,20,21]. However, they did not demonstrate detailed comparisons of accuracy among available path-finding algorithms in vessel image analysis.
In this study, we evaluated the performance of three path-finding algorithms in robustly identifying the cerebral arterial segments in the CoW. In terms of path-finding accuracy and computational time, we compared a depth-first search (DFS) based algorithm that does not involve a graph structure to a Dijkstra algorithm and an A* (A star) algorithm, both of which require a graph representation.

2. Materials and Methods

Figure 1 illustrates the flowchart of the presented vessel labeling methods. The preprocessing steps involving Otsu thresholding and 3D seeded region growing have been implemented in Matlab version 9.13 (The Mathworks Inc., Natick, MA, USA). Skeletonization, graph structure generation, path-finding methods, and visualization of the centerlines in 3D were implemented in Python version 3.10 (Python Software Foundation).

2.1. Data

We used publicly available magnetic resonance angiography (MRA) data from the IXI Dataset (https://brain-development.org/ixi-dataset) (accessed on 21 July 2023). Sixty Neuroimaging Informatics Technology Initiative (NIfTI) files were considered to evaluate the performance of path-finding algorithms. The imaging parameters are as follows: Philips Medical Systems Intera 3T, repetition time = 16.7 ms, echo time = 5.8 ms, number of phase-encoding steps = 286, reconstruction diameter = 240 × 240 mm2, acquisition matrix = 288 × 286, flip angle = 16-deg, in-plane pixel spacing = 0.4–0.6 mm, spacing between slices = 0.8 mm. We applied the bi-cubic interpolation along the slice dimension to generate iso-resolution image data with the same pixel spacing in all three directions. After the interpolation, the final voxel size was isotropic with the 0.6 mm voxel spacing.

2.2. Vessel Segmentation and Skeletonization

The user selected an axial slice that shows three cross-sections of the vessels, which are the right internal carotid artery (ICA), left ICA, and basilar artery (BA). Supplementary Figure S1 shows an example axial slice showing three vessels’ cross-sections. A seed point was manually annotated by the mouse click within each of the three cross-sections (Supplementary Figure S1). With each seed point, 3D seeded region growing was performed using the region-growing function [22] available in Matlab. The segmented arteries after the region growing included the region of the CoW. The three segmentation results were obtained in binary masks, and the union set operation was performed to combine the three segmented masks into one final binary mask of the vessels. Skeletonization was performed to find the centerlines of the vessels using the scikit-image Python library [23].

2.3. Annotation of Two Endpoints

Following the anatomical terminology provided by Dumais et al. [17], we identified 14 vessel segments by manually annotating two endpoints in each vessel segment. The vessel segments of interest were (1) AComm (anterior communicating artery), (2) R-A1 (right anterior cerebral artery A1), (3) L-A1 (left anterior cerebral artery A1), (4) R-M1 (right middle cerebral artery M1), (5) L-M1 (left middle cerebral artery M1), (6) R-ICA (right internal carotid artery), (7) L-ICA (left internal carotid artery), (8) R-PComm (right posterior communicating artery), (9) L-PComm (left posterior communicating artery), (10) R-P1 (right posterior cerebral artery P1), (11) L-P1 (left posterior cerebral artery P1), (12) R-P2 (right posterior cerebral artery P2), (13) L-P2 (left posterior cerebral artery P2), and (14) BA (basilar artery). Some endpoints were shared among the vessel segments. For example, the R-A1 and R-M1 segments share a common branch point (Figure 2). The same is true for the BA, R-P1, and L-P1 segments.
The Plotly (Plotly Technologies Inc., Montreal, QC, Canada) Python library (https://plotly.com/python) (accessed on 21 July 2023) was used to visually identify the vessel segments in the 3D space. The mouse hovering on the target endpoint shows its position information in 3D coordinates. For each vessel segment, the position information of the two endpoints was the input to our path-finding algorithms. During manual annotation, we recorded the position information of the two endpoints in vessel segments. If no centerline exists in a certain vessel segment, we recorded ‘0 0 0′ for the vessel segment. The position information was saved in a text file with .txt extension for each subject. The visualization result of the colored vessel segments was saved in a .html file, which was reserved for further investigation to check the accuracy of the vessel annotations.

2.4. Path-Finding Algorithms

In this study, we implemented three path-finding algorithms for comparisons: (Method 1) DFS-based path-finding algorithm, (Method 2) Dijkstra’s algorithm, and (Method 3) A* algorithm.
For Method 1, we implemented a maze-solving algorithm based on DFS (Supplementary Figure S2). Given two endpoints, from the start point, Method 1 attempted to find the next neighbor pixel and record the visited pixel locations using the push operation in a stack. When there is no neighboring pixel anymore, the visited locations were taken out of the stack until the algorithm found the neighboring pixel which was not visited. The algorithm was terminated when the endpoint was reached. The remaining data containing pixel locations in the stack indicate the coordinates along the path that connects the two endpoints. We implemented the method in C++ on the Microsoft Visual Studio environment. PyBind11 (https://github.com/pybind/pybind11) (accessed on 21 July 2023) was used to create Python bindings to the C++-compiled code [24]. Hence, the 3D coordinates along the path were able to be obtained and saved in a variable in Python.
Methods 2 and 3 require the generation of a graph structure from the skeleton vessel image as shown in Figure 1. To generate a graph from the skeleton image, we used the Skan version 0.10.0 Python library (https://skeleton-analysis.org/stable) (accessed on 21 July 2023) [25]. Dijkstra’s algorithm is well known as the shortest path-finding algorithm in a graph structure [26]. In our implementation, the pixels along the centerlines represented vertices (or nodes), and the connections with their neighboring pixels represented edges. The weights between the two connected nodes were calculated as the Euclidean distance of the two points. We used the Dijkstra() function provided by the SciPy Python library [27]. The Dijkstra() function is based on the Fibonacci heap implementation, which has the time complexity of O(E + V logV) (here, E is the number of edges, and V is the number of vertices) and is more time efficient than the list implementation whose time complexity is O( V 2 ) [28]. A* algorithm is another method for the shortest path search in a graph structure. In contrast to Dijkstra’s algorithm, which finds the shortest path from a starting point to all goal points, the A* algorithm only finds the shortest path from a starting point to a destination point, and it introduces heuristic cost values as well as graph weights to find the next node [29]. In our study, we used and modified the code available from the python-astar open-source library (https://github.com/jrialland/python-astar) (accessed on 21 July 2023). The heuristic cost value was calculated by computing the Euclidean distance between the nodes of interest and the destination node.

2.5. Evaluation

We compared the three methods in terms of path-finding accuracy and computational time in Python. For each vessel segment, we provided the two endpoints and let the methods automatically find a path. We counted the number of incorrect paths and the number of correct paths for each vessel segment after building consensus on the path-finding correctness. An incorrect path was defined as a path that connects the two endpoints but is not the path we expected. To measure computational time, we used the time.perf_counter() function. Since we used the already annotated endpoints, during the time measurements there were no manual annotation procedures for recording two endpoints in all the vessel segments. For Method 1, we measured the time interval of the DFS path-finding algorithm taken to find paths in all segments in each subject. For Methods 2 and 3, we first measured the time interval of the graph structure generation from the skeleton image and then measured each time interval of the path-finding algorithm (i.e., SciPy’s Dijkstra() function for Method 2 and python-astar’s find_path() function for Method 3). For each method of Methods 2 and 3, we summed the time intervals including the graph structure generation and the path finding. These time intervals were compared to the time measured from Method 1 for each subject’s data. The evaluation was performed on a Windows PC (13th Generation Intel® Core™ i7-13700K 16-Core Processor Central Processing Unit).
A Fisher’s exact test was used to analyze any differences between the path-finding methods in detecting the correct centerlines of the vessel segments. A two-sample unpaired Student’s t-test was performed to determine if the computational time between the two methods was significantly different. A p-value of <0.05 was considered statistically significant.

3. Results

Table 1 summarizes the accuracy of the three algorithms in finding the paths between two endpoints in each vessel segment. Method 1 produced 121 incorrect paths out of 735 paths, and Method 2 produced only 21 incorrect paths out of 735 paths. Method 3 produced 29 incorrect paths. The accuracy (97.1%) of Method 2 was significantly higher than that (83.5%) of Method 1. The accuracy of Method 3 was 96.1%. The AComm is highly inaccurate in Method 1 with an accuracy of 43.2%. Method 2 is superior to Method 1 in all vessel segments except for the R-ICA segment. Method 3 is comparable to Method 2 in every vessel segment, except for the AComm segment. Fisher’s exact test resulted in statistically significant associations in most vessel segments in Method 1 and Method 2. As shown in Table 1, AComm, ACA A1, PComm, L-PCA P1, L-PCA P2, and BA showed statistically significant associations, while MCA M1, R-PCA P1, R-PCA P2, and ICA were not statistically significant (p > 0.05). This implies that the choice of a path-finding method affected the correct detection of paths in AComm, ACA A1, PComm, L-PCA P1, L-PCA P2, and BA. The same was true in the association of Method 1 with Method 3. However, there was no statistically significant association between Method 2 and Method 3, except for AComm.
The seeded region-growing algorithm was advantageous because it did not segment unwanted vessels such as veins or outer intracranial vessels. The inclusion of the veins or outer intracranial vessels can obscure the visual appearances of the CoW arteries (Supplementary Figure S3) and make it challenging to locate the two endpoints when identifying a vessel segment in the CoW. The undetected paths in Table 1 resulted from a non-existent binary vessel mask possibly due to either under-segmentation of the seeded region-growing algorithm or non-enhancement of the vessel itself, and thus the path-finding algorithms were not able to be applied to the vessel segments. AComm had 16 undetected paths (i.e., 26.7% of the total subjects). R-PComm and L-PComm had 46 and 42 undetected paths (i.e., 76.7% and 70.0% of the total subjects), respectively.
Computational time was compared among the three methods in all subjects’ data (Supplementary Table S1). To isolate the path-finding algorithms from the manual annotation, we loaded the text file which contained all the manually annotated endpoints of the vessel segments, and we focused on measuring time on the path-finding procedures subject by subject. Method 1 took an average (±standard deviation) of 1489.4 (±191.4) ms. Method 2 took shorter than Method 1 with an average (±standard deviation) of 458.2 (±63.4) ms. Method 3 took an average (±standard deviation) of 458.0 (±63.4) ms, which was comparable to Method 2. The difference between Method 1 and Method 2 (or Method 3) was statistically significant (p < 0.0001), while the difference between Method 2 and Method 3 was not statistically significant (p > 0.4). In Methods 2 and 3, the conversion of the skeleton to the graph structure was the main bottleneck and took 98.95% and 99.00% of the total computational time, respectively. Moreover, it is important to note that manual annotation of the endpoints required the user to rotate and zoom in and out of the 3D skeletal vessel for correct identification of the vessel segments and took approximately 10 min per subject.

4. Discussion

Path finding is essential for annotating and quantifying vessel segments and has the potential to be useful in planning invasive procedures such as mechanical thrombectomy. This current study compared the performance of three path-finding algorithms in the total 840 arterial segments in the CoWs. Method 1 does not require a conversion of the skeleton’s centerlines to the graph representation and is based on the DFS algorithm. Since it does not construct a graph representation, it is deemed as rather simple, when compared to Methods 2 and 3, which require the construction of a graph and utilize the shortest path-finding Dijkstra and A* algorithms, respectively. Our path-finding results indicate that Method 1 is prone to errors when there is a loop in the vessel’s centerlines (Figure 3). The shortest path-finding algorithms can avoid such errors in case of the existence of a loop.
There are two possible reasons for the incorrect path-finding results. First, when the path forms a loop, Method 1 may not find the shortest path and instead can take a path that leads to the other endpoint regardless of the path’s length. The ‘loop’ path was frequent in the AComm segments, where there are alternative paths that detour via A2 segments. Second, when there are multiple paths between the two endpoints, Method 2 always finds the shortest path, which is sometimes not the correct path (Figure 4). We often detected an incorrect path by comparing it with a symmetrical vessel anatomy located in the other hemisphere of the brain. Additional routes in the main arteries may be attributed to either incorrect segmentation due to noise in the image or errors in the skeletonization process. By improving segmentation performance or denoising the gray-scale TOF-MRA images, the incorrect path-finding problem may be resolved.
Compared to the A* algorithm, Dijkstra’s algorithm finds the shortest paths for all nodes, although in our case we need only a single shortest path between the initial and destination nodes. Hence, the Dijkstra algorithm may not be the most time-efficient for path finding by nature. Since the A* algorithm only cares about finding the shortest path between the two endpoints, it is theoretically more time-efficient than the Dijkstra algorithm. However, in our evaluation, there were little or no noticeable improvements in the A* algorithm over the Dijkstra algorithm in terms of computational time, which was not expected given the inherent advantage in time efficiency in the A* algorithm.
The seeded region-growing algorithm used in this study provided under-segmentation results. This under-segmentation occurred dominantly in the PComm segments. Since we focused on the evaluation of the path-finding methods, it is a fair comparison, although there was an under-segmentation issue in several vessel segments. However, there is room for improvement in automatic segmentation. One way would be to use advanced segmentation methods based on encoder–decoder deep convolutional neural network architectures [30,31] or a multiscale image analysis approach [32]. Our study was based on MRA images, but it is possible to extract vessel segments from CTA images. Notably, a direct application of our region-growing segmentation to CTA images would be challenging because the region-growing method is based on the similarity of intensity in the blood, and contrast-enhanced arteries may not be distinguished from nearby bone structures in CTA images [33,34]. The errors in vessel segmentation are likely to occur especially in the regions where the arteries and bones are very close (e.g., internal carotid arteries, vertebral arteries, etc.). For CTA images, recent studies proposed automatic deep learning-based vessel segmentation methods to improve the accuracy of vessel segmentation [18,35].
Manual annotation of the endpoints in 3D space is tedious and can be prone to errors if not carefully checked. As such, it would be intriguing to investigate the feasibility of automatic landmark localization methods [36,37,38] to localize the vessel segments’ endpoints which we manually annotated in this study. Landmark localization of the two endpoints may also be helpful in automatically detecting occlusions of the arterial segments in ischemic stroke patients’ data.
This study has several limitations. First, the number of subjects used for the analysis is small. However, increasing the number of subjects would not significantly affect the study outcome, as there are already more than 10 segments for each subject’s MRA image data. Second, we did not perform intrarater or interrater variability in evaluating the correctness of the paths. Third, the image segmentation quality has imperfections since it showed under-segmentation results, especially in the PComm segments. However, improving the segmentation results may not significantly affect the outcome of the study because all three methods used the same segmented binary vessel masks for the analysis. Fourth, we did not consider disconnected vessel segments, which may result from the occlusions or severe stenoses of arterial segments. The development of a method for identifying disconnected vessel segments would be important for automatically detecting the vessel’s occlusion. Last but not least, we only considered arterial segments in the CoW arteries. Extending the analysis to other intracranial vessel segments such as ACA A2 and MCA M2 would be an important venue for further research.

5. Conclusions

Our study indicates that a graph representation of the centerlines of cerebral arteries in the CoW is advantageous in semi-automatically labeling the cerebral arteries when it is used with the shortest path-finding algorithms such as the Dijkstra algorithm and the A* algorithm. Among the three path-finding algorithms, the Dijkstra algorithm resulted in the highest accuracy of 97.1%, and it was slightly higher in accuracy than the A* algorithm. When the DFS approach was used, the incorrect path-finding results mostly occurred for the vessel segments with multiple paths connecting the two endpoints. Path-finding results for the Dijkstra and A* algorithms were incorrect when the paths were visually compared with the paths in the other brain hemisphere. Full automatization of vessel segmentation and landmark localization with deep learning would be desirable in order to increase time efficiency and avoid manual annotation procedures.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/tomography9040113/s1, Figure S1: A slice used for the selection of three seed points. The cross-sections of the right internal carotid artery (orange), the left internal carotid artery (blue), and the basilar artery (green) are shown. The seed locations for seeded region growing are indicated by the arrows; Figure S2: The pseudo-code of the depth first search (DFS) based maze solving algorithm. This algorithm is referred to as Method 1; Figure S3: Seeded region growing. (a) A vessel segmentation result after the Otsu’s thresholding. Disconnected vessels are present, which sometimes obscure the view of the artery of interests. (b) A vessel segmentation result after applying the seeded region growing algorithm to the Otsu-thresholded vessel. Disconnected vessels disappear, and vessel visualization is improved; Table S1: Comparison of computational time measurements.

Author Contributions

Conceptualization, Y.-C.K.; methodology, S.-O.K. and Y.-C.K.; software, S.-O.K. and Y.-C.K.; validation, S.-O.K.; formal analysis, S.-O.K.; investigation, S.-O.K.; resources, S.-O.K.; data curation, S.-O.K. and Y.-C.K.; writing—original draft preparation, Y.-C.K.; writing—review and editing, Y.-C.K.; visualization, S.-O.K. and Y.-C.K.; supervision, Y.-C.K.; project administration, Y.-C.K.; funding acquisition, Y.-C.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the “Regional Innovation Strategy (RIS)” through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (MOE) (2022RIS-005).

Institutional Review Board Statement

Ethical review and approval were waived for this study due to the use of publicly available data.

Informed Consent Statement

Patient consent was waived due to the use of publicly available data.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ritz, K.; Denswil, N.P.; Stam, O.C.; van Lieshout, J.J.; Daemen, M.J. Cause and mechanisms of intracranial atherosclerosis. Circulation 2014, 130, 1407–1414. [Google Scholar] [CrossRef] [Green Version]
  2. Hu, X.; De Silva, T.M.; Chen, J.; Faraci, F.M. Cerebral Vascular Disease and Neurovascular Injury in Ischemic Stroke. Circ. Res. 2017, 120, 449–471. [Google Scholar] [CrossRef] [PubMed]
  3. Krzyzewski, R.M.; Klis, K.M.; Kwinta, B.M.; Gackowska, M.; Gasowski, J. Increased tortuosity of ACA might be associated with increased risk of ACoA aneurysm development and less aneurysm dome size: A computer-aided analysis. Eur. Radiol. 2019, 29, 6309–6318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jeong, W.; Rhee, K. Hemodynamics of cerebral aneurysms: Computational analyses of aneurysm progress and treatment. Comput. Math. Methods Med. 2012, 2012, 782801. [Google Scholar] [CrossRef] [Green Version]
  5. Magid-Bernstein, J.; Girard, R.; Polster, S.; Srinath, A.; Romanos, S.; Awad, I.A.; Sansing, L.H. Cerebral Hemorrhage: Pathophysiology, Treatment, and Future Directions. Circ. Res. 2022, 130, 1204–1229. [Google Scholar] [CrossRef]
  6. Cirillo, M.; Scomazzoni, F.; Cirillo, L.; Cadioli, M.; Simionato, F.; Iadanza, A.; Kirchin, M.; Righi, C.; Anzalone, N. Comparison of 3D TOF-MRA and 3D CE-MRA at 3T for imaging of intracranial aneurysms. Eur. J. Radiol. 2013, 82, e853–e859. [Google Scholar] [CrossRef] [PubMed]
  7. Lell, M.M.; Anders, K.; Uder, M.; Klotz, E.; Ditt, H.; Vega-Higuera, F.; Boskamp, T.; Bautz, W.A.; Tomandl, B.F. New techniques in CT angiography. Radiographics 2006, 26 (Suppl. 1), S45–S62. [Google Scholar] [CrossRef]
  8. Han, H.C. Twisted blood vessels: Symptoms, etiology and biomechanical mechanisms. J. Vasc. Res. 2012, 49, 185–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kandil, H.; Soliman, A.; Ghazal, M.; Mahmoud, A.; Shalaby, A.; Keynton, R.; Elmaghraby, A.; Giridharan, G.; El-Baz, A. A Novel Framework for Early Detection of Hypertension using Magnetic Resonance Angiography. Sci. Rep. 2019, 9, 11105. [Google Scholar] [CrossRef] [Green Version]
  10. Kim, H.J.; Song, H.N.; Lee, J.E.; Kim, Y.C.; Baek, I.Y.; Kim, Y.S.; Chung, J.W.; Jee, T.K.; Yeon, J.Y.; Bang, O.Y.; et al. How Cerebral Vessel Tortuosity Affects Development and Recurrence of Aneurysm: Outer Curvature versus Bifurcation Type. J. Stroke 2021, 23, 213–222. [Google Scholar] [CrossRef]
  11. Klis, K.M.; Krzyzewski, R.M.; Kwinta, B.M.; Lasocha, B.; Brzegowy, P.; Stachura, K.; Popiela, T.J.; Borek, R.; Gasowski, J. Increased tortuosity of basilar artery might be associated with higher risk of aneurysm development. Eur. Radiol. 2020, 30, 5625–5632. [Google Scholar] [CrossRef] [PubMed]
  12. Sodi, A.; Guarducci, M.; Vauthier, L.; Ioannidis, A.S.; Pitz, S.; Abbruzzese, G.; Sofi, F.; Mecocci, A.; Miele, A.; Menchini, U. Computer assisted evaluation of retinal vessels tortuosity in Fabry disease. Acta Ophthalmol. 2013, 91, e113–e119. [Google Scholar] [CrossRef] [PubMed]
  13. Klis, K.M.; Krzyzewski, R.M.; Kwinta, B.M.; Stachura, K.; Gasowski, J. Tortuosity of the Internal Carotid Artery and Its Clinical Significance in the Development of Aneurysms. J. Clin. Med. 2019, 8, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bogunovic, H.; Pozo, J.M.; Cardenes, R.; San Roman, L.; Frangi, A.F. Anatomical labeling of the Circle of Willis using maximum a posteriori probability estimation. IEEE Trans. Med. Imaging 2013, 32, 1587–1599. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, L.; Mossa-Basha, M.; Balu, N.; Canton, G.; Sun, J.; Pimentel, K.; Hatsukami, T.S.; Hwang, J.N.; Yuan, C. Development of a quantitative intracranial vascular features extraction tool on 3D MRA using semiautomated open-curve active contour vessel tracing. Magn. Reson. Med. 2018, 79, 3229–3238. [Google Scholar] [CrossRef]
  16. Robben, D.; Turetken, E.; Sunaert, S.; Thijs, V.; Wilms, G.; Fua, P.; Maes, F.; Suetens, P. Simultaneous segmentation and anatomical labeling of the cerebral vasculature. Med. Image Anal. 2016, 32, 201–215. [Google Scholar] [CrossRef] [Green Version]
  17. Dumais, F.; Caceres, M.P.; Janelle, F.; Seifeldine, K.; Ares-Bruneau, N.; Gutierrez, J.; Bocti, C.; Whittingstall, K. eICAB: A novel deep learning pipeline for Circle of Willis multiclass segmentation and analysis. Neuroimage 2022, 260, 119425. [Google Scholar] [CrossRef]
  18. Nazir, A.; Cheema, M.N.; Sheng, B.; Li, H.T.; Li, P.; Yang, P.; Jung, Y.Y.; Qin, J.; Kim, J.M.; Feng, D.D. OFF-eNET: An Optimally Fused Fully End-to-End Network for Automatic Dense Volumetric 3D Intracranial Blood Vessels Segmentation. IEEE Trans. Image Process. 2020, 29, 7192–7202. [Google Scholar] [CrossRef]
  19. Suran, S.; Pattanaik, V.; Malathi, D. Discovering shortest path between points in cerebrovascular system. In Proceedings of the 6th IBM Collaborative Academia Research Exchange Conference (I-CARE) on I-CARE 2014, Bangalore, India, 9–11 October 2014; pp. 1–3. [Google Scholar]
  20. Shen, M.; Wei, J.; Fan, J.; Tan, J.; Wang, Z.; Yang, Z.; Qiao, P.; Liao, F. Automatic cerebral artery system labeling using registration and key points tracking. In Proceedings of the Knowledge Science, Engineering and Management: 13th International Conference, KSEM 2020, Hangzhou, China, 28–30 August 2020; pp. 355–367. [Google Scholar]
  21. Thamm, F.; Jurgens, M.; Taubmann, O.; Thamm, A.; Rist, L.; Ditt, H.; Maier, A. An algorithm for the labeling and interactive visualization of the cerebrovascular system of ischemic strokes. Biomed. Phys. Eng. Express 2022, 8, 065016. [Google Scholar] [CrossRef]
  22. Kroon, D.-J. Region Growing, MATLAB Central File Exchange; The MathWorks, Inc.: Natick, MA, USA, 2021. [Google Scholar]
  23. Van der Walt, S.; Schönberger, J.L.; Nunez-Iglesias, J.; Boulogne, F.; Warner, J.D.; Yager, N.; Gouillart, E.; Yu, T. Scikit-Image: Image processing in Python. PeerJ 2014, 2, e453. [Google Scholar] [CrossRef]
  24. Kim, Y.C.; Kim, K.R.; Lee, H.; Choe, Y.H. Fast calculation software for modified Look-Locker inversion recovery (MOLLI) T1 mapping. BMC Med. Imaging 2021, 21, 26. [Google Scholar] [CrossRef] [PubMed]
  25. Nunez-Iglesias, J.; Blanch, A.J.; Looker, O.; Dixon, M.W.; Tilley, L. A new Python library to analyse skeleton images confirms malaria parasite remodelling of the red blood cell membrane skeleton. PeerJ 2018, 6, e4312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lee, K.D.; Lee, K.D.; Steve Hubbard, S.H. Data Structures and Algorithms with Python; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  27. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
  28. Schrijver, A. Combinatorial Optimization: Polyhedra and Efficiency; Springer: Berlin/Heidelberg, Germany, 2003; Volume 24. [Google Scholar]
  29. Russell, S.J. Artificial Intelligence a Modern Approach; Pearson Education, Inc.: London, UK, 2010. [Google Scholar]
  30. Chen, Y.; Jin, D.; Guo, B.; Bai, X. Attention-Assisted Adversarial Model for Cerebrovascular Segmentation in 3D TOF-MRA Volumes. IEEE Trans. Med. Imaging 2022, 41, 3520–3532. [Google Scholar] [CrossRef]
  31. Livne, M.; Rieger, J.; Aydin, O.U.; Taha, A.A.; Akay, E.M.; Kossen, T.; Sobesky, J.; Kelleher, J.D.; Hildebrand, K.; Frey, D.; et al. A U-Net Deep Learning Framework for High Performance Vessel Segmentation in Patients With Cerebrovascular Disease. Front. Neurosci. 2019, 13, 97. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wozniak, T.; Strzelecki, M.; Majos, A.; Stefanczyk, L. 3D vascular tree segmentation using a multiscale vesselness function and a level set approach. Biocybern. Biomed. Eng. 2017, 37, 66–77. [Google Scholar] [CrossRef]
  33. Lell, M.M.; Ruehm, S.G.; Kramer, M.; Panknin, C.; Habibi, R.; Klotz, E.; Villablanca, P. Cranial computed tomography angiography with automated bone subtraction: A feasibility study. Investig. Radiol. 2009, 44, 38–43. [Google Scholar] [CrossRef]
  34. van Straten, M.; Venema, H.W.; Streekstra, G.J.; Majoie, C.B.; den Heeten, G.J.; Grimbergen, C.A. Removal of bone in CT angiography of the cervical arteries by piecewise matched mask bone elimination. Med. Phys. 2004, 31, 2924–2933. [Google Scholar] [CrossRef] [Green Version]
  35. Fu, F.; Wei, J.; Zhang, M.; Yu, F.; Xiao, Y.; Rong, D.; Shan, Y.; Li, Y.; Zhao, C.; Liao, F.; et al. Rapid vessel segmentation and reconstruction of head and neck angiograms using 3D convolutional neural network. Nat. Commun. 2020, 11, 4829. [Google Scholar] [CrossRef]
  36. Payer, C.; Stern, D.; Bischof, H.; Urschler, M. Integrating spatial configuration into heatmap regression based CNNs for landmark localization. Med. Image Anal. 2019, 54, 207–219. [Google Scholar] [CrossRef]
  37. Alansary, A.; Oktay, O.; Li, Y.; Folgoc, L.L.; Hou, B.; Vaillant, G.; Kamnitsas, K.; Vlontzos, A.; Glocker, B.; Kainz, B.; et al. Evaluating reinforcement learning agents for anatomical landmark detection. Med. Image Anal. 2019, 53, 156–164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Xue, H.; Artico, J.; Fontana, M.; Moon, J.C.; Davies, R.H.; Kellman, P. Landmark Detection in Cardiac MRI by Using a Convolutional Neural Network. Radiol. Artif. Intell. 2021, 3, e200197. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A flowchart of the segmental vessel-labeling algorithm. Three path-finding methods were used to label the segments in the circle of Willis. Notably, Method 1 does not require a graph structure generation process, while Methods 2 and 3 require the graph structure generation process.
Figure 1. A flowchart of the segmental vessel-labeling algorithm. Three path-finding methods were used to label the segments in the circle of Willis. Notably, Method 1 does not require a graph structure generation process, while Methods 2 and 3 require the graph structure generation process.
Tomography 09 00113 g001
Figure 2. An example of colored labeling of the arteries in the circle of Willis. The left figure is the centerlines after skeletonization, and the right figure shows the labeled segments in colors after applying a path-finding algorithm. The numbers in the middle text show the 3D coordinates of the two endpoints in each vessel segment along with the vessel segment’s name.
Figure 2. An example of colored labeling of the arteries in the circle of Willis. The left figure is the centerlines after skeletonization, and the right figure shows the labeled segments in colors after applying a path-finding algorithm. The numbers in the middle text show the 3D coordinates of the two endpoints in each vessel segment along with the vessel segment’s name.
Tomography 09 00113 g002
Figure 3. Incorrect path-finding results using Method 1 when compared with correct path-finding results using Methods 2 and 3. (a) Basilar artery (BA) segment, (b) right posterior cerebral artery P2 (R-P2) segment, and (c) anterior communicating artery (AComm) segment. Method 1 finds detoured paths in (ac). The colorful lines indicate the specific vessel segments (see Figure 2 for the correspondences between colors and vessel segments).
Figure 3. Incorrect path-finding results using Method 1 when compared with correct path-finding results using Methods 2 and 3. (a) Basilar artery (BA) segment, (b) right posterior cerebral artery P2 (R-P2) segment, and (c) anterior communicating artery (AComm) segment. Method 1 finds detoured paths in (ac). The colorful lines indicate the specific vessel segments (see Figure 2 for the correspondences between colors and vessel segments).
Tomography 09 00113 g003
Figure 4. Incorrect path-finding results using Methods 2 and 3 when compared with correct path-finding results using Method 1. (a) The R-ICA and L-ICA segments in Method 2 found the shortest paths, but the paths are incorrect, given the nature of tortuous vessels in the ICA in general. Compare the ICA vessels indicated by the purple arrows. (b) The L-ICA segment in Methods 2 and 3 found the shortest paths, which is incorrect. Compare the ICA vessels indicated by the yellow arrows.
Figure 4. Incorrect path-finding results using Methods 2 and 3 when compared with correct path-finding results using Method 1. (a) The R-ICA and L-ICA segments in Method 2 found the shortest paths, but the paths are incorrect, given the nature of tortuous vessels in the ICA in general. Compare the ICA vessels indicated by the purple arrows. (b) The L-ICA segment in Methods 2 and 3 found the shortest paths, which is incorrect. Compare the ICA vessels indicated by the yellow arrows.
Tomography 09 00113 g004
Table 1. Evaluation of the path-finding methods.
Table 1. Evaluation of the path-finding methods.
Method 1
(DFS Algorithm)
Method 2
(Dijkstra Algorithm)
Method 3
(A* Algorithm)
p-Value *p-Value †p-Value ‡No. of Undetected PathsTotal
No. of Correct PathsNo. of Incorrect PathsNo. of Correct PathsNo. of Incorrect PathsNo. of Correct PathsNo. of Incorrect Paths
AComm 1925440368<0.001<0.0010.0061660
ACA A1R5196006000.0030.0031060
L49105905900.0010.0011160
MCA M1R5556006000.0570.0571060
L5826005910.49611060
PCommR861401400.0160.01614660
L1081801800.0030.00314260
PCA P1R5285735820.2040.0951060
L5196006000.0030.0031060
PCA P2R5556006000.0570.0571060
L5196006000.0030.0031060
BA 5376006000.0130.0131060
ICAR51950105010111060
L519528528111060
Total6141217142170629 105840
* Comparison between Method 1 (DFS algorithm) and Method 2 (Dijkstra algorithm). † Comparison between Method 1 (DFS algorithm) and Method 3 (A* algorithm). ‡ Comparison between Method 2 (Dijkstra algorithm) and Method 3 (A* algorithm).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, S.-O.; Kim, Y.-C. Effects of Path-Finding Algorithms on the Labeling of the Centerlines of Circle of Willis Arteries. Tomography 2023, 9, 1423-1433. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography9040113

AMA Style

Kim S-O, Kim Y-C. Effects of Path-Finding Algorithms on the Labeling of the Centerlines of Circle of Willis Arteries. Tomography. 2023; 9(4):1423-1433. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography9040113

Chicago/Turabian Style

Kim, Se-On, and Yoon-Chul Kim. 2023. "Effects of Path-Finding Algorithms on the Labeling of the Centerlines of Circle of Willis Arteries" Tomography 9, no. 4: 1423-1433. https://0-doi-org.brum.beds.ac.uk/10.3390/tomography9040113

Article Metrics

Back to TopTop