专利名称：一种基于噪声空间特性的非线性滤波去噪方法
技术领域：
本发明属于数字图像处理技术领域，具体涉及一种基于噪声空间特性的非线性滤波去噪方法。
背景技术：
图像在形成、传输和存储过程中都往往会被噪声污染。噪声的强度一般采用信噪比(snr)或分贝数(dB)来定义。信噪比即是图像中信号和噪声的功率比值，信噪比和分贝数的转换关系为dB = IOlog10Snr一般地，噪声随机的分布在图像中，而且支撑域很小。因此传统的去噪方法是利用一种非线性的中值滤波器对图像进行滤波，去除图像中支撑域很小的噪声，从而实现对图像的去噪。对于一般情况下的弱噪声，传统的滤波去噪方法可以很好地实现图像去噪。但是， 在动平台条件下成像，成像积分时间很短，从而导致了传感器噪声很强；而某些时候，由于成像条件太差，或者目标距离成像太远，目标辐射强度很低，导致了噪声相对较强。面对这类强噪声的恶劣条件，传统的滤波去噪方法常常无能为力。其次传统的非线性滤波去噪方法(如中值滤波去噪、最大值滤波去噪、最小值滤波去噪、中点滤波去噪等方法)都没有考虑目标图像中噪声和目标的空间特征，噪声的空间特性(1)无方向性；(2)空间局部性，支撑域很小；形状目标的空间特性(1)目标有较大的支撑域，具有空间多尺度特性；( 目标的各个部件具有空间方向性。仅仅盲目的在一个模板(如3X3)内求平均值、最大值、最小值、中点值等来实现图像去噪。然而空间目标图像具有很强的结构性，去噪时其某些模板应该是求中值，一些模板是求最大值或最小值， 而其他的模板却是求中点值。利用传统的非线性滤波去噪方法对空间目标图像进行去噪， 会在去除噪声的同时丢失图像的细节，使得图像更加模糊，增加了后续处理工作(如图像校正、目标识别等)的难度。综上所述，需要研究有效的图像去噪方法，使得在图像去噪的同时保留图像的细节和图像的边缘，降低图像预处理和后续处理的难度。
发明内容
本发明的目的在于提供一种基于噪声空间特性的非线性滤波去噪方法，该方法弥补了传统非线性滤波去噪方法去噪时的盲目性，充分利用噪声空间特性和目标空间特性的差异，从而能在去噪的同时较好的保留图像的细节和图像的边缘。本发明提供的一种基于噪声空间特性的非线性滤波去噪方法，其特征在于，该方法步骤包括(1)输入红外成像探测器探测到的图像f ；(2)创建一幅与图像f大小相同，灰度值却全为0的图像g;
3
(3)对带噪声的图像f进行单点去噪若图像f中存在某一像素点(i，j)的灰度值大于其周围8个像素点的灰度值中的最大值，或者小于其周围8个像素点的灰度值中的最小值，则判定该点为噪声点；令图像g在像素点(i，j)的灰度值g(i，j)等于图像f中像素点(i，j)周围8个像素点的灰度值的平均值；(4)对带噪声的图像f进行两点去噪(4. 1)若图像f中存在相邻两像素点的灰度值大于该两点周围10个像素点的灰度值中的最大值，或小于该两点周围10个像素点的灰度值中的最小值，则判定该两点为噪声点；计算此相邻两点周围10个像素点的灰度值的平均值为Temp，如果图像g中该相邻两点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半， 否则就直接令图像g中该点的灰度值为Temp ；(4. 2)若图像f中存在相邻两像素点的灰度值大于该两点周围14个像素的灰度值中的最大值，或小于该两点周围14个像素的灰度值中的最小值，则判定该两点为噪声点； 计算此相邻两点周围14个像素灰度值的平均值为Temp，如果图像g中该相邻两点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半，否则就直接令图像g中该点的灰度值为Temp ；(5)对带噪声的图像f进行三点去噪若图像f中存在相邻3像素的灰度值都大于该3点周围13个像素的灰度值中的最大值，或小于该3点周围13个像素的灰度值中的最小值，则判定该3点均为噪声点；计算此相邻三点周围13个像素灰度值的平均值为Temp，如果图像g中该相邻三点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半，否则就直接令图像g中该点的灰度值为iTemp ；(6)令图像g中灰度值为零的像素的灰度值为图像f该点的灰度值。传统非线性滤波去噪方法虽能去除噪声但模糊了图像的细节和图像的边缘，限制了该方法的效能。同时，图像中噪声和目标的空间特性没有充分利用。本发明方法分析了噪声空间特性与形状目标空间特性的差异，提出一种单点模式的噪声形式、四种两点模式的噪声形式及四种三点模式的噪声形式。对图像中每个像素点进行单点模式、两点模式和三点模式的噪声形式的判定，如果该点判定为噪声点则作相应模式的滤波去噪处理；若不是噪声点则保留，从而实现了有效的图像去噪。本发明弥补了传统非线性滤波去噪方法去噪时的盲目性，充分利用噪声空间特性和目标空间特性的差异，从而能在去噪的同时较好的保留图像的细节和图像的边缘，便于图像校正、目标识别等后续处理工作的展开。
图1是本发明的基于噪声空间特性的非线性滤波去噪方法流程图。图2是单点形式的噪声模式示意图。图3a、图3b、图3c、图3d是两点形式的噪声模式示意图。图4a、图4b、图4c、图4d是三点形式的噪声模式示意图。图5是某空间目标在叠加信噪比为5db的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图5a是European仿真图像。
图恥是图fe叠加了信噪比为5db的高斯白噪声后的图像。图5c是中值滤波器的滤波结果。图5d是本发明算法的滤波结果。图6是某空间目标在叠加信噪比为3db的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图6a是European仿真图像。图6b是图6a叠加了信噪比为3db的高斯白噪声后的图像。图6c是中值滤波器的滤波结果。图6d是本发明算法的滤波结果。图7是某空间目标在叠加信噪比为Idb的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图7a是European仿真图像。图7b是图7a叠加了信噪比为Idb的高斯白噪声后的图像。图7c是中值滤波器的滤波结果。图7d是本发明算法的滤波结果。图8是某空间目标在叠加信噪比为5db的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图8a是HiAble仿真图像。图8b是图8a叠加了信噪比为5db的高斯白噪声后的图像。图8c是中值滤波器的滤波结果。图8d是本发明算法的滤波结果。图9是某空间目标在叠加信噪比为3db的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图9a是HiAble仿真图像。图9b是图9a叠加了信噪比为3db的高斯白噪声后的图像。图9c是中值滤波器的滤波结果。图9d是本发明的算法滤波结果。图10是某空间目标在叠加信噪比为Idb的高斯白噪声的条件下，中值滤波器和本发明提出的算法滤波的实验结果对比。图IOa是HiAble仿真图像。图IOb是图IOa叠加了信噪比为Idb的高斯白噪声后的图像。图IOc是中值滤波器的滤波结果。图IOd是本发明的算法滤波结果。
具体实施例方式下面结合附图和实例对本发明作进一步详细的说明(1)输入红外成像探测器探测到的图像f ；(2)创建一幅与图像f大小相同，灰度值却全为0的图像g;(3)首先对带噪声的图像f进行单点去噪
(3.1)若图像f在像素点(i，j)的灰度值f(i，j)大于其周围8个像素点的灰度值中的最大值，或者小于其周围8个像素点的灰度值中的最小值，则判定该点为噪声点，其示意图如图2所示，则令图像g在像素点(i，j)的灰度值g(i，j)等于图像f中像素点(i， j)周围8个像素点的灰度值的平均值；否则图像g在像素点(i，j)的灰度值不作修改；(3. 2)重复执行步骤(3. 1)直到判断完图像f的所有像素点，然后进行步骤⑷；(4)其次对带噪声的图像f进行两点去噪(4. 1)对形如图3a的两点模式进行两点去噪(4. 1. 1)若图像f在相邻两像素点(i，j)与(i，j+1) (4邻域相邻)的灰度值都大于该两点周围10个像素的灰度值中的最大值，或者都小于该两点周围10个像素的灰度值中的最小值，则判定该两点为噪声点，计算此相邻两点(i，j)与(i，j+1)周围10个像素点的灰度值的平均值为Temp，如果图像g中该相邻两像素点(i，j)与(i，j+l)中存在某点的灰度值非零(即该点的灰度值经过单点去噪时修改)，则令图像g在该点的灰度值为原来图像g在该点的灰度值和Temp求和的一半，否则直接令图像g在该点的灰度值为Temp ；如果该相邻两像素点判定不是噪声点则其灰度值都不作修改。(4. 1. 2)重复执行步骤(4. 1. 1)直到判断完图像f的所有像素点，然后进行步骤 (4. 2)；(4. 2)对形如图北的两点模式进行两点去噪对图像f中另外相邻两像素点(i， j)与(i+1，j) (4邻域相邻)，进行类似步骤(4. 1)的处理，然后进行步骤(4. 3)。(4. 3)对形如图3c的两点模式进行两点去噪(4. 3. 1)若图像f在相邻两像素点(i，j)与(i+1，j+1) (4邻域相邻)的灰度值都大于该两点周围14个像素点的灰度值中的最大值，或者都小于该两点周围14个像素点的灰度值中的最小值，则判定该两点为噪声点，计算此相邻两点(i，j)与(i+1，j+1)周围14 个像素点的灰度值的平均值为Temp，如果图像g中该相邻两点(i，j)与(i+1，j+1)中存在某点的灰度值非零(即该点的灰度值经过单点去噪修改或两点去噪的4. 1步、4. 2步修改)，则令图像g中该点的灰度值为原来图像g在该点的灰度值和Temp求和的一半，否则直接令图像g在该点的灰度值为Temp ；如果该相邻两像素点判定不是噪声点则其灰度值都不作修改。(4. 3. 2)重复执行步骤(4. 3. 1)直到判断完图像f的所有像素点，然后进行步骤 (4. 4)；(4. 4)对形如图3d的两点模式进行两点去噪对图像f中另外相邻两像素点(i， j+1)与(i+1, j) (4邻域相邻)，进行类似步骤(4. 3)的处理，然后进行步骤(5)。(5)然后对带噪声的图像f进行三点去噪(5. 1)对形如图如的三点模式进行三点去噪(5. 1. 1)若图像f中相邻三像素点(i，j)、(i，j+1)与(i+1, j+1) (4邻域相邻且 3点不在同一直线)的灰度值都大于该3点周围13个像素点的灰度值中的最大值，或者都小于该3点周围13个像素点的灰度值中的最小值，则判定该3点均为噪声点，计算此相邻三点(i，j)、(i，j+1)与(i+1, j+1)周围13个像素点的灰度值的平均值为Temp，如果图像 g中该相邻三点中存在某点的灰度值非零(即该点的灰度值经过单点去噪或两点去噪时修改)，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半，否则直接令图像g在该点的灰度值为Temp ；如果该相邻三像素点判定不是噪声点则其灰度值都不作修改。(5. 1. 2)重复执行步骤(5. 1. 1)直到判断完图像f的所有像素点，然后进行步骤 (5. 2)；(5. 2)对形如图4b的三点模式进行三点去噪对图像f中相邻三像素点(i，j)、 (i，j+1)与(i+1，j) (4邻域相邻)，进行类似步骤(5. 1)的处理，然后进行步骤(5. 3)。(5. 3)对形如图如的三点模式进行三点去噪对图像f中相邻三像素点(i，j)、 (i+1, j)与α+1，j+1) (4邻域相邻)，进行类似步骤(5. 1)的处理，然后进行步骤(5.4)。(5. 4)对形如图4d的三点模式进行三点去噪对图像f中相邻三像素点(i，j+1)、 (i+1, j)与(i+1, j+1) (4邻域相邻)，进行类似步骤(5. 1)的处理，然后进行步骤(6)。(6)令图像g中灰度值为零的像素的灰度值为图像f该点的灰度值。实例附图5至图10的所示的实验结果反映了本发明所提出的方法的有效性。在强噪声 (ldb、3db、5db等)条件下，中值滤波器虽然可以对图像进行较好的去噪，但是，在去噪的同时，也明显模糊了图像的边缘，丢失了图像的某些细节(如图5(3、6(3、7(、8(3、9(3、10(3)。而本发明提出的方法，在有效去除噪声的同时，很好地保护了图像的边缘和细节信息(如图5d、 6d、7d、8d、9d、10d)，从而验证了本发明的有效性。本发明不仅局限于上述具体实施方式
，本领域一般技术人员根据本发明公开的内容，可以采用其它多种具体实施方式
实施本发明，因此，凡是采用本发明的设计结构和思路，做一些简单的变化或更改的设计，都落入本发明保护的范围。
权利要求
1. 一种基于噪声空间特性的非线性滤波去噪方法，其特征在于，该方法步骤包括(1)输入红外成像探测器探测到的图像r；(2)创建一幅与图像f大小相同，灰度值却全为0的图像g；(3)对带噪声的图像f进行单点去噪若图像f中存在某一像素点(i，j)的灰度值大于其周围8个像素点的灰度值中的最大值，或者小于其周围8个像素点的灰度值中的最小值，则判定该点为噪声点；令图像g在像素点(i，j)的灰度值g(i，j)等于图像f中像素点 (i，j)周围8个像素点的灰度值的平均值；(4)对带噪声的图像f进行两点去噪(4. 1)若图像f中存在相邻两像素点的灰度值都大于该两点周围10个像素点的灰度值中的最大值，或小于该两点周围10个像素点的灰度值中的最小值，则判定该两点为噪声点；计算此相邻两点周围10个像素点的灰度值的平均值为Temp，如果图像g中该相邻两点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半， 否则就直接令图像g中该点的灰度值为Temp ；(4. 2)若图像f中存在相邻两像素点的灰度值都大于该两点周围14个像素的灰度值中的最大值，或小于该两点周围14个像素的灰度值中的最小值，则判定该两点为噪声点；计算此相邻两点周围14个像素灰度值的平均值为Temp，如果图像g中该相邻两点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半，否则就直接令图像g中该点的灰度值为Temp ；(5)对带噪声的图像f进行三点去噪若图像f中存在相邻3像素的灰度值都大于该3点周围13个像素的灰度值中的最大值，或小于该3点周围13个像素的灰度值中的最小值，则判定该3点均为噪声点；计算此相邻三点周围13个像素灰度值的平均值为Temp，如果图像g中该相邻三点中存在某点的灰度值非零，则令图像g中该点的灰度值为原来灰度值和Temp求和的一半，否则就直接令图像 g中该点的灰度值为iTemp ；(6)令图像g中灰度值为零的像素的灰度值为图像f该点的灰度值。
全文摘要
本发明公开了一种基于噪声空间特性的非线性滤波去噪方法，过程为①输入红外成像探测器探测到的图像f；②创建一幅与图像f大小相同，灰度值却全为0的图像g；③对带噪声的图像f进行单点去噪；④对带噪声的图像f进行两点去噪⑤对带噪声的图像f进行三点去噪；⑥令图像g中经过单点去噪、两点去噪及三点去噪修改后灰度值还为零的像素的灰度值为图像f该点的灰度值。本发明弥补了传统非线性滤波去噪方法去噪时的盲目性，充分利用噪声空间特性和目标空间特性的差异，从而能在去噪的同时较好的保留图像的细节和图像的边缘，便于图像校正、目标识别等后续处理工作的展开。
文档编号G06T5/00GK102184525SQ201110098209
公开日2011年9月14日 申请日期2011年4月19日 优先权日2011年4月19日
发明者关静, 周钢, 左芝勇, 张天序, 朱生国, 武道龙, 薛米生, 陈建冲, 陈浩 申请人:华中科技大学