Magnetic resonance angiography: From anatomical knowledge modeling to vessel segmentation

YoloCurvSeg: You only label one noisy skeleton for vessel-style curvilinear structure segmentation

Weakly-supervised learning (WSL) has been proposed to alleviate the conflict between data annotation cost and model performance through employing sparsely-grained (i.e., point-, box-, scribble-wise) supervision and has shown promising performance, particularly in the image segmentation field. However, it is still a very challenging task due to the limited supervision, especially when only a small number of labeled samples are available. Additionally, almost all existing WSL segmentation methods are designed for star-convex structures which are very different from curvilinear structures such as vessels and nerves. In this paper, we propose a novel sparsely annotated segmentation framework for curvilinear structures, named YoloCurvSeg. A very essential component of YoloCurvSeg is image synthesis. Specifically, a background generator delivers image backgrounds that closely match the real distributions through inpainting dilated skeletons. The extracted backgrounds are then combined with randomly emulated curves generated by a Space Colonization Algorithm-based foreground generator and through a multilayer patch-wise contrastive learning synthesizer. In this way, a synthetic dataset with both images and curve segmentation labels is obtained, at the cost of only one or a few noisy skeleton annotations. Finally, a segmenter is trained with the generated dataset and possibly an unlabeled dataset. The proposed YoloCurvSeg is evaluated on four publicly available datasets (OCTA500, CORN, DRIVE and CHASEDB1) and the results show that YoloCurvSeg outperforms state-of-the-art WSL segmentation methods by large margins. With only one noisy skeleton annotation (respectively 0.14%, 0.03%, 1.40%, and 0.65% of the full annotation), YoloCurvSeg achieves more than 97% of the fully-supervised performance on each dataset. Code and datasets will be released at https://github.com/llmir/YoloCurvSeg.

Human treelike tubular structure segmentation: A comprehensive review and future perspectives

Various structures in human physiology follow a treelike morphology, which often expresses complexity at very fine scales. Examples of such structures are intrathoracic airways, retinal blood vessels, and hepatic blood vessels. Large collections of 2D and 3D images have been made available by medical imaging modalities such as magnetic resonance imaging (MRI), computed tomography (CT), Optical coherence tomography (OCT) and ultrasound in which the spatial arrangement can be observed. Segmentation of these structures in medical imaging is of great importance since the analysis of the structure provides insights into disease diagnosis, treatment planning, and prognosis. Manually labelling extensive data by radiologists is often time-consuming and error-prone. As a result, automated or semi-automated computational models have become a popular research field of medical imaging in the past two decades, and many have been developed to date. In this survey, we aim to provide a comprehensive review of currently publicly available datasets, segmentation algorithms, and evaluation metrics. In addition, current challenges and future research directions are discussed.

Automatic segmentation, feature extraction and comparison of healthy and stroke cerebral vasculature

Accurate segmentation of cerebral vasculature and a quantitative assessment of its morphology is critical to various diagnostic and therapeutic purposes and is pertinent to studying brain health and disease. However, this is still a challenging task due to the complexity of the vascular imaging data. We propose an automated method for cerebral vascular segmentation without the need of any manual intervention as well as a method to skeletonize the binary segmented map to extract vascular geometric features and characterize vessel structure. We combine a Hessian-based probabilistic vessel-enhancing filtering with an active-contour-based technique to segment magnetic resonance and computed tomography angiograms (MRA and CTA) and subsequently extract the vessel centerlines and diameters to calculate the geometrical properties of the vasculature. Our method was validated using a 3D phantom of the Circle-of-Willis region, demonstrating 84% mean Dice similarity coefficient (DSC) and 85% mean Pearson's correlation coefficient (PCC) with minimal modified Hausdorff distance (MHD) error (3 surface pixels at most), and showed superior performance compared to existing segmentation algorithms upon quantitative comparison using DSC, PCC and MHD. We subsequently applied our algorithm to a dataset of 40 subjects, including 1) MRA scans of healthy subjects (n = 10, age = 30 ± 9), 2) MRA scans of stroke patients (n = 10, age = 51 ± 15), 3) CTA scans of healthy subjects (n = 10, age = 62 ± 12), and 4) CTA scans of stroke patients (n = 10, age = 68 ± 11), and obtained a quantitative comparison between the stroke and normal vasculature for both imaging modalities. The vascular network in stroke patients compared to age-adjusted healthy subjects was found to have a significantly (p < 0.05) higher tortuosity (3.24 ± 0.88 rad/cm vs. 7.17 ± 1.61 rad/cm for MRA, and 4.36 ± 1.32 rad/cm vs. 7.80 ± 0.92 rad/cm for CTA), higher fractal dimension (1.36 ± 0.28 vs. 1.71 ± 0.14 for MRA, and 1.56 ± 0.05 vs. 1.69 ± 0.20 for CTA), lower total length (3.46 ± 0.99 m vs. 2.20 ± 0.67 m for CTA), lower total volume (61.80 ± 18.79 ml vs. 34.43 ± 22.9 ml for CTA), lower average diameter (2.4 ± 0.21 mm vs. 2.18 ± 0.07 mm for CTA), and lower average branch length (4.81 ± 1.97 mm vs. 8.68 ± 2.03 mm for MRA), respectively. We additionally studied the change in vascular features with respect to aging and imaging modality. While we observed differences between features as a result of aging, statistical analysis did not show any significant differences, whereas we found that the number of branches were significantly different (p < 0.05) between the two imaging modalities (201 ± 73 for MRA vs. 189 ± 69 for CTA). Our segmentation and feature extraction algorithm can be applied on any imaging modality and can be used in the future to automatically obtain the 3D segmented vasculature for diagnosis and treatment planning as well as to study morphological changes due to stroke and other cerebrovascular diseases (CVD) in the clinic.

Evolution of Human Brain Atlases in Terms of Content, Applications, Functionality, and Availability

Human brain atlases have been evolving tremendously, propelled recently by brain big projects, and driven by sophisticated imaging techniques, advanced brain mapping methods, vast data, analytical strategies, and powerful computing. We overview here this evolution in four categories: content, applications, functionality, and availability, in contrast to other works limited mostly to content. Four atlas generations are distinguished: early cortical maps, print stereotactic atlases, early digital atlases, and advanced brain atlas platforms, and 5 avenues in electronic atlases spanning the last two generations. Content-wise, new electronic atlases are categorized into eight groups considering their scope, parcellation, modality, plurality, scale, ethnicity, abnormality, and a mixture of them. Atlas content developments in these groups are heading in 23 various directions. Application-wise, we overview atlases in neuroeducation, research, and clinics, including stereotactic and functional neurosurgery, neuroradiology, neurology, and stroke. Functionality-wise, tools and functionalities are addressed for atlas creation, navigation, individualization, enabling operations, and application-specific. Availability is discussed in media and platforms, ranging from mobile solutions to leading-edge supercomputers, with three accessibility levels. The major application-wise shift has been from research to clinical practice, particularly in stereotactic and functional neurosurgery, although clinical applications are still lagging behind the atlas content progress. Atlas functionality also has been relatively neglected until recently, as the management of brain data explosion requires powerful tools. We suggest that the future human brain atlas-related research and development activities shall be founded on and benefit from a standard framework containing the core virtual brain model cum the brain atlas platform general architecture.

Curvilinear Structure Analysis by Ranking the Orientation Responses of Path Operators

The analysis of thin curvilinear objects in 3D images is a complex and challenging task. In this article, we introduce a new, non-linear operator, called RORPO (Ranking the Orientation Responses of Path Operators). Inspired by the multidirectional paradigm currently used in linear filtering for thin structure analysis, RORPO is built upon the notion of path operator from mathematical morphology. This operator, unlike most operators commonly used for 3D curvilinear structure analysis, is discrete, non-linear and non-local. From this new operator, two main curvilinear structure characteristics can be estimated: an intensity feature, that can be assimilated to a quantitative measure of curvilinearity; and a directional feature, providing a quantitative measure of the structure's orientation. We provide a full description of the structural and algorithmic details for computing these two features from RORPO, and we discuss computational issues. We experimentally assess RORPO by comparison with three of the most popular curvilinear structure analysis filters, namely Frangi Vesselness, Optimally Oriented Flux, and Hybrid Diffusion with Continuous Switch. In particular, we show that our method provides up to 8 percent more true positive and 50 percent less false positives than the next best method, on synthetic and real 3D images.

Robust Segmentation of the Full Cerebral Vasculature in 4D CT of Suspected Stroke Patients

A robust method is presented for the segmentation of the full cerebral vasculature in 4-dimensional (4D) computed tomography (CT). The method consists of candidate vessel selection, feature extraction, random forest classification and postprocessing. Image features include among others the weighted temporal variance image and parameters, including entropy, of an intensity histogram in a local region at different scales. These histogram parameters revealed to be a strong feature in the detection of vessels regardless of shape and size. The method was trained and tested on a large database of 264 patients with suspicion of acute ischemia who underwent 4D CT in our hospital in the period January 2014 to December 2015. Five subvolumes representing different regions of the cerebral vasculature were annotated in each image in the training set by medical assistants. The evaluation was done on 242 patients. A total of 16 (<8%) patients showed severe under or over segmentation and were reported as failures. One out of five subvolumes was randomly annotated in 159 patients and was used for quantitative evaluation. Quantitative evaluation showed a Dice coefficient of 0.91 ± 0.07 and a modified Hausdorff distance of 0.23 ± 0.22 mm. Therefore, robust vessel segmentation in 4D CT is feasible with good accuracy.

A Stereotactic Probabilistic Atlas for the Major Cerebral Arteries

Improved whole brain angiographic and velocity-sensitive MRI is pushing the boundaries of noninvasively obtained cerebral vascular flow information. The complexity of the information contained in such datasets calls for automated algorithms and pipelines, thus reducing the need of manual analyses by trained radiologists. The objective of this work was to lay the foundation for such automated pipelining by constructing and evaluating a probabilistic atlas describing the shape and location of the major cerebral arteries. Specifically, we investigated how the implementation of a non-linear normalization into Montreal Neurological Institute (MNI) space improved the alignment of individual arterial branches. In a population-based cohort of 167 subjects, age 64–68 years, we performed 4D flow MRI with whole brain volumetric coverage, yielding both angiographic and anatomical data. For each subject, sixteen cerebral arteries were manually labeled to construct the atlas. Angiographic data were normalized to MNI space using both rigid-body and non-linear transformations obtained from anatomical images. The alignment of arterial branches was significantly improved by the non-linear normalization (p < 0.001). Validation of the atlas was based on its applicability in automatic arterial labeling. A leave-one-out validation scheme revealed a labeling accuracy of 96 %. Arterial labeling was also performed in a separate clinical sample (n = 10) with an accuracy of 92.5 %. In conclusion, using non-linear spatial normalization we constructed an artery-specific probabilistic atlas, useful for cerebral arterial labeling.

Digital reconstruction and morphometric analysis of human brain arterial vasculature from magnetic resonance angiography

Characterization of the complex branching architecture of cerebral arteries across a representative sample of the human population is important for diagnosing, analyzing, and predicting pathological states. Brain arterial vasculature can be visualized by magnetic resonance angiography (MRA). However, most MRA studies are limited to qualitative assessments, partial morphometric analyses, individual (or small numbers of) subjects, proprietary datasets, or combinations of the above limitations. Neuroinformatics tools, developed for neuronal arbor analysis, were used to quantify vascular morphology from 3T time-of-flight MRA high-resolution (620 μm isotropic) images collected in 61 healthy volunteers (36/25 F/M, average age=31.2 ± 10.7, range=19-64 years). We present in-depth morphometric analyses of the global and local anatomical features of these arbors. The overall structure and size of the vasculature did not significantly differ across genders, ages, or hemispheres. The total length of the three major arterial trees stemming from the circle of Willis (from smallest to largest: the posterior, anterior, and middle cerebral arteries; or PCAs, ACAs, and MCAs, respectively) followed an approximate 1:2:4 proportion. Arterial size co-varied across individuals: subjects with one artery longer than average tended to have all other arteries also longer than average. There was no net right-left difference across the population in any of the individual arteries, but ACAs were more lateralized than MCAs. MCAs, ACAs, and PCAs had similar branch-level properties such as bifurcation angles. Throughout the arterial vasculature, there were considerable differences between branch types: bifurcating branches were significantly shorter and straighter than terminating branches. Furthermore, the length and meandering of bifurcating branches increased with age and with path distance from the circle of Willis. All reconstructions are freely distributed through a public database to enable additional analyses and modeling (cng.gmu.edu/brava).

An attempt toward objective assessment of brain tumor vascularization using susceptibility weighted imaging and dedicated computer program - a preliminary study

Background:

Susceptibility weighted imaging (SWI) is a novel MRI sequence which demonstrates the susceptibility differences between adjacent tissues and it is promising to be a sequence useful in the assessment of brain tumors vascularity. The aim of our study was to demonstrate usefulness of SWI in evaluation of intratumoral vessels in comparison to CET1 sequence in a standardized, objective manner.

Material/Methods:

10 patients with supratentorial brain tumors were included in the study. All of them underwent conventional MRI examination with a 1,5 T scanner. SWI sequence was additionally performed using the following parameters: TR 49 ms,TE 40 ms. We used authors’ personal computer software – Vessels View, to assess the vessels number.

Results:

Comparison of SWI and CET1 sequences was performed using our program. Analysis of all 26 ROIs demonstrated predominance of SWI in the amount of white pixels (vessel cross-sectional) and a similar number of elongated structures (blood vessels).

Conclusions:

To conclude, the results of this study are encouraging; they confirm the added value of SWI as an appropriate and useful sequence in the process of evaluation of intratumoral vascularity. Using our program significantly improved visualization of blood vessels in cerebral tumors. The Vessel View application assists radiologists in demonstrating the vessels and facilitates distinguishing them from adjacent tissues in the image.

Quantification of the human cerebrovasculature: a 7Tesla and 320-row CT in vivo study

OBJECTIVE

Human cerebrovasculature has not been quantified in volume, length, and vascular-brain relationships. We investigated this using imaging.

METHODS

From 0.5-mm 7T and 320-row CT acquisitions, 6 arterial and 4 venous systems were reconstructed, measured, and analyzed.

RESULTS

The ratio of the volume of arterial to venous system is approximately 1:3. The ratio of the volume of dural sinuses to vasculature is 1:2. The ratio of the posterior (PCA) to anterior (ACA) to middle cerebral artery (MCA) is 1:2:4 in volume and length. Ratios of left to right vessels are 1:1 for arteries and veins. Ratios of branching frequency for the ACA, MCA, and PCA are 1:1:1. The branching frequency ratio for superficial to deep veins is 1:2. The MCA occupies 1/2 of arterial length and 1/4 of vascular length. The ratio of the length of superficial to deep veins is 1:1 and each is equal to 1/4 of the vascular length. The ratio of cerebrovasculature to brain volume is 2.5%.

CONCLUSIONS

Despite its enormous complexity, cerebrovasculature is characterized by 4 approximate proportions, 1:1, 1:2, 1:3, 1:4, and their combinations, 1:1:1 and 1:2:4.

Automatic correction of gaps in cerebrovascular segmentations extracted from 3D time-of-flight MRA datasets

OBJECTIVES

Exact cerebrovascular segmentations are required for several applications in today's clinical routine. A major drawback of typical automatic segmentation methods is the occurrence of gaps within the segmentation. These gaps are typically located at small vessel structures exhibiting low intensities. Manual correction is very time-consuming and not suitable in clinical practice. This work presents a post-processing method for the automatic detection and closing of gaps in cerebrovascular segmentations.

METHODS

In this approach, the 3D centerline is calculated from an available vessel segmentation, which enables the detection of corresponding vessel endpoints. These endpoints are then used to detect possible connections to other 3D centerline voxels with a graph-based approach. After consistency check, reasonable detected paths are expanded to the vessel boundaries using a level set approach and combined with the initial segmentation.

RESULTS

For evaluation purposes, 100 gaps were artificially inserted at non-branching vessels and bifurcations in manual cerebrovascular segmentations derived from ten Time-of-Flight magnetic resonance angiography datasets. The results show that the presented method is capable of detecting 82% of the non-branching vessel gaps and 84% of the bifurcation gaps. The level set segmentation expands the detected connections with 0.42 mm accuracy compared to the initial segmentations. A further evaluation based on 10 real automatic segmentations from the same datasets shows that the proposed method detects 35 additional connections in average per dataset, whereas 92.7% were rated as correct by a medical expert.

CONCLUSION

The presented approach can considerably improve the accuracy of cerebrovascular segmentations and of following analysis outcomes.

Comparison of magnetic resonance angiography scans on 1.5, 3, and 7 Tesla units: a quantitative study of 3-dimensional cerebrovasculature

BACKGROUND

Although multiple studies demonstrate benefits of high field imaging of cerebrovasculature, a detailed quantitative analysis of complete cerebrovascular system is unavailable. To compare quality of MR angiography (MRA) acquisitions at various field strengths, we used 3-dimensional (3D) geometric cerebrovascular models extracted from 1.5 T/3 T/7 T scans.

METHODS

The 3D cerebrovascular models were compared in volume, length, and number of branches. A relationship between the vascular length and volume was statistically derived. Acquisition performance was benchmarked against the maximum volume at infinitive length.

RESULTS

The numbers of vessels discernible on 1.5 T/3 T/7 T are 138/363/907. 3T shows 3.3(1.9) and 7 T 1.2(9.1) times more arteries (veins) than 1.5 T. The vascular lengths and volumes at 1.5 T/3 T/7 T are 3.7/12.5/22.7 m and 15.8/26.6/28.0 cm(3). For arteries: 3T-1.5 T gain is very high in length, high in volume; 7 T-3T gain is medium in length, small in volume. For veins: 3 T-1.5 T gain is moderate in length, high in volume; 7 T-3T gain is very high in length, moderate in volume. 1.5 T shows merely half of vascular volume. At 3 T 6%, while at 7 T only 1% of vascular volume is missing.

CONCLUSION

Our approach differs from standard approaches based on visual assessment and signal (contrast)-to-noise ratio. It also measures absolute acquisition performance, provides a unique length-volume relationship, and predicts length/volume for intermediate teslages.

Simulation and assessment of cerebrovascular damage in deep brain stimulation using a stereotactic atlas of vasculature and structure derived from multiple 3- and 7-tesla scans

Object

The most severe complication of deep brain stimulation (DBS) is intracranial hemorrhage. Detailed knowledge of the cerebrovasculature could reduce the rate of this disorder. Morphological scans typically acquired in stereotactic and functional neurosurgery (SFN) by using 1.5-T (or sometimes even 3-T) imaging units poorly depict the vasculature. Advanced angiographic imaging, including 3- and 7-T 3D time-of-flight and susceptibility weighted imaging as well as 320-slice CT angiography, depict the vessels in great detail. However, these acquisitions are not used in SFN clinical practice, and robust methods for their processing are not available yet. Therefore, the authors proposed the use of a detailed 3D stereotactic cerebrovascular atlas to assist in SFN planning and to potentially reduce DBS-induced hemorrhage.

Methods

A very detailed 3D cerebrovascular atlas of arteries, veins, and dural sinuses was constructed from multiple 3- and 7-T scans. The atlas contained > 900 vessels, each labeled with a name and diameter with the smallest having a 90-μm diameter. The cortical areas, ventricular system, and subcortical structures were fully segmented and labeled, including the main stereotactic target structures: subthalamic nucleus, ventral intermediate nucleus of the thalamus, and internal globus pallidus. The authors also developed a computer simulator with the embedded atlas that was able to compute the effective electrode trajectory by minimizing penetration of the cerebrovascular system and vital brain structures by a DBS electrode. The simulator provides the neurosurgeon with functions for atlas manipulation, target selection, trajectory planning and editing, 3D display and manipulation, and electrode-brain penetration calculation.

Results

This simulation demonstrated that a DBS electrode inserted in the middle frontal gyrus may intersect several arteries and veins including 1) the anteromedial frontal artery of the anterior cerebral artery as well as the prefrontal artery and the precentral sulcus artery of the middle cerebral artery (range of diameters 0.4–0.6 mm); and 2) the prefrontal, anterior caudate, and medullary veins (range of diameters 0.1–2.3 mm). This work also shows that field strength and pulse sequence have a substantial impact on vessel depiction. The numbers of 3D vascular segments are 215, 363, and 907 for 1.5-, 3-, and 7-T scans, respectively.

Conclusions

Inserting devices into the brain during microrecording and stimulation may cause microbleeds not discernible on standard scans. A small change in the location of the DBS electrode can result in a major change for the patient. The described simulation increases the neurosurgeon's awareness of this phenomenon. The simulator enables the neurosurgeon to analyze the spatial relationships between the track and the cerebrovasculature, ventricles, subcortical structures, and cortical areas, which allows the DBS electrode to be placed more effectively, and thus potentially reducing the invasiveness of the stimulation procedure for the patient.

Vascular Tree Reconstruction by Minimizing A Physiological Functional Cost

The reconstruction of complete vascular trees from medical images has many important applications. Although vessel detection has been extensively investigated, little work has been done on how connect the results to reconstruct the full trees. In this paper, we propose a novel theoretical framework for automatic vessel connection, where the automation is achieved by leveraging constraints from the physiological properties of the vascular trees. In particular, a physiological functional cost for the whole vascular tree is derived and an efficient algorithm is developed to minimize it. The method is generic and can be applied to different vessel detection/segmentation results, e.g. the classic rigid detection method as adopted in this paper. We demonstrate the effectiveness of this method on both 2D and 3D data.

A shape-based segmentation algorithm for X-ray digital subtraction angiography images

We present an efficient algorithm for vessel segmentation of digital subtraction angiography (DSA) images. Continuous DSA images (projections), obtained by X-ray fluoroscopy with contrast-media, are normally used as road maps in vessel catheterization. A more efficient technique would consist in the use of a 3D model reconstruction of the vascular tree, instead of continuous X-ray scans, as a map. By separating vessel information from the undesired background (noise and signals coming from other organs and motion artefacts), efficient segmentation can play a key role in reducing the number of projections (X-ray scans) necessary to reconstruct a 3D vascular model. In what follows, the proposed method is described and some experimental results are reported, thus illustrating the behaviour of the algorithm when compared to other segmentation methods, ideated for the same application. The automatic calculation methods for some of the parameters used are also reported and discussed.

Isolated kidney phantom for development of biothermal vascular models with application to high intensity focused ultrasound therapy

A methodology using magnetic resonance angiography (MRA) is presented for identifying thermally significant blood vessels in isolated kidneys, specifically for use in biothermal model development with application to high intensity focused ultrasound (HIFU). A combination of a proven preservation technique, newly developed MR-compatible experimental procedures and the refinement of MR pulse sequence parameters was used to determine vascular characteristics using high-resolution three-dimensional time-of-flight MRA image of flow through isolated kidneys. Results presented are twofold. First, improved vessel visibility was attained through decreasing the magnetic resonance imaging bandwidth from 150 to 30 Hz/pixel while simultaneously increasing the echo time, repetition time, and flip angle; vascular center line extraction showed an 18% improvement in the number of vessel segments detected and a 23% increase in length of the terminal segments over a base line technique without improvements. Second, the overall system was shown to be practical to determine vascular flow effects during HIFU heating; testing results from heating the kidney with HIFU are presented, showing a decrease of average kidney temperature with an increase of flow rate through the kidney with localized cooling demonstrated surrounding known vessel locations.

Vessels as 4-D curves: global minimal 4-D paths to extract 3-D tubular surfaces and centerlines

In this paper, we propose an innovative approach to the segmentation of tubular structures. This approach combines all of the benefits of minimal path techniques such as global minimizers, fast computation, and powerful incorporation of user input, while also having the capability to represent and detect vessel surfaces directly which so far has been a feature restricted to active contour and surface techniques. The key is to represent the trajectory of a tubular structure not as a 3-D curve but to go up a dimension and represent the entire structure as a 4-D curve. Then we are able to fully exploit minimal path techniques to obtain global minimizing trajectories between two user supplied endpoints in order to reconstruct tubular structures from noisy or low contrast 3-D data without the sensitivity to local minima inherent in most active surface techniques. In contrast to standard purely spatial 3-D minimal path techniques, however, we are able to represent a full tubular surface rather than just a curve which runs through its interior. Our representation also yields a natural notion of a tube's "central curve." We demonstrate and validate the utility of this approach on magnetic resonance (MR) angiography and computed tomography (CT) images of coronary arteries.

From MIP image to MRA segmentation using fuzzy set theory

The aim of this paper is to describe a semi-automatic method of segmentation in magnetic resonance angiography (MRA). This method, based on fuzzy set theory, uses the information (gray levels) contained in the maximum intensity projection (MIP) image to segment the 3D vascular structure from slices. Tests have been carried out on vascular phantom and on clinical MRA images. This 3D segmentation method has proved to be satisfactory for the detection of vascular structures even for very complex shapes. Finally, this MIP-based approach is semi-automatic and produces a robust segmentation thanks to the contrast-to-noise ratio and to the slice profile which are taken into account to determine the membership of a voxel to the vascular structure.