多角度多时段导航卫星双基地PS-InSAR三维形变反演方法
【专利摘要】本发明公开了一种多角度多时段导航卫星双基地PS?InSAR三维形变反演方法，该方法通过在待监测的场景附近放置接收机，利用全向天线接收卫星直达波信号，利用喇叭天线接收场景回波信号，利用直达波信号对回波信号进行同步、成像处理获得多角度多时段的SAR图像序列，通过PS方法实现PS点的识别与形变量估计，再利用空间、时间插值实现不同角度下的形变量空间、时间的连续化，最后利用加权最小二乘实现三维形变量的反演；本发明成本低，配置灵活，可广泛应用在形变监测领域。
【专利说明】
多角度多时段导航卫星双基地PS-1 nSAR三维形变反演方法
技术领域
[0001] 本发明属于雷达信号处理的技术领域，具体涉及一种多角度多时段导航卫星双基 地PS-InSAR三维形变反演方法。
【背景技术】
[0002] 基于导航卫星的双基地重轨差分干涉SAR利用导航卫星作为发射平台，地面配置 接收机，对需要探测形变的区域进行连续成像、差分干涉和形变反演处理可获得监测区域 的一维形变。
[0003] 然而，单颗导航卫星的重轨时间较长，同时仅能实现一维形变的反演，精度受限， 实用性较差。实际上，导航卫星系统由多颗卫星组成，对于同一地区始终有多颗卫星的信号 连续覆盖，而且随着时间的不同，卫星相对于一个特定场景的几何关系也会发生变化，即可 提供不同视角的信息。因此，该系统在提高形变监测精度方面具有很大潜力。
[0004] 因此，开发一种基于多角度多时段导航卫星照射的双基地PS-InSAR的三维连续形 变监测方法，对于高精度形变测量领域具有重要意义。
【发明内容】
[0005] 有鉴于此，本发明提供了一种多角度多时段导航卫星双基地PS-InSAR三维形变反 演方法，能够利用多卫星、多时段的回波数据实现照射场景的三维形变测量，其精度高、成 本低且时间连续。
[0006] 实现本发明的技术方案如下：
[0007] 步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号 进行同步和成像处理，得到同一场景的不同角度的双基地SAR图像序列；
[0008] 步骤二:对所述双基地SAR图像序列，利用PS方法实现各角度下的PS点识别与PS点 的一维形变量的反演；
[0009 ]步骤三:根据所述PS点的一维形变量的反演结果，利用Kr i g ing插值方法对所述双 基地SAR图像序列的每个像素的一维形变量进行反演；
[0010]步骤四：设置一个统一的时间轴，利用时间序列分析方法对所述双基地SAR图像序 列的每个像素的一维形变量的反演结果进行时间插值，得到同一场景的不同角度下的每个 像素在统一的时间轴下的一维形变量历史；
[0011]步骤五:利用加权最小二乘估计方法，对所述每个像素在统一的时间轴下的一维 型变量历史进行处理，实现每个像素的三维形变量历史的反演。
[0012] 有益效果：
[0013] 本发明对比已有技术，充分利用了导航卫星数目多、导航信号覆盖范围广、照射时 间长的特点，与监测区域附近的接收机配置构成双基地SAR系统，可对监测区域实现时间、 空间连续的三维形变测量，能够获取更加全面的形变信息。同时作为被动测量方式，本发明 成本低，配置灵活，可广泛应用在形变监测领域。
【附图说明】
[0014] 图1为本发明系统构型示意图。
[0015] 图2为本发明算法总流程图。
【具体实施方式】
[0016] 下面结合附图并举实施例，对本发明进行详细描述。
[0017] 本发明方法所利用的系统构型如图1所示，不同卫星从不同角度发射导航信号到 场景，场景反射回波到回波天线，同时，一个全向直达波天线接收卫星的直达波信号。对于 每颗卫星，采集其不同重轨时段的信号。以接收机直达波天线作为原点0,正东为X轴，正北 为Y轴，竖直向上为Z轴建立空间直角坐标系0-ΧΥΖ。如图2所示为本发明方法的流程，该方法 具体步骤如下：
[0018] 步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号 进行同步和成像处理，得到同一场景的不同角度下的双基地SAR图像序列。设第i颗卫星第j 组时段(一组时段对应同一个角度）的M^+l幅重轨图像序列为Q 1>J:
[0020] 其中，S表示卫星总数表示第i颗卫星的时段总数。qi,^)表示第i颗卫星第j组时 段的第一幅图像，表示第i颗卫星第j组时段的第二幅图像，依此类推。
[0021] 第i颗卫星第j组时段的卫星重轨时间序列Tu为：
[0022] Tj j - \tj ] (2)
[0023] 其中，表示第i颗卫星第j组时段的第一个重轨时刻，表示第i颗卫星第j 组时段的第二个重轨时刻，其他依此类推。
[0024]由于不同时段，卫星相对于观测场景的观测角度不同，因此，所有卫星所有时段的 观测角度总数L为：
[0026]步骤二:对第i颗卫星第j组时段的M^+l幅重轨图像进行PS点识别与PS点的一维形 变量反演，得到各个PS点沿着等效视线方向的一维形变反演结果10)，设为 (设q^,o为主图像）：
[0028]其中，Pi,j为第i颗卫星的第j组图像中的PS点集合。 第i颗卫星的第j组图像中PS点沿着等效视线方向的时刻相对于时刻的形变反演 结果。Δ 表示第i颗卫星的第j组图像中PS点沿着等效视线方向的 时刻相对于时刻的形变反演结果。其他依此类推。
[0029]步骤三：由于不同卫星不同角度下探测到的PS点不同，无法直接计算某个像素的 三维形变量，因此，根据不同角度下SAR图像中各个PS点的一维形变反演结果，利用Kriging 插值方法，对每个角度下整个场景的双基地SAR图像序列的每个像素进行空间插值，得到空 间插值后的每个角度下场景中每个像素的形变量矩阵：
[0031] 其中，#表示场景中的像素点集合。Δ&,·(户;為^)表示空间插值后的Δ 11;J (Pi,j;ti,j,0，ti,j,1)，Δ/，. ?;.,·」)表不;?? 间插值后的 A li,j(Pi,j;ti,j,Q，ti,j,2)，其他依 此类推。
[0032] 步骤四：由于卫星的重轨时刻未必完全相同，为实现场景中每个像素三维形变量 历史的提取，需要将所获取的不同角度下的形变量历史的时间轴进行统一。设统一时间轴 为：
[0033] t = [to，ti，."，tN-1] (6)
[0034] 其中，to为时间轴的第1个时刻，tk = to+k · T(k = 1，2，…，N-1)为该时间轴的第k+1 个时刻，其中，T表示统一时间轴的时间间隔，N表示统一时间轴的时刻总数。
[0035] 那么，对进行基于时间序列分析方法的时间插值得到空间、时间插 值后的形变量矩阵：
[0037]其中，减，/為4)表示时间插值后的形变量在to时刻的形变量，其他依此类推。 [0038] 步骤五：
[0039]设在第tm(m=0，1，…N-1)时刻，场景中的点户e #的三维形变量ΔΖΗ卢;匕,）为：
[0041 ]其中，Δχ(户;/",)、和分别为点多在第tj寸刻沿着前述坐标系 0-ΧΥΖ上的X、Y和Z轴方向的形变量。
[0042]则有：
[0044]其中，为噪声，标准差为cr,y(卢;im) , Δ?，7〇;?,")为点乡在第i颗卫星 第j组时段的形变量进行时间插值和空间插值后在时刻。的形变量，为等效视线方 向的向量，可写为：
[0046] 其中：
[0047] Ψ (Sij) = [xij,yij,Zij]T (11)
[0048] 表示卫星i在第j组时段的孔径中心时刻的三维位置，其中x小ydPzu分别表示卫 星i在第j组时段的孔径中心时刻在前述坐标系0-ΧΥΖ上的X、Y和Z轴方向的坐标；
[0049] W(R) = [xR,yR,zR]T (⑵为回波天线R的三维位置。其中xR、yR和z R分别 表示回波天线R在前述坐标系0-ΧΥΖ上的X、Y和Z轴方向的坐标；
[0051 ]为点彡的三维位置。其中·％、巧和4分别表示点多在前述坐标系0-ΧΥΖ上的Χ、Υ和 Ζ轴方向的坐标；
[0052] 那么，点#处三维形变量的加权最小二乘估计为：
[0059] 其中，diag表示对角矩阵，即：Φ的第一行第一列为1/σ?_()(兵第二行第二列 为1 / )，其他依此类推。
[0060] 综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。 凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的 保护范围之内。
【主权项】
1.多角度多时段导航卫星双基地PS-InSAR三维形变反演方法，其特征在于，包括以下 步骤： 步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号进行 同步和成像处理，得到同一场景的不同角度的双基地SAR图像序列； 步骤二:对所述双基地SAR图像序列，利用PS方法实现各角度下的PS点识别与PS点的一 维形变量的反演； 步骤三:根据所述PS点的一维变形量的反演结果，利用Kriging插值方法对所述双基地 SAR图像序列的每个像素的一维形变量进行反演。 步骤四：设置一个统一的时间轴，利用时间序列分析方法对所述双基地SAR图像序列的 每个像素的一维形变量的反演结果进行时间插值，得到同一场景的不同角度下的每个像素 在统一的时间轴下的一维形变量历史； 步骤五:利用加权最小二乘估计方法，对所述每个像素在统一的时间轴下的一维型变 量历史进行处理，实现每个像素的三维形变量历史的反演。
【文档编号】G01S13/90GK105866777SQ201610187369
【公开日】2016年8月17日
【申请日】2016年3月29日
【发明人】曾涛, 田卫明, 张天, 胡程
【申请人】北京理工大学