Publications by Year: 2011

Zhu L, Gao Y, Mohan V, Stillman A, Faber T, Tannenbaum A. Estimation of Myocardial Volume at Risk from CT Angiography. Proc SPIE Int Soc Opt Eng. 2011;7963 :79632A-79632A6.Abstract
The determination of myocardial volume at risk distal to coronary stenosis provides important information for prognosis and treatment of coronary artery disease. In this paper, we present a novel computational framework for estimating the myocardial volume at risk in computed tomography angiography (CTA) imagery. Initially, epicardial and endocardial surfaces, and coronary arteries are extracted using an active contour method. Then, the extracted coronary arteries are projected onto the epicardial surface, and each point on this surface is associated with its closest coronary artery using the geodesic distance measurement. The likely myocardial region at risk on the epicardial surface caused by a stenosis is approximated by the region in which all its inner points are associated with the sub-branches distal to the stenosis on the coronary artery tree. Finally, the likely myocardial volume at risk is approximated by the volume in between the region at risk on the epicardial surface and its projection on the endocardial surface, which is expected to yield computational savings over risk volume estimation using the entire image volume. Furthermore, we expect increased accuracy since, as compared to prior work using the Euclidean distance, we employ the geodesic distance in this work. The experimental results demonstrate the effectiveness of the proposed approach on pig heart CTA datasets.
Tokuda J, Mamata H, Gill RR, Hata N, Kikinis R, Padera RF, Lenkinski RE, Sugarbaker DJ, Hatabu H. Impact of nonrigid motion correction technique on pixel-wise pharmacokinetic analysis of free-breathing pulmonary dynamic contrast-enhanced MR imaging. J Magn Reson Imaging. 2011;33 (4) :968-73.Abstract
PURPOSE: To investigates the impact of nonrigid motion correction on pixel-wise pharmacokinetic analysis of free-breathing DCE-MRI in patients with solitary pulmonary nodules (SPNs). Misalignment of focal lesions due to respiratory motion in free-breathing dynamic contrast-enhanced MRI (DCE-MRI) precludes obtaining reliable time-intensity curves, which are crucial for pharmacokinetic analysis for tissue characterization. MATERIALS AND METHODS: Single-slice 2D DCE-MRI was obtained in 15 patients. Misalignments of SPNs were corrected using nonrigid B-spline image registration. Pixel-wise pharmacokinetic parameters K(trans) , v(e) , and k(ep) were estimated from both original and motion-corrected DCE-MRI by fitting the two-compartment pharmacokinetic model to the time-intensity curve obtained in each pixel. The "goodness-of-fit" was tested with χ(2) -test in pixel-by-pixel basis to evaluate the reliability of the parameters. The percentages of reliable pixels within the SPNs were compared between the original and motion-corrected DCE-MRI. In addition, the parameters obtained from benign and malignant SPNs were compared. RESULTS: The percentage of reliable pixels in the motion-corrected DCE-MRI was significantly larger than the original DCE-MRI (P = 4 × 10(-7) ). Both K(trans) and k(ep) derived from the motion-corrected DCE-MRI showed significant differences between benign and malignant SPNs (P = 0.024, 0.015). CONCLUSION: The study demonstrated the impact of nonrigid motion correction technique on pixel-wise pharmacokinetic analysis of free-breathing DCE-MRI in SPNs.
O'Donnell LJ, Westin C-F. An Introduction to Diffusion Tensor Image Analysis. Neurosurg Clin N Am. 2011;22 (2) :185-96.Abstract

Diffusion tensor magnetic resonance imaging (DTI) is a relatively new technology that is popular for imaging the white matter of the brain. This article provides a basic and broad overview of DTI to enable the reader to develop an intuitive understanding of these types of data, and an awareness of their strengths and weaknesses.

Risholm P, Golby AJ, Wells III WM. Multimodal Image Registration for Preoperative Planning and Image-guided Neurosurgical Procedures. Neurosurg Clin N Am. 2011;22 (2) :197-206.Abstract

Image registration is the process of transforming images acquired at different time points, or with different imaging modalities, into the same coordinate system. It is an essential part of any neurosurgical planning and navigation system because it facilitates combining images with important complementary, structural, and functional information to improve the information based on which a surgeon makes critical decisions. Brigham and Women's Hospital (BWH) has been one of the pioneers in developing intraoperative registration methods for aligning preoperative and intraoperative images of the brain. This article presents an overview of intraoperative registration and highlights some recent developments at BWH.

Lee J, Lankton S, Tannenbaum A. Object tracking and target reacquisition based on 3-D range data for moving vehicles. IEEE Trans Image Process. 2011;20 (10) :2912-24.Abstract
In this paper, we propose an approach for tracking an object of interest based on 3-D range data. We employ particle filtering and active contours to simultaneously estimate the global motion of the object and its local deformations. The proposed algorithm takes advantage of range information to deal with the challenging (but common) situation in which the tracked object disappears from the image domain entirely and reappears later. To cope with this problem, a method based on principle component analysis (PCA) of shape information is proposed. In the proposed method, if the target disappears out of frame, shape similarity energy is used to detect target candidates that match a template shape learned online from previously observed frames. Thus, we require no a priori knowledge of the target's shape. Experimental results show the practical applicability and robustness of the proposed algorithm in realistic tracking scenarios.
Menze BH, Van Leemput K, Honkela A, Konukoglu E, Weber M-A, Ayache N, Golland P. A generative approach for image-based modeling of tumor growth. Inf Process Med Imaging. 2011;22 :735-47.Abstract
Extensive imaging is routinely used in brain tumor patients to monitor the state of the disease and to evaluate therapeutic options. A large number of multi-modal and multi-temporal image volumes is acquired in standard clinical cases, requiring new approaches for comprehensive integration of information from different image sources and different time points. In this work we propose a joint generative model of tumor growth and of image observation that naturally handles multimodal and longitudinal data. We use the model for analyzing imaging data in patients with glioma. The tumor growth model is based on a reaction-diffusion framework. Model personalization relies only on a forward model for the growth process and on image likelihood. We take advantage of an adaptive sparse grid approximation for efficient inference via Markov Chain Monte Carlo sampling. The approach can be used for integrating information from different multi-modal imaging protocols and can easily be adapted to other tumor growth models.
Glasser MF, Van Essen DC. Mapping human cortical areas in vivo based on myelin content as revealed by T1- and T2-weighted MRI. J Neurosci. 2011;31 (32) :11597-616.Abstract
Noninvasively mapping the layout of cortical areas in humans is a continuing challenge for neuroscience. We present a new method of mapping cortical areas based on myelin content as revealed by T1-weighted (T1w) and T2-weighted (T2w) MRI. The method is generalizable across different 3T scanners and pulse sequences. We use the ratio of T1w/T2w image intensities to eliminate the MR-related image intensity bias and enhance the contrast to noise ratio for myelin. Data from each subject were mapped to the cortical surface and aligned across individuals using surface-based registration. The spatial gradient of the group average myelin map provides an observer-independent measure of sharp transitions in myelin content across the surface--i.e., putative cortical areal borders. We found excellent agreement between the gradients of the myelin maps and the gradients of published probabilistic cytoarchitectonically defined cortical areas that were registered to the same surface-based atlas. For other cortical regions, we used published anatomical and functional information to make putative identifications of dozens of cortical areas or candidate areas. In general, primary and early unimodal association cortices are heavily myelinated and higher, multimodal, association cortices are more lightly myelinated, but there are notable exceptions in the literature that are confirmed by our results. The overall pattern in the myelin maps also has important correlations with the developmental onset of subcortical white matter myelination, evolutionary cortical areal expansion in humans compared with macaques, postnatal cortical expansion in humans, and maps of neuronal density in non-human primates.
Risholm P, Ross J, Washko GR, Wells III WM. Probabilistic Elastography: Estimating Lung Elasticity. Inf Process Med Imaging. 2011;22 :699-710.Abstract

We formulate registration-based elastography in a probabilistic framework and apply it to study lung elasticity in the presence of emphysematous and fibrotic tissue. The elasticity calculations are based on a Finite Element discretization of a linear elastic biomechanical model. We marginalize over the boundary conditions (deformation) of the biomechanical model to determine the posterior distribution over elasticity parameters. Image similarity is included in the likelihood, an elastic prior is included to constrain the boundary conditions, while a Markov model is used to spatially smooth the inhomogeneous elasticity. We use a Markov Chain Monte Carlo (MCMC) technique to characterize the posterior distribution over elasticity from which we extract the most probable elasticity as well as the uncertainty of this estimate. Even though registration-based lung elastography with inhomogeneous elasticity is challenging due the problem's highly underdetermined nature and the sparse image information available in lung CT, we show promising preliminary results on estimating lung elasticity contrast in the presence of emphysematous and fibrotic tissue.

Wassermann D, Rathi Y, Bouix S, Kubicki M, Kikinis R, Shenton M, Westin C-F. White Matter Bundle Registration and Population Analysis Based on Gaussian Processes. Inf Process Med Imaging. 2011;22 :320-32.Abstract

This paper proposes a method for the registration of white matter tract bundles traced from diffusion images and its extension to atlas generation, Our framework is based on a Gaussian process representation of tract density maps. Such a representation avoids the need for point-to-point correspondences, is robust to tract interruptions and reconnections and seamlessly handles the comparison and combination of white matter tract bundles. Moreover, being a parametric model, this approach has the potential to be defined in the Gaussian processes' parameter space, without the need for resampling the fiber bundles during the registration process. We use the similarity measure of our Gaussian process framework, which is in fact an inner product, to drive a diffeomorphic registration algorithm between two sets of homologous bundles which is not biased by point-to-point correspondences or the parametrization of the tracts. We estimate a dense deformation of the underlying white matter using the bundles as anatomical landmarks and obtain a population atlas of those fiber bundles. Finally we test our results in several different bundles obtained from in-vivo data.

Risholm P, Balter J, Wells III WM. Estimation of delivered dose in radiotherapy: the influence of registration uncertainty. Med Image Comput Comput Assist Interv. 2011;14 (Pt 1) :548-55.Abstract

We present a probabilistic framework to estimate the accumulated radiation dose and the corresponding dose uncertainty that is delivered to important anatomical structures, e.g. the primary tumor and healthy surrounding organs, during radiotherapy. The dose uncertainty we report is a direct result of uncertainties in the estimates of the deformation which aligns the daily cone-beam CT images with the planning CT. The accumulated radiation dose is an important measure to monitor during treatment, in particular to see if it significantly deviates from the planned dose which might indicate that either the patient was not properly positioned before treatment or that the anatomy has changed due to the treatment. In the case of the latter, the treatment plan should be adaptively changed to align with the current patient anatomy. We estimate the accumulated dose distribution, and its uncertainty, retrospectively on a dataset acquired during treatment of cancer in the neck and show the dose distributions in the form of dose volume histograms.

Gooya A, Pohl KM, Bilello M, Biros G, Davatzikos C. Joint segmentation and deformable registration of brain scans guided by a tumor growth model. Med Image Comput Comput Assist Interv. 2011;14 (Pt 2) :532-40.Abstract
This paper presents an approach for joint segmentation and deformable registration of brain scans of glioma patients to a normal atlas. The proposed method is based on the Expectation Maximization (EM) algorithm that incorporates a glioma growth model for atlas seeding, a process which modifies the normal atlas into one with a tumor and edema. The modified atlas is registered into the patient space and utilized for the posterior probability estimation of various tissue labels. EM iteratively refines the estimates of the registration parameters, the posterior probabilities of tissue labels and the tumor growth model parameters. We have applied this approach to 10 glioma scans acquired with four Magnetic Resonance (MR) modalities (T1, T1-CE, T2 and FLAIR) and validated the result by comparing them to manual segmentations by clinical experts. The resulting segmentations look promising and quantitatively match well with the expert provided ground truth.
Gaonkar B, Pohl K, Davatzikos C. Pattern based morphometry. Med Image Comput Comput Assist Interv. 2011;14 (Pt 2) :459-66.Abstract
Voxel based morphometry (VBM) is widely used in the neuroimaging community to infer group differences in brain morphology. VBM is effective in quantifying group differences highly localized in space. However it is not equally effective when group differences might be based on interactions between multiple brain networks. We address this by proposing a new framework called pattern based morphometry (PBM). PBM is a data driven technique. It uses a dictionary learning algorithm to extract global patterns that characterize group differences. We test this approach on simulated and real data obtained from ADNI. In both cases PBM is able to uncover complex global patterns effectively.
Tristán-Vega A, Westin C-F. Probabilistic ODF estimation from reduced HARDI data with sparse regularization. Med Image Comput Comput Assist Interv. 2011;14 (Pt 2) :182-90.Abstract
High Angular Resolution Diffusion Imaging (HARDI) demands a higher amount of data measurements compared to Diffusion Tensor Imaging (DTI), restricting its use in practice. We propose to represent the probabilistic Orientation Distribution Function (ODF) in the frame of Spherical Wavelets (SW), where it is highly sparse. From a reduced subset of measurements (nearly four times less than the standard for HARDI), we pose the estimation as an inverse problem with sparsity regularization. This allows the fast computation of a positive, unit-mass, probabilistic ODF from 14-16 samples, as we show with both synthetic diffusion signals and real HARDI data with typical parameters.
Maddah M, Miller JV, Sullivan EV, Pfefferbaum A, Rohlfing T. Sheet-like white matter fiber tracts: representation, clustering, and quantitative analysis. Med Image Comput Comput Assist Interv. 2011;14 (Pt 2) :191-9.Abstract
We introduce an automated and probabilistic method for subject-specific segmentation of sheet-like fiber tracts. In addition to clustering of trajectories into anatomically meaningful bundles, the method provides statistics of diffusion measures by establishing point correspondences on the estimated medial representation of each bundle. We also introduce a new approach for medial surface generation of sheet-like fiber bundles in order too initialize the proposed clustering algorithm. Applying the new method to a population study of brain aging on 24 subjects demonstrates the capabilities and strengths of the algorithm in identifying and visualizing spatial patterns of group differences.
Rathi Y, Michailovich OV, Setsompop K, Bouix S, Shenton ME, Westin C-F. Sparse Multi-shell Diffusion Imaging. Med Image Comput Comput Assist Interv. 2011;14 (Pt 2) :58-65.Abstract

Diffusion magnetic resonance imaging (dMRI) is an important tool that allows non-invasive investigation of neural architecture of the brain. The data obtained from these in-vivo scans provides important information about the integrity and connectivity of neural fiber bundles in the brain. A multi-shell imaging (MSI) scan can be of great value in the study of several psychiatric and neurological disorders, yet its usability has been limited due to the long acquisition times required. A typical MSI scan involves acquiring a large number of gradient directions for the 2 (or more) spherical shells (several b-values), making the acquisition time significantly long for clinical application. In this work, we propose to use results from the theory of compressive sampling and determine the minimum number of gradient directions required to attain signal reconstruction similar to a traditional MSI scan. In particular, we propose a generalization of the single shell spherical ridgelets basis for sparse representation of multi shell signals. We demonstrate its efficacy on several synthetic and in-vivo data sets and perform quantitative comparisons with solid spherical harmonics based representation. Our preliminary results show that around 20-24 directions per shell are enough for robustly recovering the diffusion propagator.

Kikinis R, Pieper S. 3D Slicer as a tool for interactive brain tumor segmentation. Conf Proc IEEE Eng Med Biol Soc. 2011;2011 :6982-4.Abstract
User interaction is required for reliable segmentation of brain tumors in clinical practice and in clinical research. By incorporating current research tools, 3D Slicer provides a set of interactive, easy to use tools that can be efficiently used for this purpose. One of the modules of 3D Slicer is an interactive editor tool, which contains a variety of interactive segmentation effects. Use of these effects for fast and reproducible segmentation of a single glioblastoma from magnetic resonance imaging data is demonstrated. The innovation in this work lies not in the algorithm, but in the accessibility of the algorithm because of its integration into a software platform that is practical for research in a clinical setting.
Huang J, Gholami B, Agar NYR, Norton I, Haddad WM, Tannenbaum AR. Classification of astrocytomas and oligodendrogliomas from mass spectrometry data using sparse kernel machines. Conf Proc IEEE Eng Med Biol Soc. 2011;2011 :7965-8.Abstract
Glioma histologies are the primary factor in prognostic estimates and are used in determining the proper course of treatment. Furthermore, due to the sensitivity of cranial environments, real-time tumor-cell classification and boundary detection can aid in the precision and completeness of tumor resection. A recent improvement to mass spectrometry known as desorption electrospray ionization operates in an ambient environment without the application of a preparation compound. This allows for a real-time acquisition of mass spectra during surgeries and other live operations. In this paper, we present a framework using sparse kernel machines to determine a glioma sample's histopathological subtype by analyzing its chemical composition acquired by desorption electrospray ionization mass spectrometry.
Gholami B, Agar NYR, Jolesz FA, Haddad WM, Tannenbaum AR. A compressive sensing approach for glioma margin delineation using mass spectrometry. Conf Proc IEEE Eng Med Biol Soc. 2011;2011 :5682-5.Abstract
Surgery, and specifically, tumor resection, is the primary treatment for most patients suffering from brain tumors. Medical imaging techniques, and in particular, magnetic resonance imaging are currently used in diagnosis as well as image-guided surgery procedures. However, studies show that computed tomography and magnetic resonance imaging fail to accurately identify the full extent of malignant brain tumors and their microscopic infiltration. Mass spectrometry is a well-known analytical technique used to identify molecules in a given sample based on their mass. In a recent study, it is proposed to use mass spectrometry as an intraoperative tool for discriminating tumor and non-tumor tissue. Integration of mass spectrometry with the resection module allows for tumor resection and immediate molecular analysis. In this paper, we propose a framework for tumor margin delineation using compressive sensing. Specifically, we show that the spatial distribution of tumor cell concentration can be efficiently reconstructed and updated using mass spectrometry information from the resected tissue. In addition, our proposed framework is model-free, and hence, requires no prior information of spatial distribution of the tumor cell concentration.
Nakhmani A, Tannenbaum A. Particle Filtering with Region-based Matching for Tracking of Partially Occluded and Scaled Targets. SIAM J Imaging Sci. 2011;4 (1) :220-242.Abstract
Visual tracking of arbitrary targets in clutter is important for a wide range of military and civilian applications. We propose a general framework for the tracking of scaled and partially occluded targets, which do not necessarily have prominent features. The algorithm proposed in the present paper utilizes a modified normalized cross-correlation as the likelihood for a particle filter. The algorithm divides the template, selected by the user in the first video frame, into numerous patches. The matching process of these patches by particle filtering allows one to handle the target's occlusions and scaling. Experimental results with fixed rectangular templates show that the method is reliable for videos with nonstationary, noisy, and cluttered background, and provides accurate trajectories in cases of target translation, scaling, and occlusion.
Fedorov A, Li X, Pohl KM, Bouix S, Styner M, Addicott M, Wyatt C, Daunais JB, Wells WM, Kikinis R. Atlas-guided Segmentation of Vervet Monkey Brain MRI. Open Neuroimag J. 2011;5 :186-97.Abstract

The vervet monkey is an important nonhuman primate model that allows the study of isolated environmental factors in a controlled environment. Analysis of monkey MRI often suffers from lower quality images compared with human MRI because clinical equipment is typically used to image the smaller monkey brain and higher spatial resolution is required. This, together with the anatomical differences of the monkey brains, complicates the use of neuroimage analysis pipelines tuned for human MRI analysis. In this paper we developed an open source image analysis framework based on the tools available within the 3D Slicer software to support a biological study that investigates the effect of chronic ethanol exposure on brain morphometry in a longitudinally followed population of male vervets. We first developed a computerized atlas of vervet monkey brain MRI, which was used to encode the typical appearance of the individual brain structures in MRI and their spatial distribution. The atlas was then used as a spatial prior during automatic segmentation to process two longitudinal scans per subject. Our evaluation confirms the consistency and reliability of the automatic segmentation. The comparison of atlas construction strategies reveals that the use of a population-specific atlas leads to improved accuracy of the segmentation for subcortical brain structures. The contribution of this work is twofold. First, we describe an image processing workflow specifically tuned towards the analysis of vervet MRI that consists solely of the open source software tools. Second, we develop a digital atlas of vervet monkey brain MRIs to enable similar studies that rely on the vervet model.