基于SAR成像的相干干扰方法

基于sar成像的相干干扰方法
技术领域
1.本发明涉及一种相干干扰方法，具体是一种基于sar成像的相干干扰方法。


背景技术：

2.rd算法于1978年处理出第一幅机载sar数字图像，至今仍在广泛使用，它通过距离和方位上的频域操作，达到了高效的模块化处理要求，同时又具有了一维操作的简便性。该算法根据距离和方位上的大尺度时间差异，在两个一维操作之间使用距离徙动校正(rcmc)，对距离和方位进行了近似的分离处理，从而成像效果较好。现有的干扰方法从相干性角度可划分为相干干扰和非相干干扰，非相干干扰通常指的是噪声压制性干扰，相干干扰与非相干干扰噪声相比可以获得脉冲压缩处理增益，大大降低了干扰功率的需求，且更好地实现了对重要目标的保护，也让抗干扰的实现变得困难。传统相干干扰中的卷积相干干扰是利用卷积噪声的方法获得脉冲压缩增益，从而生成区域性干扰，但由于是卷积噪声生成的干扰信号，所以易被反卷积算法分离识别。通过移频或者转发式干扰的生成相干干扰会与真实目标信息更为相似，难以被分离识别。因此，对sar成像的相干干扰技术研究具有十分重大的意义。


技术实现要素：

3.本发明的目的在于提供一种基于sar成像的相干干扰方法，以克服现有干扰易被抗干扰技术识别和分离的现状。
4.实现本发明目的的技术解决方案为：一种基于sar成像的相干干扰方法，包括：
5.合成孔径雷达成像：雷达发射的是线性调频信号，对目标进行rd成像，需要对信号回波进行距离向压缩、距离徙动校正和方位向压缩，最后得到压缩数据及成像结果；
6.相干压制干扰：根据拟保护目标的位置大小、干扰位置、sar平台参数，调制一个lfm信号转发，从而在sar图像中形成一个面状压制区域；
7.欺骗干扰：根据拟保护目标的位置大小、干扰位置、sar平台参数，实现在距离向或方位向上的欺骗干扰，包括距离向上的密集假目标干扰和距离移频干扰，方位向上的转发式干扰。
8.第二方面，本发明提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现第一方面所述的方法的步骤。
9.第三方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现第一方面所述的方法的步骤。
10.第四方面，本发明提供一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现第一方面所述的方法的步骤。
11.与现有技术相比，本发明的显著优点为：相干压制干扰和移频干扰都是相干干扰技术，可以有效地降低对干扰机功率的需求；根据移频量和移频方式的不同，可以产生多种
不同的干扰效果，还可以与现存较多的运动调制干扰或卷积干扰配合使用解决欺骗目标单一的问题；本发明将不同的距离向和方位向上的干扰方法结合使用，可以达到区域型假目标欺骗干扰，扩展了干扰区域，更好地保护真实目标，大大提高了抗干扰的难度。
附图说明
12.图1为基于sar成像及干扰流程图。
13.图2为正侧视sar几何关系示意图。
14.图3为未施加干扰时目标成像示意图。
15.图4为二维失配压制干扰成像示意图。
16.图5为距离固定移频干扰示意图。
17.图6为距离随机移频干扰示意图。
18.图7为距离步进移频干扰示意图。
19.图8为距离分段移频干扰示意图。
20.图9为密集假目标干扰示意图。
21.图10为转发式干扰示意图。
具体实施方式
22.如图1所示，本发明提供一种基于sar成像的相干干扰方法，包括：
23.首先获得雷达原始数据，利用原始数据进行距离向压缩、距离徙动校正和方位向压缩，获得压缩数据和成像结果。
24.根据拟保护目标的位置大小、干扰位置、sar平台参数，可以调制一个lfm信号转发，从而在sar图像中形成一个面状压制区域，降低图像对比度，实现相干压制干扰，从而实现对目标的保护。
25.根据拟保护目标的位置大小、干扰位置、sar平台参数，可以实现在距离向或方位向上的欺骗干扰，包括距离向上的密集假目标干扰和距离移频干扰，方位向上的转发式干扰。
26.对于不同的欺骗干扰，在同一信干比情况下进行干扰效果的评估，主要评估指标包含：峰值旁瓣比、积分旁瓣比、结构相似度及相关系数。
27.对于同一干扰，在不同信干比情况下进行干扰效果的评估，并对各种干扰结果进行对比。
28.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
29.如图1所示，一种基于机动sar平台成像干扰方法，包括以下步骤：
30.步骤1，首先根据雷达信号，生成雷达原始数据，雷达发射的一般是线性调频信号，其表达式为：
[0031][0032]
式中：tr为距离快时间，fc为信号载频，t为lfm信号脉冲宽度，br为lfm信号带宽，kr为线性调频率，kr＝br/t，λ＝c/fc为载频对应的波长；
[0033]
回波信号经过去载频处理后为：
[0034][0035]
式中，ta为方位慢时间，r(ta)为目标与sar之间的慢时相关瞬时距离。
[0036]
对距离向做脉冲压缩的系统匹配函数为：
[0037]
sr(tr)＝exp(-jπkrt
r2
)
[0038]
完成距离向压缩处理后就要对方位向进行处理，首先要考虑距离徙动的影响，当正侧视工作时，只考虑距离弯曲，如图2，对最近距离rb的点目标p，其斜距与ta的关系为：
[0039][0040]
式中，v为载机速度，第二项为距离弯曲，它对合成孔径期间的包络移动影响可以忽略，但对回波相位的影响必须考虑。进行距离徙动校正后，对方位向做脉冲压缩的系统匹配函数为：
[0041]
sa(ta)＝exp(-jπkat
a2
)
[0042]
式中，多普勒调频率为
[0043][0044]
完成以上操作后，得到回波成像数据及成像结果。
[0045]
步骤2，sar发射lfm信号后，为实现相干压制干扰，干扰机对截获的sar发射信号进行调制后转发，转发信号为：
[0046][0047]
式中，t
l
为信号脉冲宽度，c为调频斜率，调制信号带宽为b
l
＝t
l
c，f(ta)为方向位慢时间ta为变量的调制信号中心频率，可表示为
[0048]
f(ta)＝m1 m2ta[0049]
其中，m1和m2为根据干扰距离可调配参数。
[0050]
假设干扰机位于真实目标附近，其距离忽略不计，且调制信号脉冲宽度与截获的sar信号脉冲宽度保持一致，则去载频后的干扰信号可表示为：
[0051][0052]
式中δt为延迟时间，主要包括调制干扰、转发干扰信号的时间，kj＝kr c为干扰信
号的调频斜率。第一个指数项主要与距离向压缩有关，第二、三个指数项主要与方位向压缩有关。根据脉冲压缩原理可知，干扰信号的调频斜率与sar发射信号的调频斜率不一致，从而导致脉冲压缩输出主瓣宽度展宽。同理，干扰信号与方位向脉冲压缩参考信号的调频斜率不一致，使得方位向脉冲压缩输出主瓣宽度展宽。干扰信号经成像处理，在sar图像中形成散焦的图像，从而形成压制区域，而二维脉压失配压制干扰正是在此基础上产生的。因此，通过调整调制信号l(tr,ta)的参数，可形成面积可调的压制干扰区域，实现对重要目标的保护。
[0053]
步骤3，基于lfm信号具有时频耦合特性，所以可以对sar进行移频欺骗干扰。
[0054]
距离固定移频干扰为：干扰机接收到雷达信号后，对信号载频进行固定的频移后转发，固定移频量为fd，此时sar接收机接收到的移频干扰信号形式为：
[0055][0056]
其中tr为距离向快时间，ta为方位向慢时间，δt为干扰机转发时延，t
l
为合成孔径时间，rj(ta)为雷达到干扰机的距离，kr为信号调频斜率。
[0057]
经雷达下变频后，得到信号为：
[0058][0059]
再经距离向和方位向压缩处理后得到：
[0060][0061]
式中表示二维信号相乘。ka为方位向调频率，t为信号脉冲宽度。
[0062]
当移频量为固定时，干扰机输出为一个距离向上相对干扰机偏移τd＝-2πfd/kr的虚假点目标。若去脉间相干干扰后，可在点偏移位置生成一条直线。
[0063]
距离随机移频干扰同固定移频干扰步骤类似，区别为其中的随机移频的干扰移频量为：
[0064]
fd(td)＝ξ(ta)
·
δ
·br
[0065]
其中ξ(ta)为方向位上的高斯随机变量，设置其范围为[1,na]，na为方向位采样点数，δ为常数，br为sar信号带宽。这种情况下的移频干扰可在sar图像上形成一定宽度的条带似噪干扰，起到保护分布式目标的作用。
[0066]
此时sar接收机处接收到的距离随机移频干扰信号可表示为：
[0067][0068]
距离随机移频干扰将在sar图像上形成沿方位向上一定宽度的条带似噪欺骗干扰，起到分布式保护目标的作用。
[0069]
距离步进移频干扰是一种移频量与方位向慢时间相关，且以固定移频量逐次递增或递减的干扰方法。
[0070]
移频量为f
step
＝n
·
δf
step
，其中δf
step
为固定的移频增量，为正时表示逐次递增，为负时表示逐次递减，且要求f
step
≤|br|，使得干扰信号形成有效干扰。
[0071]
假设t0为sar照射干扰机初始时刻，n1为n的初始值，则移频量与方位向慢时间的关系表达式为：
[0072]fstep
＝prf
·
δf
step
t
a-δf
step
(m-n1)
[0073]
其中prf为脉冲重复频率，ta为方位向慢时间，m＝t0·
prf。
[0074]
则接收到的距离随机移频干扰信号为：
[0075][0076]
距离步进移频干扰将在sar图像上形成沿距离向上一定宽度的条带欺骗干扰，其干扰效果更加密集。
[0077]
距离分段移频干扰是在合成孔径时间内的不同时段，做不同量的移频干扰，设分段数为n，分段移频量可表示为：
[0078][0079]
其中，ta为方向位慢时间，t
l
为合成孔径时间，，t
i-1
，ti为第i段起止时间和截止时间，fi为第i段移频量。
[0080]
距离移频干扰可在距离向上形成不同位移的虚假点目标，与固定移频干扰不同，可以自由选择分段数和移频量，形成需要的干扰效果。
[0081]
密集假目标干扰首先对原始脉冲时域压缩得到一个子脉冲，其表达式可以表示为：
[0082]
j1(t)＝aexp(jπμt2),μ＝n
·
kr,0≤t≤tr/n
[0083]
其中，a为子脉冲的幅度，压缩后调制频率为μ＝n
·
kr,时宽为t＝tr/n，tr为原脉冲时宽。
[0084]
将j1(t)复制n次，得到干扰信号表达式为：
[0085][0086]
转发式干扰是干扰机对接收到的sar信号按prf进行转发，sar接收到的干扰信号为：
[0087][0088]
其中tj为干扰信号到sar的时延，na为方位向的脉冲个数，t
ai
为方位向上的转发时延。可以使i从-na到na，以实现方位向上前后时延。
[0089]
步骤4，对于不同的欺骗干扰，将其信干比控制在同一情况下进行干扰效果的评估。基于图像质量评估的指标主要分为峰值旁瓣比、积分旁瓣比、结构相似度及相关系数。
[0090]
峰值旁瓣比：一般定义为点目标冲激响应的最高旁瓣峰值ps与主瓣峰值pm之比，一般用分贝表示。峰值旁瓣比的大小决定了强目标掩盖弱目标的能力。通常要求峰值旁瓣比小于一定值。
[0091]
距离向峰值旁瓣比pslrr：
[0092][0093]
式中p
sr
和p
mr
分别为距离向点目标冲激响应的主瓣峰值和最高旁瓣峰值；
[0094]
方位向峰值旁瓣比pslra：
[0095][0096]
式中p
sa
和p
ma
分别为方位向点目标冲激响应的主瓣峰值和最高旁瓣峰值。
[0097]
积分旁瓣比：一般定义为旁瓣能量与主瓣能量的比值，单位为分贝。它是图像局部对比度的衡量指标，即它定量的描述了一个局部较暗的区域被来自周围的明亮区域的能量泄露所淹没的程度。它也是通过测量点目标特性获得的。
[0098]
通常情况下，积分旁瓣比越小，图像质量越高。
[0099]
距离向积分旁瓣比lslrr：
[0100][0101]
式中e
sr
和e
mr
分别为距离向冲激响应旁瓣能量和主瓣能量；
[0102]
方位向积分旁瓣比lslra：
[0103][0104]
式中e
sa
和e
ma
分别为方位向冲激响应旁瓣能量和主瓣能量。
[0105]
结构相似度：其基本思想是，人类视觉系统的主要功能是从视觉区域提取图像的结构性信息，它比较两幅图像的相似度，而不是差值，获得图像的相似结构映像，从而得到
客观预测质量。ssim将图像分成亮度，对比度，结构，分别进行比较，然后加权乘积获得相似性，可以计算出干扰图像的质量，以此来评估sar的抗干扰能力，当ssim值越大时，图像越相似。其数学表达式如：
[0106]
ssim(x,y)＝[l(x,y)]
α
·
[c(x,y)]
β
·
[s(x,y)]
γ
[0107]
式中，x和y分别代表sar原图像和干扰之后的图像；
[0108][0109][0110][0111]
其中，l(x,y)、c(x,y)、s(x,y)分别表示亮度、对比度、结构因子的比较，可以用它们来度量sar图像受干扰的程度；其中，μ
x
和μy分别为原始图像与干扰图像中像素灰度均值，表征图像的亮度。σ
x
和σy分别为原始图像和干扰图像中像素灰度的标准差，表征图像的对比度，σ
xy
为是原始图像与干扰图像对应块的相关系数，表征结构信息的相似性；c1、c2、c3是避免分母为零或接近零的很小的正数。通常令α＝β＝γ＝1，方便结构相似度的计算。
[0112]
相关系数：根据随机信号中两随机变量相关系数的相关理论，定义sar干扰图像与无干扰图像的相关系数公式如下：
[0113][0114]
式中，f(xi,yi)和g(xi,yi)分别为同一场景sar无干扰图像和干扰图像的像素灰度值，而f和g分别为sar无干扰图像和干扰图像的灰度均值。干扰前后sar图像的相关系数ρ
fg
反映了干扰前后两幅图像之间的统计相关性，对于压制式干扰，相关系数ρ
fg
越大，表明受干扰后的sar图像与原始sar图像之间越相似，说明图像受干扰程度越小，对sar的干扰效果越不好；反之，对sar的干扰效果越好。
[0115]
步骤5，对于同一干扰，设置不一样的信干比来进行干扰评估。分析信干比大小对干扰效果的影响，并对各种干扰结果进行比对。
[0116]
下面结合实施例对本发明进行详细说明。
[0117]
实施例
[0118]
为了验证本发明方案的有效性，进行如下仿真实验。系统参数为：线性调频信号部分，载频fc＝1.5ghz，信号脉宽t＝1.5us，信号带宽br＝150mhz，调制频率kr＝b/t，距离向快时间采样频率fs＝450mhz,sar平台速度v＝100m/s，高度h＝1000m，选取距离向采样点数为1024，方位向采样点数为256。
[0119]
图3为未施加干扰时的目标rd算法成像。图4为施加二维失配压制干扰成像，二维失配压制干扰调制频率c＝-0.2kr，m1＝0.4fc，m2＝0，呈区域型压制干扰。图5为距离固定移
频干扰，固定移频量fd＝-0.3br，呈目标右拖欺骗干扰。图6为距离随机移频干扰，随机移频量fd(td)＝0.1
·
ξ(ta)
·br
，呈沿方位向上一定宽度的条带似噪欺骗干扰。图7为距离步进移频干扰，步进移频量δf
step
＝0.06br，呈沿距离向上一定宽度的条带欺骗干扰。图8为距离分段移频干扰，分段移频量呈距离向前后欺骗干扰。图9为密集假目标干扰，设置n＝10，并且在叠加后向频谱的同时叠加前向频谱，实现距离向上双向密集假目标。图10为转发式干扰，实现方位向上的欺骗干扰。
[0120]
表1为不同干扰在干信比20db时的效果评估。
[0121]
表1
[0122][0123]
可以看出在距离向上实现干扰和方位向上实现干扰的峰值旁瓣比和积分旁瓣比有明显的区别，距离向和方位向上的干扰掩盖能力越强，其对应的峰值旁瓣比就越小。距离向和方位向上的干扰对图像质量影响越大，其对应的积分旁瓣比就越大。由峰值旁瓣比可以判断干扰后信号的质量，由积分旁瓣比可以判断图像的质量。由相关系数和结构相似度可以看出欺骗干扰后的图像相似程度和结构相似性。
[0124]
表2为不同干扰在不同干信比时的相关系数效果评估。
[0125]
表2
[0126][0127]
可以看出，在假目标欺骗干扰时，随着干信比的改变，图像整体质量变化较小，相关系数变化也在一定范围中。
[0128]
表3为不同干扰在不同干信比时的结构相似度效果评估。
[0129]
表3
[0130][0131]
结合相关系数效果评估可以看出，欺骗性干扰在结构相似性和统计相关性上区分度较小，在不影响图像质量的情况下实现假目标干扰，符合欺骗性干扰特征。实验中测量结果与理论值一致，符合设计要求。
[0132]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0133]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。