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

  • Kindlmann G, Tricoche X, Westin CF. Delineating white matter structure in diffusion tensor MRI with anisotropy creases. Med Image Anal. 2007;11(5):492–502.
    Geometric models of white matter architecture play an increasing role in neuroscientific applications of diffusion tensor imaging, and the most popular method for building them is fiber tractography. For some analysis tasks, however, a compelling alternative may be found in the first and second derivatives of diffusion anisotropy. We extend to tensor fields the notion from classical computer vision of ridges and valleys, and define anisotropy creases as features of locally extremal tensor anisotropy. Mathematically, these are the loci where the gradient of anisotropy is orthogonal to one or more eigenvectors of its Hessian. We propose that anisotropy creases provide a basis for extracting a skeleton of the major white matter pathways, in that ridges of anisotropy coincide with interiors of fiber tracts, and valleys of anisotropy coincide with the interfaces between adjacent but distinctly oriented tracts. The crease extraction algorithm we present generates high-quality polygonal models of crease surfaces, which are further simplified by connected-component analysis. We demonstrate anisotropy creases on measured diffusion MRI data, and visualize them in combination with tractography to confirm their anatomic relevance.
  • Yu P, Yeo BTT, Grant E, Fischl B, Golland P. Cortical Folding Development Study based on Over-Complete Spherical Wavelets. Proc IEEE Int Conf Comput Vis. 2007;2007.
    We introduce the use of over-complete spherical wavelets for shape analysis of 2D closed surfaces. Bi-orthogonal spherical wavelets have been shown to be powerful tools in the segmentation and shape analysis of 2D closed surfaces, but unfortunately they suffer from aliasing problems and are therefore not invariant under rotations of the underlying surface parameterization. In this paper, we demonstrate the theoretical advantage of over-complete wavelets over bi-orthogonal wavelets and illustrate their utility on both synthetic and real data. In particular, we show that over-complete spherical wavelets allow us to build more stable cortical folding development models, and detect a wider array of regions of folding development in a newborn dataset.
  • DiMaio S, Kapur T, Cleary K, Aylward S, Kazanzides P, Vosburgh K, Ellis R, Duncan J, Farahani K, Lemke H, Peters T, Lorensen WB, Gobbi D, Haller J, Clarke LL, Pizer S, Taylor R, Galloway R, Fichtinger G, Hata N, Lawson K, Tempany CM, Kikinis R, Jolesz FA. Challenges in Image-guided Therapy System Design. Neuroimage. 2007;37 Suppl 1:144–51.
    System development for image-guided therapy (IGT), or image-guided interventions (IGI), continues to be an area of active interest across academic and industry groups. This is an emerging field that is growing rapidly: major academic institutions and medical device manufacturers have produced IGT technologies that are in routine clinical use, dozens of high-impact publications are published in well regarded journals each year, and several small companies have successfully commercialized sophisticated IGT systems. In meetings between IGT investigators over the last two years, a consensus has emerged that several key areas must be addressed collaboratively by the community to reach the next level of impact and efficiency in IGT research and development to improve patient care. These meetings culminated in a two-day workshop that brought together several academic and industrial leaders in the field today. The goals of the workshop were to identify gaps in the engineering infrastructure available to IGT researchers, develop the role of research funding agencies and the recently established US-based National Center for Image Guided Therapy (NCIGT), and ultimately to facilitate the transfer of technology among research centers that are sponsored by the National Institutes of Health (NIH). Workshop discussions spanned many of the current challenges in the development and deployment of new IGT systems. Key challenges were identified in a number of areas, including: validation standards; workflows, use-cases, and application requirements; component reusability; and device interface standards. This report elaborates on these key points and proposes research challenges that are to be addressed by a joint effort between academic, industry, and NIH participants.
  • Csapo I, Holland CM, Guttmann CRG. Image registration framework for large-scale longitudinal MRI data sets: strategy and validation. Magn Reson Imaging. 2007;25(6):889–93.
    Advanced magnetic resonance imaging (MRI) studies often require the transformation of large numbers of images into a common space. Calculating transformations that relate each image to every other and applying them to the images on demand are theoretically possible; however, these can be computationally prohibitive. Therefore, relating each image to only one other image, then linking those transforms together to relate any two images in the database, may be an efficient alternative. Evaluated were the feasibility and validity of image registration to bring intraindividual MR images into mutual correspondence for longitudinal analysis through the concatenation of precomputed transforms. A longitudinal data set of 10 multiple sclerosis patients with nine serial dual-echo spin-echo, 1.5-T MRI scans was used. Intrasubject registrations were performed stepwise between consecutive images and direct from each time point to the baseline. Consecutive transforms were concatenated and evaluated against direct registrations by comparing the resulting transformed images (using Pearson correlation coefficient). Confounding variables such as time between scans, brain atrophy, and change in lesion load were evaluated. We found the images resampled with the direct and the concatenated transforms to be highly correlated, and there was no significant difference between methods. Differences in brain parenchymal fraction (a measure of brain atrophy) showed significant inverse correlation with the correspondence of the resampled images. Results indicate that concatenating multiple transforms that link two images together produces near-identical results to that of direct registration; thus, this method is both useful and valid.
  • This article discusses and reviews advanced forms of serial morphometry in the context of a disease progression model in multiple sclerosis (MS). This model of disease activity distinguishes between overall disease activity and the proportion thereof that becomes permanent damage. This translates into a progression model that features a repair potential, which, when exhausted, marks the conversion or progression from relapsing to progressive disease. The level of repair capacity at a given time determines the rate of progression. Both clinical and MRI variables appear to be in support of such a model. We examine possible MRI markers for this repair capacity, particularly the short-term behavior of new MRI lesions, quantified by methods of time-series analysis—that is, capturing lesion dynamics in the form of MRI intensity change directly, rather than shape or volume change. Lower rates of individual lesion recovery may represent lower repair and greater proximity to a progressive stage. Individuals with low transient lesion turnover appear to undergo more rapid progression and atrophy. Because disease-modifying therapies aim to alter the pathophysiological chain of inflammation, demyelination, and axonal loss, a therapeutic effect may therefore be more readily apparent as a change in lesion dynamics and recovery rate and level, rather than a change in total lesion burden or enhancing lesion number.