Volume 11 Supplement 1

Nineteenth Annual Computational Neuroscience Meeting: CNS*2010

Open Access

Reconstruction of neuronal morphologies from electron microscopy images using graph cuts

BMC Neuroscience201011(Suppl 1):P142


Published: 20 July 2010

Understanding neural connectivity and structures in the brain requires detailed 3D anatomical models, and such an understanding is essential to the study of the nervous system. However, the reconstruction of 3D geometry from a large set of dense nanoscale medical images is very challenging, because of the imperfections in staining and noise in the imaging process. In this abstract, we present an automated tracking and segmentation framework to extract 2D contours and to trace them through the z direction. The image segmentation problem is posed as Maximum A Posteriori (MAP) estimation of a Markov Random Field (MRF). The MAP-MRF formulation is equivalent to minimizing the corresponding Gibbs energy function and graph cuts [1] are used to obtain the optimal solution to the energy function. The energy function carries both the information from the current slice and from the previous slice. Furthermore, to improve segmentation accuracy, interactive segmentation editing tools that allow users to refine an incorrect segmentation are developed. Segmentation evaluation metrics using χ2 distance and image intensities are applied to compute confidence values of segmentations. Confidence values are helpful to point users to segmentations with low confidence values, thus reducing the time of manual inspection. Experiments were conducted on three synthetic datasets (Dataset 3 is generated by using NeuGen [3]) and on nanoscale image sequences from the Serial Block Face Scanning Electron Microscope (SBF-SEM) [2]. Table 1 shows the quantitative Dice similarity measure on the three synthetic datasets. The average DS over all experiments is 0.9113, which shows that the proposed tracking and segmentation algorithm is robust under noisy conditions. Figure 1 presents the reconstructed 3D structures of tracked and segmented 2D contours, where part of a neuron can be seen in figure 1(c) and figure 1(d). The results indicate that our method can successfully track and segment densely packed cells in EM image stacks. We expect these reconstructed geometries to allow high-fidelity computational simulation of various kinds of neurons.
Table 1

Quantitative Dice similarity measure of the three synthetic datasets with different levels of Gaussian noise.


noise σ = 0.0447

noise σ = 0.0632

noise σ = 0.0894

Dataset 1




Dataset 2




Dataset 3




Figure 1

Reconstruction results of synthetic datasets and of an SBF-SEM image stack.



NIH/NINDS #R01-NS54252 (and supplement 03S1), NSF #0079874 (MRI), NSF #0905041 (CRCNS), Texas ARP #000512-0146-2001.

Authors’ Affiliations

Department of Computer Science and Engineering, Texas A&M University


  1. Boykov Y, Jolly M: Interactive graph cuts for optimal boundary and region segmentation of objects in N-D images. ICCV. 2001, 105-112.Google Scholar
  2. Denk W, Horstmann H: Serial block-face scanning electron microscopy to reconstruct three-dimensional tissue nanostructure. PloS Biol. 2004, 2: e329-10.1371/journal.pbio.0020329.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Eberhard JP, Wanner A, Wittum G: Neugen: A tool for the generation of realistic morphology of cortical neurons and neural networks in 3d. Neurocomputing. 2006, 327-342. 10.1016/j.neucom.2006.01.028.Google Scholar


© Yang and Choe; licensee BioMed Central Ltd. 2010

This article is published under license to BioMed Central Ltd.