[发明专利]一种全散射体FS-InSAR方法及系统

权利要求书

1.一种全散射体FS-InSAR方法，其特征在于：包括以下步骤：

选定监测区，基于第一干涉像对组合并结合所述监测区的时序SAR影像生成差分干涉图；

对所述差分干涉图进行相位解缠，得到解缠相位；

根据所述解缠相位得到每个SAR成像时刻对应的时序相位；

根据双尺度时域低通滤波，抑制所述时序相位中的大气相位和噪声相位，得到时序低频相位；

重新计算干涉像对组合，得到第二干涉像对组合，根据所述第二干涉像对组合对所述时序低频相位进行差分，得到差分低频相位；

根据所述时序SAR影像得到平均幅度，根据平均幅度提取所述监测区中非水体的全散射体；

基于所述全散射体和所述差分低频相位得到地表形变信息。

2.根据权利要求1所述的一种全散射体FS-InSAR方法，其特征在于：所述选定监测区，基于第一干涉像对组合并结合所述监测区的时序SAR影像生成差分干涉图，包括：

设置第一时间基线阈值和第一垂直基线阈值得到所述第一干涉像对组合，所述第一干涉像对组合为连通网络，所述第一时间基线阈值小于等于两个卫星重访周期。

3.根据权利要求2所述的一种全散射体FS-InSAR方法，其特征在于：所述对所述差分干涉图进行相位解缠，得到解缠相位，包括：

采用相位解缠方法对所述差分干涉图进行相位解缠，得到解缠相位；

所述相位解缠方法包括最小费用流法、枝切法和区域生长法。

4.根据权利要求3所述的一种全散射体FS-InSAR方法，其特征在于：所述根据所述解缠相位得到每个SAR成像时刻对应的时序相位，包括：

利用最小二乘方法计算每个SAR成像时刻的所述解缠相位，得到对应的时序相位，所述时序相位可表示为：

（1）

其中，为形变相位，为DEM高程误差相位，为轨道误差相位，为大气延迟相位，为噪声相位。

5.根据权利要求4所述的一种全散射体FS-InSAR方法，其特征在于：在所述根据双尺度时域低通滤波，抑制所述时序相位中的大气相位和噪声相位，得到时序低频相位前，包括：

确定双尺度时间窗口；

所述双尺度时间窗口包括第一尺度时间窗口和第二尺度时间窗口，所述第一尺度时间窗口小于第二尺度时间窗口。

6.根据权利要求5所述的一种全散射体FS-InSAR方法，其特征在于：所述根据双尺度时域低通滤波，抑制所述时序相位中的大气相位和噪声相位，得到时序低频相位，包括：

采用第一尺度时间窗口并通过低通滤波器对所述时序相位进行低通滤波，得到第一低频相位：

（2）

（3）

其中，为时序相位，为高斯函数，为第一尺度时间窗口下输出的时刻的低通滤波相位，为N幅SAR影像的成像日期，为第一尺度时间窗口；

利用所述时序相位减去所述第一低频相位，得到残余相位，所述残余相位为：

（4）

采用第二尺度时间窗口并通过低通滤波器对所述残余相位再次进行低通滤波，得到第二低频相位：

（5）

（6）

其中，为第二尺度时间窗口，为第二尺度时间窗口下输出的时刻的低通滤波相位；

将第一低频相位和第二低频相位叠加，得到时序低频相位，时序低频相位为：

（7）

其中，为所述时序低频相位。

7.根据权利要求6所述的一种全散射体FS-InSAR方法，其特征在于：所述重新计算干涉像对组合，得到第二干涉像对组合，根据所述第二干涉像对组合对所述时序低频相位进行差分，得到差分低频相位，包括：

设置第二时间基线阈值和第二垂直基线阈值得到所述第二干涉像对组合，所述第二干涉像对组合的数目大于所述第一干涉像对组合的数目；

根据所述第二干涉像对组合，对所述时序低频相位进行差分，得到差分低频相位：

（8）

其中，为主影像对应的时序低频相位，为辅影像对应的时序低频相位，,。