A probabilistic atlas based on 168 subjects for labeling of brain arteries
Tora Dunås1, Anders Wåhlin1,2, Khalid Ambarki1,3, Laleh Zarrinkoob4, Jan Malm4, and Anders Eklund1,2,3

1Department of Radiation Sciences, Umeå University, Umeå, Sweden, 2Umeå Center for Functional Brain Imaging, Umeå University, Umeå, Sweden, 3Centre for Biomedical Engineering and Physics, Umeå University, Umeå, Sweden, 4Department of Clinical Neuroscience, Umeå University, Umeå, Sweden

Synopsis

The cerebral arterial system is complex with large inter-individual spatial variations, which are potentially challenging for automatic methods. The objective of this work was to construct an artery specific probabilistic atlas of 16 large cerebral arteries, based on 168 subjects, and to investigate if the regional specificity of vascular branches was sufficient to permit atlas based arterial labeling. Voxels of the arterial centerlines was labeled according to the highest probability at the corresponding location in the atlas. The rate of correctly labeled voxels was over 80% for all arteries, which should be sufficient to permit atlas based arterial labeling.

Purpose

Automatic assessment and blood flow distribution of the major cerebral arteries have applications in various brain diseases1, for example in acute stroke care. The cerebral arterial system is complex with large inter-individual spatial variations, which are potentially challenging for automatic methods. Recently a method for automatic identification of cerebral arteries, that utilized a small sample probabilistic atlas, was presented2. However, to reach its full potential this method will require a population based probabilistic cerebrovascular atlas. The objective of this work was to construct an artery specific probabilistic atlas of 16 large cerebral arteries, based on 168 subjects, and to investigate if the regional specificity of vascular branches was sufficient to permit atlas based arterial labeling.

Methods

The atlas was based on 4D flow MRI3 from 168 subjects (65.6 ± 1.4 year old, 94 men, 74 female) collected on a 3 Tesla General Electric HD750 system and a 32-channel headcoil. FOV: 220 x 220 x 220 mm, voxel size: 0.7x0.7x0.7 mm3, scan time. 9.5 minutes, velocity encoding, 110 cm/s; TR/TE, 6.5/2.7 ms; flip angle, 8°. A temporal maximum intensity projection (tMIP) was reconstructed and used for arterial assessment. The tMIP for each subject was transformed to MNI-space using SPM8’s DARTEL4. The tMIP was smoothed with a low-pass-filter and binarized by thresholding at 18% of the maximum intensity value. The following arteries were segmented and labeled using an inhouse matlab tool: a. internal carotid arteries (ICA); b. vertebral arteries (VA); c. posterior cerebral arteries (PCA); d. proximal middle cerebral arteries (M1 segment); e. distal middle cerebral arteries (at M2-M3); f. proximal anterior cerebral arteries (A1 segment); g. distal anterior cerebral artery (at A2-A3); h. posterior communicating arteries (PCoA); i. basilar artery (BA). A probability map was constructed for each artery by averaging the binary volumes across subjects, these probability maps together formed the atlas. To investigate if the regional specificity of vascular branches was sufficient to permit atlas based arterial labeling, an internal analysis was performed. During arterial segmentation, a single-voxel-thick centerline (CL) was calculated for each artery by gradually thinning the binary image5. The average probability value along each CL was calculated, as well as the proportion of voxels in each CL that was assigned to the correct artery based on the highest probability value at its location. The rationale was that if >50% of the CL voxels were assigned to the correct artery, the arterial segment would be correctly labeled.

Results

The prevalence of arteries in the population based sample that was used for the probabilistic atlas is presented in Table 1. The highest mean probability value along the CL was obtained for ICA. High values were also seen for BA, M1, and A1, although there was a large variation among subjects (Table 1). The rate of correctly labeled CL-voxels was over 80% for all arteries, and over 90% for all arteries except PCoA.

Discussion

The main result of this study was the high accuracy of CL-voxel assignment supporting that the cerebrovascular arterial system permits atlas based arterial labeling. The average probability along the CL is related to the spatial variation across subjects. The closer to each other the arteries are located, the higher the average probability value. However, a lower mean probability along the centerlines did not equal a low rate of correct labeling. For example, the lowest probability was obtained for M2-M3, but the corresponding rate of correct labeling was among the highest because it has a more isolated spatial territory.

Conclusion

We successfully constructed an artery specific probabilistic atlas that included 16 major cerebral arteries. The atlas showed a regional specificity of vascular branches that was sufficient to permit atlas based arterial labeling. This is a first step toward medical applications such as automatic arterial vessel identification, improved vascular segmentation, automatic blood flow assessment, and characterization of cerebral arterial morphology.

Acknowledgements

This study was supported by the Swedish Research Council Grant 621-2011-5216, the Swedish Heart and Lung Foundation Grants 20110383 and 20140592 and the Swedish Brain Foundation.

References

1. Wright SN, Kochunov P, Mut F, et al. Digital reconstruction and morphometric analysis of human brain arterial vasculature from magnetic resonance angiography. Neuroimage 2013;82:170–81.

2. Dunås T, Wåhlin A, Ambarki K, et al. Automatic labeling of cerebral arteries in magnetic resonance angiography. Magn Reson Mater Phy. In press.

3. Gu T, Korosec FR, Block WF, et al. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. Am J Neuroradiol 2005;26:743–9.

4. Ashburner J. A fast diffeomorphic image registration algorithm. Neuroimage 2007;38:95–113.

5. Chen Z, Molloi S. Automatic 3D vascular tree construction in CT angiography. Comput Med Imaging Graph 2003;27:469–479.

Figures

Table 1: Prevalence of subjects where each artery could be segmented, mean probability value in percentage along the centerline of the arteries, and rate of correct labeling along the centerlines in percentage.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
4122