Publications

2015
Eriksson S, Lasič S, Nilsson M, Westin C-F, Topgaard D. NMR Diffusion-Encoding with Axial Symmetry and Variable Anisotropy: Distinguishing Between Prolate and Oblate Microscopic Diffusion Tensors with Unknown Orientation Distribution. J Chem Phys. 2015;142 (10) :104201.Abstract
We introduce a nuclear magnetic resonance method for quantifying the shape of axially symmetric microscopic diffusion tensors in terms of a new diffusion anisotropy metric, DΔ, which has unique values for oblate, spherical, and prolate tensor shapes. The pulse sequence includes a series of equal-amplitude magnetic field gradient pulse pairs, the directions of which are tailored to give an axially symmetric diffusion-encoding tensor b with variable anisotropy bΔ. Averaging of data acquired for a range of orientations of the symmetry axis of the tensor b renders the method insensitive to the orientation distribution function of the microscopic diffusion tensors. Proof-of-principle experiments are performed on water in polydomain lyotropic liquid crystals with geometries that give rise to microscopic diffusion tensors with oblate, spherical, and prolate shapes. The method could be useful for characterizing the geometry of fluid-filled compartments in porous solids, soft matter, and biological tissues.
Li M, Miller K, Joldes GR, Doyle B, Garlapati RR, Kikinis R, Wittek A. Patient-specific biomechanical model as whole-body CT image registration tool. Med Image Anal. 2015;22 (1) :22-34.Abstract
Whole-body computed tomography (CT) image registration is important for cancer diagnosis, therapy planning and treatment. Such registration requires accounting for large differences between source and target images caused by deformations of soft organs/tissues and articulated motion of skeletal structures. The registration algorithms relying solely on image processing methods exhibit deficiencies in accounting for such deformations and motion. We propose to predict the deformations and movements of body organs/tissues and skeletal structures for whole-body CT image registration using patient-specific non-linear biomechanical modelling. Unlike the conventional biomechanical modelling, our approach for building the biomechanical models does not require time-consuming segmentation of CT scans to divide the whole body into non-overlapping constituents with different material properties. Instead, a Fuzzy C-Means (FCM) algorithm is used for tissue classification to assign the constitutive properties automatically at integration points of the computation grid. We use only very simple segmentation of the spine when determining vertebrae displacements to define loading for biomechanical models. We demonstrate the feasibility and accuracy of our approach on CT images of seven patients suffering from cancer and aortic disease. The results confirm that accurate whole-body CT image registration can be achieved using a patient-specific non-linear biomechanical model constructed without time-consuming segmentation of the whole-body images.
Chen Z, Tie Y, Olubiyi O, Rigolo L, Mehrtash A, Norton I, Pasternak O, Rathi Y, Golby AJ, O'Donnell LJ. Reconstruction of the arcuate fasciculus for surgical planning in the setting of peritumoral edema using two-tensor unscented Kalman filter tractography. Neuroimage Clin. 2015;7 :815-22.Abstract
BACKGROUND: Diffusion imaging tractography is increasingly used to trace critical fiber tracts in brain tumor patients to reduce the risk of post-operative neurological deficit. However, the effects of peritumoral edema pose a challenge to conventional tractography using the standard diffusion tensor model. The aim of this study was to present a novel technique using a two-tensor unscented Kalman filter (UKF) algorithm to track the arcuate fasciculus (AF) in brain tumor patients with peritumoral edema. METHODS: Ten right-handed patients with left-sided brain tumors in the vicinity of language-related cortex and evidence of significant peritumoral edema were retrospectively selected for the study. All patients underwent 3-Tesla magnetic resonance imaging (MRI) including a diffusion-weighted dataset with 31 directions. Fiber tractography was performed using both single-tensor streamline and two-tensor UKF tractography. A two-regions-of-interest approach was applied to perform the delineation of the AF. Results from the two different tractography algorithms were compared visually and quantitatively. RESULTS: Using single-tensor streamline tractography, the AF appeared disrupted in four patients and contained few fibers in the remaining six patients. Two-tensor UKF tractography delineated an AF that traversed edematous brain areas in all patients. The volume of the AF was significantly larger on two-tensor UKF than on single-tensor streamline tractography (p < 0.01). CONCLUSIONS: Two-tensor UKF tractography provides the ability to trace a larger volume AF than single-tensor streamline tractography in the setting of peritumoral edema in brain tumor patients.
Garlapati RR, Mostayed A, Joldes GR, Wittek A, Doyle B, Miller K. Towards measuring neuroimage misalignment. Comput Biol Med. 2015;64 :12-23.Abstract
To enhance neuro-navigation, high quality pre-operative images must be registered onto intra-operative configuration of the brain. Therefore evaluation of the degree to which structures may remain misaligned after registration is critically important. We consider two Hausdorff Distance (HD)-based evaluation approaches: the edge-based HD (EBHD) metric and the Robust HD (RHD) metric as well as various commonly used intensity-based similarity metrics such as Mutual Information (MI), Normalised Mutual Information (NMI), Entropy Correlation Coefficient (ECC), Kullback-Leibler Distance (KLD) and Correlation Ratio (CR). We conducted the evaluation by applying known deformations to simple sample images and real cases of brain shift. We conclude that the intensity-based similarity metrics such as MI, NMI, ECC, KLD and CR do not correlate well with actual alignment errors, and hence are not useful for assessing misalignment. On the contrary, the EBHD and the RHD metrics correlated well with actual alignment errors; however, they have been found to underestimate the actual misalignment. We also note that it is beneficial to present HD results as a percentile-HD curve rather than a single number such as the 95-percentile HD. Percentile-HD curves present the full range of alignment errors and also facilitate the comparison of results obtained using different approaches. Furthermore, the qualities that should be possessed by an ideal evaluation metric were highlighted. Future studies could focus on developing such an evaluation metric.
Ohtani T, Bouix S, Lyall AE, Hosokawa T, Saito Y, Melonakos E, Westin C-F, Seidman LJ, Goldstein J, Mesholam-Gately R, et al. Abnormal White Matter Connections between Medial Frontal Regions Predict Symptoms in Patients with First Episode Schizophrenia. Cortex. 2015;71 :264-76.Abstract
INTRODUCTION: The medial orbitofrontal cortex (mOFC) and rostral part of anterior cingulate cortex (rACC) have been suggested to be involved in the neural network of salience and emotional processing, and associated with specific clinical symptoms in schizophrenia. Considering the schizophrenia dysconnectivity hypothesis, the connectivity abnormalities between mOFC and rACC might be associated with clinical characteristics in first episode schizophrenia patients (FESZ). METHODS: After parcellating mOFC into the anterior and posterior part, diffusion properties of the mOFC-rACC white matter connections for 21 patients with FESZ and 21 healthy controls (HCs) were examined using stochastic tractography, one of the most effective Diffusion Tensor Imaging (DTI) methods for examining tracts between adjacent gray matter (GM) regions. RESULTS: Fractional anisotropy (FA) reductions were observed in bilateral posterior, but not anterior mOFC-rACC connections (left: p < .0001; right: p < .0001) in FESZ compared to HCs. In addition, reduced FA in the left posterior mOFC-rACC connection was associated with more severe anhedonia-asociality (rho = -.633, p = .006) and total score (rho = -.520, p = .032) in the Scale for the Assessment of Negative Symptoms (SANS); reduced FA in the right posterior mOFC-rACC connection was associated with more severe affective flattening (rho = -.644, p = .005), total score (rho = -.535, p = .027) in SANS, hallucinations (rho = -.551, p = .018), delusions (rho = -.632, p = .005) and total score (rho = -.721, p = .001) in the Scale for the Assessment of Positive Symptoms (SAPS) in FESZ. CONCLUSIONS: The observed white matter abnormalities within the connections between mOFC and rACC might be associated with the psychopathology of the early stage of schizophrenia.
Ning L, Setsompop K, Michailovich O, Makris N, Westin C-F, Rathi Y. A Compressed-Sensing Approach for Super-Resolution Reconstruction of Diffusion MRI. Inf Process Med Imaging. 2015;24 :57-68.Abstract

We present an innovative framework for reconstructing high-spatial-resolution diffusion magnetic resonance imaging (dMRI) from multiple low-resolution (LR) images. Our approach combines the twin concepts of compressed sensing (CS) and classical super-resolution to reduce acquisition time while increasing spatial resolution. We use subpixel-shifted LR images with down-sampled and non-overlapping diffusion directions to reduce acquisition time. The diffusion signal in the high resolution (HR) image is represented in a sparsifying basis of spherical ridgelets to model complex fiber orientations with reduced number of measurements. The HR image is obtained as the solution of a convex optimization problem which can be solved using the proposed algorithm based on the alternating direction method of multipliers (ADMM). We qualitatively and quantitatively evaluate the performance of our method on two sets of in-vivo human brain data and show its effectiveness in accurately recovering very high resolution diffusion images.

Pujol S, Wells III WM, Pierpaoli C, Brun C, Gee J, Cheng G, Vemuri B, Commowick O, Prima S, Stamm A, et al. The DTI Challenge: Toward Standardized Evaluation of Diffusion Tensor Imaging Tractography for Neurosurgery. J Neuroimaging. 2015;25 (6) :875-82.Abstract
BACKGROUND AND PURPOSE: Diffusion tensor imaging (DTI) tractography reconstruction of white matter pathways can help guide brain tumor resection. However, DTI tracts are complex mathematical objects and the validity of tractography-derived information in clinical settings has yet to be fully established. To address this issue, we initiated the DTI Challenge, an international working group of clinicians and scientists whose goal was to provide standardized evaluation of tractography methods for neurosurgery. The purpose of this empirical study was to evaluate different tractography techniques in the first DTI Challenge workshop. METHODS: Eight international teams from leading institutions reconstructed the pyramidal tract in four neurosurgical cases presenting with a glioma near the motor cortex. Tractography methods included deterministic, probabilistic, filtered, and global approaches. Standardized evaluation of the tracts consisted in the qualitative review of the pyramidal pathways by a panel of neurosurgeons and DTI experts and the quantitative evaluation of the degree of agreement among methods. RESULTS: The evaluation of tractography reconstructions showed a great interalgorithm variability. Although most methods found projections of the pyramidal tract from the medial portion of the motor strip, only a few algorithms could trace the lateral projections from the hand, face, and tongue area. In addition, the structure of disagreement among methods was similar across hemispheres despite the anatomical distortions caused by pathological tissues. CONCLUSIONS: The DTI Challenge provides a benchmark for the standardized evaluation of tractography methods on neurosurgical data. This study suggests that there are still limitations to the clinical use of tractography for neurosurgical decision making.
Toews M, Wachinger C, Estepar RSJ, Wells III WM. A Feature-Based Approach to Big Data Analysis of Medical Images. Inf Process Med Imaging. 2015;24 :339-50.Abstract

This paper proposes an inference method well-suited to large sets of medical images. The method is based upon a framework where distinctive 3D scale-invariant features are indexed efficiently to identify approximate nearest-neighbor (NN) feature matches-in O (log N) computational complexity in the number of images N. It thus scales well to large data sets, in contrast to methods based on pair-wise image registration or feature matching requiring O(N) complexity. Our theoretical contribution is a density estimator based on a generative model that generalizes kernel density estimation and K-nearest neighbor (KNN) methods.. The estimator can be used for on-the-fly queries, without requiring explicit parametric models or an off-line training phase. The method is validated on a large multi-site data set of 95,000,000 features extracted from 19,000 lung CT scans. Subject-level classification identifies all images of the same subjects across the entire data set despite deformation due to breathing state, including unintentional duplicate scans. State-of-the-art performance is achieved in predicting chronic pulmonary obstructive disorder (COPD) severity across the 5-category GOLD clinical rating, with an accuracy of 89% if both exact and one-off predictions are considered correct.

Wachinger C, Golland P. Sampling from Determinantal Point Processes for Scalable Manifold Learning. Inf Process Med Imaging. 2015;24 :687-98.Abstract
High computational costs of manifold learning prohibit its application for large datasets. A common strategy to overcome this problem is to perform dimensionality reduction on selected landmarks and to successively embed the entire dataset with the Nyström method. The two main challenges that arise are: (i) the landmarks selected in non-Euclidean geometries must result in a low reconstruction error, (ii) the graph constructed from sparsely sampled landmarks must approximate the manifold well. We propose to sample the landmarks from determinantal distributions on non-Euclidean spaces. Since current determinantal sampling algorithms have the same complexity as those for manifold learning, we present an efficient approximation with linear complexity. Further, we recover the local geometry after the sparsification by assigning each landmark a local covariance matrix, estimated from the original point set. The resulting neighborhood selection .based on the Bhattacharyya distance improves the embedding of sparsely sampled manifolds. Our experiments show a significant performance improvement compared to state-of-the-art landmark selection techniques on synthetic and medical data.
Parisot S, Arslan S, Passerat-Palmbach J, Wells III WM, Rueckert D. Tractography-Driven Groupwise Multi-scale Parcellation of the Cortex. Inf Process Med Imaging. 2015;24 :600-12.Abstract

The analysis of the connectome of the human brain provides key insight into the brain's organisation and function, and its evolution in disease or ageing. Parcellation of the cortical surface into distinct regions in terms of structural connectivity is an essential step that can enable such analysis. The estimation of a stable connectome across a population of healthy subjects requires the estimation of a groupwise parcellation that can capture the variability of the connectome across the population. This problem has solely been addressed in the literature via averaging of connectivity profiles or finding correspondences between individual parcellations a posteriori. In this paper, we propose a groupwise parcellation method of the cortex based on diffusion MR images (dMRI). We borrow ideas from the area of cosegmentation in computer vision and directly estimate a consistent parcellation across different subjects and scales through a spectral clustering approach. The parcellation is driven by the tractography connectivity profiles, and information between subjects and across scales. Promising qualitative and quantitative results on a sizeable data-set demonstrate the strong potential of the method.

Vaughan T, Lasso A, Ungi T, Fichtinger G. Hole filling with oriented sticks in ultrasound volume reconstruction. J Med Imaging (Bellingham). 2015;2 (3) :034002.Abstract
Volumes reconstructed from tracked planar ultrasound images often contain regions where no information was recorded. Existing interpolation methods introduce image artifacts and tend to be slow in filling large missing regions. Our goal was to develop a computationally efficient method that fills missing regions while adequately preserving image features. We use directional sticks to interpolate between pairs of known opposing voxels in nearby images. We tested our method on 30 volumetric ultrasound scans acquired from human subjects, and compared its performance to that of other published hole-filling methods. Reconstruction accuracy, fidelity, and time were improved compared with other methods.
Gao Y, Zhu L, Cates J, MacLeod RS, Bouix S, Tannenbaum A. A Kalman Filtering Perspective for Multiatlas Segmentation. SIAM J Imaging Sci. 2015;8 (2) :1007-29.Abstract
In multiatlas segmentation, one typically registers several atlases to the novel image, and their respective segmented label images are transformed and fused to form the final segmentation. In this work, we provide a new dynamical system perspective for multiatlas segmentation, inspired by the following fact: The transformation that aligns the current atlas to the novel image can be not only computed by direct registration but also inferred from the transformation that aligns the previous atlas to the image together with the transformation between the two atlases. This process is similar to the global positioning system on a vehicle, which gets position by inquiring from the satellite and by employing the previous location and velocity-neither answer in isolation being perfect. To solve this problem, a dynamical system scheme is crucial to combine the two pieces of information; for example, a Kalman filtering scheme is used. Accordingly, in this work, a Kalman multiatlas segmentation is proposed to stabilize the global/affine registration step. The contributions of this work are twofold. First, it provides a new dynamical systematic perspective for standard independent multiatlas registrations, and it is solved by Kalman filtering. Second, with very little extra computation, it can be combined with most existing multiatlas segmentation schemes for better registration/segmentation accuracy.
Wang D, Buckner RL, Fox MD, Holt DJ, Holmes AJ, Stoecklein S, Langs G, Pan R, Qian T, Li K, et al. Parcellating Cortical Functional Networks in Individuals. Nat Neurosci. 2015;18 (12) :1853-60.Abstract

The capacity to identify the unique functional architecture of an individual's brain is a crucial step toward personalized medicine and understanding the neural basis of variation in human cognition and behavior. Here we developed a cortical parcellation approach to accurately map functional organization at the individual level using resting-state functional magnetic resonance imaging (fMRI). A population-based functional atlas and a map of inter-individual variability were employed to guide the iterative search for functional networks in individual subjects. Functional networks mapped by this approach were highly reproducible within subjects and effectively captured the variability across subjects, including individual differences in brain lateralization. The algorithm performed well across different subject populations and data types, including task fMRI data. The approach was then validated by invasive cortical stimulation mapping in surgical patients, suggesting potential for use in clinical applications.

Ning L, Laun F, Gur Y, DiBella EVR, Deslauriers-Gauthier S, Megherbi T, Ghosh A, Zucchelli M, Menegaz G, Fick R, et al. Sparse Reconstruction Challenge for diffusion MRI: Validation on a physical phantom to determine which acquisition scheme and analysis method to use?. Med Image Anal. 2015;26 (1) :316-31.Abstract
Diffusion magnetic resonance imaging (dMRI) is the modality of choice for investigating in-vivo white matter connectivity and neural tissue architecture of the brain. The diffusion-weighted signal in dMRI reflects the diffusivity of water molecules in brain tissue and can be utilized to produce image-based biomarkers for clinical research. Due to the constraints on scanning time, a limited number of measurements can be acquired within a clinically feasible scan time. In order to reconstruct the dMRI signal from a discrete set of measurements, a large number of algorithms have been proposed in recent years in conjunction with varying sampling schemes, i.e., with varying b-values and gradient directions. Thus, it is imperative to compare the performance of these reconstruction methods on a single data set to provide appropriate guidelines to neuroscientists on making an informed decision while designing their acquisition protocols. For this purpose, the SPArse Reconstruction Challenge (SPARC) was held along with the workshop on Computational Diffusion MRI (at MICCAI 2014) to validate the performance of multiple reconstruction methods using data acquired from a physical phantom. A total of 16 reconstruction algorithms (9 teams) participated in this community challenge. The goal was to reconstruct single b-value and/or multiple b-value data from a sparse set of measurements. In particular, the aim was to determine an appropriate acquisition protocol (in terms of the number of measurements, b-values) and the analysis method to use for a neuroimaging study. The challenge did not delve on the accuracy of these methods in estimating model specific measures such as fractional anisotropy (FA) or mean diffusivity, but on the accuracy of these methods to fit the data. This paper presents several quantitative results pertaining to each reconstruction algorithm. The conclusions in this paper provide a valuable guideline for choosing a suitable algorithm and the corresponding data-sampling scheme for clinical neuroscience applications.
Zhu L, Kolesov I, Karasev P, Tannenbaum A. Interactive Image Segmentation Framework Based On Control Theory. Proc SPIE Int Soc Opt Eng. 2015;9413.Abstract
Segmentation of anatomical structures in medical imagery is a key step in a variety of clinical applications. Designing a generic, automated method that works for various structures and imaging modalities is a daunting task. Instead of proposing a new specific segmentation algorithm, in this paper, we present a general design principle on how to integrate user interactions from the perspective of control theory. In this formulation, Lyapunov stability analysis is employed to design and analyze an interactive segmentation system. The effectiveness and robustness of the proposed method are demonstrated.
Chen GH, Shah D, Golland P. A Latent Source Model for Patch-Based Image Segmentation. Med Image Comput Comput Assist Interv. 2015;9351 :140-8.Abstract

Despite the popularity and empirical success of patch-based nearest-neighbor and weighted majority voting approaches to medical image segmentation, there has been no theoretical development on when, why, and how well these nonparametric methods work. We bridge this gap by providing a theoretical performance guarantee for nearest-neighbor and weighted majority voting segmentation under a new probabilistic model for patch-based image segmentation. Our analysis relies on a new local property for how similar nearby patches are, and fuses existing lines of work on modeling natural imagery patches and theory for nonparametric classification. We use the model to derive a new patch-based segmentation algorithm that iterates between inferring local label patches and merging these local segmentations to produce a globally consistent image segmentation. Many existing patch-based algorithms arise as special cases of the new algorithm.

Liu S, Cai W, Liu S, Zhang F, Fulham M, Feng D, Pujol S, Kikinis R. Multimodal Neuroimaging Computing: A Review of the Applications in Neuropsychiatric Disorders. Brain Inform. 2015;2 :167-80.Abstract
Multimodal neuroimaging is increasingly used in neuroscience research, as it overcomes the limitations of individual modalities. One of the most important applications of multimodal neuroimaging is the provision of vital diagnostic data for neuropsychiatric disorders. Multimodal neuroimaging computing enables the visualization and quantitative analysis of the alterations in brain structure and function, and has reshaped how neuroscience research is carried out. Research in this area is growing exponentially, and so it is an appropriate time to review the current and future development of this emerging area. Hence, in this paper, we review the recent advances in multimodal neuroimaging (MRI, PET) and electrophysiological (EEG, MEG) technologies, and their applications to the neuropsychiatric disorders. We also outline some future directions for multimodal neuroimaging where researchers will design more advanced methods and models for neuropsychiatric research.
Liu S, Cai W, Liu S, Zhang F, Fulham M, Feng D, Pujol S, Kikinis R. Multimodal Neuroimaging Computing: The Workflows, Methods, and Platforms. Brain Inform. 2015;2 :181-95.Abstract
The last two decades have witnessed the explosive growth in the development and use of noninvasive neuroimaging technologies that advance the research on human brain under normal and pathological conditions. Multimodal neuroimaging has become a major driver of current neuroimaging research due to the recognition of the clinical benefits of multimodal data, and the better access to hybrid devices. Multimodal neuroimaging computing is very challenging, and requires sophisticated computing to address the variations in spatiotemporal resolution and merge the biophysical/biochemical information. We review the current workflows and methods for multimodal neuroimaging computing, and also demonstrate how to conduct research using the established neuroimaging computing packages and platforms.
Ratner V, Zhu L, Kolesov I, Nedergaard M, Benveniste H, Tannenbaum A. Optimal-mass-transfer-based estimation of glymphatic transport in living brain. Proc SPIE Int Soc Opt Eng. 2015;9413.Abstract
It was recently shown that the brain-wide cerebrospinal fluid (CSF) and interstitial fluid exchange system designated the 'glymphatic pathway' plays a key role in removing waste products from the brain, similarly to the lymphatic system in other body organs(1,2). It is therefore important to study the flow patterns of glymphatic transport through the live brain in order to better understand its functionality in normal and pathological states. Unlike blood, the CSF does not flow rapidly through a network of dedicated vessels, but rather through para-vascular channels and brain parenchyma in a slower time-domain, and thus conventional fMRI or other blood-flow sensitive MRI sequences do not provide much useful information about the desired flow patterns. We have accordingly analyzed a series of MRI images, taken at different times, of the brain of a live rat, which was injected with a paramagnetic tracer into the CSF via the lumbar intrathecal space of the spine. Our goal is twofold: (a) find glymphatic (tracer) flow directions in the live rodent brain; and (b) provide a model of a (healthy) brain that will allow the prediction of tracer concentrations given initial conditions. We model the liquid flow through the brain by the diffusion equation. We then use the Optimal Mass Transfer (OMT) approach(3) to derive the glymphatic flow vector field, and estimate the diffusion tensors by analyzing the (changes in the) flow. Simulations show that the resulting model successfully reproduces the dominant features of the experimental data.
Langs G, Golland P, Ghosh SS. Predicting Activation Across Individuals with Resting-State Functional Connectivity Based Multi-Atlas Label Fusion. Med Image Comput Comput Assist Interv. 2015;9350 :313-20.Abstract

The alignment of brain imaging data for functional neuroimaging studies is challenging due to the discrepancy between correspondence of morphology, and equivalence of functional role. In this paper we map functional activation areas across individuals by a multi-atlas label fusion algorithm in a functional space. We learn the manifold of resting-state fMRI signals in each individual, and perform manifold alignment in an embedding space. We then transfer activation predictions from a source population to a target subject via multi-atlas label fusion. The cost function is derived from the aligned manifolds, so that the resulting correspondences are derived based on the similarity of intrinsic connectivity architecture. Experiments show that the resulting label fusion predicts activation evoked by various experiment conditions with higher accuracy than relying on morphological alignment. Interestingly, the distribution of this gain is distributed heterogeneously across the cortex, and across tasks. This offers insights into the relationship between intrinsic connectivity, morphology and task activation. Practically, the mechanism can serve as prior, and provides an avenue to infer task-related activation in individuals for whom only resting data is available.

Pages