【active contours without edges】level set水平集方法实现图像分割 论文详解+代码实现

1. 传统依赖边界的active contour 算法

传统的levelset方法依赖于图像边界的梯度，需要引入一个单调递减函数edge-detector $g(z)$定义levelset方法的迭代速度。
$$st:\underset{z\to \infty }{\lim }g(z)=0$$

例如：$g(\nabla u_0(x,y))=\frac{1}{1+|\nabla G_{\sigma }\ast u_0(x,y){|}^p}$
$u_0(x,y)=z$
原图像看作曲面$\Sigma (x,y,u_0(x,y))=0$

level set方法定义了一系列等高线
$\varphi (x,y)=c_i$
我们不妨用曲线$C$:$\varphi (x,y)=0$作为初始曲线，在根据需要对进行迭代
此时曲线$C_{t_i}$ :$\varphi (t_i,x,y)=0$在迭代的过程中不断向图像的边界逼近，最后稳定下来即得到边界。

Remark： $\varphi (t,x,y)$其中$\varphi (0,x,y)={\varphi }_0(x,y)$，根据定义的$\frac{\partial \varphi }{\partial t}=V|\nabla \varphi |$利用迭代法求解PDE

level set方法有以下结论：
单位切向量$T$,$n=-\frac{\nabla \varphi }{|\nabla \varphi |}$ , $\varphi (x,y)=c$确定曲线$C(x,y)$

1. $n\perp T$
   证明：由IFT（隐函数定理）
   $\varphi (x,y)$