一种高低轨双基地合成孔径雷达频域快速成像方法

1.本发明属于雷达技术领域，特别涉及一种合成孔径雷达成像技术。


背景技术：

2.合成孔径雷达(synthetic aperture radar，sar)能够实现全天时、全天候的成像观测能力，在地形测绘、资源勘测、洋流及水文观测、灾害监测、植被分析、军事侦察等领域发挥着愈来愈重要的作用。双基sar相对于单基sar而言，因其收发分置的构型特点，具有强隐蔽性、抗干扰强、生存能力高等优势，能够实现接收站的轻量化、小型化、低成本和多样化，并且具备前视、侧视、后视全方位成像能力。
3.高低轨构型的双基地sar，在具备常规双基sar优势的基础上，充分发挥了收发双站作用距离远的特性，具有波束覆盖范围广、重访周期短的优点，可以满足长期大范围观测的应用场景需求。在近年，这种特殊构型的双基sar受到国内外各高校及科研单位的重视和研究。
4.目前高低轨双基sar的研究主要针对最佳分辨性能的轨道设计，如文献“f y.wang,z.lu,z.suo,y.zhu,z.li,and q.zhang,
‘
optimal configuration of spaceborne bistatic sar with geo transmitter and leo receiver,’iet radar,sonar navigat.,vol.13,no.2,pp.229
–
235,feb.2019”，针对高低轨双基sar复杂的双基构型以及差异的轨道特性，推导了成像性能评估参数，并构建了多元目标函数，通过基于模拟退火算法的优化方法设计出了最优高低轨双基构型。文献“g.krieger and a.moreira,“spaceborne bi-and multistatic sar:potential and challenges,”ieee proc.radar,sonar navigat.,vol.153,no.3,pp.184
–
198,2006.”分析了星载双多基sar的成像性能，阐述了星载双多基sar的不同应用，以及探讨了星载双多基sar存在的问题。
5.目前针对高低轨双基sar的成像方法研究主要为时域成像方法，因其不受双基sar构型的影响，文献：y.wang,y.liu,z.li,z.suo,c.fang,and j.chen,“high-resolution wide-swath imaging of spaceborne multichannel bistatic sar with inclined geosynchronous illuminator,”ieee geosci.remote sens.lett.,vol.14,no.12,pp.2380
–
2384,dec.2017，提出了非等效相位补偿的加权后向投影算法，通过将各接收通道的子图像加权求和得到最终图像，实现无方位模糊、无相位失真的图像聚焦。但是此方法的算法复杂度为n3(其中，n为数据方阵的大小)，高于频域成像方法的复杂度，不能实现快速成像。为了适应高低轨双基sar高精度成像的应用需求，需要研究高低轨双基地sar频域成像方法。
6.然而针对高低轨双基地sar的频域成像算法研究，存在三个方面的问题：首先传统高低轨双基sar平台的距离模型精度较差，在滑动聚束模式下，孔径时间较长，传统斜视距离模型无法满足高精度成像要求；滑动聚束模式会引起多普勒质心的空变，导致方位频谱展宽，从而引起方位模糊；对于高低轨移变模式下，由于接收机和发射机的速度大小和方向不同，导致相同双基斜距和的目标具有不同的距离单元徙动(rcm)和多普勒参数，因而高低
轨双基sar回波具有二维空变性。


技术实现要素：

7.为解决解决高低轨双基地sar频域成像方法受限于精确频谱未知和现有双基地sar时域成像方法运算量大，从而无法实现精确聚焦的问题，本发明提出一种高低轨双基地合成孔径雷达频域快速成像方法，能够有效实现滑动聚束模式下高低轨双基sar高精度频域成像。
8.本发明采用的技术方案为：一种高低轨双基地合成孔径雷达频域快速成像方法，包括：
9.a1、基于圆轨假设，构建高、低轨距离历史的圆轨模型；
10.a2、引入“非停走停”时延误差，将圆轨模型进行转化；
11.a3、将转化后的圆轨模型等效为圆轨轨迹距离模型；
12.a4、根据步骤a3的圆轨轨迹距离模型，得到高低轨双基sar回波模型；
13.a5、基于步骤a4的高低轨双基sar回波模型，进行高低轨双基sar成像。
14.步骤a3所述圆轨轨迹距离模型表达式为：
[0015][0016]
其中，ωe、re、θe为等效参数，k
t4s
表示高低轨双基sar的距离历史的四阶泰勒展开系数，
△kerr
表示圆轨轨迹距离模型四阶项的误差项。
[0017]
步骤a4所述高低轨双基sar回波模型表达式为：
[0018][0019]
其中，t
ac
表示目标点的方位中心照射时刻，wr(g)表示距离窗，wa(g)表示方位窗，kr表示距离向调频斜率，f0表示载频。
[0020]
步骤a5的实现过程包括以下步骤：
[0021]
a51、对高低轨双基sar回波模型进行距离向傅里叶变换，得到s0(f
τ
,t)：
[0022][0023]
其中，f
τ
表示距离向频率变量。
[0024]
a52、根据高低轨双站位置信息与“非停走停”误差项
△
τd，得到高低轨双基sar的多普勒参数f
dc
：
[0025]
[0026]
对s0(f
τ
,t)进行去质心处理，将原始方位频谱搬移到方位零频位置，得到s1(f
τ
,t)：
[0027]
s1(f
τ
,t)＝s0(f
τ
,t)exp(-j2πf
dc
t)；
[0028]
a53、将s1(f
τ
,t)乘以去斜函数h1，得到s2(f
τ
,t)：
[0029]
s2(f
τ
,t)＝s1(f
τ
,t)h1[0030]
a54、将回波数据s2(f
τ
,t)按照vet
′
进行方位向重采样，得到s2(f
τ
,t
′
)：
[0031][0032]
其中，ve为参考点位置的卫星运动速度，t
′
为新的方位时间变量，
[0033]
a55、将s2(f
τ
,t
′
)方位向傅里叶变换，并乘以补偿相位得到信号s3(f
τ
,t
′
)；
[0034][0035]
a56、将s3(f
τ
,t
′
)方位向傅里叶变换，乘以卷积补偿相位得到信号s4(f
τ
,f
t
)；
[0036][0037]
其中，f
t
表示方位向频率。
[0038]
a57、将s4(f
τ
,f
t
)乘以相位h2，得到信号s5(f
τ
,f
t
)；
[0039]
s5(f
τ
,f
t
)＝s4(f
τ
,f
t
)h2[0040]
其中，
[0041]
a58、将
△rt0s
与
△
t1按照y与
△
r1进行泰勒展开，
△rt0s
＝p1y+q1△
r1，
△
t1＝p2y+q2△
r1，得到信号s6(f
τ
,f
t
)；
[0042][0043]
其中，p1、q1、p2、q2是
△rt0s
与
△
t1关于空间变量y、
△
r1一阶展开系数；
[0044]
a59、对s6(f
τ
,f
t
)进行二维stolt插值，得到信号s7(fy,f
△
r1
)，通过对s7(fy,f
△
r1
)进行二维逆傅里叶变换，实现最终成像；
[0045][0046]
其中，
[0047]
本发明的有益效果：本发明的方法具体采用圆形轨道模型来构建精确的高低轨双基sar距离模型，通过计算“非停走停”的时延误差，以及方位重采样方法，将距离模型等效变换为一站固定式双基斜视模型，从而有利于求解高低轨双基sar的二维频谱模型，为频域成像方法奠定理论基础。成像处理部分通过方位预滤波方法消除雷达波束调整引入的方位频谱混叠的问题，适用于滑动聚束模式，后通过ω-k方位实现精确聚焦成像；
[0048]
本发明的方法创新点在于摒弃了传统的星载sar匀速直线运动模型，提出了高精度的等效圆轨轨迹距离模型，结合非停走停假设和方位重采样技术，得到了精度高且利于二维频谱推导的高低轨双基地sar距离模型，并提出方位预滤波方法和改进的ω-k成像方法，实现了回波数据的精确聚焦；
[0049]
本发明的方法可以应用于地球遥感、区域监视等领域。
附图说明
[0050]
图1是本发明提供方法的流程框图。
[0051]
图2是本发明具体实施采用的双基sar构型图。
[0052]
图3是本发明具体实施例中采用的目标场景布置图。
[0053]
图4是本发明具体实施例中对图3中9个点目标进行成像的结果图。
具体实施方式
[0054]
为便于本领域技术人员理解本发明的技术内容，首先对以下技术术语进行解释：
[0055]
术语1：高低轨双基地sar(bisar)
[0056]
双基地sar是指系统发射站和接收站分置于不同平台上的sar系统，其中至少有一个平台为运动平台，在概念上属于双基地雷达。而高低轨双基sar指发射平台为高轨sar卫星，接收平台为低轨sar卫星的双基sar系统。
[0057]
本发明实施例的几何结构图如图2所示，下面就具体实施方式对本发明的方法作进一步的详细描述。
[0058]
实施例1
[0059]
如图1所示，本发明的方法包括以下步骤：
[0060]
步骤一：根据设置的高低轨卫星轨道参数生成高低轨运动轨迹，产生回波信号s0(τ,t)；
[0061]
将高低轨平台的运动参数转换为目标成像坐标系，xoy平面与地球表面相切，切点为成像场景中心点o，也是坐标原点，z轴与地心指向切点的连线方向一致，y轴为低轨sar卫星方位零时刻的速度方向，x轴与y轴垂直，满足右手定则。
[0062]
将方位时间向量记为：t＝{-pri
·
na/2,-pri
·
(na/2-1),
…
,pri
·
(na/2-1)}
t
，
pri为脉冲重复间隔，na为目标回波方位点数。
[0063]
将距离时间向量记为：τ＝{-1/fs
·
nr/2,-1/fs
·
(nr/2-1),
…
,1/fs
·
(nr/2-1)}
t
，fs为脉冲重复间隔，nr为目标回波距离点数。
[0064]
发射站与场景中心点的斜距历史记为r
t
(t)，方位零时刻距离为r
t0
，接收站与场景中心点的斜距历史记为rr(t)，方位零时刻距离为r
r0
，任意目标点的多普勒中心时刻记为t0。
[0065]
基于圆轨假设将高轨发射站的距离历史表示为如下形式，
[0066]rt
(t)＝r
t0
+k
t1
t+k
t2
t2+k
t3
t3+k
t4
t4[0067]
其中，r
t0
、k
t1
、k
t2
、k
t3
、k
t4
是r
t
(t)的各阶泰勒展开系数。
[0068]
基于圆轨假设将低轨接收站的距离历史表示为如下形式，
[0069][0070]
其中ωs为等效角速度，θs为等效斜视角，r为低轨接收站的等效半径。
[0071]
计算“非停走停”假设的高低轨双基sar总时延τd，
[0072]
τd＝τ1+τ2[0073]
其中，高轨sar在t时刻发射信号，经过发射站斜距r
t
(t)到达目标p0，对应延时为τ1，此时，接收站运动至rr(t+τ1)，经过时间τ2，接收站接收到信号，所以
[0074][0075]
其中，c为光速。
[0076]
建立“非停走停”时延方程，
[0077]rbi
(t)＝cτd＝r
t
(t)+rr(t+τd)
[0078]
首先将rr(t)对t展开到四阶，
[0079]rr
(t)≈r
r0
+k
r1
t+k
r2
t2+k
r3
t3+k
r4
t4[0080]
其中，r
r0
、k
r1
、k
r2
、k
r3
、k
r4
是rr(t)的各阶泰勒展开系数。
[0081]
则可以得到
[0082]rr
(t+τd)≈r
r0
+k
r1
(t+τd)+k
r2
(t+τd)2+k
r3
(t+τd)3+k
r4
(t+τd)4[0083]
由于k
r3
和k
r4
很小，可忽略τd对三阶和四阶的影响。于是可以得到“非停走停”时延方程
[0084][0085]
求解上述一元二次方程可得
[0086][0087]
将τd对t展开，从而可以求出“非停走停”假设的误差如下
[0088]
[0089]
其中，
△
k0、
△
k1、
△
k2、
△
k3、
△
k4是的各阶泰勒展开系数。
[0090]
因此，高低轨双基sar的距离历史表示为
[0091][0092]
其中，r
t0s
和k
tis
为发射站距离历史和停走停误差的各阶展开系数和。
[0093]
为了得到简洁的高低轨双基sar距离模型，便于后续成像算法推倒，将上式等效为如下距离模型
[0094][0095]
其中
△kerr
表示距离模型四阶项的误差项，ωe，re和θe为等效参数，可以通过将上式的一阶、二阶和三阶展开系数与等效前距离历史的各阶展开系数相等求得，其表达式如下
[0096][0097]
其中，
[0098]kb1
＝-sinθe·
ωere[0099][0100][0101]
因此，高低轨双基sar回波模型为
[0102]
[0103]
其中，t
ac
表示目标点的方位中心照射时刻，wr(g)表示距离窗，wa(g)表示方位窗，kr表示距离向调频斜率，f0表示载频。
[0104]
步骤二：距离向傅里叶变换，将回波信号变为s0(f
τ
,t)；
[0105][0106]
其中，f
τ
表示距离向频率。上式中忽略了距离以及方位窗函数，只考察回波相位。
[0107]
步骤三：方位去质心，乘以相位exp(-j2πf
dc
t)，得到信号s1(f
τ
,t)；
[0108]
根据高低轨双站位置信息，以及步骤一推导的“非停走停”误差项
△
τd，解算高低轨双基sar的多普勒参数：
[0109][0110]
对回波数据进行去质心处理，将原始方位频谱搬移到方位零频位置，得到信号s1(f
τ
,t)：
[0111]
s1(f
τ
,t)＝s0(f
τ
,t)exp(-j2πf
dc
t)
[0112]
步骤四：乘以去斜函数h1；
[0113]
s2(f
τ
,t)＝s1(f
τ
,t)h1[0114]
其中，r
rot
是滑动聚束模式天线转动中心斜距，vr是leo卫星的等效速度，vb是leo卫星的地面波束移动速度，
△kt1s
是k
t1s
沿方位向的一阶展开系数。
[0115][0116]
其中，r是目标点斜距。
[0117]
步骤五：将回波数据s2(f
τ
,t)按照进行方位向重采样，得到信号s2(f
τ
,t
′
)，ve＝ωer
eref
为参考点位置的卫星运动速度。t
′
为新的方位时间变量；
[0118][0119]
其中，
[0120]
步骤六：将s2(f
τ
,t
′
)方位向傅里叶变换，并乘以补偿相位得到信号s3(f
τ
,t
′
)；
[0121]
[0122]
步骤七：将s3(f
τ
,t
′
)方位向傅里叶变换，乘以卷积补偿相位得到信号s4(f
τ
,f
t
)；
[0123][0124]
其中，f
t
表示方位向频率。
[0125]
步骤八：参考函数匹配滤波，即乘以相位h2，得到信号s5(f
τ
,f
t
)；
[0126]
由驻定相位原理得到s4(f
τ
,f
t
)二维频谱表达式：
[0127][0128]
其中，φ
err4
是四阶补偿因子
[0129][0130]
其中，t
′
posp
为驻定相位点
[0131][0132]
参考点补偿相位表达式为
[0133][0134]
则s5(f
τ
,f
t
)表达式为：
[0135][0136]
其中，
[0137]
[0138][0139][0140]
步骤九：二维stolt插值，将
△rt0s
与
△
t1按照y与
△
r1进行泰勒展开，
△rt0s
＝p1y+q1△
r1，
△
t1＝p2y+q2△
r1，得到信号s6(f
τ
,f
t
)；
[0141]
p1、q1、p2、q2是
△rt0s
与
△
t1关于空间变量y、
△
r1一阶展开系数，计算方式如下：
[0142][0143][0144]
其中，
[0145][0146][0147][0148]
其中，
[0149][0150][0151][0152][0153]
对于q1，可以表示为
[0154][0155]
其中，可以通过与相同的方式求出，能够表示为
[0156][0157]
其中，与可以表示为
[0158][0159][0160]
其中，
[0161][0162]
其中，
[0163][0164][0165]
其中与数值较小，可以忽略。而和与发射站接收站各自的多普勒参数有关。以下列出与的计算方法，同理可以推出和
[0166]kt1
＝-λf
dct
[0167][0168]
其中，
[0169][0170][0171]
r2＝(r
s-rr)
·
(r
s-rr)
[0172]
上述公式里的r、v和a表示位置、速度和加速度矢量，下标t和s分别表示发射站和
接收站。
[0173]
p2、q2可以表示为
[0174][0175][0176]
上式中的可以通过与相同的方式计算得到。
[0177]
将
△rt0s
和
△
t1表达式代入s5(f
τ
,f
t
)可以得到以下表示：
[0178][0179]
按照与对距离和方位频率进行变标操作，即对回波数据进行二维插值，得到信号s7(fy,f
△
r1
)；
[0180][0181]
步骤十：对进行二维逆傅里叶变换，实现最终成像。
[0182]
实施例2
[0183]
本实施例基于的目标场景如图3所示，采用的仿真参数如表1所示，图4为图3中9个点目标进行成像的结果图。
[0184]
表1双基sar仿真参数
[0185][0186]
图4是对图3中9个点目标进行成像的结果，从图4中可以看出，本发明提供的方法可以很好的实现高低轨双基地sar成像处理，能够实现对高低轨双基地sar回波的精确聚焦。
[0187]
上述步骤描述了本发明的优选实例，显然本领域技术人员通过参考本发明的优选实例和附图可以对本发明做出各种修改和替换，这些修改和替换都应落入本发明的保护范围之内。本发明所涉及的数学符号均为本领域常用符号。
[0188]
通过本发明具体实施方式可以看出，本发明可以实现对高低轨双基sar的精确聚焦。采用圆形轨道模型来构建精确的高低轨双基sar距离模型，通过计算“非停走停”的时延误差，以及方位重采样方法，将距离模型等效变换为一站固定式双基斜视模型，从而有利于求解高低轨双基sar的二维频谱模型，为频域成像方法奠定理论基础。并且成像处理部分通过方位预滤波方法消除雷达波束调整引入的方位频谱混叠的问题，适用于滑动聚束模式，后通过ω-k方位实现精确聚焦成像。本发明的方法可以用于地球遥感、区域监视等领域。
[0189]
本发明的方法构建了有利于二维频谱推导的高精度高低轨双基sar距离模型，并推导了相应的成像处理方法，能够解决高低轨双基sar波束调整引入的方位频谱混叠的问题和实现精确聚焦。具体的：
[0190]
通过将传统匀速、匀加速双基斜视距离模型改进为圆轨模型，精确模拟了高低轨的轨道运动特性，并且引入了高低轨双基构型超远距长基线高速度导致的“非停走停”时延误差，从而能够长时间高精度地拟合实际距离历史。一方面，高精度保证了后续成像方法聚焦效果和成像景深；另一方面，通过将该双基距离模型等效为一站固定式斜视距离模型，有利于二维频谱以及后续成像方法推导。
[0191]
由于高低轨不等速基线长，成像过程通常采用滑动聚束模式，天线调整会引入多普勒质心的方位空变，扩展了整体成像区域的频谱宽度，导致方位混叠。因此本成像方法结合方位预滤波方法，解决回波方位频谱混叠的问题，从而可以使用二维频域方法对回波进
行处理。后续使用步骤八至步骤十的ω-k成像方法实现了回波信号的精确聚焦，先通过参考函数相乘实现了距离向的匹配滤波，回波信号的一致rcmc(range cell migration correction，距离徙动校正)，以及方位向信号的一致聚焦；再通过stolt插值操作去除了方位向回波信号的二维空变，实现了补余rcmc以及补偿了方位向信号的残余相位。本发明方法填补了高低轨双基sar频域处理方法的空缺，相较于高低轨双基sar的时域处理方法，具有更高效的处理速度。
[0192]
本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。