水平集图像处理入门.doc

水平集图像处理入门 (A TUTORIAL OF LEVEL SET FOR IMAGE PROCESSING) 冯向军(FENG XIANGJUN) 05/29/2006 [内容提要]:本文简明扼要地阐明了水平集技术的基本概念、基本思想、基本方法和基本技术。通过展示水平集去噪的机理和实效，使读者不难对水平集技术实用于图像处理有个真切的体会。本文可供图像处理专业的大学生、研究生参考。附录给出了水平集技术的关键细节：曲率和梯度计算的一个MATLAB程序。 Abstract ---In this report, the basic concepts, basic thought, basic methodology and basic technology of level set technology for image processing were briefly introduced. Undergraduate students and graduate students related with image processing may take this report as a reference. 一、水平集的定义 [1]与实数c对应的可微函数的水平集是实点集 {(x1, x2, ...,xn) | f(x1, x2,...,xn) = c} 称可微函数f为水平集函数。 [举例] 函数 对应于常数c的水平集是 以(0，0，0)为球心，sqrt(c) 为半径的球面。 当 n=2, 称水平集为水平曲线(LEVEL CURVE)。 当 n=3, 称水平集为水平曲面(LEVEL SURFACE)。 二、水平集图像处理的核心思想 水平集图像处理的核心思想是把n维描述视为高一维(n+1)维的水平集,或者说是把n维描述视为有n维变量的水平集函数f的水平集.这样一来就把求解n维描述的演化过程转化为求解关于有n维变量的水平集函数f的演化所导致的水平集的演化过程。其要害是通过这种转化，引入了变中的相对不变：水平集函数f的水平c不变。我们把这种变中的相对不变叫做泛对称。引入了泛对称，就引入了规律，而引入了规律就能推演出水平集在此规律下依各种具体条件而演化的具体演化方程。也即是说，引进了泛对称这一规律，我们就有了从一般到特殊的演绎过程的出发点和依据。这种思想方法的实质是以关系来决定对象。 三、水平集图像处理的奠基性工作 水平集图像处理的奠基性工作包含在下列论文之中 Osher Sethian (Fronts propogating with curvature-dependent speed:Algorithms based on Hamilton-Jacobi formulations, Journal of Computational Physics 79, 12-49, 1988) 以下是开创人之一STANLEY OSHER的样子 [4]四、水平集图像处理的基本方程 考虑零水平集x(t)所对应的水平集函数 ，则有  (4-1) 对方程(4-1)两边求关于时间的偏导数,有  (4-2) 假设F为外法向方向的速度，那么 这其中 因此，我们便得到基本方程式  (4-3) 除此基本方程式以外，还有其它一些方程式，它们都能从方程(4-1)推导出来。 基本方程式式(4-3)是水平集函数及相应的水平集在法向力F的推动下的演化方程。 国内外有很多人介绍水平集图像处理的思想时一般都太偏向技术层面，没有看到这思想的实质：以极为普遍的关系:泛对称来决定对象。通过上述简单推导，读者诸君不难看出这水平集图像处理的精神实质。 五、一般性算法 水平集图像处理的这种思想方法直接形成了一种一般性算法 (1)设定水平集函数的初态； (2)确定动力F的形式； (3)按基本方程推演水平集函数的各状态； (4)对于每一水平集函数的状态求解零水平集。 尽管具体技术细节千变万化，水平集图像处理的算法万变不离其宗，都是围绕上述 一般性算法而展开的。  六、水平集技术去噪的机理[2] 所谓水平集技术，就是利用水平集的概念和将求解n维描述的演化转化为求解由于具有n个实变量的水平集函数的演化而导致的水平集的演化这么一种思想方法来达到某种实用目的的技术。水平集技术可用于图像去噪。其机理如下 图6.1 简单曲线 图6.2 不太简单的曲线 图6.3 不简单的曲线 图6.1-6.3给出了曲线的三种类型：简单曲线、不太简单的曲线，不简单的曲线。 我们所感兴趣的是简单曲线。通俗地说，所谓简单曲线即是无交叉点的曲线。 简单曲线在曲率力的