Abstract
Brain surface profile is essential for brain studies, including registration, segmentation of brain structure and drawing neuronal circuits. Recent advances in high-throughput imaging techniques enable imaging whole mouse brain at micron spatial resolution and provide a basis for more fine quantitative studies in neuroscience. However, reconstructing micron resolution brain surface from newly produced neuronal dataset still faces challenges. Most current methods apply global analysis, which are neither applicable to a large imaging dataset nor to a brain surface with an inhomogeneous signal intensity. Here, we proposed a resampling-based variational model for this purpose. In this model, the movement directions of the initial boundary elements are fixed, the final positions of the initial boundary elements that form the brain surface are determined by the local signal intensity. These features assure an effective reconstruction of the brain surface from a new brain dataset. Compared with conventional typical methods, such as level set based method and active contour method, our method significantly increases the recall and precision rates above 97% and is approximately hundreds-fold faster. We demonstrated a fast reconstruction at micron level of the whole brain surface from a large dataset of hundreds of GB in size within 6 hours.