一种SAR图像超像素滑窗CFAR检测方法技术

本发明专利技术公开了一种SAR图像超像素滑窗CFAR检测方法，该方法设置待检测超像素和背景超像素构成超像素滑窗，采用自适应阈值对背景超像素进行杂波截断处理，剔除影响杂波建模精度的异质像素，采用伽玛分布的截断形式进行杂波参数估计，根据给定的虚警概率求解CFAR检测阈值，对待检测超像素内的像素点进行目标判别，实现基于截断伽玛杂波统计特性的CFAR检测。该方法可有效提升多目标干扰等复杂环境下的目标检测精度和检测实时性。标检测精度和检测实时性。标检测精度和检测实时性。

【技术实现步骤摘要】
一种SAR图像超像素滑窗CFAR检测方法


[0001]本专利技术属于合成孔径雷达SAR图像目标检测
，特别涉及一种基于截断伽玛杂波的SAR图像超像素滑窗CFAR检测方法。

技术介绍

[0002]作为一种主动微波传感器，合成孔径雷达SAR系统不受光照、天气等条件的限制，具有全天候、全天时的观测能力，利用SAR图像进行海面和地面目标检测是目前雷达探测领域的研究热点。
[0003]在现有的SAR图像目标检测方法中，双参数恒虚警CFAR检测方法应用广泛。该检测算法基于背景杂波服从高斯或对数正态分布的假设，设置由目标窗口、保护窗口和背景窗口三部分组成的滑动窗口，利用该滑动窗口遍历整幅SAR图像中的像素。传统CFAR方法中，保护窗口是为了防止复杂背景下目标的部分像素泄露到背景窗口、影响杂波参数估计的准确性，然而在港口等航运繁忙区域，仅采用保护窗口不能保证完全排除邻近目标像素、方位向模糊等异质像素的干扰。因此，双参数CFAR方法的缺点主要在于：多目标等复杂场景下背景杂波建模不准确，造成目标检测准确率下降；目标检测采用像素点滑动窗口实现，对每个像素点进行检测判决时，均需对背景杂波参数进行估计，导致算法速度偏慢。
[0004]目前针对多目标等复杂环境下检测率降低的问题，提出了基于样本筛选的CFAR检测方法，通过对迭代筛选出的杂波进行参数估计和统计建模，有效提升了目标检测精度，然而采用固定阈值进行杂波筛选，使得真实杂波样本也被丢弃，导致杂波参数估计精度降低，并且像素点滑窗迭代筛选的方式计算效率低。
专利技术内容
[0005]针对上述已有技术的不足，本专利技术提出了一种基于截断伽玛杂波的SAR图像超像素滑窗CFAR检测方法，该方法采用超像素作为CFAR检测的基本单元，构建包含待检测超像素和背景超像素的超像素滑窗，采用自适应阈值截断筛选背景超像素中的真实杂波样本，利用截断式伽玛分布进行精确的杂波建模，根据设定的虚警概率计算CFAR检测阈值得到目标检测结果。
[0006]一种基于截断伽玛杂波的SAR图像超像素滑窗CFAR检测方法，包括如下步骤：
[0007]步骤1：设定超像素尺寸s，采用简单线性迭代聚类SLIC方法对待检测的SAR图像进行超像素分割，得到超像素{S
m
},m＝1,2,
…
,M，其中M为超像素的个数；
[0008]步骤2：设置由待检测超像素和背景超像素组成的局部滑动窗口，通过背景超像素中的杂波统计概率直方图自适应地计算杂波截断阈值，去除泄露到背景超像素中的目标、方位模糊的异质像素点，保留真实的杂波样本；
[0009]优选地，超像素局部滑动窗口中的背景超像素采用如下方法确定：
[0010]以待检测超像素的质心为圆心，以超像素尺寸s为半径画圆，选取圆内覆盖的所有超像素作为背景超像素。
[0011]优选地，自适应的杂波截断阈值计算方法如下：
[0012]首先统计背景超像素中像素点灰度值的概率直方图P，P满足0≤P(i)≤1，其中H为灰度级数，P(i)表示第i个灰度级的概率值，搜索概率直方图P最大值所在的灰度级H
max
，计算杂波截断阈值为t＝(H
max
+H)/2H，设背景超像素中某一像素灰度值为I
B
，则截断规则为I
B
≤t，即将灰度值小于t的杂波保留认定为真实海杂波。
[0013]步骤3：采用截断伽玛分布对背景超像素中自适应阈值截断后的杂波进行建模，根据截断矩估计方法得到杂波参数，具体如下：
[0014]采用如下的伽玛分布对采用自适应阈值t截断后的杂波样本的灰度概率密度进行建模：
[0015][0016]式中，L和μ分别为伽玛分布的形状参数和均值参数，通过计算截断杂波样本的一阶矩和二阶矩建立如下的联立方程组：
[0017][0018]其中，截断杂波样本的一阶矩和二阶矩通过和近似表示，n为截断杂波样本数，以作为初始点迭代求解得到数值解其中
[0019]步骤4：根据估计方法得到的杂波参数，给定虚警概率，自适应地求解截断伽玛分布CFAR检测阈值，对待检测超像素中的待检测像素点进行目标判别，对所有M个待检测超像素中的像素点完成目标判别，从而实现SAR图像目标检测。
[0020]具体如下：
[0021]设第m个待检测超像素S
m
的杂波参数估计结果为给定虚警概率P
fa
，根据下式计算CFAR检测阈值T
m
：
[0022][0023]采用T
m
对待检测超像素S
m
中的K个像素点的灰度值{I
m1
,I
m2
,
…
,I
mk
,
…
,I
mK
}进行判决，则CFAR检测目标判决规则为I
mk
≥T
m
。
[0024]相比于现有技术，本专利技术的显著优点在于：
[0025]1.本专利技术方法能够自适应地计算截断阈值，减少固定阈值带来的杂波拟合偏差；
[0026]2.本专利技术方法采用伽玛分布对截断杂波统计建模，相比于截断对数正态分布，提升了杂波拟合精度；
[0027]3.本专利技术方法利用超像素滑窗实现目标检测，大大提高了CFAR检测的计算效率。
附图说明
[0028]图1为超像素滑动窗口结构图；
[0029]图2为本专利技术中的截断式伽玛模型和对比方法中的截断对数正态模型对截断筛选后的杂波灰度直方图拟合对比结果；
[0030]图3为待检测的原始SAR图像，其中方框标记的为目标真值；
[0031]图4为本专利技术方法与现有方法的目标检测结果对比图；
[0032]图5为本专利技术方法与现有方法的ROC曲线对比图。
具体实施方式
[0033]下面结合附图对本专利技术的实施步骤和效果作进一步的说明。
[0034]步骤1：设定超像素尺寸s，采用简单线性迭代聚类SLIC方法对待检测的SAR图像进行超像素分割，得到超像素{S
m
},m＝1,2,
…
,M，其中M为超像素的个数；
[0035]步骤2：如附图1所示，设置由待检测超像素和背景超像素组成的局部滑动窗口，通过背景超像素中的杂波统计直方图自适应地计算杂波截断阈值，去除泄露到背景超像素中的目标、方位模糊等异质像素点，保留真实的杂波样本。
[0036]具体地，超像素局部滑动窗口中的背景超像素采用如下方法确定：
[0037]以待检测超像素S
m
的质心C
m
,m＝1,2,
…
,M为圆心，以超像素尺寸s为半径画圆，选取圆内覆盖的所有超像素为背景超像素，具体地，选取规则为：所选的背景超像素集合为{S
k
}，满足并且k∈[1,M],k≠m，其中C
k
为超像素S本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种SAR图像超像素滑窗CFAR检测方法，其特征在于：包括如下步骤：步骤1：设定超像素尺寸s，采用简单线性迭代聚类SLIC方法对待检测的SAR图像进行超像素分割，得到超像素{S
m
},m＝1,2,
…
,M，其中M为超像素的个数；步骤2：设置由待检测超像素和背景超像素组成的局部滑动窗口，通过背景超像素中的杂波统计概率直方图自适应地计算杂波截断阈值，去除泄露到背景超像素中的目标、方位模糊的异质像素点，保留真实的杂波样本；步骤3：采用截断伽玛分布对背景超像素中自适应阈值截断后的杂波进行建模，根据截断矩估计方法得到杂波参数；步骤4：根据估计方法得到的杂波参数，给定虚警概率，自适应地求解截断伽玛分布CFAR检测阈值，对待检测超像素中的待检测像素点进行目标判别，对所有M个待检测超像素中的像素点完成目标判别，从而实现SAR图像目标检测。2.根据权利要求1所述的一种SAR图像超像素滑窗CFAR检测方法，其特征在于：步骤2中所述的背景超像素中的根据杂波统计直方图计算自适应杂波截断阈值t的计算方法如下：首先统计背景超像素中像素点灰度值的概率直方图P，P满足0≤P(i)≤1，其中H为灰度级数，P(i)表示第i个灰度级的概率值，搜索概率直方图P最大值所在的灰度级H
max
，计算杂波截断阈值为t＝(H
max
+H)/2H，设背景超像素中某一像素灰度值为I
B
，则截断规则为I
B
≤t，即将灰度值小于t的杂波保留认定为真实海杂波。3.根据权利要求1所述的一种SAR图像超像素滑窗CFAR检测方法，其特征在于：步骤3中所述的采...

【专利技术属性】
技术研发人员：李焘，吕俊阳，彭冬亮，陈志坤，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：