Multi-domain, higher order level set scheme for 3D image segmentation on the GPU

Level set method based segmentation provides an efficient tool for topological and geometrical shape handling. Conventional level set surfaces are only C⁰ continuous since the level set evolution involves linear interpolation to compute derivatives. Bajaj et al. present a higher order method to evaluate level set surfaces that are C² continuous, but are slow due to high computational burden. In this paper, we provide a higher order GPU based solver for fast and efficient segmentation of large volumetric images. We also extend the higher order method to multi-domain segmentation. Our streaming solver is efficient in memory usage.