Statistical-Based Approach for Extracting 3D Blood Vessels from TOF-MyRA Data

Cerebrovascular segmentation from TOF using stochastic models

In this paper, we present an automatic statistical approach for extracting 3D blood vessels from time-of-flight (TOF) magnetic resonance angiography (MRA) data. The voxels of the dataset are classified as either blood vessels or background noise. The observed volume data is modeled by two stochastic processes. The low level process characterizes the intensity distribution of the data, while the high level process characterizes their statistical dependence among neighboring voxels. The low level process of the background signal is modeled by a finite mixture of one Rayleigh and two normal distributions, while the blood vessels are modeled by one normal distribution. The parameters of the low level process are estimated using the expectation maximization (EM) algorithm. Since the convergence of the EM is sensitive to the initial estimate of the model parameters, an automatic method for parameter initialization, based on histogram analysis, is provided. To improve the quality of segmentation achieved by the proposed low level model especially in the regions of significantly vascular signal loss, the high level process is modeled as a Markov random field (MRF). Since MRF is sensitive to edges and the intracranial vessels represent roughly 5% of the intracranial volume, 2D MRF will destroy most of the small and medium sized vessels. Therefore, to reduce this limitation, we employed 3D MRF, whose parameters are estimated using the maximum pseudo likelihood estimator (MPLE), which converges to the true likelihood under large lattice. Our proposed model exhibits a good fit to the clinical data and is extensively tested on different synthetic vessel phantoms and several 2D/3D TOF datasets acquired from two different MRI scanners. Experimental results showed that the proposed model provides good quality of segmentation and is capable of delineating vessels down to 3 voxel diameters.

Magnetic resonance angiography: From anatomical knowledge modeling to vessel segmentation

Magnetic resonance angiography (MRA) has become a common way to study cerebral vascular structures. Indeed, it enables to obtain information on flowing blood in a totally non-invasive and non-irradiant fashion. MRA exams are generally performed for three main applications: detection of vascular pathologies, neurosurgery planning, and vascular landmark detection for brain functional analysis. This large field of applications justifies the necessity to provide efficient vessel segmentation tools. Several methods have been proposed during the last fifteen years. However, the obtained results are still not fully satisfying. A solution to improve brain vessel segmentation from MRA data could consist in integrating high-level a priori knowledge in the segmentation process. A preliminary attempt to integrate such knowledge is proposed here. It is composed of two methods devoted to phase contrast MRA (PC MRA) data. The first method is a cerebral vascular atlas creation process, composed of three steps: knowledge extraction, registration, and data fusion. Knowledge extraction is performed using a vessel size determination algorithm based on skeletonization, while a topology preserving non-rigid registration method is used to fuse the information into the atlas. The second method is a segmentation process involving adaptive sets of gray-level hit-or-miss operators. It uses anatomical knowledge modeled by the cerebral vascular atlas to adapt the parameters of these operators (number, size, and orientation) to the searched vascular structures. These two methods have been tested by creating an atlas from a 18 MRA database, and by using it to segment 30 MRA images, comparing the results to those obtained from a region-growing segmentation method.