Publications

Diffusion MR imaging has received increasing attention in the neuroimaging community, as it yields new insights into the microstructural organization of white matter that are not available with conventional MRI techniques. While the technology has enormous potential, diffusion MRI suffers from a unique and complex set of image quality problems, limiting the sensitivity of studies and reducing the accuracy of findings. Furthermore, the acquisition time for diffusion MRI is longer than conventional MRI due to the need for multiple acquisitions to obtain directionally encoded Diffusion Weighted Images (DWI). This leads to increased motion artifacts, reduced signal-to-noise ratio (SNR), and increased proneness to a wide variety of artifacts, including eddy-current and motion artifacts, "venetian blind" artifacts, as well as slice-wise and gradient-wise inconsistencies. Such artifacts mandate stringent Quality Control (QC) schemes in the processing of diffusion MRI data. Most existing QC procedures are conducted in the DWI domain and/or on a voxel level, but our own experiments show that these methods often do not fully detect and eliminate certain types of artifacts, often only visible when investigating groups of DWI’s or a derived diffusion model, such as the most-employed diffusion tensor imaging (DTI). Here, we propose a novel regional QC measure in the DTI domain that employs the entropy of the regional distribution of the principal directions (PD). The PD entropy quantifies the scattering and spread of the principal diffusion directions and is invariant to the patient’s position in the scanner. High entropy value indicates that the PDs are distributed relatively uniformly, while low entropy value indicates the presence of clusters in the PD distribution. The novel QC measure is intended to complement the existing set of QC procedures by detecting and correcting residual artifacts. Such residual artifacts cause directional bias in the measured PD and here called dominant direction artifacts. Experiments show that our automatic method can reliably detect and potentially correct such artifacts, especially the ones caused by the vibrations of the scanner table during the scan. The results further indicate the usefulness of this method for general quality assessment in DTI studies.

PURPOSE: Performing lobe-based quantitative analysis of the lung in computed tomography (CT) scans can assist in efforts to better characterize complex diseases such as chronic obstructive pulmonary disease (COPD). While airways and vessels can help to indicate the location of lobe boundaries, segmentations of these structures are not always available, so methods to define the lobes in the absence of these structures are desirable. METHODS: The authors present a fully automatic lung lobe segmentation algorithm that is effective in volumetric inspiratory and expiratory computed tomography (CT) datasets. The authors rely on ridge surface image features indicating fissure locations and a novel approach to modeling shape variation in the surfaces defining the lobe boundaries. The authors employ a particle system that efficiently samples ridge surfaces in the image domain and provides a set of candidate fissure locations based on the Hessian matrix. Following this, lobe boundary shape models generated from principal component analysis (PCA) are fit to the particles data to discriminate between fissure and nonfissure candidates. The resulting set of particle points are used to fit thin plate spline (TPS) interpolating surfaces to form the final boundaries between the lung lobes. RESULTS: The authors tested algorithm performance on 50 inspiratory and 50 expiratory CT scans taken from the COPDGene study. Results indicate that the authors’ algorithm performs comparably to pulmonologist-generated lung lobe segmentations and can produce good results in cases with accessory fissures, incomplete fissures, advanced emphysema, and low dose acquisition protocols. Dice scores indicate that only 29 out of 500 (5.85%) lobes showed Dice scores lower than 0.9. Two different approaches for evaluating lobe boundary surface discrepancies were applied and indicate that algorithm boundary identification is most accurate in the vicinity of fissures detectable on CT. CONCLUSIONS: The proposed algorithm is effective for lung lobe segmentation in absence of auxiliary structures such as vessels and airways. The most challenging cases are those with mostly incomplete, absent, or near-absent fissures and in cases with poorly revealed fissures due to high image noise. However, the authors observe good performance even in the majority of these cases.

The main contribution of this work is the careful syntactical definition of major white matter tracts in the human brain based on a neuroanatomist’s expert knowledge. We present a technique to formally describe white matter tracts and to automatically extract them from diffusion MRI data. The framework is based on a novel query language with a near-to-English textual syntax. This query language allows us to construct a dictionary of anatomical definitions describing white matter tracts. The definitions include adjacent gray and white matter regions, and rules for spatial relations. This enables automated coherent labeling of white matter anatomy across subjects. We use our method to encode anatomical knowledge in human white matter describing 10 association and 8 projection tracts per hemisphere and 7 commissural tracts. The technique is shown to be comparable in accuracy to manual labeling. We present results applying this framework to create a white matter atlas from 77 healthy subjects, and we use this atlas in a proof-of-concept study to detect tract changes specific to schizophrenia.

Estimation of the diffusion propagator from a sparse set of diffusion MRI (dMRI) measurements is a field of active research. Sparse reconstruction methods propose to reduce scan time and are particularly suitable for scanning un-coperative patients. Recent work on reconstructing the diffusion signal from very few measurements using compressed sensing based techniques has focussed on propagator (or signal) estimation at each voxel independently. However, the goal of many neuroscience studies is to use tractography to study the pathology in white matter fiber tracts. Thus, in this work, we propose a joint framework for robust estimation of the diffusion propagator from sparse measurements while simultaneously tracing the white matter tracts. We propose to use a novel multi-tensor model of diffusion which incorporates the biexponential radial decay of the signal. Our preliminary results on in-vivo data show that the proposed method produces consistent and reliable fiber tracts from very few gradient directions while simultaneously estimating the bi-exponential decay of the diffusion propagator.

This paper presents a novel pipeline for the registration of diffusion tensor images (DTI) with large pathological variations to normal controls based on the use of a novel feature map derived from white matter (WM) fiber tracts. The research presented aims towards an atlas based DTI analysis of subjects with considerable brain pathologies such as tumors or hydrocephalus. In this paper, we propose a novel feature map that is robust against variations in WM fiber tract integrity and use these feature maps to determine a landmark correspondence using a 3D point correspondence algorithm. This correspondence drives a deformation field computed using Gaussian radial basis functions(RBF). This field is employed as an initialization to a standard deformable registration method like demons. We present early preliminary results on the registration of a normal control dataset to a dataset with abnormally enlarged lateral ventricles affected by fatal demyelinating Krabbe disease. The results are analyzed based on a regional tensor matching criterion and a visual assessment of overlap of major WM fiber tracts. While further evaluation and improvements are necessary, the results presented in this paper highlight the potential of our method in handling registration of subjects with severe WM pathology.

In this article we investigate the suitability of a manifold learning technique to classify different types of emphysema based on embedded Probabilistic PCA (PPCA). Our approach finds the most discriminant linear space for each emphysema pattern against the remaining patterns where lung CT image patches can be embedded. In this embedded space, we train a PPCA model for each pattern. The main novelty of our technique is that it is possible to compute the class membership posterior probability for each emphysema pattern rather than a hard assignment as it is typically done by other approaches. We tested our algorithm with six emphysema patterns using a data set of 1337 CT training patches. Using a 10-fold cross validation experiment, an average recall rate of 69% is achieved when the posterior probability is greater than 75%. A quantitative comparison with a texture-based approach based on Local Binary Patterns and with an approach based on local intensity distributions shows that our method is competitive. The analysis of full lungs using our approach shows a good visual agreement with the underlying emphysema types and a smooth spatial relation.

BACKGROUND: A common promoter polymorphism (rs35705950) in MUC5B, the gene encoding mucin 5B, is associated with idiopathic pulmonary fibrosis. It is not known whether this polymorphism is associated with interstitial lung disease in the general population. METHODS: We performed a blinded assessment of interstitial lung abnormalities detected in 2633 participants in the Framingham Heart Study by means of volumetric chest computed tomography (CT). We evaluated the relationship between the abnormalities and the genotype at the rs35705950 locus. RESULTS: Of the 2633 chest CT scans that were evaluated, interstitial lung abnormalities were present in 177 (7%). Participants with such abnormalities were more likely to have shortness of breath and chronic cough and reduced measures of total lung and diffusion capacity, as compared with participants without such abnormalities. After adjustment for covariates, for each copy of the minor rs35705950 allele, the odds of interstitial lung abnormalities were 2.8 times greater (95% confidence interval [CI], 2.0 to 3.9; P

We present a probability model for lung airways in computed tomography (CT) images. Lung airways are tubular structures that display specific features, such as low intensity and proximity to vessels and bronchial walls. From these features, the posterior probability for the airway feature space was computed using a Bayesian model based on 20 CT images from subjects with different degrees of Chronic Obstructive Pulmonary Disease (COPD). The likelihood probability was modeled using both a Gaussian distribution and a nonparametric kernel density estimation method. After exhaustive feature selection, good specificity and sensitivity were achieved in a cross-validation study for both the Gaussian (0.83, 0.87) and the nonparametric method (0.79, 0.89). The model generalizes well when trained using images from a late stage COPD group. This probability model may facilitate airway extraction and quantitative assessment of lung diseases, which is useful in many clinical and research settings.

RATIONALE: Angiographic investigation suggests that pulmonary vascular remodeling in smokers is characterized by distal pruning of the blood vessels. OBJECTIVES: Using volumetric computed tomography scans of the chest we sought to quantitatively evaluate this process and assess its clinical associations. METHODS: Pulmonary vessels were automatically identified, segmented, and measured. Total blood vessel volume (TBV) and the aggregate vessel volume for vessels less than 5 mm(2) (BV5) were calculated for all lobes. The lobe-specific BV5 measures were normalized to the TBV of that lobe and the nonvascular tissue volume (BV5/T(issue)V) to calculate lobe-specific BV5/TBV and BV5/T(issue)V ratios. Densitometric measures of emphysema were obtained using a Hounsfield unit threshold of -950 (%LAA-950). Measures of chronic obstructive pulmonary disease severity included single breath measures of diffusing capacity of carbon monoxide, oxygen saturation, the 6-minute-walk distance, St George’s Respiratory Questionnaire total score (SGRQ), and the body mass index, airflow obstruction, dyspnea, and exercise capacity (BODE) index. MEASUREMENTS AND MAIN RESULTS: The %LAA-950 was inversely related to all calculated vascular ratios. In multivariate models including age, sex, and %LAA-950, lobe-specific measurements of BV5/TBV were directly related to resting oxygen saturation and inversely associated with both the SGRQ and BODE scores. In similar multivariate adjustment lobe-specific BV5/T(issue)V ratios were inversely related to resting oxygen saturation, diffusing capacity of carbon monoxide, 6-minute-walk distance, and directly related to the SGRQ and BODE. CONCLUSIONS: Smoking-related chronic obstructive pulmonary disease is characterized by distal pruning of the small blood vessels (

BACKGROUND: In CT scans of smokers with COPD, the subsegmental airway wall area percent (WA%) is greater and more strongly correlated with FEV1 % predicted than WA% obtained in the segmental airways. Because emphysema is linked to loss of airway tethering and may limit airway expansion, increases in WA% may be related to emphysema and not solely to remodeling. We aimed to first determine whether the stronger association of subsegmental vs segmental WA% with FEV1 % predicted is mitigated by emphysema and, second, to assess the relationships among emphysema, WA%, and total bronchial area (TBA).