一种非平行轨迹分布式星载SAR空间同步方法技术

本发明专利技术提供一种非平行轨迹分布式星载SAR空间同步方法，该空间同步方法可首先结合主从星的速度和观测场景精确计算主星孔径中心时间，并得到主星天线方位角和俯仰角；进一步，基于令主辅图像的几何相关性最优的准则，在非平行轨迹条件下计算得到双基地观测的最优孔径中心时刻；最后，基于该时刻得到从星天线方位角和俯仰角，完成空间同步。本发明专利技术结合了主从星各自的速度和场景高程信息，相比于传统方法可有效提高主辅图像的相关性。可有效提高主辅图像的相关性。可有效提高主辅图像的相关性。

【技术实现步骤摘要】
一种非平行轨迹分布式星载SAR空间同步方法


[0001]本专利技术属于合成孔径雷达
，尤其涉及一种非平行轨迹分布式星载SAR空间同步方法。

技术介绍

[0002]分布式星载合成孔径雷达(SAR)利用双/多星形成近距离编队，以同时/近同时获取场景的多幅图像，实现对场景的单航过干涉测高功能。由于各卫星同时/近同时观测同一个场景，各星获得的图像之间的场景去相关可忽略不计，因此分布式星载SAR可获得精度远优于重轨星载SAR的地面数字高程模型(DEM)。由于上述优势，分布式星载SAR成为近年来国内外的研究热点。
[0003]分布式星载SAR由于收发卫星分置，需解决时间、频率和空间同步问题。其中，空间同步指的是通过实时有效地控制发射和接收平台的天线指向，使得发射和接收波束同时、精确地指向地面同一个目标，以尽可能提高系统产品的性能(如干涉测高精度等)。
[0004]目前，在已有的分布式星载SAR(TanDEM
‑
X、陆探一号等)中，采用的空间同步方法都是使主星和从星都独立、正侧视地指向场景，使图像具有相同的多普勒中心。与主从星共同指向场景中心的方式相比，这种空间同步方法尽管会略微降低辅图像的信噪比，但是能够保证主辅图像具有相同的多普勒中心，从而具有最高的多普勒相关性。
[0005]然而，上述已有分布式星载SAR空间同步方法都是基于主从星具有平行轨迹的假设前提下。当卫星之间的轨迹不平行时，即使令各主辅图像多普勒中心相同，也无法保证主辅图像的多普勒相关性最高。对于分布式星载SAR而言，图像几何相关性是决定高程反演的首要因素，因此，迫切需要突破不平行轨迹下的分布式星载SAR空间同步方法，以提高单航过测高的精度。这在现有的分布式星载SAR空间同步技术中未有提及。

技术实现思路

[0006]为解决上述问题，本专利技术提供一种非平行轨迹分布式星载SAR空间同步方法，该空间同步方法可结合主从星的速度和观测场景精确计算主星孔径中心时间、主星天线方位角和俯仰角、从星天线方位角和俯仰角，完成空间同步，提高主辅图像的相关性。
[0007]本专利技术的非平行轨迹分布式星载SAR空间同步方法包括：
[0008]S1、获取卫星的空间参数，包括主星的可开机时间段Ω＝[t
beg
,t
end
]、Ω时间内地心固定(ECF)坐标系下主星的平均速度v1、从星的平均速度v2、主星的起始位置s
1,0
、从星的起始位置s
2,0
，以及地心惯性坐标系(ECI)坐标系下主星的平均速度v
1,ECI
、从星的平均速度v
2,ECI
、主星的起始位置s
1,0,ECI
、从星的起始位置s
2,0,ECI
；获取待成像场景信息，包括可开机时间段中间时刻ECI和ECF坐标系下场景中心的坐标p
ECI
和p、场景的DEM、t
beg
时刻场景中心在ECI坐标系下由地球自转引起的瞬时运动速度v
p
；获取卫星的载荷参数，包括波长λ、带宽B、合成孔径时间T
s
；获取干涉所需的最短垂直基线b
min
。
[0009]S2、分别计算可开机时间段中间时刻，ECI和ECF坐标系下场景的平均法向量n
ECI
和
n。在ECI坐标系下，根据场景DEM构建场景中各点的三维坐标，形成三维曲面并得到各个网格内曲面的法线。对所有法线方向取均值并归一化，得到场景的平均法向量n
ECI
；类似地，在ECF坐标系下执行此操作，得到n。
[0010]S3、计算可开机时间段中间时刻主从星的最优孔径中心位置s
1,c
、s
2,c
，具体步骤包括：
[0011]S31、计算可开机时间段中间时刻，计算方式为t
1,c
＝(t
beg
+t
end
)/2；进一步，获得此时主星的位置，计算方式为s
1,c
＝s
1,0
+v1·
(t
1,c
‑
t
beg
)。
[0012]S32、计算地面等距离方向u
g
⊥
。首先计算视线方向矢量，计算方式为u＝(s
1,c
‑
p)/||s
1,c
‑
p||；然后地距方向u
g
，计算方式为u
g
＝P
g
u，其中P
g
＝I
‑
nn
T
表示以n为法向量的投影矩阵，I为单位矩阵，上标T表示转置操作；进一步，根据下式计算地面等距离方向
[0013]u
g
⊥
＝u
g
×
n/||u
g
×
n||
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0014]其中
×
为向量的叉乘运算。
[0015]S33、计算与可开机时间段中间时刻相应的最优从星孔径中心时刻t
2,c
，计算方式为
[0016][0017]其中P＝I
‑
uu
T
表示以视线方向矢量为法向量的投影矩阵。
[0018]S34、计算开可开机时间段中间时刻从星的最优孔径中心位置s
2,c
，计算方式为
[0019]s
2,c
＝s
2,0
+v2·
(t
2,c
‑
t
beg
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0020]S4、计算实际垂直基线小于最短垂直基线b
min
的时间段Φ＝[t1+t
1,c
,t2+t
1,c
]，具体步骤包括：
[0021]S41、计算开机时间段中任意时刻t对应的主从星星间垂直基线，计算方式为
[0022][0023]其中表示高程方向。Δs
c
＝s
1,c
‑
s
2,c
表示t
1,c
时刻的星间基线，Δv＝v1‑
v2表示主从星速度偏差。
[0024]S42、求解方程b
⊥
≤b
min
，得到时间段Φ＝[t1+t
1,c
,t2+t
1,c
]，其中
[0025][0026]S5、根据可开机时间段、垂直基线小于最短垂直基线的时间段，计算主星的孔径中心时刻t
c
，计算方式为
[0027][0028]其中Ω
‑
Φ表示集合Ω与Φ的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种非平行轨迹分布式星载SAR空间同步方法，其特征在于，包括以下步骤：S1、获取卫星的空间参数；S2、分别计算可开机时间段中间时刻，ECI和ECF坐标系下场景的平均法向量n
ECI
和n。在ECI坐标系下，根据场景DEM构建场景中各点的三维坐标，形成三维曲面并得到各个网格内曲面的法线；S3、计算可开机时间段中间时刻主从星的最优孔径中心位置s
1,c
、s
2,c
；S4、计算实际垂直基线小于最短垂直基线b
min
的时间段Φ＝[t1+t
1,c
,t2+t
1,c
]；S5、根据可开机时间段、垂直基线小于最短垂直基线的时间段，计算主星的孔径中心时刻t
c
；S6、确定主星开机时的天线方位向角度α1和俯仰向角度β1；S7、确定从星开机时的天线方位向角度α2和俯仰向角度β2。2.如权利要求1所...

【专利技术属性】
技术研发人员：陈志扬，李元昊，胡程，
申请(专利权)人：北京理工大学前沿技术研究院，
类型：发明
国别省市：