Fusing speed and phase information for vascular segmentation of phase contrast MR angiograms

Vascular segmentation of head phase-contrast magnetic resonance angiograms using grayscale and shape features

BACKGROUND AND OBJECTIVE

In neurosurgery planning, vascular structures must be predetermined, which can guarantee the security of the operation carried out in the case of avoiding blood vessels. In this paper, an automatic algorithm of vascular segmentation, which combined the grayscale and shape features of the blood vessels, is proposed to extract 3D vascular structures from head phase-contrast magnetic resonance angiography dataset.

METHODS

First, a cost function of mis-segmentation is introduced on the basis of traditional Bayesian statistical classification, and the blood vessel of weak grayscale that tended to be misclassified into background will be preserved. Second, enhanced vesselness image is obtained according to the shape-based multiscale vascular enhancement filter. Third, a new reconstructed vascular image is established according to the fusion of vascular grayscale and shape features using Dempster-Shafer evidence theory; subsequently, the corresponding segmentation structures are obtained. Finally, according to the noise distribution characteristic of the data, segmentation ratio coefficient, which increased linearly from top to bottom, is proposed to control the segmentation result, thereby preventing over-segmentation.

RESULTS

Experiment results show that, through the proposed method, vascular structures can be detected not only when both grayscale and shape features are strong, but also when either of them is strong. Compared with traditional grayscale feature- and shape feature-based methods, it is better in the evaluation of testing in segmentation accuracy, and over-segmentation and under-segmentation ratios.

CONCLUSIONS

The proposed grayscale and shape features combined vascular segmentation is not only effective but also accurate. It may be used for diagnosis of vascular diseases and planning of neurosurgery.

Segmentation of intracranial vessels and aneurysms in phase contrast magnetic resonance angiography using multirange filters and local variances

Segmentation of intensity varying and low-contrast structures is an extremely challenging and rewarding task. In computer-aided diagnosis of intracranial aneurysms, segmenting the high-intensity major vessels along with the attached low-contrast aneurysms is essential to the recognition of this lethal vascular disease. It is particularly helpful in performing early and noninvasive diagnosis of intracranial aneurysms using phase contrast magnetic resonance angiographic (PC-MRA) images. The major challenges of developing a PC-MRA-based segmentation method are the significantly varying voxel intensity inside vessels with different flow velocities and the signal loss in the aneurysmal regions where turbulent flows occur. This paper proposes a novel intensity-based algorithm to segment intracranial vessels and the attached aneurysms. The proposed method can handle intensity varying vasculatures and also the low-contrast aneurysmal regions affected by turbulent flows. It is grounded on the use of multirange filters and local variances to extract intensity-based image features for identifying contrast varying vasculatures. The extremely low-intensity region affected by turbulent flows is detected according to the topology of the structure detected by multirange filters and local variances. The proposed method is evaluated using a phantom image volume with an aneurysm and four clinical cases. It achieves 0.80 dice score in the phantom case. In addition, different components of the proposed method-the multirange filters, local variances, and topology-based detection-are evaluated in the comparison between the proposed method and its lower complexity variants. Owing to the analogy between these variants and existing vascular segmentation methods, this comparison also exemplifies the advantage of the proposed method over the existing approaches. It analyzes the weaknesses of these existing approaches and justifies the use of every component involved in the proposed method. It is shown that the proposed method is capable of segmenting blood vessels and the attached aneurysms on PC-MRA images.

A fast and fully automatic method for cerebrovascular segmentation on time-of-flight (TOF) MRA image

The precise three-dimensional (3-D) segmentation of cerebral vessels from magnetic resonance angiography (MRA) images is essential for the detection of cerebrovascular diseases (e.g., occlusion, aneurysm). The complex 3-D structure of cerebral vessels and the low contrast of thin vessels in MRA images make precise segmentation difficult. We present a fast, fully automatic segmentation algorithm based on statistical model analysis and improved curve evolution for extracting the 3-D cerebral vessels from a time-of-flight (TOF) MRA dataset. Cerebral vessels and other tissue (brain tissue, CSF, and bone) in TOF MRA dataset are modeled by Gaussian distribution and combination of Rayleigh with several Gaussian distributions separately. The region distribution combined with gradient information is used in edge-strength of curve evolution as one novel mode. This edge-strength function is able to determine the boundary of thin vessels with low contrast around brain tissue accurately and robustly. Moreover, a fast level set method is developed to implement the curve evolution to assure high efficiency of the cerebrovascular segmentation. Quantitative comparisons with 10 sets of manual segmentation results showed that the average volume sensitivity, the average branch sensitivity, and average mean absolute distance error are 93.6%, 95.98%, and 0.333 mm, respectively. By applying the algorithm to 200 clinical datasets from three hospitals, it is demonstrated that the proposed algorithm can provide good quality segmentation capable of extracting a vessel with a one-voxel diameter in less than 2 min. Its accuracy and speed make this novel algorithm more suitable for a clinical computer-aided diagnosis system.

Blood Vessel Segmentation in Complex-Valued Magnetic Resonance Images with Snake Active Contour Model

Accurate blood vessel segmentation plays a crucial role in non-invasive blood flow velocity measurement based on complex-valued magnetic resonance images. We propose a specific snake active contour model-based blood vessel segmentation framework for complex-valued magnetic resonance images. The proposed framework combines both magnitude and phase information from a complex-valued image representation to obtain an optimum segmentation result. Magnitude information of the complex-valued image provides a structural localization of the target object, while phase information identifies the existence of flowing matters within the object. Snake active contour model, which models the segmentation procedure as a force-balancing physical system, is being adopted as a framework for this work due to its interactive, dynamic, and customizable characteristics. Two snake-based segmentation models are developed to produce a more accurate segmentation result, namely the Model-constrained Gradient Vector Flow-snake (MC GVF-snake) and Stochastic-snake. MC GVF-snake elaborates a prior knowledge on common physical structure of the target object to restrict and guide the segmentation mechanism, while Stochastic-snake implements the simulated annealing stochastic procedure to produce improved segmentation accuracy. The developed segmentation framework has been evaluated on actual complex-valued MRI images, both in noise-free and noisy simulated conditions. Evaluation results indicate that both of the developed algorithms give an improved segmentation performance as well as increased robustness, in comparison to the conventional snake algorithm.

Optimum fuzzy filters for phase-contrast magnetic resonance imaging segmentation

PURPOSE

To develop and validate a multidimensional segmentation and filtering methodology for accurate blood flow velocity field reconstruction from phase-contrast magnetic resonance imaging (PC MRI).

MATERIALS AND METHODS

The proposed technique consists of two steps: (1) the boundary of the vessel is automatically segmented using the active contour approach; and (2) the noise embedded within the segmented vector field is selectively removed using a novel fuzzy adaptive vector median filtering (FAVMF) technique. This two-step segmentation process was tested and validated on 111 synthetically generated PC MRI slices and on 10 patients with congenital heart disease.

RESULTS

The active contour technique was effective for segmenting blood vessels having a sensitivity and specificity of 93.1% and 92.1% using manual segmentation as a reference standard. FAVMF was the superior technique in filtering out noise vectors, when compared with other commonly used filters in PC MRI (P < 0.05). The peak wall shear rate calculated from the PC MRI data (248 +/- 39 sec(-1)), was significantly decreased to (146 +/- 26 sec(-1)) after the filtering process.

CONCLUSION

The proposed two-step segmentation and filtering methodology is more accurate compared to a single-step segmentation process for post-processing of PC MRI data.

A new adaptive probabilistic model of blood vessels for segmenting MRA images

A new physically justified adaptive probabilistic model of blood vessels on magnetic resonance angiography (MRA) images is proposed. The model accounts for both laminar (for normal subjects) and turbulent blood flow (in abnormal cases like anemia or stenosis) and results in a fast algorithm for extracting a 3D cerebrovascular system from the MRA data. Experiments with synthetic and 50 real data sets confirm the high accuracy of the proposed approach.

State-space analysis of cardiac motion with biomechanical constraints

Quantitative estimation of nonrigid motion from image sequences has important technical and practical significance. State-space analysis provides powerful and convenient ways to construct and incorporate the physically meaningful system dynamics of an object, the image-derived observations, and the process and measurement noise disturbances. In this paper, we present a biomechanical-model constrained state-space analysis framework for the multiframe estimation of the periodic cardiac motion and deformation. The physical constraints take the roles as spatial regulator of the myocardial behavior and spatial filter/interpolator of the data measurements, while techniques from statistical filtering theory impose spatiotemporal constraints to facilitate the incorporation of multiframe information to generate optimal estimates of the heart kinematics. Physiologically meaningful results have been achieved from estimated displacement fields and strain maps using in vivo left ventricular magnetic resonance tagging and phase contrast image sequences, which provide the tag-tag and tag-boundary displacement inputs, and the mid-wall instantaneous velocity information and boundary displacement measures, respectively.

Coronary vessel trees from 3D imagery: a topological approach

We propose a simple method for reconstructing vascular trees from 3D images. Our algorithm extracts persistent maxima of the intensity on all axis-aligned 2D slices of the input image. The maxima concentrate along 1D intensity ridges, in particular along blood vessels. We build a forest connecting the persistent maxima with short edges. The forest tends to approximate the blood vessels present in the image, but also contains numerous spurious features and often fails to connect segments belonging to one vessel in low contrast areas. We improve the forest by applying simple geometric filters that trim short branches, fill gaps in blood vessels and remove spurious branches from the vascular tree to be extracted. Experiments show that our technique can be applied to extract coronary trees from heart CT scans.

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.

Statistical cerebrovascular segmentation in three-dimensional rotational angiography based on maximum intensity projections

Segmentation of three-dimensional rotational angiography (3D-RA) can provide quantitative 3D morphological information of vasculature. The expectation maximization-(EM-) based segmentation techniques have been widely used in the medical image processing community, because of the implementation simplicity, and computational efficiency of the approach. In a brain 3D-RA, vascular regions usually occupy a very small proportion (around 1%) inside an entire image volume. This severe imbalance between the intensity distributions of vessels and background can lead to inaccurate statistical modeling in the EM-based segmentation methods, and thus adversely affect the segmentation quality for 3D-RA. In this paper we present a new method for the extraction of vasculature in 3D-RA images. The new method is fully automatic and computationally efficient. As compared with the original 3D-RA volume, there is a larger proportion (around 20%) of vessels in its corresponding maximum intensity projection (MIP) image. The proposed method exploits this property to increase the accuracy of statistical modeling with the EM algorithm. The algorithm takes an iterative approach to compiling the 3D vascular segmentation progressively with the segmentation of MIP images along the three principal axes, and use a winner-takes-all strategy to combine the results obtained along individual axes. Experimental results on 12 3D-RA clinical datasets indicate that the segmentations obtained by the new method exhibit a high degree of agreement to the ground truth segmentations and are comparable to those produced by the manual optimal global thresholding method.

Bayesian image segmentation using local iso-intensity structural orientation

Image segmentation is a fundamental problem in early computer vision. In segmentation of flat shaded, nontextured objects in real-world images, objects are usually assumed to be piecewise homogeneous. This assumption, however, is not always valid with images such as medical images. As a result, any techniques based on this assumption may produce less-than-satisfactory image segmentation. In this work, we relax the piecewise homogeneous assumption. By assuming that the intensity nonuniformity is smooth in the imaged objects, a novel algorithm that exploits the coherence in the intensity profile to segment objects is proposed. The algorithm uses a novel smoothness prior to improve the quality of image segmentation. The formulation of the prior is based on the coherence of the local structural orientation in the image. The segmentation process is performed in a Bayesian framework. Local structural orientation estimation is obtained with an orientation tensor. Comparisons between the conventional Hessian matrix and the orientation tensor have been conducted. The experimental results on the synthetic images and the real-world images have indicated that our novel segmentation algorithm produces better segmentations than both the global thresholding with the maximum likelihood estimation and the algorithm with the multilevel logistic MRF model.

A comparison of a similarity-based and a feature-based 2-D-3-D registration method for neurointerventional use

Two-dimensional (2-D)-to-three-dimensional (3-D) registration can improve visualization which may aid minimally invasive neurointerventions. Using clinical and phantom studies, two state-of-the-art approaches to rigid registration are compared quantitatively: an intensity-based algorithm using the gradient difference similarity measure; and an iterative closest point (ICP)-based algorithm. The gradient difference approach was found to be more accurate, with an average registration accuracy of 1.7 mm for clinical data, compared to the ICP-based algorithm with an average accuracy of 2.8 mm. In phantom studies, the ICP-based algorithm proved more reliable, but with more complicated clinical data, the gradient difference algorithm was more robust. Average computation time for the ICP-based algorithm was 20 s per registration, compared with 14 min and 50 s for the gradient difference algorithm.

Vascular segmentation of phase contrast magnetic resonance angiograms based on statistical mixture modeling and local phase coherence

In this paper, we present an approach to segmenting the brain vasculature in phase contrast magnetic resonance angiography (PC-MRA). According to our prior work, we can describe the overall probability density function of a PC-MRA speed image as either a Maxwell-uniform (MU) or Maxwell-Gaussian-uniform (MGU) mixture model. An automatic mechanism based on Kullback-Leibler divergence is proposed for selecting between the MGU and MU models given a speed image volume. A coherence measure, namely local phase coherence (LPC), which incorporates information about the spatial relationships between neighboring flow vectors, is defined and shown to be more robust to noise than previously described coherence measures. A statistical measure from the speed images and the LPC measure from the phase images are combined in a probabilistic framework, based on the maximum a posteriori method and Markov random fields, to estimate the posterior probabilities of vessel and background for classification. It is shown that segmentation based on both measures gives a more accurate segmentation than using either speed or flow coherence information alone. The proposed method is tested on synthetic, flow phantom and clinical datasets. The results show that the method can segment normal vessels and vascular regions with relatively low flow rate and low signal-to-noise ratio, e.g., aneurysms and veins.

Stochastic finite element framework for simultaneous estimation of cardiac kinematic functions and material parameters

A stochastic finite element framework is presented for the simultaneous estimation of the cardiac kinematic functions and material model parameters from periodic medical image sequences. While existing biomechanics studies of the myocardial material constitutive laws have assumed known tissue kinematic measurements, and image analysis efforts on cardiac kinematic functions have relied on fixed constraining models of mathematical or mechanical nature, we illustrate through synthetic data that a probabilistic joint estimation strategy is needed to achieve more robust and accurate analysis of the kinematic functions and material parameters at the same time. For a particular a priori constraining material model with uncertain subject-dependent parameters and a posteriori noisy imaging based observations, our strategy combines the stochastic differential equations of the myocardial dynamics with the finite element method, and the material parameters and the imaging data are treated as random variables with known prior statistics. After the conversion to state space representation, the extended Kalman filtering procedures are adopted to linearize the equations and to provide the joint estimates in an approximate optimal sense. The estimation bias and convergence issues are addressed, and we conclude experimentally that it is possible to adopt this biomechanical model based multiframe estimation approach to achieve converged estimates because of the periodic nature of the cardiac dynamics. The effort is validated using synthetic data sequence with known kinematics and material parameters. Further, under linear elastic material model, estimation results using canine magnetic resonance phase contrast image sequences are presented, which are in very good agreement with histological tissue staining results, the current gold standards.

Intensity-based 2-D-3-D registration of cerebral angiograms

We propose a new method for aligning three-dimensional (3-D) magnetic resonance angiography (MRA) with 2-D X-ray digital subtraction angiograms (DSA). Our method is developed from our algorithm to register computed tomography volumes to X-ray images based on intensity matching of digitally reconstructed radiographs (DRRs). To make the DSA and DRR more similar, we transform the MRA images to images of the vasculature and set to zero the contralateral side of the MRA to that imaged with DSA. We initialize the search for a match on a user defined circular region of interest. We have tested six similarity measures using both unsegmented MRA and three segmentation variants of the MRA. Registrations were carried out on images of a physical neuro-vascular phantom and images obtained during four neuro-vascular interventions. The most accurate and robust registrations were obtained using the pattern intensity, gradient difference, and gradient correlation similarity measures, when used in conjunction with the most sophisticated MRA segmentations. Using these measures, 95% of the phantom start positions and 82% of the clinical start positions were successfully registered. The lowest root mean square reprojection errors were 1.3 mm (standard deviation 0.6) for the phantom and 1.5 mm (standard deviation 0.9) for the clinical data sets. Finally, we present a novel method for the comparison of similarity measure performance using a technique borrowed from receiver operator characteristic analysis.