Markus Nilsson, Samo Lasič, Ivana Drobnjak, Daniel Topgaard, and Carl-Fredrik Westin. 7/2017. “Resolution Limit of Cylinder Diameter Estimation by Diffusion MRI: The Impact of Gradient Waveform and Orientation Dispersion.” NMR Biomed, 30, 7.Abstract
Diffusion MRI has been proposed as a non-invasive technique for axonal diameter mapping. However, accurate estimation of small diameters requires strong gradients, which is a challenge for the transition of the technique from preclinical to clinical MRI scanners, since these have weaker gradients. In this work, we develop a framework to estimate the lower bound for accurate diameter estimation, which we refer to as the resolution limit. We analyse only the contribution from the intra-axonal space and assume that axons can be represented by impermeable cylinders. To address the growing interest in using techniques for diffusion encoding that go beyond the conventional single diffusion encoding (SDE) sequence, we present a generalised analysis capable of predicting the resolution limit regardless of the gradient waveform. Using this framework, waveforms were optimised to minimise the resolution limit. The results show that, for parallel cylinders, the SDE experiment is optimal in terms of yielding the lowest possible resolution limit. In the presence of orientation dispersion, diffusion encoding sequences with square-wave oscillating gradients were optimal. The resolution limit for standard clinical MRI scanners (maximum gradient strength 60-80 mT/m) was found to be between 4 and 8 μm, depending on the noise levels and the level of orientation dispersion. For scanners with a maximum gradient strength of 300 mT/m, the limit was reduced to between 2 and 5 μm.
M Zhang, R Liao, Adrian V Dalca, E Turk, J Luo, E Grant, and Polina Golland. 6/2017. “Frequency Diffeomorphisms for Efficient Image Registration.” Inf Process Med Imaging., 10265, Pp. 559-70.
Adrian V Dalca, K. L. Bouman, William T. Freeman, Natalia S Rost, Mert R Sabuncu, and Polina Golland. 6/2017. “Population Based Image Imputation.” Inf Process Med Imaging., 10265, 659-71.
Lena Maier-Hein, Swaroop Vedula, Stefanis Speidel, Nassir Navab, Ron Kikinis, Adrian Park, Matthias Eisenman, Hubertus Feussner, and Germain Forestier. 6/2017. “Surgical Data Science: Enabling Next-generation Surgery.” Nature Biomedical Engineering. Maier-Hein-NBE2017.pdf
Stephen SF Yip, Chintan Parmar, Daniel Blezek, Raul San Jose Estepar, Steve Pieper, John Kim, and Hugo JWL Aerts. 6/2017. “Application of the 3D Slicer Chest Imaging Platform Segmentation Algorithm for Large Lung Nodule Delineation.” PLoS One, 12, 6, Pp. e0178944.Abstract
PURPOSE: Accurate segmentation of lung nodules is crucial in the development of imaging biomarkers for predicting malignancy of the nodules. Manual segmentation is time consuming and affected by inter-observer variability. We evaluated the robustness and accuracy of a publically available semiautomatic segmentation algorithm that is implemented in the 3D Slicer Chest Imaging Platform (CIP) and compared it with the performance of manual segmentation. METHODS: CT images of 354 manually segmented nodules were downloaded from the LIDC database. Four radiologists performed the manual segmentation and assessed various nodule characteristics. The semiautomatic CIP segmentation was initialized using the centroid of the manual segmentations, thereby generating four contours for each nodule. The robustness of both segmentation methods was assessed using the region of uncertainty (δ) and Dice similarity index (DSI). The robustness of the segmentation methods was compared using the Wilcoxon-signed rank test (pWilcoxon<0.05). The Dice similarity index (DSIAgree) between the manual and CIP segmentations was computed to estimate the accuracy of the semiautomatic contours. RESULTS: The median computational time of the CIP segmentation was 10 s. The median CIP and manually segmented volumes were 477 ml and 309 ml, respectively. CIP segmentations were significantly more robust than manual segmentations (median δCIP = 14ml, median dsiCIP = 99% vs. median δmanual = 222ml, median dsimanual = 82%) with pWilcoxon~10-16. The agreement between CIP and manual segmentations had a median DSIAgree of 60%. While 13% (47/354) of the nodules did not require any manual adjustment, minor to substantial manual adjustments were needed for 87% (305/354) of the nodules. CIP segmentations were observed to perform poorly (median DSIAgree≈50%) for non-/sub-solid nodules with subtle appearances and poorly defined boundaries. CONCLUSION: Semi-automatic CIP segmentation can potentially reduce the physician workload for 13% of nodules owing to its computational efficiency and superior stability compared to manual segmentation. Although manual adjustment is needed for many cases, CIP segmentation provides a preliminary contour for physicians as a starting point.
Ariel D Stock, Sivan Gelb, Ofer Pasternak, Ayal Ben-Zvi, and Chaim Putterman. 6/2017. “The Blood Brain Barrier and Neuropsychiatric Lupus: New Perspectives in Light of Advances in Understanding the Neuroimmune Interface.” Autoimmun Rev, 16, 6, Pp. 612-9.Abstract
Experts have previously postulated a linkage between lupus associated vascular pathology and abnormal brain barriers in the immunopathogenesis of neuropsychiatric lupus. Nevertheless, there are some discrepancies between the experimental evidence, or its interpretation, and the working hypotheses prevalent in this field; specifically, that a primary contributor to neuropsychiatric disease in lupus is permeabilization of the blood brain barrier. In this commonly held view, any contribution of the other known brain barriers, including the blood-cerebrospinal fluid and meningeal barriers, is mostly excluded from the discussion. In this review we will shed light on some of the blood brain barrier hypotheses and try to trace their roots. In addition, we will suggest new research directions to allow for confirmation of alternative interpretations of the experimental evidence linking the pathology of intra-cerebral vasculature to the pathogenesis of neuropsychiatric lupus.
Christian Wachinger, Matthew Brennan, Greg C Sharp, and Polina Golland. 6/2017. “Efficient Descriptor-Based Segmentation of Parotid Glands With Nonlocal Means.” IEEE Trans Biomed Eng, 64, 7, Pp. 1492-1502.Abstract
OBJECTIVE: We introduce descriptor-based segmentation that extends existing patch-based methods by combining intensities, features, and location information. Since it is unclear which image features are best suited for patch selection, we perform a broad empirical study on a multitude of different features. METHODS: We extend nonlocal means segmentation by including image features and location information. We search larger windows with an efficient nearest neighbor search based on kd-trees. We compare a large number of image features. RESULTS: The best results were obtained for entropy image features, which have not yet been used for patch-based segmentation. We further show that searching larger image regions with an approximate nearest neighbor search and location information yields a significant improvement over the bounded nearest neighbor search traditionally employed in patch-based segmentation methods. CONCLUSION: Features and location information significantly increase the segmentation accuracy. The best features highlight boundaries in the image. SIGNIFICANCE: Our detailed analysis of several aspects of nonlocal means-based segmentation yields new insights about patch and neighborhood sizes together with the inclusion of location information. The presented approach advances the state-of-the-art in the segmentation of parotid glands for radiation therapy planning.
Miaomiao Zhang, Ruizhi Liao, Adrian V Dalca, Esra A Turk, Jie Luo, Ellen P Grant, and Polina Golland. 6/2017. “Frequency Diffeomorphisms for Efficient Image Registration.” Inf Process Med Imaging, 10265, Pp. 559-570.Abstract
This paper presents an efficient algorithm for large deformation diffeomorphic metric mapping (LDDMM) with geodesic shooting for image registration. We introduce a novel finite dimensional Fourier representation of diffeomorphic deformations based on the key fact that the high frequency components of a diffeomorphism remain stationary throughout the integration process when computing the deformation associated with smooth velocity fields. We show that manipulating high dimensional diffeomorphisms can be carried out entirely in the bandlimited space by integrating the nonstationary low frequency components of the displacement field. This insight substantially reduces the computational cost of the registration problem. Experimental results show that our method is significantly faster than the state-of-the-art diffeomorphic image registration methods while producing equally accurate alignment. We demonstrate our algorithm in two different applications of image registration: neuroimaging and in-utero imaging.
Vincent Koppelmans, Ofer Pasternak, Jacob J Bloomberg, Yiri De E Dios, Scott J Wood, Roy Riascos, Patricia A Reuter-Lorenz, Igor S Kofman, Ajitkumar P Mulavara, and Rachael D Seidler. 6/2017. “Intracranial Fluid Redistribution but No White Matter Microstructural Changes During a Spaceflight Analog.” Sci Rep, 7, 1, Pp. 3154.Abstract
The neural correlates of spaceflight-induced sensorimotor impairments are unknown. Head down-tilt bed rest (HDBR) serves as a microgravity analog because it mimics the headward fluid shift and axial body unloading of spaceflight. We investigated focal brain white matter (WM) changes and fluid shifts during 70 days of 6° HDBR in 16 subjects who were assessed pre (2x), during (3x), and post-HDBR (2x). Changes over time were compared to those in control subjects (n = 12) assessed four times over 90 days. Diffusion MRI was used to assess WM microstructure and fluid shifts. Free-Water Imaging was used to quantify distribution of intracranial extracellular free water (FW). Additionally, we tested whether WM and FW changes correlated with changes in functional mobility and balance measures. HDBR resulted in FW increases in fronto-temporal regions and decreases in posterior-parietal regions that largely recovered by two weeks post-HDBR. WM microstructure was unaffected by HDBR. FW decreases in the post-central gyrus and precuneus correlated negatively with balance changes. We previously reported that gray matter increases in these regions were associated with less HDBR-induced balance impairment, suggesting adaptive structural neuroplasticity. Future studies are warranted to determine causality and underlying mechanisms.
Jenna Schabdach, William M Wells, Michael Cho, and Kayhan N Batmanghelich. 6/2017. “A Likelihood-Free Approach for Characterizing Heterogeneous Diseases in Large-Scale Studies.” Inf Process Med Imaging, 10265, Pp. 170-183.Abstract
We propose a non-parametric approach for characterizing heterogeneous diseases in large-scale studies. We target diseases where multiple types of pathology present simultaneously in each subject and a more severe disease manifests as a higher level of tissue destruction. For each subject, we model theof local image descriptors as samples generated by an unknown subject-specific probability density. Instead of approximating the probability density via a parametric family, we propose to side step the parametric inference by directly estimating the divergence between subject densities. Our method maps the collection of local image descriptors to a signaturethat is used to predict a clinical measurement. We are able to interpret the prediction of the clinical variable in the population and individual levels by carefully studying the divergences. We illustrate an application this method on simulated data as well as on a large-scale lung CT study of Chronic Obstructive Pulmonary Disease (COPD). Our approach outperforms classical methods on both simulated and COPD data and demonstrates the state-of-the-art prediction on an important physiologic measure of airflow (the forced respiratory volume in one second, FEV1).
Ruizhi Liao, Lipeng Ning, Zhenrui Chen, Laura Rigolo, Shun Gong, Ofer Pasternak, Alexandra J Golby, Yogesh Rathi, and Lauren J O'Donnell. 6/2017. “Performance of Unscented Kalman Filter Tractography in Edema: Analysis of the Two-tensor Model.” Neuroimage Clin, 15, Pp. 819-31.Abstract
Diffusion MRI tractography is increasingly used in pre-operative neurosurgical planning to visualize critical fiber tracts. However, a major challenge for conventional tractography, especially in patients with brain tumors, is tracing fiber tracts that are affected by vasogenic edema, which increases water content in the tissue and lowers diffusion anisotropy. One strategy for improving fiber tracking is to use a tractography method that is more sensitive than the traditional single-tensor streamline tractography. We performed experiments to assess the performance of two-tensor unscented Kalman filter (UKF) tractography in edema. UKF tractography fits a diffusion model to the data during fiber tracking, taking advantage of prior information from the previous step along the fiber. We studied UKF performance in a synthetic diffusion MRI digital phantom with simulated edema and in retrospective data from two neurosurgical patients with edema affecting the arcuate fasciculus and corticospinal tracts. We compared the performance of several tractography methods including traditional streamline, UKF single-tensor, and UKF two-tensor. To provide practical guidance on how the UKF method could be employed, we evaluated the impact of using various seed regions both inside and outside the edematous regions, as well as the impact of parameter settings on the tractography sensitivity. We quantified the sensitivity of different methods by measuring the percentage of the patient-specific fMRI activation that was reached by the tractography. We expected that diffusion anisotropy threshold parameters, as well as the inclusion of a free water model, would significantly influence the reconstruction of edematous WM fiber tracts, because edema increases water content in the tissue and lowers anisotropy. Contrary to our initial expectations, varying the fractional anisotropy threshold and including a free water model did not affect the UKF two-tensor tractography output appreciably in these two patient datasets. The most effective parameter for increasing tracking sensitivity was the generalized anisotropy (GA) threshold, which increased the length of tracked fibers when reduced to 0.075. In addition, the most effective seeding strategy was seeding in the whole brain or in a large region outside of the edema. Overall, the main contribution of this study is to provide insight into how UKF tractography can work, using a two-tensor model, to begin to address the challenge of fiber tract reconstruction in edematous regions near brain tumors.
Adrian V Dalca, Katherine L Bouman, William T. Freeman, Natalia S Rost, Mert R Sabuncu, and Polina Golland. 6/2017. “Population Based Image Imputation.” Inf Process Med Imaging, 10265, Pp. 659-671.Abstract
We present an algorithm for creating high resolution anatomically plausible images consistent with acquired clinical brain MRI scans with large inter-slice spacing. Although large databases of clinical images contain a wealth of information, medical acquisition constraints result in sparse scans that miss much of the anatomy. These characteristics often render computational analysis impractical as standard processing algorithms tend to fail when applied to such images. Highly specialized or application-specific algorithms that explicitly handle sparse slice spacing do not generalize well across problem domains. In contrast, our goal is to enable application of existing algorithms that were originally developed for high resolution research scans to significantly undersampled scans. We introduce a model that captures fine-scale anatomical similarity across subjects in clinical image collections and use it to fill in the missing data in scans with large slice spacing. Our experimental results demonstrate that the proposed method outperforms current upsampling methods and promises to facilitate subsequent analysis not previously possible with scans of this quality.
Yangming Ou, Lilla Zöllei, Kallirroi Retzepi, Victor Castro, Sara V Bates, Steve Pieper, Katherine P Andriole, Shawn N Murphy, Randy L Gollub, and Patricia Ellen Grant. 6/2017. “Using Clinically Acquired MRI to Construct Age-specific Atlases: Quantifying Spatiotemporal ADC Changes from Birth to 6-year Old.” Hum Brain Mapp, 38, 6, Pp. 3052-68.Abstract
Diffusion imaging is critical for detecting acute brain injury. However, normal apparent diffusion coefficient (ADC) maps change rapidly in early childhood, making abnormality detection difficult. In this article, we explored clinical PACS and electronic healthcare records (EHR) to create age-specific ADC atlases for clinical radiology reference. Using the EHR and three rounds of multiexpert reviews, we found ADC maps from 201 children 0-6 years of age scanned between 2006 and 2013 who had brain MRIs with no reported abnormalities and normal clinical evaluations 2+ years later. These images were grouped in 10 age bins, densely sampling the first 1 year of life (5 bins, including neonates and 4 quarters) and representing the 1-6 year age range (an age bin per year). Unbiased group-wise registration was used to construct ADC atlases for 10 age bins. We used the atlases to quantify (a) cross-sectional normative ADC variations; (b) spatiotemporal heterogeneous ADC changes; and (c) spatiotemporal heterogeneous volumetric changes. The quantified age-specific whole-brain and region-wise ADC values were compared to those from age-matched individual subjects in our study and in multiple existing independent studies. The significance of this study is that we have shown that clinically acquired images can be used to construct normative age-specific atlases. These first of their kind age-specific normative ADC atlases quantitatively characterize changes of myelination-related water diffusion in the first 6 years of life. The quantified voxel-wise spatiotemporal ADC variations provide standard references to assist radiologists toward more objective interpretation of abnormalities in clinical images. Our atlases are available at
Walid I Essayed, Fan Zhang, Prashin Unadkat, Rees G Cosgrove, Alexandra J Golby, and Lauren J O'Donnell. 6/2017. “White Matter Tractography for Neurosurgical Planning: A Topography-based Review of the Current State of the Art.” Neuroimage Clin, 15, Pp. 659-72.Abstract
We perform a review of the literature in the field of white matter tractography for neurosurgical planning, focusing on those works where tractography was correlated with clinical information such as patient outcome, clinical functional testing, or electro-cortical stimulation. We organize the review by anatomical location in the brain and by surgical procedure, including both supratentorial and infratentorial pathologies, and excluding spinal cord applications. Where possible, we discuss implications of tractography for clinical care, as well as clinically relevant technical considerations regarding the tractography methods. We find that tractography is a valuable tool in variable situations in modern neurosurgery. Our survey of recent reports demonstrates multiple potentially successful applications of white matter tractography in neurosurgery, with progress towards overcoming clinical challenges of standardization and interpretation.
Vadim Ratner, Yi Gao, Hedok Lee, Rena Elkin, Maiken Nedergaard, Helene Benveniste, and Allen Tannenbaum. 5/2017. “Cerebrospinal and Interstitial Fluid Transport via the Glymphatic Pathway Modeled by Optimal Mass Transport.” Neuroimage, 152, Pp. 530-7.Abstract
The glymphatic pathway is a system which facilitates continuous cerebrospinal fluid (CSF) and interstitial fluid (ISF) exchange and plays a key role in removing waste products from the rodent brain. Dysfunction of the glymphatic pathway may be implicated in the pathophysiology of Alzheimer's disease. Intriguingly, the glymphatic system is most active during deep wave sleep general anesthesia. By using paramagnetic tracers administered into CSF of rodents, we previously showed the utility of MRI in characterizing a macroscopic whole brain view of glymphatic transport but we have yet to define and visualize the specific flow patterns. Here we have applied an alternative mathematical analysis approach to a dynamic time series of MRI images acquired every 4min over ∼3h in anesthetized rats, following administration of a small molecular weight paramagnetic tracer into the CSF reservoir of the cisterna magna. We use Optimal Mass Transport (OMT) to model the glymphatic flow vector field, and then analyze the flow to find the network of CSF-ISF flow channels. We use 3D visualization computational tools to visualize the OMT defined network of CSF-ISF flow channels in relation to anatomical and vascular key landmarks from the live rodent brain. The resulting OMT model of the glymphatic transport network agrees largely with the current understanding of the glymphatic transport patterns defined by dynamic contrast-enhanced MRI revealing key CSF transport pathways along the ventral surface of the brain with a trajectory towards the pineal gland, cerebellum, hypothalamus and olfactory bulb. In addition, the OMT analysis also revealed some interesting previously unnoticed behaviors regarding CSF transport involving parenchymal streamlines moving from ventral reservoirs towards the surface of the brain, olfactory bulb and large central veins.
Jorge L Bernal-Rusiel, Nicolas Rannou, Randy L Gollub, Steve Pieper, Shawn Murphy, Richard Robertson, Patricia E Grant, and Rudolph Pienaar. 5/2017. “Reusable Client-Side JavaScript Modules for Immersive Web-Based Real-Time Collaborative Neuroimage Visualization.” Front Neuroinform, 11, Pp. 32.Abstract
In this paper we present a web-based software solution to the problem of implementing real-time collaborative neuroimage visualization. In both clinical and research settings, simple and powerful access to imaging technologies across multiple devices is becoming increasingly useful. Prior technical solutions have used a server-side rendering and push-to-client model wherein only the server has the full image dataset. We propose a rich client solution in which each client has all the data and uses the Google Drive Realtime API for state synchronization. We have developed a small set of reusable client-side object-oriented JavaScript modules that make use of the XTK toolkit, a popular open-source JavaScript library also developed by our team, for the in-browser rendering and visualization of brain image volumes. Efficient realtime communication among the remote instances is achieved by using just a small JSON object, comprising a representation of the XTK image renderers' state, as the Google Drive Realtime collaborative data model. The developed open-source JavaScript modules have already been instantiated in a web-app called MedView, a distributed collaborative neuroimage visualization application that is delivered to the users over the web without requiring the installation of any extra software or browser plugin. This responsive application allows multiple physically distant physicians or researchers to cooperate in real time to reach a diagnosis or scientific conclusion. It also serves as a proof of concept for the capabilities of the presented technological solution.
Inbar Seroussi, Denis S Grebenkov, Ofer Pasternak, and Nir Sochen. 4/2017. “Microscopic Interpretation and Generalization of the Bloch-Torrey Equation for Diffusion Magnetic Resonance.” J Magn Reson, 277, Pp. 95-103.Abstract
In order to bridge microscopic molecular motion with macroscopic diffusion MR signal in complex structures, we propose a general stochastic model for molecular motion in a magnetic field. The Fokker-Planck equation of this model governs the probability density function describing the diffusion-magnetization propagator. From the propagator we derive a generalized version of the Bloch-Torrey equation and the relation to the random phase approach. This derivation does not require assumptions such as a spatially constant diffusion coefficient, or ad hoc selection of a propagator. In particular, the boundary conditions that implicitly incorporate the microstructure into the diffusion MR signal can now be included explicitly through a spatially varying diffusion coefficient. While our generalization is reduced to the conventional Bloch-Torrey equation for piecewise constant diffusion coefficients, it also predicts scenarios in which an additional term to the equation is required to fully describe the MR signal.
Chantal MW Tax, Carl-Fredrik Westin, Tom Dela Haije, Andrea Fuster, Max A Viergever, Evan Calabrese, Luc Florack, and Alexander Leemans. 4/2017. “Quantifying the Brain's Sheet Structure with Normalized Convolution.” Med Image Anal, 39, Pp. 162-77.Abstract
The hypothesis that brain pathways form 2D sheet-like structures layered in 3D as "pages of a book" has been a topic of debate in the recent literature. This hypothesis was mainly supported by a qualitative evaluation of "path neighborhoods" reconstructed with diffusion MRI (dMRI) tractography. Notwithstanding the potentially important implications of the sheet structure hypothesis for our understanding of brain structure and development, it is still considered controversial by many for lack of quantitative analysis. A means to quantify sheet structure is therefore necessary to reliably investigate its occurrence in the brain. Previous work has proposed the Lie bracket as a quantitative indicator of sheet structure, which could be computed by reconstructing path neighborhoods from the peak orientations of dMRI orientation density functions. Robust estimation of the Lie bracket, however, is challenging due to high noise levels and missing peak orientations. We propose a novel method to estimate the Lie bracket that does not involve the reconstruction of path neighborhoods with tractography. This method requires the computation of derivatives of the fiber peak orientations, for which we adopt an approach called normalized convolution. With simulations and experimental data we show that the new approach is more robust with respect to missing peaks and noise. We also demonstrate that the method is able to quantify to what extent sheet structure is supported for dMRI data of different species, acquired with different scanners, diffusion weightings, dMRI sampling schemes, and spatial resolutions. The proposed method can also be used with directional data derived from other techniques than dMRI, which will facilitate further validation of the existence of sheet structure.
Tina Kapur and Clare M. Tempany. 3/2017. “Proceedings of the 9th Image Guided Therapy Workshop” 9, Pp. 1-54. 2017 IGT Workshop Proceedings
Xiaojun Chen, Lu Xu, Huixiang Wang, Fang Wang, Qiugen Wang, and Ron Kikinis. 3/2017. “Development of a Surgical Navigation System Based on 3D Slicer for Intraoperative Implant Placement Surgery.” Med Eng Phys, 41, Pp. 81-9.Abstract

Implant placement has been widely used in various kinds of surgery. However, accurate intraoperative drilling performance is essential to avoid injury to adjacent structures. Although some commercially-available surgical navigation systems have been approved for clinical applications, these systems are expensive and the source code is not available to researchers. 3D Slicer is a free, open source software platform for the research community of computer-aided surgery. In this study, a loadable module based on Slicer has been developed and validated to support surgical navigation. This research module allows reliable calibration of the surgical drill, point-based registration and surface matching registration, so that the position and orientation of the surgical drill can be tracked and displayed on the computer screen in real time, aiming at reducing risks. In accuracy verification experiments, the mean target registration error (TRE) for point-based and surface-based registration were 0.31±0.06mm and 1.01±0.06mm respectively, which should meet clinical requirements. Both phantom and cadaver experiments demonstrated the feasibility of our surgical navigation software module.