[发明专利]基于核回归的SAR图像相干斑抑制

说明书

技术领域

本发明属于遥感图像数据处理技术领域，具体地说是一种图像相干斑抑制方法，该方法可用于SAR图像的相干斑抑制。

背景技术

合成孔径雷达SAR是一种高分辨率成像雷达。SAR不受气候、昼夜的影响，具备全天时、全天候的特点，在对地观测系统中，装备高分辨率的合成孔径雷达的航空系统机动性好、响应速度快，在灾害应急处理方面优势显著。作为卫星数据的有效补充，雷达在地形测绘、灾害监测评估、环境与国土资源监测等领域具有广泛的应用价值。特别是在灾害应急处理方面，雷达系统能及时为灾害应急处理提供信息丰富、清晰准确的图像，对及时开展救援、进行受灾评估、减轻灾害损失具有重要意义。但由于SAR系统采用的是相干微波源进行成像，回波的振幅在相干叠加时会相互叠加或抵消，得到的总回波强度与子回波平均强度之间会存在一定的偏差，在SAR图像中这种随机偏差以强噪声的形式呈现，这种噪声通常被称为相干斑。因为原理的缺陷性，SAR图像中相干斑影响了后续的图像解译，因此如何抑制SAR图像中的相干斑，提高图像的解读能力以及获得更多的信息成为一个重要的问题。SAR图像相干斑抑制的目标是在去除图像中的相干斑及尽可能的保持图像的细节信息。

SAR图像的相干斑是一种复杂的乘性噪声模型，对于斑点噪声这种特殊的性质，在过去的近二十年中，人们提出了很多经典的SAR图像相干斑抑制方法。如Lee滤波，增强Lee滤波，Kuan滤波等。这些方法是用一个已定义好的滤波器窗来估计局域相干斑的方差，并进行滤波处理，其结果通常过分的平滑了图像的边缘细节信息，在一定程度上取得了较好的效果，但存在模糊了边缘和点目标的缺陷。除了在空域进行图像相干斑抑制外，1995年，美国学者Donoho把小波理论引入到图像去噪中，提出了小波软阈值方法。小波软阈值法开创了变换域进行图像去噪的先河，之后涌现了许多变换域去噪的优秀成果，包括多尺度变换的图像去噪方法。小波软阈值方法是一种非线性的算法，依然存在破坏图像细节信息的问题，对图像的辐射特性保持也不好。

核回归方法本质上是一种空域局部平均的去噪方法，不同的是它依赖于被处理数据本身来指定这个模型的结构，这个模型是基于泰勒局部展开，采用的是非参数化方法，模型建立在图像的数据特征上。核回归分析应用于图像去噪最早是由H.Takeda和S.Farsiu等人于2005年提出，随后又得到了进一步的发展。其滤波效果与目前非常优秀的非局部均值方法NLM和三维块匹配的联合滤波BM3D方法相差无几，但是以上核回归算法目前只应用于噪声模型为加性的自然图像中，不能用到噪声模型为乘性的SAR图像中。

发明内容

本发明的目的在于提出了一种基于核回归的SAR图像相干斑抑制方法，将性能优异的核回归去噪方法扩展到噪声模型为乘性的SAR图像中，实现对SAR图像充分滤波的同时保持点目标和边缘细节信息。

为实现上述目的，本发明技术方案包括如下步骤：

(1)读入一幅SAR图像，并进行边界镜像对称扩展，扩展半径为N，取图像中的一个像素点x，以该像素点为中心取(2×N+1)大小的邻域窗w;

(2)用索贝尔（Sobel）算子计算邻域窗内w每一个像素点的水平梯度与垂直梯度；

(3)根据得到的水平梯度与垂直梯度计算以像素点x为中心的邻域窗w自适应协方差矩阵C；

(4)依据步骤(4)得到的局部自适应协方差矩阵C，计算以像素点x为中心的邻域窗w平滑矩阵H；

(5)依据计算得到的平滑矩阵H，利用权重影响，求得该像素点步骤1中的像素点x的相干斑抑制结果；

(6)对SAR图像每一个像素点按步骤(1)—(5)进行处理，得到每一个像素点的相干斑抑制结果，从而得到所述SAR图像的相干斑抑制结果。

(7)将步骤6得到的结果，按照步骤(1)—步骤(6)进行迭代处理，得到最终的所述SAR图像的相干斑抑制结果。

本发明与现有技术相比具有以下优点：

1、本发明是基于SAR图像的乘性噪声模型而提出的，更加符合SAR图像的相干斑抑制，对同质区域的滤波更加平滑；

2、本发明充分利用了图像的结构信息，克服了增强Lee的滤波等方法无法表示图像大量丰富的局部结构；

3、本发明用索贝尔Sobel算子计算图像的初始水平梯度与垂直梯度，减少计算量，能够快速有效的处理大幅面的SAR图像；

附图说明

图1是本发明的流程图；

图2是本发明仿真使用的barbara自然图像模拟的SAR图像；

图3是本发明仿真使用的光学图像模拟的SAR图像；