Publications

Advances in high-resolution computed tomography (CT) imaging are making a full three-dimensional analysis of the lungs feasible. In particular, airway morphology can be studied in vivo and quantitative metrics of airway size and shape can be extracted. The thickening process associated with the inflammatory response in the diseased lung can be quantified by means of image processing techniques that extract the airway lumen and airway wall. In this article, we give an overview of these imaging techniques and their diverse nature. We also offer a comprehensive view of the analysis pipeline for three-dimensional airway trees and a validation framework that is needed to compare different techniques.

This paper presents a tract-oriented analysis of diffusion tensor (DT) images of the human brain. We demonstrate that unlike the commonly used ROI-based methods for population studies, our technique is sensitive to the local variation of diffusivity parameters along the fiber tracts. We show the strength of the proposed approach in identifying the differences in schizophrenic data compared to controls. Statistically significant drops in fractional anisotropy are observed along the genu and bilaterally in the splenium, as well as an increase in principal eigenvalue in uncinate fasciculus. This is the first tract-oriented clinical study in which an anatomical atlas is used to guide the algorithm.

We introduce a versatile framework for characterizing and extracting salient structures in three-dimensional symmetric second-order tensor fields. The key insight is that degenerate lines in tensor fields, as defined by the standard topological approach, are exactly crease (ridge and valley) lines of a particular tensor invariant called mode. This reformulation allows us to apply well-studied approaches from scientific visualization or computer vision to the extraction of topological lines in tensor fields. More generally, this main result suggests that other tensor invariants, such as anisotropy measures like fractional anisotropy (FA), can be used in the same framework in lieu of mode to identify important structural properties in tensor fields. Our implementation addresses the specific challenge posed by the non-linearity of the considered scalar measures and by the smoothness requirement of the crease manifold computation. We use a combination of smooth reconstruction kernels and adaptive refinement strategy that automatically adjust the resolution of the analysis to the spatial variation of the considered quantities. Together, these improvements allow for the robust application of existing ridge line extraction algorithms in the tensor context of our problem. Results are proposed for a diffusion tensor MRI dataset, and for a benchmark stress tensor field used in engineering research.

Several recent studies explored the use of unsupervised segmentation methods for segmenting thalamic nuclei from diffusion tensor images. These methods provide a plausible segmentation on individual subjects; however, they do not address the problem of consistently identifying the same functional areas in a population. The lack of correspondence between the segmented nuclei make it more difficult to use the results from the unsupervised segmentation tools for morphometry. In this paper we present a novel segmentation algorithm to automatically segment the gray matter nuclei while ensuring consistency between subjects in a population. This new algorithm, referred to as Consistency Clustering, finds correspondence between the nuclei as the segmentation is achieved through a single model for the whole population, similar to the brain atlases experts use to identify thalamic nuclei.

We propose a Bayesian approach to incorporate anatomical information in the clustering of fiber trajectories. An expectation-maximization (EM) algorithm is used to cluster the trajectories, in which an atlas serves as the prior on the labels. The atlas guides the clustering algorithm and makes the resulting bundles anatomically meaningful. In addition, it provides the seed points for the tractography and initial settings of the EM algorithm. The proposed approach provides a robust and automated tool for tract-oriented analysis both in a single subject and over a population.

This paper introduces and analyzes a linear minimum mean square error (LMMSE) estimator using a Rician noise model and its recursive version (RLMMSE) for the restoration of diffusion weighted images. A method to estimate the noise level based on local estimations of mean or variance is used to automatically parametrize the estimator. The restoration performance is evaluated using quality indexes and compared to alternative estimation schemes. The overall scheme is simple, robust, fast, and improves estimations. Filtering diffusion weighted magnetic resonance imaging (DW-MRI) with the proposed methodology leads to more accurate tensor estimations. Real and synthetic datasets are analyzed.

A new method for noise filtering in images that follow a Rician model-with particular attention to magnetic resonance imaging-is proposed. To that end, we have derived a (novel) closed-form solution of the linear minimum mean square error (LMMSE) estimator for this distribution. Additionally, a set of methods that automatically estimate the noise power are developed. These methods use information of the sample distribution of local statistics of the image, such as the local variance, the local mean, and the local mean square value. Accordingly, the dynamic estimation of noise leads to a recursive version of the LMMSE, which shows a good performance in both noise cleaning and feature preservation. This paper also includes the derivation of the probability density function of several local sample statistics for the Rayleigh and Rician model, upon which the estimators are built.

This paper presents a segmentation method that extends geodesic active region methods by the incorporation of a statistical classifier trained using feature selection. The classifier provides class probability maps based on class representative local features, and the geodesic active region formulation enables the partitioning of the image according to the region information. We demonstrate automatic segmentation results of the myocardium in cardiac late gadolinium-enhanced magnetic resonance imaging (CE-MRI) data using coupled level set curve evolutions, in which the classifier is incorporated both from a region term and from a shape term from particle filtering. The results show potential for clinical studies of scar tissue in late CE-MRI data.

We introduce a fibre tract segmentation algorithm based on the geometric coherence of fibre orientations as indicated by a streamline flow model. The inference of local flow approximations motivates a pairwise consistency measure between fibre ODF maxima. We use this measure in a recursive algorithm to cluster consistent ODF maxima, leading to the segmentation of white matter pathways. The method requires minimal seeding compared to streamline tractography-based methods, and allows multiple tracts to pass through the same voxels. We illustrate the approach with a segmentation of the corpus callosum and one of the cortico-spinal tract, with each example seeded at a single voxel.

OBJECTIVE: Impairment of white matter connecting frontal and temporal cortices has been reported in schizophrenia. Yet, not much is known about the effects of age on fibers connecting these brain regions. Using diffusion tensor imaging tractography, we investigated the relationship between age and fiber integrity in patients with schizophrenia vs. healthy adults. METHODS: DTI tractography was used to create 3D reconstructions of the cingulum, uncinate and inferior occipito-frontal fasciculi in 27 patients with schizophrenia and 34 healthy volunteers (23-56 years of age, group-matched on age). Fractional anisotropy (FA), describing fiber integrity, was then calculated along the entire length of these tracts, and correlated with subjects’ age. RESULTS: Patients revealed a significant decline in FA with age in both the cingulum and uncinate, but not in the inferior occipito-frontal fasciculi. No statistically significant correlations were found in these fiber bundles in controls. CONCLUSIONS: These results suggest an age-associated reduction of frontal-temporal connectivity in schizophrenia, but not in healthy controls.