Anatomy of large biological specimens is often reconstructed from serially sectioned volumes imaged by high-resolution microscopy. We developed a method to reassemble a continuous volume from such large section series that explicitly minimizes artificial deformation by applying a global elastic constraint. We demonstrate our method on a series of transmission electron microscopy sections covering the entire 558-cell Caenorhabditis elegans embryo and a segment of the Drosophila melanogaster larval ventral nerve cord. DOI: https://doi.org/10. 1038/nmeth.2072 Posted at the Zurich Open Repository and Archive, University of Zurich ZORA URL: https://doi.org/10.5167/uzh-127852 Accepted Version Originally published at: Saalfeld, Stephan; Fetter, Richard; Cardona, Albert; Tomancak, Pavel (2012). Elastic volume reconstruction from series of ultra-thin microscopy sections. Nature Methods, 9(7):717-720. DOI: https://doi.org/10.1038/nmeth.2072
Elastic Volume Reconstruction from Series of Ultra-thin Microscopy SectionsStephan The downside of the method is that by physically cutting a block into sections the continuity between sections is lost and individual sections are deformed. To recover the imaged volume and extract biologically interesting information such as the reconstruction of neuronal circuits 2,3,5 , sections need to be aligned and distortion must be removed.Alignment can be achieved by maximizing the overlap of similar image content between adjacent sections. However, there are two unknowns that change image content across the section series: specimen's shape and independent section distortion introduced during preparation. Naively warping one section into another would compensate for the shape of the specimen and by that introduce artificial deformation. Our method exploits the fact that the biological specimen's shape typically changes smoothly across sections whereas the independent distortion in each section is random and uncorrelated with neighboring sections. We align all sections not only to their direct neighbors in the series but to all sections in a local neighborhood by modelling sections as two-dimensional elastic sheets that penalize non-rigid deformation ( Fig. 1a and Supplementary Fig. 1).All sections are treated as moving targets in a template-free global alignment. The elastic constraint is implemented as a spring-connected particle system where each section is represented as a triangular spring-mesh ( Fig. 1b and Supplementary Video 1 and Online methods).For each vertex of the spring-mesh, we search for the corresponding location in other sections by pairwise block-matching using normalized cross-correlation (NCC). To that end, we explore all translation vectors in an immediate neighborhood which requires sections to be in approximate alignment (Fig. 1c,d). We estimate this approximate alignment using automatically extracted landmark correspondences from invariant local image features as described previously 7,8 . Originally proposed for robust object recognition under partial occlus...