基于核回归的sar图像相干斑抑制的制作方法
【专利摘要】本发明公开了一种基于核回归的SAR图像相干斑抑制方法，主要是解决现有的核回归去噪方法不能应用于噪声模型为乘性的SAR图像中。其实现过程是：对SAR图像选取邻域窗；用索贝尔sobel算子计算图像块的水平梯度与垂直梯度；依据得到的梯度求得局部自适应协方差矩阵，计算平滑矩阵；通过权影响估计该像素点的理想值得到基本的去噪结果；在基本的滤波结果上重复该过程，得到满意的滤波效果。本发明在抑制相干斑的同时保持点目标和边缘纹理细节信息，提高了SAR图像的相干斑抑制效果，可用于SAR图像目标识别和地物分类。
【专利说明】基于核回归的SAR图像相干斑抑制
【技术领域】
[0001]本发明属于遥感图像数据处理【技术领域】，具体地说是一种图像相干斑抑制方法，该方法可用于SAR图像的相干斑抑制。
【背景技术】
[0002]合成孔径雷达SAR是一种高分辨率成像雷达。SAR不受气候、昼夜的影响，具备全天时、全天候的特点，在对地观测系统中，装备高分辨率的合成孔径雷达的航空系统机动性好、响应速度快，在灾害应急处理方面优势显著。作为卫星数据的有效补充，雷达在地形测绘、灾害监测评估、环境与国土资源监测等领域具有广泛的应用价值。特别是在灾害应急处理方面，雷达系统能及时为灾害应急处理提供信息丰富、清晰准确的图像，对及时开展救援、进行受灾评估、减轻灾害损失具有重要意义。但由于SAR系统采用的是相干微波源进行成像，回波的振幅在相干叠加时会相互叠加或抵消，得到的总回波强度与子回波平均强度之间会存在一定的偏差，在SAR图像中这种随机偏差以强噪声的形式呈现，这种噪声通常被称为相干斑。因为原理的缺陷性，SAR图像中相干斑影响了后续的图像解译，因此如何抑制SAR图像中的相干斑，提高图像的解读能力以及获得更多的信息成为一个重要的问题。SAR图像相干斑抑制的目标是在去除图像中的相干斑及尽可能的保持图像的细节信息。
[0003]SAR图像的相干斑是一种复杂的乘性噪声模型，对于斑点噪声这种特殊的性质，在过去的近二十年中，人们提出了很多经典的SAR图像相干斑抑制方法。如Lee滤波，增强Lee滤波，Kuan滤波等。这些方法是用一个已定义好的滤波器窗来估计局域相干斑的方差，并进行滤波处理，其结果通常过分的平滑了图像的边缘细节信息，在一定程度上取得了较好的效果，但存在模糊了边缘和点目标的缺陷。除了在空域进行图像相干斑抑制外，1995年，美国学者Donoho把小波理论引入到图像去噪中，提出了小波软阈值方法。小波软阈值法开创了变换域进行图像去噪的先河，之后涌现了许多变换域去噪的优秀成果，包括多尺度变换的图像去噪方法。小波软阈值方法是一种非线性的算法，依然存在破坏图像细节信息的问题，对图像的辐射特性保持也不好。
[0004]核回归方法本质上是一种空域局部平均的去噪方法，不同的是它依赖于被处理数据本身来指定这个模型的结构，这个模型是基于泰勒局部展开，采用的是非参数化方法，模型建立在图像的数据特征上。核回归分析应用于图像去噪最早是由H.Takeda和S.Farsiu等人于2005年提出，随后又得到了进一步的发展。其滤波效果与目前非常优秀的非局部均值方法NLM和三维块匹配的联合滤波BM3D方法相差无几，但是以上核回归算法目前只应用于噪声模型为加性的自然图像中，不能用到噪声模型为乘性的SAR图像中。
【发明内容】
[0005]本发明的目的在于提出了一种基于核回归的SAR图像相干斑抑制方法，将性能优异的核回归去噪方法扩展到噪声模型为乘性的SAR图像中，实现对SAR图像充分滤波的同时保持点目标和边缘细节信息。[0006]为实现上述目的，本发明技术方案包括如下步骤:
[0007](I)读入一幅SAR图像，并进行边界镜像对称扩展，扩展半径为N，取图像中的一个像素点X，以该像素点为中心取(2XN+1)大小的邻域窗W;
[0008](2)用索贝尔(Sobel)算子计算邻域窗内w每一个像素点的水平梯度与垂直梯度；
[0009](3)根据得到的水平梯度与垂直梯度计算以像素点X为中心的邻域窗w自适应协方差矩阵C ；
[0010](4)依据步骤⑷得到的局部自适应协方差矩阵C，计算以像素点X为中心的邻域窗w平滑矩阵H;
[0011](5)依据计算得到的平滑矩阵H，利用权重影响，求得该像素点步骤I中的像素点X的相干斑抑制结果；
[0012](6)对SAR图像每一个像素点按步骤⑴一(5)进行处理，得到每一个像素点的相干斑抑制结果，从而得到所述SAR图像的相干斑抑制结果。
[0013](7)将步骤6得到的结果，按照步骤⑴一步骤(6)进行迭代处理，得到最终的所述SAR图像的相干斑抑制结果。
[0014]本发明与现有技术相比具有以下优点:
[0015]1、本发明是基于SAR图像的乘性噪声模型而提出的，更加符合SAR图像的相干斑抑制，对同质区域的滤波更加平滑；
[0016]2、本发明充分利用了图像的结构信息，克服了增强Lee的滤波等方法无法表示图像大量丰富的局部结构；
[0017]3、本发明用索贝尔Sobel算子计算图像的初始水平梯度与垂直梯度，减少计算量，能够快速有效的处理大幅面的SAR图像；
【专利附图】
【附图说明】
[0018]图1是本发明的流程图；
[0019]图2是本发明仿真使用的barbara自然图像模拟的SAR图像；
[0020]图3是本发明仿真使用的光学图像模拟的SAR图像；
[0021]图4是本发明仿真使用的真实的幅度SAR图像；
[0022]图5是本发明仿真使用的真实的强度SAR图像；
[0023]图6是用现有方法及本发明对barbara自然图像模拟的SAR图像滤波后的结果图；
[0024]图7是用现有方法及本发明对光学图像模拟的SAR图像滤波后的结果图；
[0025]图8是用现有方法及本发明对真实的幅度SAR图像滤波后的结果图；
[0026]图9是用现有方法及本发明对真实的强度SAR图像滤波后的结果图。
[0027]具体实施步骤
[0028]参照图1，本发明的具体实现步骤如下:
[0029]步骤1，过二幅SAR图像，进行边界镜像对称扩展，扩展半径为N，取图像中的一个像素点X，以该像素点为中心取(2XN+1)大小的邻域窗W。
[0030]邻域窗太小会出现过平滑现象，太大会导致降斑效果不明显。优选地，N=3。[0031]步骤2，用索贝尔(Sobel)算子计算邻域窗内每一个像素点的水平梯度与垂直梯度；
[0032]具体过程如下:
[0033]2a)设置水平梯度算子模板与垂直梯度算子模板分别用来检测图像的垂直边缘及水平边缘，模板内的数字为模板系数，水平梯度算子模板的方向与垂直边缘的方向相垂直，垂直梯度算子模板的方向与水平边缘的方向相垂直。
【权利要求】
1.基于核回归的SAR图像相干斑抑制方法，包括如下步骤: 步骤1，对一幅SAR图像，进行边界镜像对称扩展，扩展半径为N，取图像中的一个像素点X，以该像素点为中心取(2XN+1)大小的邻域窗w ; 步骤2，用索贝尔(Sobel)算子计算邻域窗内每一个像素点的水平梯度与垂直梯度；步骤3，根据步骤2得到的水平梯度与垂直梯度计算以像素点X为中心的邻域窗w的自适应协方差矩阵C ； 步骤4，依据步骤3得到的局部自适应协方差矩阵C，计算以像素点X为中心的邻域窗w平滑矩阵H ； 步骤5，依据步骤4得到的平滑矩阵H，求得步骤I中的像素点X的相干斑抑制结果；步骤6，对步骤I中所述的SAR图像每一个像素点按步骤1-5进行处理，得到每一个像素点的相干斑抑制结果，从而得到所述SAR图像的相干斑抑制结果； 步骤7，将步骤6得到的结果，按照步骤1-步骤6进行迭代处理，得到最终的所述SAR图像的相干斑抑制结果。
2.如权利要求1所述的SAR图像相干斑抑制方法，其中，N=3。
3.如权利要求1所述的SAR图像相干斑抑制方法，其中，用索贝尔(Sobel)算子计算邻域窗内每一个像素点的水平梯度与垂直梯度具体包括: 2a)设置水平梯度算子模板与垂直梯度算子模板分别用来检测图像的垂直边缘及水平边缘，模板内的数字为模板系数，水平梯度算子模板的方向与垂直边缘的方向相垂直，垂直梯度算子模板的方向与水平边缘的方向相垂直；
4.如权利要求1所述的SAR图像相干斑抑制方法，其中，
5.如权利要求1所述的SAR图像相干斑抑制方法，其中，
6.如权利要求1所述的SAR图像相干斑抑制方法，其中，步骤5求得的像素点X的相干斑抑制结果为
7.如权利要求1所述的SAR图像相干斑抑制方法，其中，对于给定的强度格式的SAR图像，该SAR图像的相干斑方差为
8.如权利要求1所述的SAR图像相干斑抑制方法，其中，迭代次数为[8，10]次。
9.如权利要求1所述的SAR图像相干斑抑制方法，其中，h= 4.3。
【文档编号】G06T5/00GK103871031SQ201410081817
【公开日】2014年6月18日 申请日期:2014年3月7日 优先权日:2014年3月7日
【发明者】王爽, 焦李成, 于佳平, 霍丽娜, 侯彪, 马文萍, 马晶晶 申请人:西安电子科技大学