基于偏微分方程的图像平滑方法比较研究.docVIP

基于偏微分方程的图像平滑方法比较研究 　　摘要:近年来,基于偏微分方程的图像平滑技术受到越来越多的关注。该文详细介绍了几种基于偏微分方程的图像平滑方法,进而根据理论分析和实验结果对几种方法进行比较。 　　关键词:偏微分方程;图像平滑;比较研究 　　中图分类号:TP3文献标识码:A文章编号:1009-3044(2009)14-3804-02 　　The Comparison Research of Image Smoothing Method Based on PDE 　　WANG Zhu-xia1,2, ZANG Shun-quan1,2 　　(1.Xian Institute of Post and Telecommunications, Xian 710121, China; 2.Xidian University; Xian 710071, China) 　　Abstract: Recently, the image smoothing method based on PDE attract attention of more and more. In this paper,several image smoothing models based on PDE are introduced in detail. Image smoothing models arecompared according to theory analysis and experiment results. 　　Key words: PDE; image smoothing; comparison research 　　1 引言 　　由于获取和传输过程往往会对图像引入噪声,使得图像的视觉效果受到严重的影响,对图像后继处理和分析带来很大不便。 图像平滑作为图像的预处理过程,平滑质量的好坏直接影响到后继处理和分析的结果。传统的图像平滑算法有均值滤波、中值滤波和高斯滤波等,这些方法在去噪的同时会破坏图像的重要特征,如边缘、纹理等。 　　基于偏微分方程(PDE)的图像处理方法已经在图像处理领域得到了广泛应用。偏微分方程方法与通常的图像处理算法相比,虽然计算量较大,耗时较长,但由于其灵活的拓扑学结构、广阔的应用领域逐渐受到人们的重视。 　　图像平滑作为图像处理重要的环节,基于偏微分方程的图像平滑方法得到了迅速的发展。本文在介绍经典的基于偏微分方程的图像平滑方法基础上,详细分析各种模型的优缺点,以便为图像平滑方法的改进提供一些启示。 　　2 整体变分(TV )平滑模型[1-2] 　　在基于偏微分方程方法的平滑模型中,整体变分(TV,Total Variation)模型表现出色。设Ω是Rn中的有界开子集,μ为局部可积函数,其整体变分定义为 　　(1) 　　有界变差函数空间可定义为 　　(2) 　　设μ0为观测到的图像,于是μ0=μ+n, 则: 　　定义图像的能量函数: 　　其中λ为拉格朗日乘数。 　　由变分法知识,能量泛函取得极小值得必要条件是满足Euler-Lagrange方程,即 　　(4) 　　由梯度下降法,则得到TV 平滑模型: 　　(5) 　　对|?μ|进行正则化,可得到正则化模型: 　　(6) 　　由有限差分法得到离散格式: 　　(7) 　　3 四阶偏微分方程平滑模型 　　TV模型的偏微分方程是一个二阶微分方程,为了解决二阶偏微分方程处理时对图像分块产生的阶梯效应,可以采用四阶偏微分方程。 　　定义图像的正则能量为: 　　 (8) 　　由变分法知识,得欧拉公式为: 　　(9) 　　则梯度下降模型为: 　　 (10) 　　定义g(s)=f(s)/s并规定: 　　则(10)可化为: 　　 (11) 　　其离散过程为: 　　1)计算u的二阶差分: 　　(12) 　　2)计算函数g: 　　(13) 　　3)计算函数g的二阶差分: 　　(14) 　　4)离散格式为: 　　(15) 　　4 基于各向异性扩散方程的图像平滑 　　4.1 P-M模型 　　为了达到去噪同时保护边缘,可以采用扩散过程中的传导系数依赖于图像的局部特征。具体来说,在图像比较平坦的区域,传导系数能自动增大。这可使平坦区域中较小的起伏被平滑;而在图像的边缘附近,传导系数能自动减小,可使边缘几乎不受影响。Perona和Malik于1990年在热传导方程的基础上得到了以下的各向异性扩散模型: 　　(16) 　　其中,g(|?μ|)规定了扩散程度。 　　若g(0)=1当s-∞时g(s) -0。在平滑过程中,在边缘处|?μ|很大,g(s)很小,扩散程度变小;在平坦区域|?μ|很小,g(s)很大