CN107358616B

CN107358616B - 基于各向异性形态学方向比率的sar图像边缘检测方法

Info

Publication number: CN107358616B
Application number: CN201710523319.6A
Authority: CN
Inventors: 水鹏朗; 聂源盈; 樊书辰; 许述文
Original assignee: Xidian University
Current assignee: Xidian University
Priority date: 2017-06-30
Filing date: 2017-06-30
Publication date: 2020-04-14
Anticipated expiration: 2037-06-30
Also published as: CN107358616A

Abstract

本发明公开了一种基于各向异性形态学方向比率的SAR图像边缘检测方法，主要解决现有技术边缘检测零星强散射点鲁棒性差的问题，其技术方案是：根据输入的SAR图像X设置高斯伽马形双窗的参数；逆时针旋转高斯伽马形双窗，得到多方向的高斯伽马形双窗，并对其进行加权中值滤波；利用加权中值计算多方向边缘响应并对其进行增强；根据边缘响应计算边缘强度映射；从边缘强度映射中抽取极大值点，作为候选边缘点；对候选边缘点进行双阈值门限决策，得到最终图像的边缘检测结果。本发明增强了SAR图像边缘检测对零星强散射点的鲁棒性，满足实时大场景的SAR图像边缘检测要求，可用于含零星强散射点的大场景下提取SAR图像边缘。

Description

基于各向异性形态学方向比率的SAR图像边缘检测方法

技术领域

本发明属于图像处理技术领域，具体涉及一种SAR图像边缘检测方法，可用于含零星强散射点的大场景下的图像分割。

背景技术

近年来，随着合成孔径雷达SAR传感器的不断发展，SAR图像处理逐渐成为SAR领域研究的热点，SAR图像背景下的边缘检测技术是SAR图像分割应用技术中一个至关重要的研究方向。与传统的光学图像相比，SAR图像的人眼视觉可读性较差，而且容易受到相干斑及诸如迎坡缩短等形变的影响，这使得SAR图像的信息处理非常困难，最突出的困难在于SAR图像难以被解译，截至目前大多数的数据还是靠专业人士进行人工解译。图像的分割和分类一直以来都是图像解译的关键步骤，图像分割技术有助于对图像中所含分布式目标区域的信息进行分析和解译，但SAR图像中存在强烈的乘性相干斑噪声，使传统的光学图像处理算法直接用于SAR图像处理时，无法达到预期效果。大量的SAR图像数据不能得到有效的解读与应用，严重制约了SAR技术的应用。

SAR图像处理通常是要从图像中提取一些有意义的特征，如区域分割、边缘提取和目标检测等。边缘检测技术在SAR图像分割研究中占据着重要位置，它位于系统的最底层，为其他模块所依赖，它标记的图像中场景及目标的轮廓表征了图像的基本特征，其性能在很大程度上影响着一系列处理的整体效果。目前SAR图像边缘检测主要的方法分为四类：基于单边缘模型的方法、基于多边缘模型的方法、多分辨率的方法以及其他方法。其中基于单边缘模型的检测方法，按检测方法的相似性又可分为以均值比ROA检测器为代表的基于局部均值差异的检测方法和以广义似然比GLR检测器为代表的假设检验方法。

ROA检测器主要是基于局部均值对比度的差异，在目标像素点两侧设置两个窗口，计算两个窗口内像素均值之比，该检测器克服了经典梯度边缘检测器对乘性噪声非常敏感的缺点，同时考虑到了边缘的方向性，是一种恒虚警边缘检测器。人们根据不同的应用需要发展出一些新的基于ROA的边缘检测器，如基于高斯伽马双窗的ROA检测器，它比传统的ROA检测器在边缘附近虚警少且边缘连接度好。但高分辨SAR图像中呈现出更多的强散射点，这些强散射点虽然提供了更为精细的目标结构的信息，却不利于ROA检测器检测图像边缘。由于强散射点的出现增加了SAR图像的不均匀度，其统计特性发生改变，使得ROA检测器在处理含强散射点的SAR图像时很难得到较好的检测结果。

发明内容

本发明的目的在于针对上述现有技术的不足，提出一种基于各向异性形态学方向比率的SAR图像边缘检测方法，以提高含强散射点SAR图像的处理能力，精准定位边缘方向，降低边缘检测的错检率。

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

(1)根据输入的SAR图像X设置像素(x,y)的高斯伽马形双窗的参数，得到像素(x,y)的高斯伽马形上窗W_U和下窗W_L；

(2)将像素(x,y)的高斯伽马形上窗W_U和下窗W_L沿中心线逆时针旋转角度θ_n，得到像素(x,y)多方向的高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)，其中θ_n为[0，2π)内取的16个方向的角度，即θ_n＝n*π/8(0≤n＜16)，n是整数；

(3)分别利用上述的上窗W_U(θ_n)和下窗W_L(θ_n)对像素(x,y)进行加权中值滤波：

3a)从输入的SAR图像X中得到每一个像素(x,y)的灰度值I；

3b)在每个方向上分别将像素(x,y)的高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)内所有像素的灰度值I按升序排列，像素的灰度值I对应权值的顺序也被重新排列；

3c)分别将高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)重新排列后的权值进行累加，当累计权值达到1/2时，该权值对应的像素灰度值就是高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)在这个方向上的加权中值m_U(θ_n)和m_L(θ_n)；

(4)根据上述的加权中值m_U(θ_n)和m_L(θ_n)，计算像素(x,y)的多方向边缘响应R(θ_n)：

(5)增强像素(x,y)的多方向边缘响应R(θ_n)：

5a)像素(x,y)处的高斯伽马形双窗函数在当前方向θ处的理想阶跃边缘响应，对其进行对数变换并归一化得到方向匹配滤波器h(Δθ)：

其中Δθ表示其他方向与当前方向θ的差值；

5b)将方向匹配滤波器h(Δθ)作用于像素(x,y)的多方向边缘响应R(θ_n)，利用连乘和比率的形式对像素(x,y)的多方向边缘响应R(θ_n)进行增强，得到增强后像素(x,y)的边缘响应R_e(θ)：

其中∏表示连乘形式；

(6)根据增强后像素(x,y)的边缘响应R_e(θ)得到像素(x,y)的边缘强度映射E并提取边缘方向映射D：

D＝argmax{R_e(θ)}

(7)从像素(x,y)的边缘强度映射E和边缘方向映射D中抽取极大值点，作为候选边缘点；

(8)对候选边缘点进行双阈值门限决策，得到最终图像的边缘检测结果。

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

1)本发明由于利用改进的加权中值滤波，增强了在大场景边缘检测中边缘检测器对于零星强散射点的鲁棒性，减少了出现细碎边缘的可能性；

2)本发明由于使用了高斯伽马形双窗进行加权中值滤波，增强了边缘的连续性与平滑性，减少了错检的边缘像素个数；

3)本发明由于对多方向边缘强度响应进行了增强，使得边缘方向对应的边缘强度映射更大，因此与一般形态学方法相比，本发明可以更加精准地确定边缘方向，从而有利于后续对单像素边缘的提取。

附图说明

图1为本发明的实现流程图；

图2为用本发明和现有基于高斯伽马形双窗的ROA边缘检测方法对仿真SAR图像进行边缘检测的对比图。

具体实施方式

下面结合附图对本发明的实施率和效果作进一步详细说明：

参照图1，本发明的实现步骤如下：

步骤1，根据输入的SAR图像X设置高斯伽马形双窗的参数。

高斯伽马形双窗是沿着中心线对称分布的一种双窗函数，分为上窗W_U和下窗W_L，其沿着中心线方向服从高斯分布，垂直于中心线方向服从伽马分布，其参数包括x方向的高斯参数σ_x，y方向的伽马参数α和β；

本实例中，在假定高斯伽马形双窗中心线在水平轴时，SAR图像X中的像素(x,y)的高斯伽马形双窗函数为：

其中Γ(α)表示伽马函数，σ_x用于调节窗口长度，β用于控制窗口宽度，α用于控制双窗之间的间距，且σ_x、α>1，β>0，取

β＝1.8，α＝3.5。

步骤2，逆时针旋转高斯伽马形双窗，得到多方向的高斯伽马形双窗。

将像素(x,y)的高斯伽马形上窗W_U和下窗W_L沿中心线逆时针旋转角度θ_n，得到像素(x,y)多方向的高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)：

其中θ_n为[0，2π)内取的16个方向，即θ_n＝n*π/8,0≤n＜16，n是整数。

步骤3，利用多方向的高斯伽马形双窗进行加权中值滤波，得到双窗的加权中值。

这里的滤波方法可采用均值滤波，中值滤波，加权中值滤波，本实例选择的是改进的加权中值滤波，其步骤如下：

(3.1)从输入的SAR图像X中得到每一个像素(x,y)的灰度值I；

(3.2)在每个方向上分别将像素(x,y)的高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)内所有像素的灰度值I按升序排列，像素的灰度值I对应权值的顺序也被重新排列；

(3.3)分别将高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)重新排列后的权值进行累加，当累计权值达到1/2时，该权值对应的像素灰度值就是高斯伽马形上窗W_U(θ_n)和下窗W_L(θ_n)在这个方向上的加权中值m_U(θ_n)和m_L(θ_n)。

步骤4，根据高斯伽马形双窗的加权中值，计算多方向边缘响应R(θ_n)。

步骤5，增强多方向边缘响应R(θ_n)得到增强后的边缘响应R_e(θ)。

(5.1)当像素(x,y)的高斯伽马形双窗函数旋转到方向θ_n时，对当前方向θ_n处的理想阶跃边缘响应进行对数变换并归一化得到方向匹配滤波器h(Δθ)：

其中Δθ表示其他方向与当前方向θ_n的差值；

(5.2)将方向匹配滤波器h(Δθ)作用于像素(x,y)的多方向边缘响应R(θ_n)，利用连乘和比率的形式对像素(x,y)的多方向边缘响应R(θ_n)进行增强，得到增强后像素(x,y)的边缘响应R_e(θ)：

其中∏表示连乘形式。

步骤6，根据边缘响应R_e(θ)计算边缘强度映射E并提取边缘方向映射D：

步骤7，从边缘强度映射E和边缘方向映射D中抽取极大值点，作为候选边缘点。

从像素(x,y)的边缘强度映射E和边缘方向映射D中提取极大值点，是在边缘方向映射D的垂直方向上，提取像素(x,y)的3×3邻域内边缘强度映射最大的像素。

步骤8，对候选边缘点进行双阈值门限决策，得到最终图像的边缘检测结果。

(8.1)设定双阈值门限，包括高门限H和低门限L；

(8.2)将候选边缘点的边缘强度映射E与高门限H进行比较，若E≥H，则认为该点为是边缘点，否则，执行(8.3)；

(8.3)将候选边缘点的边缘强度映射E与低门限L进行比较，若E≤L，则该点不是边缘点，若E＞L，则认为该点是弱边缘，否则，执行(8.4)；

(8.4)判断弱边缘是否和边缘点是八邻域连通的，若是，则认为该弱边缘是边缘点，否则，不是边缘点。

基于步骤1到步骤8，实现了基于各向异性形态学方向比率的SAR图像边缘检测方法。

下面结合仿真实验对本发明的效果做进一步说明。

1.仿真参数

仿真实验中采用的数据是6视仿真幅度格式SAR图像。

2.仿真实验内容

仿真实验中分别采用本发明方法和基于高斯伽马双窗的ROA边缘检测方法得到的SAR图像边缘检测结果，通过边缘检测结果图比较两种边缘检测方法的质量。

仿真实验，对6视幅度格式SAR图像人为加入15个零星强散射点且散射点的像素值在同质区像素值的前10％，分别采用本发明方法和基于高斯伽马双窗的ROA边缘检测方法对SAR图像进行边缘检测结果，结果如图2所示，其中：

图2(a)表示6视仿真幅度格式SAR图像；

图2(b)表示采用本发明得到的6视仿真幅度格式SAR图像的边缘检测结果；

图2(c)表示采用基于高斯伽马双窗的ROA边缘检测方法得到的6视仿真SAR图像边缘检测结果。

从图2中可以看出，采用本发明方法得到的边缘检测结果明显优于现有方法得到的边缘检测结果。

综上所述，本发明提出的基于各向异性形态学方向比率的SAR图像边缘检测方法，可以提高对零星强散射点的鲁棒性，更加精准地确定边缘方向，满足实际雷达系统的图像处理要求，有利于后续操作对于单像素边缘的提取。

Claims

1.一种基于各向异性形态学方向比率的SAR图像边缘检测方法，包括：

3a)从输入的SAR图像X中得到每一个像素(x,y)的灰度值I；

(5)增强像素(x,y)的多方向边缘响应R(θ_n)：

5a)当像素(x,y)的高斯伽马形双窗函数旋转到方向θ_n时，对当前方向θ_n处的理想阶跃边缘响应进行对数变换并归一化得到方向匹配滤波器h(Δθ)：

其中Δθ表示其他方向与当前方向θ_n的差值；

其中Π表示连乘形式；

D＝arg max{R_e(θ)}

2.如权利要求1所述的方法，其中步骤(1)中根据SAR图像X设置像素(x,y)处的高斯伽马形双窗的参数，包括x方向的高斯参数σ_x，y方向的伽马参数α和β，且

β＝1.8，α＝3.5。

3.如权利要求1所述的方法，其中步骤(7)中从像素(x,y)的边缘强度映射E和边缘方向映射D中提取极大值点，是在边缘方向映射D的垂直方向上，提取像素(x,y)的3×3邻域内边缘强度映射最大的像素。

4.如权利要求1所述的方法，其中步骤(8)中对候选边缘点进行双阈值门限决策，按如下步骤进行：

8a)设定双阈值门限，包括高门限H和低门限L；

8b)将候选边缘点的边缘强度映射E与高门限H进行比较，若E≥H，则认为该点为边缘点，否则，执行8c)；

8c)将候选边缘点的边缘强度映射E与低门限L进行比较，若E≤L，则该点不是边缘点，若E＞L，则认为该点是弱边缘，否则，执行8d)；

8d)判断弱边缘是否和边缘点是八邻域连通的，若是，则认为该弱边缘是边缘点，否则，不是边缘点。