分数阶偏微分方程在图像处理的应用.ppt

分数阶偏微分方程在图像处理的应用 刘路逍 数字图像处理简介 图像去噪一般可以分为两类：空域去噪方法和变换域去噪方法 空域去噪方法是直接对图像的像素进行处理，如均值滤波，中值滤波和偏微分方程滤波方法 变换域去噪方法主要是利用有用信号和噪声信号在变换域中表现出的不同特征来有效的去除噪声，如傅里叶变换，小波变换滤波方法 数字图像处理(Digital Image Processing)是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。在航天航空、生物医学、通讯工程等多个方面有广泛的应用。 偏微分方程图像去噪方法 偏微分方程（Partial Differential Equations?）图像处理一般是采用某一能量泛函,通过变分法,得到欧拉-拉格朗日方程，并用梯度下降法求得相应的解。 初次出现的PDE滤波模型是线性的热扩散方程，该模型的扩散行为是朝四周各个方向的，不可避免地会破坏图像的边缘等特色。为了克服这种缺陷，许多研究者从各种角度提出各种方法来避免这种“同向扩散（Isotropic Diffusion）”行为 ，于是就诞生了各种整数阶PDE滤波模型，如P-M模型，ROF模型等等。 分数阶偏微分方程图像处理的优点 从数学性质上讲,对纹理结构的本身特性而言,纹理是具有弱导数(即分数阶导数)特性的信息,整数阶微分算子并不适合于处理这类具有弱导数的信息。 上图为分数阶微分算子在阶次不同时的幅频特性。从图中可以清楚地得出：分数阶微分算子对信号都有加强的作用，并且加强的幅度随频率和微分阶次的增加而非线性地急剧增强，在信号的低频部分(ω1)，分数阶微分算子对信号的幅值进行了一定的提升，并且提升的幅度稍大于1阶次和2阶次微分； 在信号的高频部分(ω1)，分数阶微分算子对信号的幅值进行一定的提升，但是提升的幅度也明显小于1阶次和2阶次微分。 上述性质表明，分数阶微分算子在加强信号中高频成分的同时，对信号的甚低频分量进行了非线性保留。所以，分数阶微分可以大幅提升高频成分，增强中频成分，非线性保留低频成分。所以采用分数阶微分进行图像去噪时，不仅能够较好地保持图像边缘特征，还能较好地保留图像平滑区域内灰度变化不大的纹理细节信息。 分数阶导数的定义 分数阶微积分的定义主要分为空域中的定义和频域中的定义两大类,空域中的定义主要包括Grumwald-Letnikov定义、Riemann-Liouville定义和Caputo定义,频域中的定义主要包括在Fourier变换域、LaPlace变换域中的定义形式。 下面只介绍本文中使用的Grumwald-Letnikov分数阶导数 : 设 在闭区间 上连续，则 关于 的 阶导数可以表示为, 其中， 称为广义二项式系数， 为伽玛函数， 为取整运算 全变分（TVD）去噪模型 令u（x , y）为原始的清晰图像，f（x，y）为受到噪声污染的图像。即 n（x，y）为具有零均值、方差为 的随机噪声，则TV去噪模型可表示为： 其中：右边第一项为图像的正则项，它在极小化过程中可以起到抑制噪声的作用；第二项为保真项，它主要起保持图像边缘特征和降低图像失真度的作用。 为梯度算子； 为图像的定义域，像素点(x , y) ∈Ω。 该模型存在的不足：一方面，由于该模型中正则项仅包含|▽u|，很多研究结果表明该模型本身存在固有缺陷，即在处理平滑区域时，有可能将噪声当成边缘，从而容易产生“阶梯”效应；另一方面，虽然TV图像去噪方法进行图像去噪有利于保持图像边缘信息，但对于图像的纹理细节的刻画却不够理想。因此有必要在此模型上进行改进。 新分数阶全变分图像去噪模型 全变分模型的一般形式： 受分数阶微分理论和TV模型的启发，本文将上式中的整数阶微分 运算用 阶微分 来替换： 其中 是一个 的对角矩阵， 代表第i个像素的正则化参数， ，其中 代表第i个元素微分的分数阶次，且 为离散化分数阶梯度算子， 这样替换有两个好处: 1) 全变分模型去噪效果好与保持图像边缘细节特征的优点得到了很好地继承; 2)分数阶微分可以大幅提升图像高频成分、增强图像中频成分、非线性保留图像低频成分的特性, 较好地保留了更多图像平滑区域中灰度变化不大的纹理细节信息。这样一来, 本文建立的新模型较好地解决在去除噪声的同时, 保持图像边缘特征这一矛盾,并且还能保留了图像更多的纹理细节信息。 分别代表x和y方向的分数阶差分算子。 为了求解上式，我们应用交替最小化的步骤，即，对 ，依次求解 运用半二次正则化方法与分离变量法将问题转化， 前两式的每一步迭代都有确切的解 , 第三个最小化的解 可以