RESUMEN
We present a novel semi-automatic method for segmenting neural processes in large, highly anisotropic EM (electron microscopy) image stacks. Our method takes advantage of sparse scribble annotations provided by the user to guide a 3D variational segmentation model, thereby allowing our method to globally optimally enforce 3D geometric constraints on the segmentation. Moreover, we leverage a novel algorithm for propagating segmentation constraints through the image stack via optimal volumetric pathways, thereby allowing our method to compute highly accurate 3D segmentations from very sparse user input. We evaluate our method by reconstructing 16 neural processes in a 1024 x 1024 x 50 nanometer-scale EM image stack of a mouse hippocampus. We demonstrate that, on average, our method is 68% more accurate than previous state-of-the-art semi-automatic methods.
Asunto(s)
Hipocampo/metabolismo , Procesamiento de Imagen Asistido por Computador/métodos , Imagenología Tridimensional/métodos , Microscopía Electrónica/métodos , Neuronas/patología , Algoritmos , Animales , Anisotropía , Automatización , Humanos , Cadenas de Markov , Ratones , Modelos Neurológicos , Programas InformáticosRESUMEN
OBJECTIVE: Subsampling of radially encoded MRI acquisitions in combination with sparsity promoting methods opened a door to significantly increased imaging speed, which is crucial for many important clinical applications. In particular, it has been shown recently that total variation (TV) regularization efficiently reduces undersampling artifacts. The drawback of the method is the long reconstruction time which makes it impossible to use in daily clinical practice, especially if the TV optimization problem has to be solved repeatedly to select a proper regularization parameter. MATERIALS AND METHODS: The goal of this work was to show that for the case of MR Angiography, TV filtering can be performed as a post-processing step, in contrast to the common approach of integrating TV penalties in the image reconstruction process. With this approach, it is possible to use TV algorithms with data fidelity terms in image space, which can be implemented very efficiently on graphic processing units (GPUs). The combination of a special radial sampling trajectory and a full 3D formulation of the TV minimization problem is crucial for the effectiveness of the artifact elimination process. RESULTS AND CONCLUSION: The computation times of GPU-TV show that interactive elimination of undersampling artifacts is possible even for large volume data sets, in particular allowing the interactive determination of the regularization parameter. Results from phantom measurements and in vivo angiography data sets show that 3D TV, together with the proposed sampling trajectory, leads to pronounced improvements in image quality. However, while artifact removal was very efficient for angiography data sets in this work, it cannot be expected that the proposed method of TV post-processing will work for arbitrary types of scans.