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

  • Limperopoulos C, Soul JS, Haidar H, Hüppi PS, Bassan H, Warfield SK, Robertson RL, Moore M, Akins P, Volpe JJ, Plessis A\ e J du. Impaired trophic interactions between the cerebellum and the cerebrum among preterm infants. Pediatrics 2005;116(4):844-50.
    BACKGROUND: Advanced neuroimaging techniques have brought increasing recognition of cerebellar injury among premature infants. The developmental relationship between early brain injury and effects on the cerebrum and cerebellum remains unclear. OBJECTIVES: To examine whether cerebral parenchymal brain lesions among preterm infants are associated with subsequent decreases in cerebellar volume and, conversely, whether primary cerebellar injury is associated with decreased cerebral brain volumes, with advanced, 3-dimensional, volumetric MRI at term gestational age equivalent. METHODS: Total cerebellar volumes and cerebellar gray and myelinated white matter volumes were determined through manual outlining for 74 preterm infants with unilateral periventricular hemorrhagic infarction (14 infants), bilateral diffuse periventricular leukomalacia (20 infants), cerebellar hemorrhage (10 infants), or normal term gestational age equivalent MRI findings (30 infants). Total brain and right/left cerebral and cerebellar hemispheric volumes were calculated. RESULTS: Unilateral cerebral brain injury was associated with significantly decreased volume of the contralateral cerebellar hemisphere. Conversely, unilateral primary cerebellar injury was associated with a contralateral decrease in supratentorial brain volume. Cerebellar gray matter and myelinated white matter volumes were reduced significantly not only among preterm infants with primary cerebellar hemorrhage but also among infants with cerebral parenchymal brain injury. CONCLUSIONS: These data suggest strongly that both reduction in contralateral cerebellar volume with unilateral cerebral parenchymal injury and reduction in total cerebellar volume with bilateral cerebral lesions are related to trophic transsynaptic effects. Early-life cerebellar injury may contribute importantly to the high rates of cognitive, behavioral, and motor deficits reported for premature infants.
  • Tsai A, Wells WM, Warfield SK, Willsky AS. An EM Algorithm for Shape Classification Based on Level Sets. Med Image Anal 2005;9(5):491-502.
    In this paper, we propose an expectation-maximization (EM) approach to separate a shape database into different shape classes, while simultaneously estimating the shape contours that best exemplify each of the different shape classes. We begin our formulation by employing the level set function as the shape descriptor. Next, for each shape class we assume that there exists an unknown underlying level set function whose zero level set describes the contour that best represents the shapes within that shape class. The level set function for each example shape in the database is modeled as a noisy measurement of the appropriate shape class’s unknown underlying level set function. Based on this measurement model and the judicious introduction of the class labels as the hidden data, our EM formulation calculates the labels for shape classification and estimates the shape contours that best typify the different shape classes. This resulting iterative algorithm is computationally efficient, simple, and accurate. We demonstrate the utility and performance of this algorithm by applying it to two medical applications.
  • Kubicki M, Park H, Westin CF, Nestor PG, Mulkern RV, Maier SE, Niznikiewicz M, Connor E, Levitt JJ, Frumin M, Kikinis R, Jolesz FA, McCarley RW, Shenton ME. DTI and MTR Abnormalities in Schizophrenia: Analysis of White Matter Integrity. Neuroimage 2005;26(4):1109-18.
    Diffusion tensor imaging (DTI) studies in schizophrenia demonstrate lower anisotropic diffusion within white matter due either to loss of coherence of white matter fiber tracts, to changes in the number and/or density of interconnecting fiber tracts, or to changes in myelination, although methodology as well as localization of such changes differ between studies. The aim of this study is to localize and to specify further DTI abnormalities in schizophrenia by combining DTI with magnetization transfer imaging (MTI), a technique sensitive to myelin and axonal alterations in order to increase specificity of DTI findings. 21 chronic schizophrenics and 26 controls were scanned using Line-Scan-Diffusion-Imaging and T1-weighted techniques with and without a saturation pulse (MT). Diffusion information was used to normalize co-registered maps of fractional anisotropy (FA) and magnetization transfer ratio (MTR) to a study-specific template, using the multi-channel daemon algorithm, designed specifically to deal with multidirectional tensor information. Diffusion anisotropy was decreased in schizophrenia in the following brain regions: the fornix, the corpus callosum, bilaterally in the cingulum bundle, bilaterally in the superior occipito-frontal fasciculus, bilaterally in the internal capsule, in the right inferior occipito-frontal fasciculus and the left arcuate fasciculus. MTR maps demonstrated changes in the corpus callosum, fornix, right internal capsule, and the superior occipito-frontal fasciculus bilaterally; however, no changes were noted in the anterior cingulum bundle, the left internal capsule, the arcuate fasciculus, or inferior occipito-frontal fasciculus. In addition, the right posterior cingulum bundle showed MTR but not FA changes in schizophrenia. These findings suggest that, while some of the diffusion abnormalities in schizophrenia are likely due to abnormal coherence, or organization of the fiber tracts, some of these abnormalities may, in fact, be attributed to or coincide with myelin/axonal disruption.
  • Aische A du B d, De Craene M, Geets X, Gregoire V, Macq B, Warfield SK. Efficient multi-modal dense field non-rigid registration: alignment of histological and section images. Med Image Anal 2005;9(6):538-46.
    We describe a new algorithm for non-rigid registration capable of estimating a constrained dense displacement field from multi-modal image data. We applied this algorithm to capture non-rigid deformation between digital images of histological slides and digital flat-bed scanned images of cryotomed sections of the larynx, and carried out validation experiments to measure the effectiveness of the algorithm. The implementation was carried out by extending the open-source Insight ToolKit software. In diagnostic imaging of cancer of the larynx, imaging modalities sensitive to both anatomy (such as MRI and CT) and function (PET) are valuable. However, these modalities differ in their capability to discriminate the margins of tumor. Gold standard tumor margins can be obtained from histological images from cryotomed sections of the larynx. Unfortunately, the process of freezing, fixation, cryotoming and staining the tissue to create histological images introduces non-rigid deformations and significant contrast changes. We demonstrate that the non-rigid registration algorithm we present is able to capture these deformations and the algorithm allows us to align histological images with scanned images of the larynx. Our non-rigid registration algorithm constructs a deformation field to warp one image onto another. The algorithm measures image similarity using a mutual information similarity criterion, and avoids spurious deformations due to noise by constraining the estimated deformation field with a linear elastic regularization term. The finite element method is used to represent the deformation field, and our implementation enables us to assign inhomogeneous material characteristics so that hard regions resist internal deformation whereas soft regions are more pliant. A gradient descent optimization strategy is used and this has enabled rapid and accurate convergence to the desired estimate of the deformation field. A further acceleration in speed without cost of accuracy is achieved by using an adaptive mesh refinement strategy.
  • Brehmer A, Lindig TM, Schrödl F, Neuhuber W, Ditterich D, Rexer M, Rupprecht H. Morphology of enkephalin-immunoreactive myenteric neurons in the human gut. Histochem Cell Biol 2005;123(2):131-8.
    The aim of this study was the morphological and further chemical characterisation of neurons immunoreactive for leu-enkephalin (leuENK). Ten wholemounts of small and large intestinal segments from nine patients were immunohistochemically triple-stained for leuENK/neurofilament 200 (NF)/substance P (SP). Based on their simultaneous NF-reactivity and 3D reconstruction of single NF-reactive cells, 97.5% of leuENK-positive neurons displayed the appearance of stubby neurons: small somata; short, stubby dendrites and one axon. Of these leuENK-reactive stubby neurons, 91.3% did not display co-reactivity for SP whereas 8.7% were SP-co-reactive. As to their axonal projection pattern, 50.4% of the recorded leuENK stubby neurons had axons running orally whereas in 29.4% they ran anally; the directions of the remaining 20.2% could not be determined. No axons were seen to enter into secondary strands of the myenteric plexus. Somal area measurements revealed clearly smaller somata of leuENK-reactive stubby neurons (between 259+/-47 microm(2) and 487+/-113 microm(2)) than those of putative sensory type II neurons (between 700+/-217 microm(2) and 1,164+/-396 microm(2)). The ratio dendritic field area per somal area of leuENK-reactive stubby neurons was between 2.0 and 2.8 reflecting their short dendrites. Additionally, we estimated the proportion of leuENK-positive neurons in comparison to the putative whole myenteric neuron population in four leuENK/anti-Hu doublestained wholemounts. This proportion ranged between 5.9% and 8.3%. We suggest leuENK-reactive stubby neurons to be muscle motor neurons and/or ascending interneurons. Furthermore, we explain why we do not use the term "Dogiel type I neurons" for this population.