基于海面场景感知的雷达目标检测方法、系统及存储介质

1.本发明涉及雷达检测技术领域，尤其涉及一种基于海面场景感知的雷达目标检测方法、系统及存储介质。


背景技术：

2.高分辨预警雷达海杂波表现为显著的“三非”特性，即非高斯特性、时域非平稳特性和空域非均匀性。除岛礁地杂波带来的干扰外，海洋气候和地理环境差异也会导致海杂波在空域上的非均匀性。国内外学者采用混沌特性理论、分形理论、时频分析和恒虚警率(constant false alarm rate,cfar)等理论检测海面目标，已经取得了较大进展，这些方法往往是通过原始雷达回波在子空间上的特征值划分门限来检测目标。然而，复杂海面场景的雷达目标检测仍面临以下难题：1)混杂于海杂波中的岛礁地杂波干扰因其较大幅度，产生大量虚警；2)海杂波各区域起伏程度差异较大，杂波强弱交界处参考单元内的非独立同分布数据将导致参数估计准确度降低，导致虚警增加甚至弱目标漏检。
3.为了应对复杂场景下的非均匀回波带来的性能下降问题，认知雷达利用目标和环境的先验知识，优化发射波束、波形以及接收处理过程。有文献利用检测环境的先验信息，结合杂波分布参数的nlr(nonlocal low-rank)正则化表示，来降低非均匀杂波背景下的虚警率。但先验知识获取难度大，且其准确性差或使用不当时，甚至会造成性能恶化。同时学者们对地、海杂波的分类方法也进行了大量研究，如相位相关度、回波相似度、直流分量占比、能量、相关系数等。准确地分离岛海杂波，有助于提取出纯海面进而专注于海域内目标检测，降低强地杂波干扰。因此，自主感知海面杂波场景信息是提升雷达目标检测性能的必然途径。此外，有文献通过估计共形阵列的四维最小方差无失真响应谱，对近程非平稳杂波区域的杂波数据进行补偿。还有文献通过构建ard(range-doppler processed，等距离多普勒)单元分布参数的稀疏二维曲面来反映待测单元的空间信息，从而避免非均匀背景下的参考单元选择问题。还有文献则通过均分俯视角进行区域划分，却未考虑实际空间环境变化对回波的影响，预设的分区数目对分区效果也有一定影响。此外基于相关时间的小场景分区法常被用于提高协方差估计准确性，空变k分布模型等方法也被用于场景分区过程中。然而，这些仅适用于空间范围较小的场景，而对于距离范围广而观测时间较短的回波数据来说性能较差。


技术实现要素：

4.本发明提供基于海面场景感知的雷达目标检测方法、系统及存储介质，解决的技术问题在于：如何有效检测海面场景中的目标，而且大幅度提升非平稳海杂波下弱目标检测能力。
5.为解决以上技术问题，本发明提供基于海面场景感知的雷达目标检测方法，其关键在于，包括步骤：
6.s1、接收雷达回波数据；
7.s2、依据空域回波幅度均值将所述雷达回波数据分离为区域间非平稳、区域内平稳的多个回波区域；
8.s3、对每个所述回波区域进行杂波识别，分辨出不同波动程度的岛屿区域和海杂波区域；
9.s4、确定恒虚警门限，并根据该恒虚警门限计算每个所述海杂波区域的幅度门限；
10.s5、将幅度值大于其对应的幅度门限的距离门判定为可能存在目标的可疑距离门；
11.s6、采用检测器计算每个所述可疑距离门的子空间特征，并重新计算该检测器的特征门限值；
12.s7、将子空间特征大于该特征门限值的可疑距离门判定为该处存在目标的目标距离门。
13.进一步地，所述步骤s2具体包括步骤：
14.s21、求取m
×
n维的雷达回波数据r(m,n)在每个距离门上的幅度均值得到x(m)，而第m个距离门的幅度均值x(m)由下式计算：
[0015][0016]
r(m,i)表示第m个距离门的第i个脉冲的回波数据，m为距离门数，n为脉冲维数；
[0017]
s22、基于衡量相邻区域的幅度均值是否一致的检验统计量对x(m)进行粗划分，得到多个第一回波区域；
[0018]
s23、基于检验统计量显著性值对每个所述第一回波区域进行细划分，得到多个第二回波区域。
[0019]
进一步地，步骤s22具体包括步骤：
[0020]
s221、确定起始区域长度l0；
[0021]
s222、每次增加m0个距离门作为新增区域，计算衡量起始区域与新增区域的幅度均值是否一致的检验统计量t(m)；
[0022]
s223、若t(m)小于其门限t0，则将起始区域和新增区域合并为新的起始区域，返回至步骤s222进行重新计算，直至t(m)≥t0；若t(m)≥t0，则当前区域划分结束，确定起始区域的前端边界并返回至步骤s221进行迭代，直至遍历完每个距离门，得到多个第一回波区域。
[0023]
进一步地，所述步骤s23具体包括步骤：
[0024]
s231、遍历每个所述第一回波区域，计算每个距离门的t(m)值并确定其t(m)的最大值t
max
，
[0025]
s232、计算每个所述第一回波区域的检验统计量显著性值p(t
max
)即该第一回波区域不会被分割的概率；
[0026]
s233、p(t
max
)大于其门限p0时，在取得t
max
的距离门处将该第一回波区域分割为两段，获得两段新的子回波区域，返回至步骤s221进行迭代，直至p(t
max
)≤p0时停止分割，得到多个子回波区域。
[0027]
进一步地，第m个距离门的检验统计量t(m)由下式计算：
[0028][0029]
其中，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的平均值，sd(m)为第m个距离门处的合并偏差；sd(m)由下式计算：
[0030][0031]
其中，m1、m2分别为第m个距离门的左、右子序列的距离门个数，s1(m)、s2(m)分别为第m个距离门的左、右子序列的幅度平均值，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的标准差；
[0032]
p(t
max
)由下式计算：
[0033][0034]
其中，参数γ＝4.19lnm-11.54，参数δ＝0.40，参数λ＝m-2，i
x
(a,b)为非完全β函数。
[0035]
进一步地，所述步骤s3具体包括步骤：
[0036]
s31、计算每个所述回波区域中每个距离门的直流分量占比，其中第m个距离门的直流分量占比w(m)由下式计算：
[0037][0038]
其中，r(m,i)表示第m个距离门的第i个脉冲的回波数据；
[0039]
s32、将直流分量占比w(m)大于其设定阈值的距离门所对应的回波区域判定为岛屿区域，其他区域则为海杂波区域。
[0040]
进一步地，在步骤s4中，每个所述海杂波区域的幅度门限t根据下式计算：
[0041][0042]
其中pf为恒虚警门限，其中a为k-分布的尺度参数，v为k-分布的形状参数，γ(
·
)为伽马函数，k
v-1
(x/a)为v-1阶第二类修正贝塞尔函数，x表示x(m)。
[0043]
进一步地，在步骤s6中，检测器的特征门限值t
fea
由下式重新计算：
[0044][0045]
其中，t
fea_0
为检测器的原特征门限值，k为用于调节形状参数v和尺度参数a对于t
fea
影响大小的常数。
[0046]
本发明还提供一种基于海面场景感知的雷达目标检测系统，其关键在于，包括信
号接收模块、区域划分模块、杂波识别模块、幅度门限计算模块、可疑目标判定模块、检测模块和目标判定模块，分别用于上述方法中所述的步骤s1、s2、s3、s4、s5、s6、s7。
[0047]
本发明还提供一种存储介质，其关键在于：其上存储有计算机程序，执行该计算机程序可实现上述方法中所述的步骤s1、s2、s3、s4、s5、s6、s7中的至少一个步骤。
[0048]
本发明提供的基于海面场景感知的雷达目标检测方法、系统及存储介质，首先，依据回波幅度均值将整个距离门分离为多个区域间非平稳、区域内平稳的回波区域；之后，感知各平稳区域回波波动程度，并结合地杂波直流分量高的特点，识别出岛屿区域和海杂波区域；然后，对纯海杂波区域进行k分布拟合计算出其幅度门限；之后，根据幅度门限提取可疑目标的信号特征；最后，采用检测器进行最后的检测。基于雷达学报公开数据测试试验表明，本发明不仅能有效识别海面场景中的目标，而且能大幅度提升非平稳海杂波下弱目标检测能力。
附图说明
[0049]
图1是本发明实施例提供的基于海面场景感知的雷达目标检测方法的流程图；
[0050]
图2是本发明实施例提供的0614号数据部分脉冲的空间幅度均值分布图；
[0051]
图3是本发明实施例提供的0614号数据部分脉冲的分区效果图；
[0052]
图4是本发明实施例提供的0614号数据部分脉冲的场景感知结果图；
[0053]
图5是本发明实施例提供的实测数据检测的roc曲线图；
[0054]
图6是本发明实施例提供的pf＝0.01时检测率曲线图；
[0055]
图7是本发明实施例提供的纯海杂波数据roc曲线图；
[0056]
图8是本发明实施例提供的pf＝0.01时纯海杂波数据检测率曲线图。
具体实施方式
[0057]
下面结合附图具体阐明本发明的实施方式，实施例的给出仅仅是为了说明目的，并不能理解为对本发明的限定，包括附图仅供参考和说明使用，不构成对本发明专利保护范围的限制，因为在不脱离本发明精神和范围基础上，可以对本发明进行许多改变。
[0058]
本发明实施例首先提供一种基于海面场景感知的雷达目标检测方法，如图1所示，包括步骤：
[0059]
s1、接收雷达回波数据；
[0060]
s2、依据空域回波幅度均值将所述雷达回波数据分离为区域间非平稳、区域内平稳的多个回波区域；
[0061]
s3、对每个所述回波区域进行杂波识别，分辨出不同波动程度的岛屿区域和海杂波区域；
[0062]
s4、确定恒虚警门限，并根据该恒虚警门限计算每个所述海杂波区域的幅度门限；
[0063]
s5、将幅度值大于其对应的幅度门限的距离门判定为可能存在目标的可疑距离门；
[0064]
s6、采用检测器计算每个所述可疑距离门的子空间特征，并重新计算该检测器的特征门限值；
[0065]
s7、将子空间特征大于该特征门限值的可疑距离门判定为该处存在目标的目标距
离门。
[0066]
下面对步骤进行详细说明。
[0067]
(1)步骤s1
[0068]
s1、接收雷达回波数据。
[0069]
本例假设一个cpi(相参处理时间)的雷达回波数据r(m,n)大小为m
×
n，m为距离门数，n为脉冲维数。
[0070]
(2)步骤s2
[0071]
s2、依据空域回波幅度均值将所述雷达回波数据分离为区域间非平稳、区域内平稳的多个回波区域。
[0072]
该步骤s2具体包括步骤：
[0073]
s21、求取m
×
n维的雷达回波数据r(m,n)在每个距离门上的幅度均值得到x(m)，而第m个距离门的幅度均值x(m)由下式计算：
[0074][0075]
r(m,i)表示第m个距离门的第i个脉冲的回波数据，m为距离门数，n为脉冲维数。
[0076]
s22、基于衡量相邻区域的幅度均值是否一致的检验统计量对x(m)进行粗划分，得到多个第一回波区域；
[0077]
s23、基于检验统计量显著性值对每个所述第一回波区域进行细划分，得到多个第二回波区域。
[0078]
图2为雷达学报公布数据集0614号数据部分脉冲的空间幅度均值分布。在第1100距离门以内的回波波动较大且幅度均值大，这是由于靠岸海水更浅，水流撞击海底礁石增加了海杂波混乱程度；在第2050到第2800距离门，回波幅度较大且相邻距离门的均值变化不大，位于岛礁区；其它距离门的海杂波幅度均值较小，且变化相对不算明显。由此可见，海杂波幅度均值随空间距离变化，难以服从某一种既定的分布规律。
[0079]
事实上，复杂场景的非平稳海杂波特性参数不仅随着距离与方位变化，而且在时间尺度上也有一定变化。与之相比，特性参数随着距离和方位变化更为明显。海洋表面气象条件的变化一般以小时为单位，在一个相参处理时间(coherent processing interval,cpi)内海杂波变化可以忽略不计。因此，空间幅度均值体现了雷达回波在空间的分布特性。
[0080]
此外，海尖峰持续时间短，使用幅度均值作为区域分离依据，可以有效避免海尖峰在短时内带来的影响，体现“背景”的概念，即回波幅度在空间上的不同分布情况对应不同的杂波状态。因此，本例选择回波幅度均值作为场景区域分离的依据。
[0081]
对于步骤s22，其具体包括步骤：
[0082]
s221、确定起始区域长度l0；
[0083]
s222、每次增加m0个距离门作为新增区域，计算衡量起始区域与新增区域的幅度均值是否一致的检验统计量t(m)；
[0084]
s223、若t(m)小于其门限t0，则将起始区域和新增区域合并为新的起始区域，返回至步骤s222进行重新计算，直至t(m)≥t0；若t(m)≥t0，则当前区域划分结束，确定起始区域的前端边界并返回至步骤s221进行迭代，直至遍历完每个距离门，得到多个第一回波区域。
[0085]
对于步骤s23，其具体包括步骤：
[0086]
s231、遍历每个所述第一回波区域，计算每个距离门的t(m)值并确定其t(m)的最大值t
max
，
[0087]
s232、计算每个所述第一回波区域的检验统计量显著性值p(t
max
)即该第一回波区域不会被分割的概率；
[0088]
s233、p(t
max
)大于其门限p0时，在取得t
max
的距离门处将该第一回波区域分割为两段，获得两段新的子回波区域，返回至步骤s221进行迭代，直至p(t
max
)≤p0时停止分割，得到多个子回波区域。
[0089]
第m个距离门的检验统计量t(m)具体由下式计算：
[0090][0091]
其中，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的平均值，sd(m)为第m个距离门处的合并偏差；sd(m)由下式计算：
[0092][0093]
其中，m1、m2分别为第m个距离门的左、右子序列的距离门个数，s1(m)、s2(m)分别为第m个距离门的左、右子序列的幅度平均值，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的标准差。
[0094]
当t(m)越大时，表明左右子序列的幅值差异越大。对空间幅度均值x(m)的每个距离门分别求得相应的检验统计量t(m)，并计算t(m)的最大值t
max
所对应的显著性值p(t
max
)。
[0095]
p(t
max
)由下式计算：
[0096][0097]
其中，参数γ＝4.19lnm-11.54，参数δ＝0.40，参数λ＝m-2，i
x
(a,b)为非完全β函数。
[0098]
经过上述粗划分和细划分后，雷达回波数据在空间上将会被分为若干区域，每个相邻区域的幅度均值具有明显差异，而区域内的幅度具有一定的平稳性。图3是雷达学报0614号数据的回波均值分区示意图，其中红色折线代表了均值大小。整个场景4346个距离门被分为了7个区域，每个区域内回波幅度围绕在该区域红线附近，且其与相邻区域的红线幅度具有一定差异。此外，本算法在幅度均值变化较为明显的距离门处的分界点划分非常准确，如第2040距离门处。
[0099]
(3)步骤s3
[0100]
s3、对每个所述回波区域进行杂波识别，分辨出不同波动程度的岛屿区域和海杂波区域。
[0101]
海面回波会因为潮汐、风速、洋流等环境因素产生随机变化，从而在相邻脉冲时间内产生随机相位变化，多普勒维具有较大的频谱宽度。相反，岛礁回波在相邻脉冲内相位变化较小，其多普勒能量集中在零频附近，表现为直流分量占比大。因此可以通过直流分量占
比来区分陆地和海杂波区域，计算公式如下：
[0102][0103]
其中r(m,i)表示第m个距离门的第i个脉冲的回波数据，s(m)和p(m)分别为第m个距离门的直流分量与平均功率。
[0104]
步骤s3具体包括步骤：
[0105]
s31、计算每个所述回波区域中每个距离门的直流分量占比，其中第m个距离门的直流分量占比表示为w(m)；
[0106]
s32、将直流分量占比w(m)大于其设定阈值的距离门所对应的回波区域判定为岛屿区域，其他区域则为海杂波区域。
[0107]
(4)步骤s4
[0108]
s4、确定恒虚警门限，并根据该恒虚警门限计算每个所述海杂波区域的幅度门限。
[0109]
在场景分离并识别岛屿区域后，再对非岛礁区域的海杂波提取其波动程度参数，以此指导检测器选取阈值。考虑到纯海杂波幅度主要服从k分布，本例则采用该分布对海杂波暂平稳区域进行拟合，提取拟合参数来衡量杂波起伏程度。各平稳分区分别获取对应的海杂波波动参数，可以有效避免参考窗口内非均匀数据对参数估计的影响。
[0110]
用x表示x(m)，k-分布的概率密度函数(pdf)为：
[0111][0112]
其中a为k-分布的尺度参数，由实际杂波的平均功率决定；v为形状参数，描述k-分布的倾斜度，决定分布形状。γ(
·
)为伽马函数，k
v-1
(x/a)为v-1阶第二类修正贝塞尔函数。当形状参数v不变时，随着尺度参数a的增加，pdf的波峰逐渐随着幅度变大的方向移动，拖尾也逐渐增大。随着a的增大，信号波峰右移，构成信号的主要幅度值增大，且信号的分布范围更广，即尺度参数更大的数据具有更大的波动性。同样地，尺度参数a一定时，形状参数v也具有同样的规律。
[0113]
在场景感知的基础上，对每个海杂波区域采取k-分布恒虚警门限，从幅度上判断是否存在目标。本例将采用xln(x)的一、二阶矩估计k分布参数。对每个海杂波区域内的x(m)归一化后，令z＝ln(x)，则z的一阶原点矩为：
[0114][0115]
其中ψ(
·
)为psi函数，ω为常数。z的二阶原点矩为：
[0116][0117]
将(1)代入(2)中消除ω，得到形状参数的估计值的表达式：
[0118][0119]
其中ψ(1,
·
)为一阶psi函数。再通过式(3)得到尺度参数的估计值的表达式：
[0120][0121]
接着根据恒虚警门限pf和幅度门限的关系式，确定该海杂波区域的幅度门限t。
[0122][0123]
(5)步骤s5
[0124]
s5、将幅度值大于其对应的幅度门限的距离门判定为可能存在目标的可疑距离门。
[0125]
每个海杂波区域内根据恒虚警门限进行判决，大于门限的距离门被判定为可疑距离门(可疑目标)。相比k分布的ca-cfar来说，本例提出的方法避免了因参考单元内数据非均匀而引起的参数估计不准确问题。
[0126]
(6)步骤s6
[0127]
s6、采用检测器计算每个所述可疑距离门的子空间特征，并重新计算该检测器的特征门限值。
[0128]
对通过幅度检测的可疑目标进一步提取该距离门回波数据多特征来确认目标。本例采用当前应用较广的anmf(adaptive normalized matched filter)，sm(s-method)，rft(radon-fourier transform)和sact(strong auto-correlation time)检测器。
[0129]
与此同时，相应区域的幅度波动程度将会用于调节检测方法在该区域的特征值门限t
fea
，从而自适应调整门限，以避免目标在波动较强的区域被淹没。衡量某区的波动程度的参数为k分布的形状参数v和尺度参数a，其检测器特征门限t
fea
为：
[0130][0131]
其中t
fea_0
为anmf等检测器的原特征门限值；k为常数，用于调节形状参数v和尺度参数a对于门限的影响大小。
[0132]
(7)步骤s7
[0133]
s7、将子空间特征大于该特征门限值的可疑距离门判定为该处存在目标的目标距离门(目标位于该处)。
[0134]
本实施例还提供了一种基于海面场景感知的雷达目标检测系统，包括信号接收模块、区域划分模块、杂波识别模块、幅度门限计算模块、可疑目标判定模块、检测模块和目标判定模块，分别用于执行上述方法中的步骤s1、s2、s3、s4、s5、s6、s7。
[0135]
本实施例还提供了一种存储介质，其上存储有计算机程序，执行该计算机程序可实现上述方法中的步骤s1、s2、s3、s4、s5、s6、s7中的至少一个步骤。
[0136]
下面进行试验及分析。
[0137]
本例试验采用的数据为雷达学报2020年公布的第1期x波段数据，如表1所示，雷达数据距离门数为4346。雷达架设地点为烟台第1海水浴场，架高约为80m，凝视观测模式。极
化方式为hh，脉宽为3us，距离分辨率为2.5m，采样起始距离为距离岸边0.345km。雷达观测场景同时包含海杂波和岛礁地杂波，以及船只强目标和浮标弱目标，弱目标信杂比甚至会低于-10。
[0138]
表1数据具体情况
[0139][0140]
本例分别对岛礁地杂波和海杂波估计感知能力。岛礁抑制率isr衡量岛礁地杂波感知性能，即
[0141][0142]
其中m
land
(1)和m
land
(2)分别表示感知到的岛礁两侧边缘距离门，m'
land
(1)和m'
land
(2)分别表示实际岛礁两侧边缘距离门。图4示出了感知到的岛屿位置和标记的岛屿位置，两者吻合程度(平均isr)达到83.94％。
[0143]
同时，引入海杂波抑制率sr和感知率pr分别来描述海杂波感知性能和可疑目标提取效果。
[0144][0145][0146]
其中，ms为未被标记为可疑目标位置的距离门数；tn(i)表示处理第i个数据时正确感知到的目标数，tn'(i)表示第i个数据实际的目标数量，l为检测数据总个数。经过测试，平均sr为98.65％，pr高达99.22％，在保证几乎不丢失目标的高感知率情况下，本例提出的方法可以获得较高的海杂波抑制能力。
[0147]
实验数据为20210106160919和20210106172511数据，其中数据20210106160919有强弱目标各一个，数据20210106172511则有2个强目标。对提取到的可疑目标进行自适应目标检测，试验结果如图5和图6所示。主要用到anmf、sm、sact和rft检测器。图5为实测数据检测的roc曲线，图6为相同虚警率pf下，检测率pd随着信噪比snr(signal noise ratio)的变化曲线。从试验结果可以看出，对于非均匀岛礁场景下的强目标检测，无论是在相同pf或者是相同snr的条件下，结合场景感知的目标检测算法的检测率均高于对应不采用场景感知的原检测算法。
[0148]
为进一步验证非平稳海杂波平稳化分区处理对弱小目标检测性能提升的影响，本
试验对删除岛礁回波的纯海杂波数据进行测试。选用的数据是航道浮标弱目标数据，即20210106150614。类似地，对比测试直接使用anmf、sm、sact和rft检测器，以及在场景感知基础上使用这四个检测器，试验结果如图7和图8。实验结果表明，排除岛礁杂波带来的虚警影响后，在非均匀纯海杂波背景下，本例算法相较对比算法仍然保持了更低的虚警率，且虚警率相同时获得了更高的检测率。
[0149]
综上，本发明实施例提供的基于海面场景感知的雷达目标检测方法、系统及存储介质，首先，依据回波幅度均值将整个距离门分离为多个区域间非平稳、区域内平稳的回波区域；之后，感知各平稳区域回波波动程度，并结合地杂波直流分量高的特点，识别出岛屿区域和海杂波区域；然后，对纯海杂波区域进行k分布拟合计算出其幅度门限；之后，根据幅度门限提取可疑目标的信号特征；最后，采用检测器进行最后的检测。基于雷达学报公开数据测试试验表明，本发明不仅能有效识别海面场景中的目标，而且能大幅度提升非平稳海杂波下弱目标检测能力，有效降低虚警率，有效解决了非平稳、非均匀雷达回波数据降低参数估计准确性，恶化弱小目标检测的难题。
[0150]
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

技术特征：
1.基于海面场景感知的雷达目标检测方法，其特征在于，包括步骤：s1、接收雷达回波数据；s2、依据空域回波幅度均值将所述雷达回波数据分离为区域间非平稳、区域内平稳的多个回波区域；s3、对每个所述回波区域进行杂波识别，分辨出不同波动程度的岛屿区域和海杂波区域；s4、确定恒虚警门限，并根据该恒虚警门限计算每个所述海杂波区域的幅度门限；s5、将幅度值大于其对应的幅度门限的距离门判定为可能存在目标的可疑距离门；s6、采用检测器计算每个所述可疑距离门的子空间特征，并重新计算该检测器的特征门限值；s7、将子空间特征大于该特征门限值的可疑距离门判定为该处存在目标的目标距离门。2.根据权利要求1所述的基于海面场景感知的雷达目标检测方法，其特征在于，所述步骤s2具体包括步骤：s21、求取m
×
n维的雷达回波数据r(m,n)在每个距离门上的幅度均值得到x(m)，而第m个距离门的幅度均值x(m)由下式计算：r(m,i)表示第m个距离门的第i个脉冲的回波数据，m为距离门数，n为脉冲维数；s22、基于衡量相邻区域的幅度均值是否一致的检验统计量对x(m)进行粗划分，得到多个第一回波区域；s23、基于检验统计量显著性值对每个所述第一回波区域进行细划分，得到多个第二回波区域。3.根据权利要求2所述的基于海面场景感知的雷达目标检测方法，其特征在于，步骤s22具体包括步骤：s221、确定起始区域长度l0；s222、每次增加m0个距离门作为新增区域，计算衡量起始区域与新增区域的幅度均值是否一致的检验统计量t(m)；s223、若t(m)小于其门限t0，则将起始区域和新增区域合并为新的起始区域，返回至步骤s222进行重新计算，直至t(m)≥t0；若t(m)≥t0，则当前区域划分结束，确定起始区域的前端边界并返回至步骤s221进行迭代，直至遍历完每个距离门，得到多个第一回波区域。4.根据权利要求3所述的基于海面场景感知的雷达目标检测方法，其特征在于，所述步骤s23具体包括步骤：s231、遍历每个所述第一回波区域，计算每个距离门的t(m)值并确定其t(m)的最大值t
max
，s232、计算每个所述第一回波区域的检验统计量显著性值p(t
max
)即该第一回波区域不会被分割的概率；s233、p(t
max
)大于其门限p0时，在取得t
max
的距离门处将该第一回波区域分割为两段，获
得两段新的子回波区域，返回至步骤s221进行迭代，直至p(t
max
)≤p0时停止分割，得到多个子回波区域。5.根据权利要求4所述的基于海面场景感知的雷达目标检测方法，其特征在于，第m个距离门的检验统计量t(m)由下式计算：其中，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的平均值，s
d
(m)为第m个距离门处的合并偏差；s
d
(m)由下式计算：其中，m1、m2分别为第m个距离门的左、右子序列的距离门个数，s1(m)、s2(m)分别为第m个距离门的左、右子序列的幅度平均值，μ1(m)、μ2(m)分别为第m个距离门的左、右子序列的标准差；p(t
max
)由下式计算：其中，参数γ＝4.19lnm-11.54，参数δ＝0.40，参数λ＝m-2，i
x
(a,b)为非完全β函数。6.根据权利要求1所述的基于海面场景感知的雷达目标检测方法，其特征在于，所述步骤s3具体包括步骤：s31、计算每个所述回波区域中每个距离门的直流分量占比，其中第m个距离门的直流分量占比w(m)由下式计算：其中，r(m,i)表示第m个距离门的第i个脉冲的回波数据；s32、将直流分量占比w(m)大于其设定阈值的距离门所对应的回波区域判定为岛屿区域，其他区域则为海杂波区域。7.根据权利要求1所述的基于海面场景感知的雷达目标检测方法，其特征在于，在步骤s4中，每个所述海杂波区域的幅度门限t根据下式计算：其中p
f
为恒虚警门限，其中a为k-分布的尺度参数，v为k-分布的形状参数，γ(
·
)为伽马函数，k
v-1
(x/a)为v-1阶第二类修正贝塞尔函数，x表示x(m)。8.根据权利要求7所述的基于海面场景感知的雷达目标检测方法，其特征在于，在步骤
s6中，检测器的特征门限值t
fea
由下式重新计算：其中，t
fea_0
为检测器的原特征门限值，k为用于调节形状参数v和尺度参数a对于t
fea
影响大小的常数。9.基于海面场景感知的雷达目标检测系统，其特征在于，包括信号接收模块、区域划分模块、杂波识别模块、幅度门限计算模块、可疑目标判定模块、检测模块和目标判定模块，分别用于执行权利要求1～8任意一项所述的步骤s1、s2、s3、s4、s5、s6、s7。10.存储介质，其特征在于：其上存储有计算机程序，执行该计算机程序可实现权利要求1～8任意一项所述的步骤s1、s2、s3、s4、s5、s6、s7中的至少一个步骤。

技术总结
本发明涉及雷达检测技术领域，具体公开了一种基于海面场景感知的雷达目标检测方法、系统及存储介质，首先，依据回波幅度均值将整个距离门分离为多个区域间非平稳、区域内平稳的回波区域；之后，感知各平稳区域回波波动程度，并结合地杂波直流分量高的特点，识别出岛屿区域和海杂波区域；然后，对纯海杂波区域进行K分布拟合计算出其幅度门限；之后，根据幅度门限提取可疑目标的信号特征；最后，采用检测器进行最后的检测。基于雷达学报公开数据测试试验表明，本发明不仅能有效识别海面场景中的目标，而且能大幅度提升非平稳海杂波下弱目标检测能力。测能力。测能力。


技术研发人员：邓稼麒 李正周
受保护的技术使用者：重庆大学
技术研发日：2023.03.17
技术公布日：2023/7/26