Slide 1
An Anatomically Curated Fiber Clustering White Matter Atlas for Consistent White Matter Tract Parcellation across the Lifespan
Slide 2
An Immersive Virtual Reality Environment for Diagnostic Imaging
Slide 3
Inter-site and Inter-scanner Diffusion MRI Data Harmonization
Slide 4
The Open Anatomy Browser: A Collaborative Web-Based Viewer for Interoperable Anatomy Atlases
Slide 5
Unsupervised Discovery of Emphysema Subtypes in a Large Clinical Cohort
Slide 6
Identifying Shared Brain Networks in Individuals by Decoupling Functional and Anatomical Variability
Slide 7
Supra-Threshold Fiber Cluster Statistics for Data-Driven Whole Brain Tractography Analysis
Slide 8
Free Water Modeling of Peritumoral Edema using Multi-fiber Tractography
Slide 8
Estimation of Bounded and Unbounded Trajectories in Diffusion MRI
Slide 9
Principal Gradient of Macroscale Cortical Organization
Slide 10
Evolution of a Simultaneous Segmentation and Atlas Registration
Slide 11
Multi-modality MRI-based Atlas of the Brain
Slide 12
Intracranial Fluid Redistribution
Slide 13
Corticospinal Tract Modeling for Neurosurgical Planning by Tracking through Regions of Peritumoral Edema and Crossing Fibers
Slide 14
Automated White Matter Fiber Tract Identification in Patients with Brain Tumors
Slide 15
State-space Models of Mental Processes from fMRI
Slide 16
Robust Initialization of Active Shape Models for Lung Segmentation in CT Scans: A Feature-Based Atlas Approach
Slide 17
Tractography-driven Groupwise Multi-Scale Parcellation of the Cortex
Slide 18
Gray Matter Alterations in Early Aging
Slide 19
Statistical Shape Analysis: From Landmarks to Diffeomorphisms
Slide 20
A Generative Probabilistic Model and Discriminative Extensions for Brain Lesion Segmentation
Slide 21
Joint Modeling of Imaging and Genetic Variability
Slide 22
MR-Ultrasound Fusion for Neurosurgery
Slide 23
Diffusion MRI and Tumor Heterogeneity
Slide 24
SlicerDMRI: Open Source Diffusion MRI Software for Brain Cancer Research

Neuroimage Analysis Center


The Neuroimaging Analysis Center is a research and technology center with the mission of advancing the role of neuroimaging in health care. The ability to access huge cohorts of patient medical records and radiology data, the emergence of ever-more detailed imaging modalities, and the availability of unprecedented computer processing power marks the possibility for a new era in neuroimaging, disease understanding, and patient treatment. We are excited to present a national resource center with the goal of finding new ways of extracting disease characteristics from advanced imaging and computation, and to make these methods available to the larger medical community through a proven methodology of world-class research, open-source software, and extensive collaboration.

Our Sponsor



The NAC is a Biomedical Technology Resource Center supported by the National Institute of Biomedical Imaging and Bioengineering (NIBIB) (P41 EB015902). It was supported by the National Center for Research Resources (NCRR) (P41 RR13218) through December 2011.

Contact the Center Directors



Carl-Fredrik Westin, PhD
Laboratory of Mathematics in Imaging
Brigham and Women's Hospital
1249 Boylston St., Room 240
Boston, MA 02215
Phone: +1 617 525-6209
E-mail: westin at

Ron Kikinis

Ron Kikinis, MD
Surgical Planning Laboratory 
Brigham and Women's Hospital 
75 Francis St, L1 Room 050
Boston, MA 02115
Phone: +1 617 732-7389
E-mail: kikinis at


Recent Publications

  • Ou W, Golland P, Hämäläinen M. Sources of variability in MEG. Med Image Comput Comput Assist Interv 2007;10(Pt 2):751-9.
    This paper investigates and characterizes sources of variability in MEG signals in multi-site, multi-subject studies. Understanding these sources will help to develop efficient strategies for comparing and pooling data across repetitions of an experiment, across subjects, and across sites. In this work, we investigated somatosensory MEG data collected at three different sites and applied variance component analysis and nonparametric KL divergence analysis in order to characterize the sources of variability. Our analysis showed that inter-subject differences are the biggest factor in the signal variability. We demonstrated that the timing of the deflections is very consistent in the early somatosensory response, which justifies a direct comparison of deflection peak times acquired from different visits, subjects, and systems. Compared with deflection peak times, deflection magnitudes have larger variation across sites; modeling of this variability is necessary for data pooling.
  • Bhalerao A, Westin C-F. Hyperspherical von Mises-Fisher mixture (HvMF) modelling of high angular resolution diffusion MRI. Med Image Comput Comput Assist Interv 2007;10(Pt 1):236-43.
    A mapping of unit vectors onto a 5D hypersphere is used to model and partition ODFs from HARDI data. This mapping has a number of useful and interesting properties and we make a link to interpretation of the second order spherical harmonic decompositions of HARDI data. The paper presents the working theory and experiments of using a von Mises-Fisher mixture model for directional samples. The MLE of the second moment of the HvMF pdf can also be related to fractional anisotropy. We perform error analysis of the estimation scheme in single and multi-fibre regions and then show how a penalised-likelihood model selection method can be employed to differentiate single and multiple fibre regions.
  • Kindlmann G, Estepar RSJ, Niethammer M, Haker S, Westin C-F. Geodesic-loxodromes for diffusion tensor interpolation and difference measurement. Med Image Comput Comput Assist Interv 2007;10(Pt 1):1-9.
    In algorithms for processing diffusion tensor images, two common ingredients are interpolating tensors, and measuring the distance between them. We propose a new class of interpolation paths for tensors, termed geodesic-loxodromes, which explicitly preserve clinically important tensor attributes, such as mean diffusivity or fractional anisotropy, while using basic differential geometry to interpolate tensor orientation. This contrasts with previous Riemannian and Log-Euclidean methods that preserve the determinant. Path integrals of tangents of geodesic-loxodromes generate novel measures of over-all difference between two tensors, and of difference in shape and in orientation.
  • Golland P, Golland Y, Malach R. Detection of spatial activation patterns as unsupervised segmentation of fMRI data. Med Image Comput Comput Assist Interv 2007;10(Pt 1):110-8.
    In functional connectivity analysis, networks of interest are defined based on correlation with the mean time course of a user-selected ’seed’ region. In this work we propose to simultaneously estimate the optimal representative time courses that summarize the fMRI data well and the partition of the volume into a set of disjoint regions that are best explained by these representative time courses. Our approach offers two advantages. First, is removes the sensitivity of the analysis to the details of the seed selection. Second, it substantially simplifies group analysis by eliminating the need for a subject-specific threshold at which correlation values are deemed significant. This unsupervised technique generalizes connectivity analysis to situations where candidate seeds are difficult to identify reliably or are unknown. Our experimental results indicate that the functional segmentation provides a robust, anatomically meaningful and consistent model for functional connectivity in fMRI.
  • Yeo BTT, Sabuncu MR, Desikan R, Fischl B, Golland P. Effects of registration regularization and atlas sharpness on segmentation accuracy. Med Image Comput Comput Assist Interv 2007;10(Pt 1):683-91.
    In this paper, we propose a unified framework for computing atlases from manually labeled data at various degrees of "sharpness" and the joint registration-segmentation of a new brain with these atlases. In non-rigid registration, the tradeoff between warp regularization and image fidelity is typically set empirically. In segmentation, this leads to a probabilistic atlas of arbitrary "sharpness": weak regularization results in well-aligned training images and a "sharp" atlas; strong regularization yields a "blurry" atlas. We study the effects of this tradeoff in the context of cortical surface parcellation by comparing three special cases of our framework, namely: progressive registration-segmentation of a new brain to increasingly "sharp" atlases with increasingly flexible warps; secondly, progressive registration to a single atlas with increasingly flexible warps; and thirdly, registration to a single atlas with fixed constrained warps. The optimal parcellation in all three cases corresponds to a unique balance of atlas "sharpness" and warp regularization that yield statistically significant improvements over the previously demonstrated parcellation results.