基于非局部双边滤波的极化sar相干斑噪声抑制方法
【专利摘要】本发明公开了一种基于非局部双边滤波的极化SAR相干斑噪声的抑制方法，主要解决现有技术边缘纹理细节丢失的问题。其实现过程是：1.对极化SAR图像寻找相似的像素点，组成相似像素点集，得到相似性权重；2.以每一个相似像素点为中心取邻域窗，对邻域窗进行双边滤波，得到滤波结果；3.对相似像素点的所有滤波结果进行非局部权重平均，得到初始滤波后的协方差矩阵；4.对初始滤波后的协方差矩阵进行滤波，得到最终滤波后的协方差矩阵；5.将最终滤波后的协方差矩阵合成伪彩图，并将其转换为灰度图像。本发明能在抑制相干斑的同时保持边缘纹理细节信息，提高了极化SAR图像的相干斑抑制效果，可用于极化SAR图像地物分类和目标识别。
【专利说明】基于非局部双边滤波的极化SAR相干斑噪声抑制方法
【技术领域】
[0001]本发明属于图像数据处理【技术领域】，具体地说是一种相干斑抑制方法，可用于极化合成孔径雷达SAR图像数据的去噪，提高图像的后续解译能力。
【背景技术】
[0002]随着雷达技术的发展，极化SAR已成为SAR的发展趋势，极化SAR能够得到更丰富的目标信息，在提高目标检测，辨别和分类能力等方面，能够体现极化SAR系统的优势。但和SAR图像一样，极化SAR图像中也不可避免的会出现相干斑，直接影响了对数据信息的提取与解译。因此，相干斑噪声的抑制已经成为极化SAR数据信息处理的首要问题，也是国内外雷达图像处理与应用的热点之一。对于极化SAR数据，相干斑抑制的目的在于能够在抑制相干斑的同时保持住数据的极化特性，边缘细节以及纹理信息。现有对极化SAR数据相干斑抑制的方法主要涉及以下几种:
[0003]一、极化白化滤波PWF。这种方法是最早的一个滤波方法，该方法主要是针对Span数据进行相干斑抑制，而其余极化SAR数据的各元素并没有进行相干斑的抑制。
[0004]二、精致极化Lee滤波。它是在SAR图像Lee滤波的基础上针对极化SAR数据的特点改进得到的最具代表性的一种滤波方法。通过使用边缘窗口进行滤波，使得滤波后的数据在边缘的特性保持方面效果显著，但是，在纹理细节信息的保持上，其滤波效果并非特别理想，因此在相干斑的抑制中，数据原始的一些特性无法很好的保留。
[0005]三、双边滤波。这种方法在处理相邻各像素的相关信息时，既考虑了几何上的邻近关系又考虑了亮度上的相似性，因此，用于自然图像去噪时，取得了比较好的效果。依据极化SAR数据噪声模型及复矩阵形式，最近有人将双边滤波推广至极化SAR数据处理，该方法可以直接对极化协方差矩阵进行处理，在对相干斑具有较好的抑制能力的同时能够很好地保持极化信息。但由于此方法是基于局部邻域窗，会产生局部块效应。
[0006]四、非局部均值滤波。这种方法用在自然图像的噪声抑制上取得了非常优异的效果。将非局部的方法扩展到极化SAR相干斑抑制，仍然保留了自然图像噪声抑制的诸多优点，比如大大抑制了相干斑，同时边缘、纹理和点目标保持良好，极化SAR的极化特性也得到了保持。但是由于有参数的选择，造成保留细节、点目标和最大抑制噪声方面的矛盾，例如邻域窗选取过大时，可以更好的抑制噪声，但可能会丢失点目标和细节；邻域窗选取过小时，可以较好的保留点目标和细节，但噪声抑制能力相对减弱。
【发明内容】
[0007]本发明的目的在于针对上述已有技术的缺点，提出一种基于非局部双边滤波的极化SAR相干斑噪声抑制方法，以实现在抑制相干斑的同时保持亮目标和边缘纹理细节信息，提高极化SAR图像数据的相干斑抑制效果。
[0008]实现本发明的技术思路是:将非局部滤波与双边滤波相结合应用到极化SAR图像数据相干斑抑制中，其具体步骤包括如下:
[0009](I)读取极化SAR图像数据，进行边界镜像对称扩展，扩展半径为7 ；
[0010](2)利用非局部滤波方法，寻找极化SAR图像数据中每个像素点X的相似像素点，组成相似像素点集S，得到第j个相似像素点的相似性权重wml (j)，I彡j彡J，J为相似像素点集S中像素点的个数；
[0011 ] (3)从像素点X的相似像素点集S中选取第j个相似像素点y」，以y」为中心取5 X 5大小的邻域窗Yp并以像素点X为参考相素，对邻域窗Yj进行双边滤波，得到滤波后的协方差矩阵C\r,:h
[0012](4)对滤波后的协方差矩阵进行权重wMl (j)平均，得到滤波后的协方差矩阵 ?χ);
[0013](5)对极化SAR图像数据的每个像素点，进行步骤(3)—步骤⑷处理，得到整个极化SAR图像数据初始滤波后的协方差矩阵(
[0014](6)对步骤⑵和步骤(5)进行2次迭代处理，得到最终的滤波后的协方差矩阵
Λ.Λ
C.
[0015](7)用Pauli向量法将最终的滤波后的协方差矩阵^合成伪彩图，再将伪彩图转换为灰度图像。
[0016]本发明与现有技术相比具有以下优点:
[0017]1.本发明首先基于图像块的相似性找到相似像素点，然后以相似像素点为中心取邻域窗，并以要滤波像素点为参考对邻域窗进行双边滤波；在考虑了图像的结构信息的同时考虑了单个像素点的相似性，能够找到更好的相似像素点，不仅能平滑同质区域，而且对边缘和纹理细节的保留效果理想；
[0018]2.本发明基于非局部的思想，在非局部区域寻找相似像素点，克服了双边滤波在局部区域选择相似像素点，导致块效应，纹理和边缘的保持也很难让人满意的缺点，同时基于双边滤波的思想，处理找到的相似像素点的邻域区域，既考虑了结构信息又考虑了单像素点信息，克服了非局部均值滤波依赖于邻域窗大小的选择，邻域窗太大会丢失点目标，太小会模糊边缘和纹理的缺陷；
[0019]3.本发明采用了像素的协方差矩阵信息来计算两个像素点的相似距离，保证了选出的相似像素点的合理性，保持了极化信息。
【专利附图】
【附图说明】
[0020]图1是本发明的实现流程图；
[0021]图2是本发明仿真使用的模拟极化SAR图像和两幅真实的极化SAR图像数据图；
[0022]图3是用现有方法及本发明对模拟SAR图像的滤波结果图；
[0023]图4是用现有方法及本发明对第一幅真实极化SAR图像的滤波结果图；
[0024]图5是用现有方法及本发明对第一幅真实极化SAR图像图像块A的滤波放大后结果图；
[0025]图6是用现有方法及本发明对第一幅真实极化SAR图像图像块B的滤波放大后结果图；
[0026]图7是用现有方法及本发明对第二幅真实极化SAR图像的滤波结果图。
【具体实施方式】
[0027]参照图1，本发明的实现步骤如下:
[0028]步骤I，读取极化SAR图像数据。
[0029]读取的极化SAR图像数据的每个像素点为含有9个元素的3X3协方差矩阵C，协方差矩阵C定义如下:
【权利要求】
1.一种基于非局部双边滤波的极化SAR相干斑噪声抑制方法，包括如下步骤: (1)读取极化SAR图像数据，并进行边界镜像对称扩展，扩展半径为7； (2)利用非局部滤波方法，寻找极化SAR图像数据中每个像素点X的相似像素点，组成相似像素点集S，得到第j个相似像素点的相似性权重Wml (j), I≤j ≤ J, J为相似像素点集S中像素点的个数； (3)从像素点X的相似像素点集S中选取第j个相似像素点y」，以y」为中心取5X5大小的邻域窗Yj，并以像素点X为参考相素，对邻域窗Yj进行双边滤波，得到滤波后的协方差矩阵; (4)对滤波后的协方差矩阵("(&)进行权重(j)平均，得到滤波后的协方差矩阵C(-V)； (5)对极化SAR图像数据的每个像素点，进行步骤(3)—步骤(4)处理，得到整个极化SAR图像数据初始滤波后的协方差矩阵 (6)对步骤(2)和步骤(5)进行2次迭代处理，得到最终的滤波后的协方差矩阵 (7)用Pauli向量法将最终的滤波后的协方差矩阵^合成伪彩图，再将伪彩图转换为灰度图像。
2.根据权利要求1所述的基于非局部双边滤波的极化SAR相干斑噪声抑制方法，其中步骤(2)所述的利用非局部滤波方法，寻找像素点X的相似像素点，组成相似像素点集S，得到相似性权重Wntjl (j)，按如下步骤进行: 2a)以像素点X为中心取5X5大小的邻域窗，记为图像块X，以像素点x为中心取15X15大小的邻域窗，记为搜索窗Ω，在搜索窗Ω内以各像素点为中心取5X5大小的邻域窗，记为图像块八￥21~4~^"』为搜索窗Ω内的像素点数；计算图像块X和搜索窗内各图像块Y1，Y2，".Υη…，YN2间的相似距离:
式中，k为图像块内像素点的个数，取值为25，?/Λ., =61n2 + ln|^(| + ln|^|-21n|x>}f I
表示两个像素点之间的相似距离，IXiI为图像块X的第i个像素点矩阵的行列式值，IYinI为第η个图像块Yn的第i个像素点矩阵的行列式值； 2b)对步骤2a)得到的相似距离七y与给定的一个阈值T做比较，如果>『，则将第η个图像块Yn所对应的中心像素点放入相似像素点集S，并得到相似性权重Wml (j)，其中阈值T和权重wnaL (j)的定义分别如下:
其中K是一个调节参数，设定为K = 20，I为极化SAR图像数据的视数，取值为4。
3.根据权利要求1所述的基于非局部双边的极化SAR相干斑噪声抑制方法，其中步骤(3)所述的以&为中心取5X5大小的邻域窗Yp对邻域窗I进行双边滤波，按如下步骤进行: 3a)设置空间平滑参数σ d = 2.9，相似性散度平滑参数σ h = 3.1 ; 3b)对邻域窗Yp根据公式
求解空间距离度量的权重值fs，其中为邻域窗Y」的中心像素点，&为邻域窗Y」中的任意点，||*||为2范数，exp (.)为指数函数； 3c)以像素点X为参考相素点，求解邻域图像块Yj内的相似性散度距离度量的权重值fr:
其中，C(x)为像素点X处的协方差矩阵，
为距离度量，用来度量极化SAR数据中两像素点矩阵的相似性，det (.)表示协方差矩阵C行列式的值； 3d)根据步骤3b)和3c)，计算双边滤波的权重Wbf (ji)和滤波后的协方差矩?.ν, ):
其中，fs( Yj-Yji I)为图像块Yj的空间距离度量的权重值，fr Irf(COc),CX.r,.,.))j为图像块I的相似性散度距离度量的权重值。
4.根据权利要求1所述的基于非局部双边的极化SAR相干斑噪声抑制方法，其中步骤(4)所述的对滤波后的协方差矩阵(\^)进行权重Wntjl(j)平均，得到滤波后的协方差矩阵Π.Ο.按如下公式进行:
其中
为非局部归一化滤波权重，Yj为相似像素点集S中的第j个相似像素点。
【文档编号】G06T5/00GK104182945SQ201410437634
【公开日】2014年12月3日 申请日期:2014年8月30日 优先权日:2014年8月30日
【发明者】王爽, 马文萍, 雷晓珍, 霍丽娜, 侯彪, 刘坤, 谢慧明 申请人:西安电子科技大学