Fast and Efficient Numerical Finite Difference Method for Multiphase Image Segmentation

We present a simple numerical solution algorithm for a gradient flow for the Modica-Mortola functional and numerically investigate its dynamics. The proposed numerical algorithm involves both the operator splitting and the explicit Euler methods. A time step formula is derived from the stability analysis, and the goodness of fit of transition width is tested. We perform various numerical experiments to investigate the property of the gradient flow equation, to verify the characteristics of our method in the image segmentation application, and to analyze the effect of parameters. In particular, we propose an initialization process based on target objects. Furthermore, we conduct comparison tests in order to check the performance of our proposed method.