基于偏微分方程滤波器的图像去噪方法技术

本发明专利技术公开了一种基于偏微分方程滤波器的图像去噪方法，主要解决传统去噪方法去噪效果性能较差的问题。其实现过程包括：（１）输入噪声图像ｕ，计算该图像ｕ的偏导数和（２）计算噪声图像ｕ的梯度模值｜▽ｕ｜；（３）根据梯度▽ｕ和梯度模值｜▽ｕ｜，建立偏微分方程；（４）计算偏微分方程中的扩散系数和ψ；（５）利用系数和ψ，求解偏微分方程得到滤波图像；（６）计算滤波图像的峰值信噪比ＰＳＮＲ；（７）重复步骤１到步骤６，当某一次迭代输出的滤波图像的ＰＳＮＲ值小于上一次迭代输出的滤波图像的ＰＳＮＲ值时，终止迭代，输出上一次迭代的滤波图像。本发明专利技术利用图像的具体结构特征进行滤波，计算简单，运行速度快，能够平滑噪声的同时更好的保持图像纹理细节。可用于自然图像的去噪处理。

【技术实现步骤摘要】

本专利技术属于图像处理
，涉及图像去噪方法，适用于SAR图像以及自然图 像的噪声去除。
技术介绍
图像去噪旨在通过算法对被噪声污染的图像进行某种处理以降低噪声对原始有 用信息的影响尽可能的还原出更接近理想化的图像，它是开展森林资源调查、土地利用、覆 盖变化研究、环境灾害评估、城市规划、国防军情监控、医学影像和天文学影像等领域图像 处理中经常会用到的预处理技术，具有迫切的需要和广泛的应用前景。合成孔径雷达SAR 图像和自然图像都会需要去噪处理，研究SAR图像和自然图像去噪处理技术有着非常广阔 的应用前景。为了满足对图像去噪应用的迫切需要，目前涌现出了非常多的去噪方法，如小波 方法、beamlelshearlelcountlet、非局部均值方法等等。这些去噪方法虽然能够比较好 的实现去噪的目的，但是对于图像细节非常丰富的图像，这些方法的去噪结果都不够理想， 不能达到去噪效果的特殊要求。为了解决上述方法的问题，基于偏微分方程滤波器的去噪方法已成为图像去噪领 域研究的热点问题，许多学者从不同的角度对现有的去噪方法进行了分类，分析和改进，但 是对于图像不同的特征区域，去噪的效果不是十分理想。特别是对于一些实时性要求比较 高的图像，去噪的速度很慢，影响后续系统的正常工作。
技术实现思路
本专利技术的目的在于克服上述已有技术的缺点，提出了一种基于偏微分方程滤波器 的图像去噪方法，以自动实现基于SAR图像和自然图像特征的图像去噪，提高去噪的效果 和速度。实现本专利技术的技术方案是先计算待处理噪声图像的梯度模值和偏导数 du/dx,加/办，再将噪声图像代入基于图像特征的偏微分方程方法去噪模型中，通过求解 模型中的偏微分方程计算出各个像素点的恢复值，从而得到最终的滤波结果图像。具体步 骤包括如下(1)输入大小为MXN单幅噪声图像U，计算该图像u在1方向的偏导数和 y方向的偏导数彻/办；(2)利用梯度公式计算噪声图像u的梯度模值I Vu|; 其中，Vu为噪声图像u的梯度，Ux表示彻/Sx，uy表示加/办；(3)根据步骤(2)中计算出的梯度Vu和梯度模值I Vu|，建立偏微分方程如下 其中，洳/沩表示噪声图像U关于时间t的偏导数；φ为平坦区域的主扩散系数；ψ为边缘区域的主扩散系数；div(·)为散度;g(| Vu|)为扩散调节函数，用于控制噪声图像U在某个方向上的扩散程度，忿(|Vw|) = 1/(1 +(I)2)，其中k为用于判断某个像素点是图像边缘还是平坦区域的阈值，k= I^Tt，其中Ictl为初始值，t为迭代时间，t= At(n-l)，其中At为迭代步长，η为迭代次 数；u(0)表示零时刻的图像，u(0) = Utl表示零时刻的初始输入图像为U。；(4)计算偏微分方程中的平坦区域的主扩散系数识和边缘区域的主扩散系数Ψ 其中，h为经验值，取0.5 0.9 ；(5)利用计算出的偏微分方程中的系数识和Ψ，通过求解步骤(3)中的偏微分方程 得到每个像素点的灰度值，这些像素点组成滤波图像； (6)计算滤波图像的峰值信噪比PSNR = 201og10(255/RMSE)，其中，255是最大灰度级， 为步骤(5)中得到的滤波图像的像素灰度值，F(i，j)为步骤⑴中输入的噪声图像u的像素灰度值， i和j为图像中的像素坐标；(7)重复步骤1到步骤6，当某一次迭代输出的滤波图像的PSNR值小于上一次迭 代输出的滤波图像的PSNR值时，终止迭代，输出上一次迭代的滤波图像即去噪结果。本专利技术与现有的技术相比具有以下优点1.本专利技术由于提出的模型利用了扩散调节系数识和Ψ，针对不同的细节区域采取 不同的扩散平滑策略，具有较强的自适应性；2.本专利技术采用的迭代终止条件是一旦出现效果评价指数PSNR下降，去噪过程立 即停止，这样提高了去噪的运行速度；3.本专利技术采用随迭代次数不断更新的阈值K，使得去噪的精确性得到了很大的提 升，提高了去噪效果。附图说明 图1是本专利技术的实现流程图2是本专利技术中对Iena图的去噪结果对比图；图3是本专利技术中对Barbara图的去噪结果对比图；图4是本专利技术中对camera图的去噪结果对比图。具体实施例方式参照图1，本专利技术的具体实施步骤如下步骤1.输入大小为MXN单幅噪声图像U，计算该图像u在χ方向的偏导数^//δχ 和y方向的偏导数彻/办。步骤2.根据计算出的偏导数和 ，利用梯度公式计算噪声图像U的梯 度Vu和梯度模值I Vul。梯度▽U = (ux，Uy)，梯度模值 ，其中，Ux 轰示 ，Uy 表示du丨Qy .在图像去噪过程中，要求尽可能多的保留图像本身的细节信息并去除噪声的影 响。在图像细节较多的区域梯度模值较大，进行较少的扩散平滑；在细节较少的平坦区域梯 度模值较小，进行较多的扩散平滑。步骤3.根据步骤2中计算出的梯度Vu和梯度模值I Vu|，建立偏微分方程。建立偏微分方程如下 其中，彻/沩表示噪声图像U关于时间t的偏导数；<p为平坦区域的主扩散系数…值较大时，平坦区域主要按照偏微分方程中的前半 部分div(Vu/| VuD进行扩散；ψ为边缘区域的主扩散系数，Ψ值较大时，边缘区域主要按照偏微分方程中的后 半部分div(g(| Vu|) · Vu)来进行扩散；div(·)为散度；g(| Vu|)为扩散调节函数，用于控制噪声图像U在某个方向上的扩散程度， )，其中k为用于判断某个像素点是图像边缘还是平坦区域的阈值，k= I^Tt，其中Ictl为初始值，t为迭代时间，t= At(n-l)，其中At为迭代步长，η为迭代次 数；u(0)表示零时刻的图像，U(O) = U0表示零时刻的初始输入图像为U。；在图像u的平滑过程中，图像u的梯度I Vu|随着迭代演化不断的变化，所以判 断某个像素点是图像边缘还是平坦区域的阈值k不能再设置为常数，而应该是一个随迭代 次数不断变化的递减函数k = IvTt,其中e为指数常数，阈值不断的更新，去噪的效果就能 得到显著的提高。步骤4.计算偏微分方程中的平坦区域的主扩散系数φ和边缘区域的主扩散系数Ψ O平坦区域的主扩散系数口 其中，h为经验值，取0.5 0.9 ；在平坦区域，梯度模值I Vul较小，故φ值较小，则Ψ值较大，在div(g(| Vu|) .Vu)的扩散系数 Ψ较大的情况下，图像主要按照div(g(| VuD . Vu)来进行扩散平滑，可以比较好的去 除高梯度噪声；边缘区域的主扩散系数t/A 二 2(1-φ),在图像细节较丰富的区域，梯度模值| Vu|较 大，故炉值较大，则Ψ值较小，在div(VU/| VuD的扩散系数《P较大的情况下，图像主要按 照div(Vu/| VuD来进行扩散平滑，在去除部分噪声的同时能够较好的保留图像的原始 细节信息。步骤5.根据计算出的偏微分方程中的系数φ和Ψ，通过求解步骤3中的偏微分方 程得到每个像素点的灰度值，这些像素点组成滤波图像。偏微分方程的扩散调节系数φ和Ψ是针对不同的细节区域采取不同的扩散平滑 策略，将系数(？和Ψ带入偏微分方程中。在求解偏微分方程的过程中，可以自适应的调节扩 散方向去平滑图像的每个区域，在图像边缘较多的区域进行较少的扩散平滑作用，在边缘 较少的平坦区域进行较多的扩散平滑作用。步骤6.计算滤波图像的峰值信噪比。峰值信噪本文档来自技高网...

【技术保护点】
一种基于偏微分方程滤波器的图像去噪方法，包括如下步骤：（１）输入大小为Ｍ×Ｎ单幅噪声图像ｕ，计算该图像ｕ在ｘ方向的偏导数＊ｕ／＊ｘ和ｙ方向的偏导数＊ｕ／＊ｙ；（２）利用梯度公式计算噪声图像ｕ的梯度模值｜▽ｕ｜：▽ｕ＝（ｕ↓［ｘ］，ｕ↓［ｙ］），＊＊＊其中，▽ｕ为噪声图像ｕ的梯度，ｕ↓［ｘ］表示＊ｕ／＊ｘ，ｕ↓［ｙ］表示＊ｕ／＊ｙ；（３）根据步骤（２）中计算出的梯度▽ｕ和梯度模值｜▽ｕ｜，建立偏微分方程如下：＊＊＊其中，＊ｕ／＊ｔ表示噪声图像ｕ关于时间ｔ的偏导数；φ为平坦区域的主扩散系数；ψ为边缘区域的主扩散系数；ｄｉｖ（.）为散度；ｇ（｜▽ｕ｜）为扩散调节函数，用于控制噪声图像ｕ在某个方向上的扩散程度，ｇ（｜▽ｕ｜）＝１／（１＋▽ｕ／ｋ）↑［２］），其中ｋ为用于判断某个像素点是图像边缘还是平坦区域的阈值，ｋ＝ｋ↓［０］ｅ↑［－ｔ］，其中ｋ↓［０］为初始值，ｔ为迭代时间，ｔ＝Δｔ（ｎ－１），其中Δｔ为迭代步长，ｎ为迭代次数；ｕ（０）表示零时刻的图像，ｕ（０）＝ｕ↓［０］表示零时刻的初始输入图像为ｕ↓［０］；（４）计算偏微分方程中的平坦区域的主扩散系数φ和边缘区域的主扩散系数ψ：φ＝ｅ↑［１／｜▽ｕ｜］－ｈ，ψ＝２（１－φ），其中，ｈ为经验值，取０．５～０．９；（５）利用计算出的偏微分方程中的系数φ和ψ，通过求解步骤（３）中的偏微分方程得到每个像素点的灰度值，这些像素点组成滤波图像；（６）计算滤波图像的峰值信噪比：ＰＳＮＲ＝２０ｌｏｇ↓［１０］（２５５／ＲＭＳＥ），其中，２５５是最大灰度级，＊＊＊，ｆ（ｉ，ｊ）为步骤（５）中得到的滤波图像的像素灰度值，Ｆ（ｉ，ｊ）为步骤（１）中输入的噪声图像ｕ的像素灰度值，ｉ和ｊ为图像中的像素坐标；（７）重复步骤１到步骤６，当某一次迭代输出的滤波图像的ＰＳＮＲ值小于上一次迭代输出的滤波图像的ＰＳＮＲ值时，终止迭代，输出上一次迭代的滤波图像即去噪结果。...

【技术特征摘要】

【专利技术属性】
技术研发人员：张小华，焦李成，王然，王爽，侯彪，马文萍，尚荣华，盖超，张强，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：87[中国|西安]