一种SAR滑动聚束模式方位去斜方法[发明专利]

(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 202010874366.7
(22)申请日 2020.08.26
(71)申请人 中国科学院空天信息创新研究院
地址 100190 北京市海淀区北四环西路19
号
(72)发明人 黄丽佳　孟大地　李莹莹　胡玉新　
仇晓兰　韩冰　仲利华　尹迪　
(74)专利代理机构 中科专利商标代理有限责任
公司 11021
代理人 孙蕾
(51)Int.Cl.
G01S 13/90(2006.01)
(54)发明名称
一种SAR滑动聚束模式方位去斜方法
(57)摘要
一种SAR滑动聚束模式方位去斜方法。

方法
包括：根据成像时间内的每个SAR方位时刻的轨
道六根数、卫星姿态角以及天线扫描角，确定每
个方位时刻对应的多普勒频率；根据每个方位时
刻对应的多普勒频率，确定成像时间内的多普勒
历程；根据多普勒历程，确定成像时间内的相对
距离历程；根据相对距离历程，确定方位去斜函
数；将回波信号与方位去斜函数相乘，以对回波
信号进行方位去斜。

本发明提供的方位去斜方法
不依赖虚拟转动点，可以解决大方位扫描角滑动
聚束模式虚拟转动点模型失效的问题，同时可以
精确去除波束方位扫描引起的多普勒混叠。

权利要求书4页 说明书6页 附图3页CN 111983612 A 2020.11.24
C N 111983612
A
1.一种SAR滑动聚束模式方位去斜方法，其特征在于，包括：
根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角，确定每个方位时刻对应的多普勒频率；
根据所述每个SAR方位时刻对应的多普勒频率，确定所述成像时间内的多普勒历程；
根据所述多普勒历程，确定所述成像时间内的相对距离历程；
根据所述相对距离历程，确定方位去斜函数；
将回波信号与所述方位去斜函数相乘，以对所述回波信号进行方位去斜。

2.根据权利要求1所述方法，其特征在于，所述根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态以及天线扫描角，确定每个方位时刻对应的多普勒频率，包括以下步骤：
根据所述成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角，计算在地球惯性坐标系下，卫星到波束中心地表交点的单位位置矢量、地心到卫星的位置矢量以及卫星的速度矢量；
根据所述单位位置矢量以及所述地心到卫星的位置矢量，计算在所述地球惯性坐标系下，卫星到波束中心地表交点的位置矢量；
根据所述卫星到波束中心地表交点的位置矢量以及所述地心到卫星的位置矢量，计算在所述地球惯性坐标系下，地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量；
根据所述地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量，以及所述地心到卫星的位置矢量与卫星的速度矢量，计算所述方位时刻对应的多普勒频率。

3.根据权利要求2所述方法，其特征在于，所述地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量，通过以下公式计算得到：
其中，
为地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量，为天线坐标系下的卫星到波束中心地表交点的单位位置矢量，A ov为地球惯性坐标系到轨道坐标系的转换矩阵的逆矩阵，A ve为轨道坐标系到本体坐标系的转换矩阵的逆矩阵，A ea为本体坐标系到天线坐标系的转换矩阵的逆矩阵，θsq为天线斜视角。

4.根据权利要求2所述方法，其特征在于，在所述地球惯性坐标系下的所述地心到卫星的位置矢量以及卫星的速度矢量，通过以下公式计算得到：
其中，
为所述地球惯性坐标系下的所述地心到卫星的位置矢量，A ov为地球惯性坐标系到轨道坐标系的转换矩阵的逆矩阵，为轨道坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的卫星的速度矢量，为轨道坐标系下的卫星的速度矢量，a为轨道的半长轴，e为轨道的离心率，u为卫星纬度角，w为近地点辐角，μ为地球引力常数。

5.根据权利要求2所述方法，其特征在于，所述地球惯性坐标系下，卫星到波束中心地表交点的位置矢量，通过以下公式计算得到：
其中，
为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量，为卫星到
波束中心地表交点的位置矢量的模值，为所述地球惯性坐标系下的卫星到波束中心
地表交点的单位位置矢量，以及分别为在地球惯性坐标系的x 轴、y轴、z轴上的投影，为地球惯性坐标系下的所述地心到卫星的位置矢量，E a为地球的
长半轴，E b为地球的短半轴，以及分别为在地球惯性坐标系的x轴、y轴、z轴上的投影。

6.根据权利要求2所述方法，其特征在于，在所述地球惯性坐标系下，地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量，通过以下公式计算得到：
为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量，为所述地球惯性坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量，为所述地球惯性坐标系下的波束中心地表交点的速度矢量，ωe为地球的自转角速度。

7.根据权利要求2所述方法，其特征在于，所述方位时刻对应的多普勒频率通过以下公式计算得到：
其中，f a为所述多普勒频率，为所述地球惯性坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量，为所述地球惯性坐标系下的卫星的速度矢量，为所述地球惯性坐标系下的波束中心地表交点的速度矢量，λ为SAR信号波长。

8.根据权利要求1所述方法，其特征在于，通过多项式拟合，计算多普勒历程的各阶系数，以得到所述多普勒历程：
其中，f a(t a)为所述多普勒历程，f dc、f dr1、f dr2、f dr3以及f dr4分别为各阶系数，t a为方位时刻。

9.根据权利要求1所述方法，其特征在于，所述相对距离历程通过以下公式计算：
其中，ΔR a(t a)为所述相对距离历程，f dc、f dr1、f dr2、F dr3以及f dr4分别为所述多普勒历程中的各阶系数，t a为SAR方位时刻，λ为SAR信号波长。

10.根据权利要求1所述方法，其特征在于，所述方位去斜函数为距离时域的方位去斜函数，或距离频域的方位去斜函数；其中，
距离时域的方位去斜函数通过以下公式计算得到：
距离频域的方位去斜函数通过以下公式计算得到：
H dechirp(t a)为距离时域的方位去斜函数，t a为SAR方位时刻，λ为SAR信号波长，ΔR a(t a)为所述相对距离历程，H dechirp(fτ，t a)为距离频域的方位去斜函数，fτ为距离频率，c为光速，
f c为SAR信号载频。

一种SAR滑动聚束模式方位去斜方法
技术领域
[0001]本发明涉及合成孔径雷达(SAR)数据处理技术领域，具体涉及一种SAR滑动聚束模式方位去斜方法。

背景技术
[0002]滑动聚束模式是先进体制高分宽幅SAR最为常用的一种成像模式，在滑动聚束模式下，天线波束中心的交点在地平面下方的某一点处，这一点被称为虚拟转动点。

现有的频域成像算法的第一步通常是对回波信号进行方位去斜处理，采用的方位去斜方法通常先计算虚拟转动点到卫星的距离历程，再利用该距离历程实现方位去斜。

可以看出，现有的方位去斜方法依赖虚拟转动点。

[0003]但是，随着SAR方位幅宽和方位分辨率的提高，天线波束方位扫描角范围迅速增大，匀速扫描控制下波束中心会明显偏离虚拟转动点，导致虚拟转动点模型失效，此时，无法继续利用虚拟转动点进行方位去斜，使得基于方位去斜的频域成像算法也难以继续适用。

[0004]因此，亟需一种不依赖虚拟转动点的方位去斜方法，以解决大方位扫描角滑动聚束模式下的虚拟转动点模型失效的问题。

发明内容
[0005]本发明的目的是提供一种SAR滑动聚束模式方位去斜方法，以解决大方位扫描角滑动聚束模式下的虚拟转动点模型失效的问题。

[0006]为达到上述目的，本发明实施例提供一种SAR滑动聚束模式方位去斜方法，方法包括：
[0007]根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角，确定每个方位时刻对应的多普勒频率；
[0008]根据所述每个方位时刻对应的多普勒频率，确定所述成像时间内的多普勒历程；[0009]根据所述多普勒历程，确定所述成像时间内的相对距离历程；
[0010]根据所述相对距离历程，确定方位去斜函数；
[0011]将回波信号与所述方位去斜函数相乘，以对所述回波信号进行方位去斜。

[0012]由以上本发明实施例提供的技术方案可见，本发明提供的方位去斜方法不再依赖虚拟转动点，通过波束中心地表交点计算瞬时多普勒历程，再积分得到相对距离历程，就可以构建方位去斜函数，从而实现方位去斜处理。

此外，基于本发明提供的方位去斜方法，SAR 系统设计也不再需要严格按照虚拟转动点模型，方位波束可以采用更加简单的匀速扫描，使得系统设计和实现变得更加灵活和简单，并且本发明提供的方位去斜方法还能够精确地去除波束方位扫描引起的多普勒混叠。

附图说明
[0013]为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0014]图1是依赖虚拟转动点的滑动聚束模式示意图；
[0015]图2是不依赖虚拟转动点的滑动聚束模式示意图；
[0016]图3是本公开实施例提供的方位去斜方法的流程图；
[0017]图4是本公开实施例提供的方位去斜方法的算法流程框图。

具体实施方式
[0018]下面将结合附图和具体实施方式，对本发明的技术方案作详细说明，应理解这些实施方式仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落入本发明所附权利要求所限定的范围内。

[0019]除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。

本文中在本发明的说明书中所使用的术语只是为了描述具体的实施方式的目的，不是旨在于限制本发明。

[0020]参考图1所示，为依赖虚拟转动点的滑动聚束模式示意图。

参考图2所示，为不依赖虚拟转动点的滑动聚束模式示意图。

现有的方位去斜方法都是针对图1所示的依赖虚拟转动点的滑动聚束模式，为了使SAR数据处理不再依赖虚拟转动点，从而可以在图2所示的滑动聚束模式下进行SAR系统设计与数据处理，本公开实施例提供了一种SAR滑动聚束模式方位去斜方法，参考图3所示，该方法可以包括：
[0021]S1：根据成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角，确定每个方位时刻对应的多普勒频率。

[0022]其中，参考图4所示，步骤S1还可以包括子步骤S101至S104：
[0023]S101：根据所述成像时间内的每个SAR方位时刻的轨道六根数、卫星姿态角以及天线扫描角，计算在地球惯性坐标系下，卫星到波束中心地表交点的单位位置矢量、地心到卫星的位置矢量以及卫星的速度矢量。

[0024]具体的，地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量，通过以下公式计算得到：
[0025]
[0026]其中，
[0027]
[0028]
[0029]
[0030]
[0031]u1＝-u-π/2
[0032]i1＝i-π/2
[0033]上述公式中的为地球惯性坐标系下的卫星到波束中心地表交点的单位位置矢量，为天线坐标系下的卫星到波束中心地表交点的单位位置矢量，A ov为地球惯性坐
标系到轨道坐标系的转换矩阵的逆矩阵，A ve为轨道坐标系到本体坐标系的转换矩阵的逆矩阵，A ea为本体坐标系到天线坐标系的转换矩阵的逆矩阵，θlook为天线视角，θsq为天线斜视角，q为升交点赤径，θyaw为卫星的偏航角，θpitch为卫星的俯仰角，θroll为卫星的横滚角，u为卫星纬度角，i为轨道的倾角。

[0034]具体的，地球惯性坐标系下的地心到卫星的位置矢量以及卫星的速度矢量，通过以下公式计算得到：
[0035]
[0036]
[0037]其中，
[0038]
[0039]
[0040]上述公式中的为所述地球惯性坐标系下的所述地心到卫星的位置矢量，A ov为地球惯性坐标系到轨道坐标系的转换矩阵的逆矩阵，为轨道坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的卫星的速度矢量，为轨道坐标系下的卫星的速度矢量，a为轨道的半长轴，e为轨道的离心率，u为卫星纬度角，w为近地点辐角，μ为
地球引力常数。

[0041]S102：根据所述单位位置矢量以及所述地心到卫星的位置矢量，计算在所述地球惯性坐标系下，卫星到波束中心地表交点的位置矢量。

[0042]具体的，地球惯性坐标系下，卫星到波束中心地表交点的位置矢量，可以通过以下公式计算得到：
[0043]
[0044]其中，
[0045]
[0046]
[0047]
[0048]
[0049]上述公式中的为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量，为卫星到波束中心地表交点的位置矢量的模值，为所述地球惯性坐标系
下的卫星到波束中心地表交点的单位位置矢量，以及分别为在地球惯性坐标系的x轴、y轴、z轴上的投影，为地球惯性坐标系下的所述地心到卫星的
位置矢量，E a为地球的长半轴，E b为地球的短半轴，以及分别为在地球惯性坐标系的x轴、y轴、z轴上的投影。

[0050]S103：根据所述卫星到波束中心地表交点的位置矢量以及所述地心到卫星的位置矢量，计算在所述地球惯性坐标系下，地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量。

[0051]具体的，在地球惯性坐标系下，地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量，通过以下公式计算得到：
[0052]
[0053]
[0054]上述公式中的为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量，为所述地球惯性坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的卫星到波束中心地表交点的位置矢量，为所述地球惯性坐标系下的波束中
心地表交点的速度矢量，ωe为地球的自转角速度。

[0055]S104：根据所述地心到波束中心地表交点的位置矢量与波束中心地表交点的速度矢量，以及所述地心到卫星的位置矢量与卫星的速度矢量，计算所述方位时刻对应的多普勒频率。

[0056]具体的，方位时刻对应的多普勒频率通过以下公式计算得到：
[0057]
[0058]上述公式中的f a为所述多普勒频率，为所述地球惯性坐标系下的所述地心到卫星的位置矢量，为所述地球惯性坐标系下的地心到波束中心地表交点的位置矢量，
为所述地球惯性坐标系下的卫星的速度矢量，为所述地球惯性坐标系下的波束中心地表交点的速度矢量，λ为SAR信号波长。

[0059]S2：根据所述每个方位时刻对应的多普勒频率，确定所述成像时间内的多普勒历程。

[0060]具体的，通过多项式拟合，计算多普勒历程的各阶系数，以得到所述多普勒历程：
[0061]
[0062]上述公式中的f a(t a)为所述多普勒历程，f dc、f dr1、f dr2、f dr3以及f dr4分别为各阶系数，t a为SAR方位时刻。

[0063]S3：根据所述多普勒历程，确定所述成像时间内的相对距离历程。

[0064]具体的，根据相对距离历程和多普勒历程的关系表达式：
[0065]f a(t a)＝-2ΔR′a(t a)/λ
[0066]通过积分可以计算得到相对距离历程：
[0067]
[0068]上述公式中的ΔR a(t a)为所述相对距离历程，f dc、f dr1、f dr2、f dr3以及f dr4分别为所述多普勒历程中的各阶系数，t a为SAR方位时刻，λ为SAR信号波长。

[0069]S4：根据所述相对距离历程，确定方位去斜函数。

[0070]具体的，方位去斜函数为距离时域的方位去斜函数，或距离频域的方位去斜函数；其中，距离时域的方位去斜函数通过以下公式计算得到：
[0071]
[0072]距离频域的方位去斜函数通过以下公式计算得到：
[0073]
[0074]上述公式中的H dechirp(t a)为距离时域的方位去斜函数，t a为SAR方位时刻，λ为SAR 信号波长，ΔR a(t a)为所述相对距离历程，H dechirp(fτ，t a)为距离频域的方位去斜函数，fτ为距离频率，c为光速，f c为SAR信号载频。

[0075]S5：将回波信号与所述方位去斜函数相乘，以对所述回波信号进行方位去斜。

[0076]由以上本发明实施例提供的技术方案可见，本发明提供的方位去斜方法不再依赖虚拟转动点，通过波束中心地表交点计算瞬时多普勒历程，再积分得到相对距离历程，就可以构建方位去斜函数，从而实现方位去斜处理。

此外，基于本发明提供的方位去斜方法，SAR 系统设计也不再需要严格按照虚拟转动点模型，方位波束可以采用更加简单的匀速扫描，使得系统设计和实现变得更加灵活和简单。

并且本发明提供的方位去斜方法还能够精确地去除波束方位扫描引起的多普勒混叠。

[0077]本说明书中的上述各个实施方式均采用递进的方式描述，各个实施方式之间相同相似部分相互参照即可，每个实施方式重点说明的都是与其他实施方式不同之处。

[0078]以上所述仅为本发明的几个实施方式，虽然本发明所揭露的实施方式如上，但所述内容只是为了便于理解本发明的技术方案而采用的实施方式，并非用于限定本发明。

任何本发明所属技术领域的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施方式的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附权利要求书所界定的范围为准。

图1
图2
图3
图4。