最小十字窗中值与均值迭代滤波的图像噪声滤除方法
【技术领域】
[0001] 本发明属于图像处理技术领域，特别设及一种可用于航空航天、工业、农业、医学， 军事等领域的数字图像处理的最小十字窗中值与均值迭代滤波的图像噪声滤除方法。
【背景技术】
[0002] 图像在景物成像、空间采样和量化等过程中经常受外界各种噪声干扰，使图像质 量下降。而图像噪声是影响人眼观察最直接、最有害和最关键的问题。为尽可能减小噪声影 响，对噪声污染的图像必须进行去噪处理。1971年，著名学者化k巧在他的开拓性论文中，提 出了一种非线性滤波器一中值滤波器。研究表明，中值滤波在一定的条件下可W克服线性 滤波所带来的图像模糊，而且是消除图像脉冲噪声的最佳办法。其原理是一个含有奇数像 素的滑动窗口，窗口中屯、的像素灰度值用窗口内像素灰度值的中值代替。然而正如众多文 献指出，中值滤波从一维信号扩展到多维信号处理，尽管中值滤波具有良好的窄脉干扰抑 审IJ能力，但在二维图像中的点、细线、尖锐顶角，经中值滤波后，可能丢失；另外其去噪效果 与窗口大小有关，窗口小，去噪效果差一点，保护细节好一点，运行速度快一点，反之，窗口 大，去噪效果好些，运行速度慢些，保护细节差些。因为标准中值滤波的缺陷，各种中值滤波 的改进算法应运而生，例如自适应中值滤波算法、开关中值滤波算法、加权中值滤波算法、 =态中值滤波算法、软切换中值滤波等等，运些算法在去噪和保护细节方面都有较好的改 进效果，但处理速度并不理想，不能满足实时要求。如传统的中值滤波算法对每个窗口内的 像素进行排序求中值，其复杂度（W比较次数来衡量)均正比于排序对象个数之平方，且当 窗口尺寸增大时，计算量将按四次方增大，很难满足数字图像处理的实时性要求。基于此， 张新明等人提出了一种基于十字滑动窗口的快速自适应图像中值滤波方法(张新明，党留 群，徐久成.基于十字滑动窗口的快速自适应图像中值滤波.计算机工程与应用，2007， 43(27) :37-39.),提高了细节保护能力和运行速度，但运种方法由于图像非噪声点的信息 利用率不高，且采用是非裁剪中值滤波等原因，所W去噪效果有限。
【发明内容】
[0003] 本发明的目的在于克服上述已有技术的不足，提出了一种最小十字窗中值与均值 迭代滤波的图像噪声滤除方法。本发明不仅是因为最小十字窗口是最少像素排序窗口，能 提高滤波速度，而且能很好地保护图像细节。
[0004] 本发明的目的是运样实现的： 一种最小十字窗中值与均值迭代滤波的图像噪声滤除方法，包括如下步骤： 步骤1:输入大小为mXn，灰度级为0至化之间的含噪声图像I，其中最小灰度级为0，最大 灰度级为L，它常取255; 步骤2 :用极值方法判断噪声像素点，产生噪声0-1二值映射矩阵Ni，即：
步骤4:比较噪声密度P与参数C的大小，依据比较结果采用最小十字窗口迭代中值滤波 或者最小十字窗口迭代均值滤波，其中，C取O . 525。设R为恢复图像，其对应的噪声0-1二值 映射矩阵妒，令R=I ,N^=Ni; 步骤5:如果P小于或等于C，采用最小十字窗口迭代中值滤波，具体步骤如下： 步骤5.1:对于R中的每一个像素(x,y)，采用最小十字窗口中值滤波，最小十字窗口对 应的点为：.换"扣知父;古一，…&父f如料; Ri: W上采用了开关滤波技术
，即对噪声点(Ni(x，y)=0)作中值滤波，非噪声点(Ni(x，y) = 1)保持不变;也采用了裁剪滤波技术，即只对Sxy中为1的点（非噪声点)排序取中值，W下类 同， 如果Sxy为零矩阵，则RiU，y)=0，其噪声0-1二值映射矩阵为：
賣法结束，输出去噪结果； 步骤5.2:对于Ri中的每一个像素(X, y )，采用最小十字窗口中值滤波，最小十字窗口对 应的点为：二妍…1災…化妒妒n、沪皆夺!.，>')}，获得滤 波结果R2:
W上采用递归滤波技术，即对前一次滤波结果Ri再作中值滤波，W下类同。
[0005]如果Sxy为零矩阵，贝化2(x，y)=0,其噪声0-1二值映射矩阵为：
如身
，则R=I?2，算法结束，输出去噪结果； 步骤5.3:对于R2中的每一个像素 (X,y)，采用最小十字窗口中值滤波，最小十字窗口对 应的点为:?-{災&咐-押炎瑪碱絮…耗狄孩終冻采满竭战；V-f觀祕换寺1，键}，获得 滤波结果R3:
如果Sxy为零矩阵,则Rs (X, y )=0; 其中，步骤5.巧1设骤5.3中的1116(11日]1表示排序后取中值，从上义=1,2，。'，111，7=1,2， …，n，W下类同； 步骤5.4:用R3替换R，并更新护和妒，即：
:，否则转到步 骤5.1进行最小十字窗口迭代中值滤波； 步骤6:如果P大于C，采用最小十字窗口迭代均值滤波，具体步骤如下： 步骤6.1:对于R中的每一个像素(x,y)，采用最小十字窗口均值滤波，最小十字窗口对 应的点为， 获得滤波结呆Ki:
如果Sxy为零矩阵，则Ri (X，y)=0，其噪声0-1二值映射矩阵为：
结果； 步骤6.2:对于Ri中的每一个像素(x，y)，采用最小十字窗口均值滤波，最小十字窗口对 应的点为：
，获得 滤波结果R2:
， 束，输出去噪结果； 步骤6.3:对于化中的每一个像素(x，y)，采用最小十字窗口均值滤波，最小十字窗口对 应的点为：
^?，获得 滤波结果R3:
如果Sxy为零矩阵，则R3 (X，y )=0; 其中，步骤6.1到步骤6.3中的mean表示取均值； 步骤6.4:用R3替换R，并更新护和妒，即：
I结果R，否则转到 步骤6.1进行最小十字窗口迭代均值滤波； 步骤7如果噪声已处理完，则输出去噪结果。
[0006]本发明具有W下优点： 第一、本发明首次使用最小十字模板将开关滤波、递归滤波、裁剪滤波、中值滤波和均 值滤波等技术融合在一起，大幅度提高去噪性能和图像细节保持度。由于采用开关滤波，即 仅对噪声点作均值滤波或中值滤波，而非噪声点保留，因此能避免细节模糊； 第二、由于采用递归滤波，即下一次滤波是针对上次滤波的结果进行，本发明中用1次 普通滤波和2次递归滤波的值取代噪声点的值，因此能提高了去噪效果； 第=、由于采用裁剪滤波，即仅采用非噪声点值进行均值计算或排序取中值，故既能提 高去噪效果，又能降低计算复杂度； 第四、由于最小十字窗口是最少像素窗口，不仅能提高滤波速度，而且它能很好地保护 图像细节，尤其可应用到高保真恢复图像的场合，另外也无需考虑设置窗口大小的问题； 第五、中值滤波在有效数据较多的情况下效果好(在低噪声密度情况下，去噪效果好）， 但排序操作增加计算复杂度;而均值滤波在有效数据较少的情况下有效(在高噪声密度情 况下，去噪效果好），且计算复杂度低。本发明由于将中值滤波与均值滤波有机结合，即在低 密度噪声时采用最小十字窗中值滤波，而在高密度噪声时，采用最小十字窗均值滤波，不仅 能发挥中值滤波与均值滤波的各自优势，又能弥补二者不足，故能够提高去噪效率； 第六、本发明首次采用最小十字模板迭代中值滤波和迭代均值滤波，由于采用最小十 字模板的少数像素点而不是采用整个窗口中像素点，降低了计算复杂度，所W滤除噪声的 运行速度快，可应用到实时处理的场合。
【附图说明】
[0007] 图1是本发明提供的最小十字窗中值与均值迭代滤波的图像噪声滤除方法的流 程图。
[0008] 图2是本发明提供的最小十字窗中值与均值迭代滤波的图像噪声滤除方法中采用 的最小十字窗（黑色菱形组成的图案)和测试使用的自然图像。
[0009] 图3是本发明提供的最小十字窗中值与均值迭代滤波的图像噪声滤除方法测试使 用的含噪声的图像，其中，图3(a)对图2(b)添加了强度为0.1的椒盐噪声图像，图3(b)对图2 (b)添加了强度为0.5的椒盐噪声图像，图3(c)对图2(b)添加了强度为0.9的椒盐噪声图像。
[0010] 图4是采用本发明提供的最小十字窗中值与均值迭代滤波的图像噪声滤除方法针 对图3中的图像处理后得到的去噪结果图，其中图4(a)是对图3(a)的去噪结果，图4(b)是对 图3(b)的去噪结果，图4(c)是对图3(c)的去噪结果。
[0011 ] 图5是MNMF方法、RAMF方法、SMMF1方法、IMF方法及本发明方法在Lena图像上的去 噪结果PSNR/IEF值对比图表。
【具体实施方式】
[0012] 本发明的核屯、是提出了一种最小十字窗中值与均值迭代滤波的图像噪声滤除方 法。
[0013] 下面结合附图对本发明的内容作进一步说明： 一种最小十字窗中值与均值迭代滤波的图像噪声滤除方法，如图1所示，包括如下步 骤： 步骤1:输入大小为mXn，灰度级为0至化之间的含噪声图像I，其中最小灰度级为0，最大 灰度级为L，它常取255; 步骤2:用极值方法判断噪声，产生噪声0-1二值映射矩阵Ni，即：
步骤4:比较噪声密度P与参数C的大小，依据比较结果采用最小十字窗口迭代中值滤波 或者最小十字窗口迭代均值滤波。其中，参数C取0.525。设R为恢复图像，其对应的噪声0-1 二值映射矩阵妒，令R=I，N^=Ni; 步骤5:如果P小于或等于C，采用最小十字窗口迭代中值滤波，具体步骤如下： 步骤5.1:对于R中的每一个像素(