图像平滑作为图像预处理基本步骤之一，为后继图像处理带来很大方便，常见图像平滑方法是线性高斯滤波。实质上是一种信号的滤波器，其用途是信号的平滑处理，我们知道数字图像用于后期应用，其噪声是的问题，由于误差会累计传递等原因，其实编程运算的话就是一个模板运算，拿图像的八连通区域来说，中间点的像素值就等于八连通区的像素值的均值，这样达到平滑的效果，若使用理想滤波器，会在图像中产生振铃现象。KOENDERINK指出图像I0（x,y）与不同尺度的高斯核卷积所得到的平滑图像等价于传导系数为常数的热扩散方程的解，此解属于各向同性扩散，在平滑同时造成图像特征的模糊化。PERONA等提出了如下各向异性扩散模型：

　　就被认为该点处于区域内部并对其进行平滑，否则认为该点处于边缘对其进行增强。由于P-M各向异性滤波具有选择性平滑机制，在平滑图像区域内部同时又保持原始图像边缘结构，因此达到了优于线性滤波器的效果。滤波器阈值参数决定滤波器的性能，该值起到一种边缘检测作用。“阶梯”效应一直是困扰着P-M方法的严重问题，这是由它的偏微分方程所决定的，是所有低阶非线性扩散方法所固有的。

　　基于上述问题，在图像平滑时提出了基于自适应的高阶偏微分方程图像平滑方法，避免了传统P-M方法图像平滑方法的缺陷，获取比较好的视觉效果。

　　1 高阶方法的提出

　　P-M方法处理结果是分段恒定的，容易导致结果图像“阶梯”状分布，视觉效果不理想。而高阶方法的处理结果是分段线性的，在视觉感知上明显优于P-M方法。

　　2 自适应阈值的参数估计

　　参数k的选取是非线性扩散方法的一个主要难题。如何确定扩散的范围及扩散的程度，使降噪和强化顺利进行的同时，图像信息又不致因过度平滑而大量损失，是应用时需要解决的关键问题。传统确定参数的方法是人为指定一个固定常数。一般来说，不同图像需要设置不同参数值，后来也有学者提出自动估计梯度阈值的方法，如提出自动估计梯度阈值的公式Sapiro：

　　由于图像不同区域边缘强度分布不一致，噪声也不同，而且不同尺度空间边缘强度和噪声也不一样，因此对整幅图像使用同一个全局固定的参数是不合适的。本文提出了自适应阈值参数选择方法，使得该算法阈值参数可以完全自动确定，真正实现自适应阈值参数非线性滤波。滤波是将信号中特定波段频率滤除的操作，是抑制和防止干扰的一项重要措施。根据观察某一随机过程的结果，对另一与之有关的随机过程进行估计的概率理论与方法。滤波一词起源于通信理论，它是从含有干扰的接收信号中提取有用信号的一种技术。“接收信号”相当于被观测的随机过程，“有用信号”相当于被估计的随机过程。历史上早考虑的是维纳滤波，后来R.E.卡尔曼和R.S.布西于20世纪60年代提出了卡尔曼滤波。现对一般的非线性滤波问题的研究相当活跃。非线性滤波可归结为求条件期望 的问题。对于有限多个观测值的情形，条件期望原则上可以用贝叶斯公式来计算。但即使在比较简单的场合，这样得出的结果也是相当繁杂的，无论对实际应用或理论研究都很不方便。

　　上述分析完全是一种理想化假设，对于实际图像处理存在如下两大困难：

　　（1）由于邻域内所取样本有限，从样本直方图不能直接得到完全规则分布的曲线，很难正确地选取阈值参数k（x）值。

　　（2）现实图像中，即使在同一区域内，平均灰度值往往是变化的，而不是严格满足分段常数模型。

　　为克服以上问题，采取Mean Shift核平滑方法对算子直方图进行处理，使其尽可能准确地刻画多峰分布情况。

　　直方图平滑效果示例如图2所示。

　　为了降低计算复杂度， 在估计参数k（x）时，每隔r/2距离计算（r为邻域半径），然后对每个像素点选择近邻的值近似为该点阈值参数，并且每迭代之后更新参数k（x）值，也可设置每隔若干迭代次数更新参数k（x）值，这样可保证一定的前提下，大大提高效率。由于基于直方图算法复杂度较低，阈值参数估计计算相对于图像平滑本身计算开销较小。

　　3 实验结果及分析

　　为了比较自适应参数的高阶方法和P-M方法平滑的性能，进行不同实验并加以分析。

　　（1）选用若干幅标准图像进行测试，图3显示没有加入噪声的Grid图像平滑效果，从图中可看出当迭代次数比较多时，P-M方法平滑的部分边缘明显被模糊了。

　　（2）图4是对加入噪声的图像去噪的实验结果，图5给出了对图4各边缘的提取结果，其中图5（b）采用经典的Canny边缘提取算子。anny边缘检测器是使用函数edge的有效边缘检测器。该方法总结如下：1、图像使用带有指定标准偏差σ的高斯滤波器来平滑，从而可以减少噪声。2、在每一点处计算局部梯度g（x,y）=[G2x+G2y]1/2 和边缘方向α（x,y）=arctan（Gy/Gx）。边缘点定义为梯度方向上其强度局部的点。脊像素使用两个阈值T1和T2做阈值处理，其中T1<T2。值大于T2的脊像素称为强边缘像素，T1和T2之间的脊像素称为弱边缘像素。从图4和图5中可看出，P-M方法抑制孤立噪声点和保持图像边缘的效果并不十分理想。图4（b）和图5（b）中，平滑后天空和地面存在很多孤立噪声点，建筑物边缘也很模糊，图像中细节信息也有较大程度损失。从图4（c）和图5（c）来看，本文方法很好地实现了图像平滑和保持边缘的折衷，保证了图像视觉连续性，其视觉效果明显优于P-M方法。

　　对于图4（b）和图4（c）而言，在迭代次数相同的情况下，自适应参数高阶方法所能达到峰值信噪比比P-M方法所能达到峰值信噪比要高得多，并且平均耗时比P-M方法增加不到20%。表1为对于Cameraman图自适应参数高阶方法与P-M方法运行时间和峰值信噪比的对比结果。

　　本文首先对经典P-M方程存在两大问题进行分析，提出了自适应参数高阶偏微分方法，有效解决了P-M方法“阶梯”效应及其阈值参数选取问题，图像平滑效果比较好，在耗时相对不长的情况下所能达到的峰值信噪比也P-M方法高。