GPU-Accelerated RSF Level Set Evolution for Large-Scale Microvascular Segmentation

Meher Niger,Helya Goharbavang,Taeyong Ahn,Emily K. Alley,Joshua D. Wythe,Guoning Chen,David Mayerich
2024-04-03
Abstract:Microvascular networks are challenging to model because these structures are currently near the diffraction limit for most advanced three-dimensional imaging modalities, including confocal and light sheet microscopy. This makes semantic segmentation difficult, because individual components of these networks fluctuate within the confines of individual pixels. Level set methods are ideally suited to solve this problem by providing surface and topological constraints on the resulting model, however these active contour techniques are extremely time intensive and impractical for terabyte-scale images. We propose a reformulation and implementation of the region-scalable fitting (RSF) level set model that makes it amenable to three-dimensional evaluation using both single-instruction multiple data (SIMD) and single-program multiple-data (SPMD) parallel processing. This enables evaluation of the level set equation on independent regions of the data set using graphics processing units (GPUs), making large-scale segmentation of high-resolution networks practical and inexpensive.
Image and Video Processing,Computer Vision and Pattern Recognition
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the large - scale segmentation of microvascular networks. Specifically: 1. **Background Challenges**: - The structure of the microvascular network is very complex, and currently, the resolution of most advanced three - dimensional imaging techniques (such as confocal microscopy and light - sheet microscopy) is close to the diffraction limit, which makes semantic segmentation difficult. - Although the level - set method is suitable for solving this problem because it can provide surface and topological constraints, these active contour techniques are extremely time - consuming and impractical when dealing with large - scale images. 2. **Research Objectives**: - Propose a reformulation and implementation of the Region - Scalable Fitting (RSF) level - set model to make it suitable for Single Instruction Multiple Data (SIMD) and Single Program Multiple Data (SPMD) parallel processing. - Independently evaluate different regions of the data set by the Graphics Processing Unit (GPU) to make the large - scale segmentation of high - resolution networks practical and economical. 3. **Specific Methods**: - **Theoretical Basis**: Based on the RSF model, propose a parallel implementation method suitable for three - dimensional images. - **Parallel Computing**: Utilize the parallel computing power of the GPU to efficiently segment large - scale images. - **Initialization Method**: Propose a seed - surface method to initialize the level - set model. - **Verification Method**: Use the Monte Carlo method for verification, compare the results with other segmentation methods, and show the processing speed improvement brought by parallel hardware through strict performance analysis. 4. **Application and Testing**: - This method has been tested on multiple data sets, which use state - of - the - art imaging techniques (such as micro - CT, light - sheet fluorescence microscopy, and milling microscopy) to obtain microvascular data. - Verify the performance and accuracy of the RSF model through the Monte Carlo method, and provide a detailed performance analysis to show the advantages of parallel hardware. In conclusion, this paper aims to solve the computational efficiency and memory requirements in large - scale microvascular network segmentation by improving and parallelizing the RSF level - set model, thereby providing a more powerful tool for medical and biological research.