In this paper, we present a novel approach for the segmentation of white matter tracts based on Finsler active contours. This technique provides an optimal measure of connectivity, explicitly segments the connecting fiber bundle, and is equipped with a metric which is able to utilize the directional information of high angular resolution data. We demonstrate the effectiveness of the algorithm for segmenting the cingulum bundle.
We present software engineering methods to provide free open-source software for MR-guided therapy. We report that graphical representation of the surgical tools, interconnectively with the tracking device, patient-to-image registration, and MRI-based thermal mapping are crucial components of MR-guided therapy in sharing such software. Software process includes a network-based distribution mechanism by multi-platform compiling tool CMake, CVS, quality assurance software DART. We developed six procedures in four separate clinical sites using proposed software engineering and process, and found the proposed method is feasible to facilitate multicenter clinical trial of MR-guided therapies. Our future studies include use of the software in non-MR-guided therapies.
BACKGROUND AND PURPOSE: Formation of lesions in multiple sclerosis (MS) shows pronounced short-term fluctuation of MR imaging hyperintensity and size, a qualitatively known but poorly characterized phenomenon. With the use of time-series modeling of MR imaging intensity, our study relates the short-term dynamics of new T2 lesion formation to those of contrast enhancement and markers of long-term progression of disease. MATERIALS AND METHODS: We analyzed 915 examinations from weekly to monthly MR imaging in 40 patients with MS using a time-series model, emulating 2 opposing processes of T2 prolongation and shortening, respectively. Patterns of activity, duration, and residual hyperintensity within new T2 lesions were measured and evaluated for relationships to disability, atrophy, and clinical phenotype in long-term follow-up. RESULTS: Significant T2 activity was observed for 8 to 10 weeks beyond contrast enhancement, which suggests that T2 MR imaging is sensitive to noninflammatory processes such as degeneration and repair. Larger lesions showed longer subacute phases but disproportionally more recovery. Patients with smaller average peak lesion size showed trends toward greater disability and proportional residual damage. Higher rates of disability or atrophy were associated with subjects whose lesions showed greater residual hyperintensity. CONCLUSION: Smaller lesions appeared disproportionally more damaging than larger lesions, with lesions in progressive MS smaller and of shorter activity than in relapsing-remitting MS. Associations of lesion dynamics with rates of atrophy and disability and clinical subtype suggest that changes in lesion dynamics may represent a shift from inflammatory toward degenerative disease activity and greater proximity to a progressive stage, possibly allowing staging of the progression of MS earlier, before atrophy or disability develops.
A method to estimate the magnitude MR data from several noisy samples is presented. It is based on the Linear Minimum Mean Squared Error (LMMSE) estimator for the Rician noise model when several scanning repetitions are available. This method gives a closed-form analytical solution that takes into account the probability distribution of the data as well as the existing level of noise, showing a better performance than methods such as the average or the median.
Donnell LJO, Westin CF, Golby AJ. Tract-based morphometry. Med Image Comput Comput Assist Interv. 2007;10(Pt 2):161–8.
Multisubject statistical analyses of diffusion tensor images in regions of specific white matter tracts have commonly measured only the mean value of a scalar invariant such as the fractional anisotropy (FA), ignoring the spatial variation of FA along the length of fiber tracts. We propose to instead perform tract-based morphometry (TBM), or the statistical analysis of diffusion MRI data in an anatomical tract-based coordinate system. We present a method for automatic generation of white matter tract arc length parameterizations, based on learning a fiber bundle model from tractography from multiple subjects. Our tract-based coordinate system enables TBM for the detection of white matter differences in groups of subjects. We present example TBM results from a study of interhemispheric differences in FA.
This paper addresses the problem of image segmentation by means of active contours, whose evolution is driven by the gradient flow derived from an energy functional that is based on the Bhattacharyya distance. In particular, given the values of a photometric variable (or of a set thereof), which is to be used for classifying the image pixels, the active contours are designed to converge to the shape that results in maximal discrepancy between the empirical distributions of the photometric variable inside and outside of the contours. The above discrepancy is measured by means of the Bhattacharyya distance that proves to be an extremely useful tool for solving the problem at hand. The proposed methodology can be viewed as a generalization of the segmentation methods, in which active contours maximize the difference between a finite number of empirical moments of the "inside" and "outside" distributions. Furthermore, it is shown that the proposed methodology is very versatile and flexible in the sense that it allows one to easily accommodate a diversity of the image features based on which the segmentation should be performed. As an additional contribution, a method for automatically adjusting the smoothness properties of the empirical distributions is proposed. Such a procedure is crucial in situations when the number of data samples (supporting a certain segmentation class) varies considerably in the course of the evolution of the active contour. In this case, the smoothness properties of the empirical distributions have to be properly adjusted to avoid either over- or underestimation artifacts. Finally, a number of relevant segmentation results are demonstrated and some further research directions are discussed.
We present an algorithm to generate samples from probability distributions on the space of curves. We view a traditional curve evolution energy functional as a negative log probability distribution and sample from it using a Markov chain Monte Carlo (MCMC) algorithm. We define a proposal distribution by generating smooth perturbations to the normal of the curve and show how to compute the transition probabilities to ensure that the samples come from the posterior distribution. We demonstrate some advantages of sampling methods such as robustness to local minima, better characterization of multi-modal distributions, access to some measures of estimation error, and ability to easily incorporate constraints on the curve.
We argue that registration should be thought of as a means to an end, and not as a goal by itself. In particular, we consider the problem of predicting the locations of hidden labels of a test image using observable features, given a training set with both the hidden labels and observable features. For example, the hidden labels could be segmentation labels or activation regions in fMRI, while the observable features could be sulcal geometry or MR intensity. We analyze a probabilistic framework for computing an optimal atlas, and the subsequent registration of a new subject using only the observable features to optimize the hidden label alignment to the training set. We compare two approaches for co-registering training images for the atlas construction: the traditional approach of only using observable features and a novel approach of only using hidden labels. We argue that the alternative approach is superior particularly when the relationship between the hidden labels and observable features is complex and unknown. As an application, we consider the task of registering cortical folds to optimize Brodmann area localization. We show that the alignment of the Brodmann areas improves by up to 25% when using the alternative atlas compared with the traditional atlas. To the best of our knowledge, these are the most accurate Brodmann area localization results (achieved via cortical fold registration) reported to date.
We propose a new white matter atlas creation method that learns a model of the common white matter structures present in a group of subjects. We demonstrate that our atlas creation method, which is based on group spectral clustering of tractography, discovers structures corresponding to expected white matter anatomy such as the corpus callosum, uncinate fasciculus, cingulum bundles, arcuate fasciculus, and corona radiata. The white matter clusters are augmented with expert anatomical labels and stored in a new type of atlas that we call a high-dimensional white matter atlas. We then show how to perform automatic segmentation of tractography from novel subjects by extending the spectral clustering solution, stored in the atlas, using the Nystrom method. We present results regarding the stability of our method and parameter choices. Finally we give results from an atlas creation and automatic segmentation experiment. We demonstrate that our automatic tractography segmentation identifies corresponding white matter regions across hemispheres and across subjects, enabling group comparison of white matter anatomy.
Organ motion compensation in image-guided therapy is an active area of research. However, there has been little research on motion tracking and compensation in magnetic resonance imaging (MRI)-guided therapy. In this paper, we present a method to track a moving organ in MRI and control an active mechanical device for motion compensation. The method proposed is based on MRI navigator echo tracking enhanced by Kalman filtering for noise robustness. We also developed an extrapolation scheme to resolve any discrepancies between tracking and device control sampling rates. The algorithm was tested in a simulation study using a phantom and an active mechanical tool holder. We found that the method is feasible to use in a clinical MRI scanner with sufficient accuracy (0.36 mm to 1.51 mm depending on the range of phantom motion) and is robust to noise. The method proposed may be useful in MRI-guided targeted therapy, such as focused ultrasound therapy for a moving organ.