Median filter method for mean curvature flow using a random Jacobi algorithm

Anton Ullrich,Tim Laux
2024-10-10
Abstract:We present an efficient scheme for level set mean curvature flow using a domain discretization and median filters. For this scheme, we show convergence in $L^\infty$-norm under mild assumptions on the number of points in the discretization. In addition, we strengthen the weak convergence result for the MBO thresholding scheme applied to data clustering of Lelmi and one of the authors. This is done through a strong convergence of the discretized heat flow in the optimal regime. Different boundary conditions are also discussed.
Analysis of PDEs,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop an efficient algorithm to approximate the level set mean curvature flow and prove the convergence of this algorithm. Specifically, the paper introduces an algorithm based on the median filter. This algorithm can operate on all level sets in the discrete domain simultaneously and can be easily generalized to high - dimensional spaces. The core contributions of the paper are as follows: 1. **Algorithm Design**: An evolutionary algorithm based on the median filter is proposed. In each time step, this algorithm updates the value of each point by calculating the median in the local neighborhood. This algorithm is essentially equivalent to the application of the MBO (Merriman - Bence - Osher) scheme on any level set. 2. **Convergence Proof**: It is proved that, under the assumption that there are enough points in the discretized domain, the algorithm converges almost surely in the strong \( L^\infty \)-topology. Specifically, Theorem 3.14 elaborates on this point. 3. **Heat Flow Convergence**: The strong convergence of the discretized heat flow related to the MBO scheme is further proved, thereby strengthening the previous results on the weak convergence of the MBO scheme. This is achieved through the \( L^2 \)-convergence of the mean scheme of the fully discretized heat equation. See Theorem 4.2 and Theorem 4.3. 4. **Γ - Convergence of Energy**: The Γ - convergence of the energy associated with the evolution under the non - uniform Young angle condition is proved. The main innovation of the paper lies in applying the median filter to approximate the level set mean curvature flow, thereby providing a simple and efficient numerical method, especially suitable for problems in high - dimensional spaces. In addition, the paper also discusses the application of different boundary conditions, such as Neumann conditions and Dirichlet conditions, and proposes corresponding improvement schemes. These results are not only of great theoretical significance but also provide strong support for practical applications.