Reconstruction of neuronal morphologies from electron microscopy images using graph cuts

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 χ² 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.

	noise σ = 0.0447	noise σ = 0.0632	noise σ = 0.0894
Dataset 1	0.9342	0.9315	0.9297
Dataset 2	0.8948	0.8903	0.8868
Dataset 3	0.9192	0.9112	0.9040