[发明专利]时间差分基线集时序干涉SAR处理方法

权利要求书

1.一种时间差分基线集时序干涉SAR处理方法，其特征在于包括：

对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对；

利用外部高程模型确定处理区域地形相位和平地相位，用于去除处理区域的地形和平地相位；

根据已形成的干涉对确定空间基线；

通过一维回归方式确定高程校正:

对空间上相近的高相干散射点CS做二阶差分，得到如下相位：

将上式中除高程误差相位外的误差相位统一记为残余相位：

则有：

建立如下的目标函数：

对目标函数优化后可得到任意两点之间的高程误差值Δε_i,j，然后在空间上构网进行平差处理得到任意一点的高程误差值Δε_i；

通过一维回归方式确定形变；

获取大气相位和非线性形变，经高程校正和形变解算，确定线性形变、非线性形变、大气相位和高程校正信号量。

2.根据权利要求1所述的方法，其特征在于：

N+1景已配准的SAR影像确定方式包括：

在N+1景SAR影像中选取一景影像作为主影像，其他N景影像作为从影像分别与主影像进行配准，最终得到N+1景已经配准的SAR影像；

通过对N+1景已经配准的SAR影像进行处理，获取处理区域的高相干散射点。

3.根据权利要求1所述的方法，其特征在于，所述外部高程模型通过以下至少一种方式获取：

引入外部DEM，通过地理编码转换到SAR坐标系；

任选SAR影像中两景影像数据，通过InSAR处理获取。

4.根据权利要求2所述的方法，其特征在于，对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对，包括：

对SAR影像进行二轨差分处理，确定处理区域的形变信息。

5.根据权利要求1所述的方法，其特征在于，对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对，包括：

将N+1幅已经配准的影像按时间顺序排列为：

{φ₀,φ₁,φ₂,...,φ_N}；

相邻时间的影像组合进行干涉，生成N幅干涉图：

{δφ₁,δφ₂,...,δφ_N}；

确定相干散射点的干涉相位，表示如下：

δφ_k＝φ_k-φ_k-1,k＝1,2,...,N；

确定相位信号，表示如下：

其中，δφ_k,i表示第k幅干涉图的第i个点的干涉相位，注意它是缠绕的，范围在-π到π之间，{·}_2π表示缠绕操作，表示高程误差相位，由于不精确的外部数字高程模型引起，表示形变信号相位，表示轨道误差相位，表示大气延迟相位，表示随机噪声相位。

6.根据权利要求5所述的方法，其特征在于，所述高程误差相位为：

其中，表示垂直基线，λ表示传输信号波长，R_k,i表示卫星天线到散射点的距离，θ_k,i表示信号的入射角，Δε_i表示高程误差；

所述形变相位包括线性形变项与非线性形变项，表示如下：

其中，δT_k表示第k幅干涉图的时间基线，表示沿视线方向的线性形变速率，表示非线性形变相位。

7.根据权利要求1所述的方法，其特征在于：一维回归方式确定形变，包括：

获取线性形变速率，经过高程校正后，得到包含形变信息的残余相位积分后如下式：

其中：

于是，满足：

建立如下形变的目标函数：

对于形变的目标函数优化，得到任意两点间的线性形变速率然后在空间上构网进行平差处理得到任意一点的线性形变速率