[发明专利]基于分布式散射体时序干涉SAR技术的遗产地形变监测方法

权利要求书

1.基于分布式散射体时序干涉SAR技术的遗产地形变监测方法，其特征在于，包括如下步骤：

(一)确定遗产地的监测区域及资料收集

(11)确定遗产地的监测范围，包括遗产地和所述遗产地周边的缓冲带；

(12)调查所述遗产地的背景知识，并收集所述遗产地的时间序列SAR单视复数影像；

(二)SAR数据处理与参数反演

(21)对所述步骤(12)收集的时间序列SAR单视复数影像进行互配准、辐射校正和影像裁剪预处理；

(22)提取永久散射体目标和分布式散射体目标；

(23)去除分布式散射体目标的地形与平地相位，采用最大似然估计和LBFGS技术进行优化相位估计；

(24)融合初始相位的永久散射体目标和经过优化相位的分布式散射体目标，设置短时间/空间基线阈值生成小基线集干涉对，对所述小基线集干涉对进行相位解缠并去除地形、大气延迟造成的误差；

(25)求解平均年形变速率及形变历史变化值。

2.根据权利要求1所述的基于分布式散射体时序干涉SAR技术的遗产地形变监测方法，其特征在于，所述方法还包括结果分析与解译步骤，具体为：利用野外实地调查测量结果对步骤(25)所得结果进行精度评定；并结合野外实地调查测量结果与步骤(25)所得结果共同分析所述遗产地的形变特征及引发不稳定的因素。

3.根据权利要求1所述的基于分布式散射体时序干涉SAR技术的遗产地形变监测方法，其特征在于，所述步骤(12)中所述遗产地的背景知识包括所处的环境、遭受的威胁、陆地景观、地表特征、空间尺度、遗产单体的几何形状和材料属性。

4.根据权利要求1所述的基于分布式散射体时序干涉SAR技术的遗产地形变监测方法，其特征在于，所述步骤(22)中提取永久散射体目标的具体方法为：根据单幅影像的光谱特征与多幅影像的振幅离差阈值法提取。

5.根据权利要求1所述的基于分布式散射体时序干涉SAR技术的遗产地形变监测方法，其特征在于，所述步骤(22)中提取分布式散射体目标的具体方法包括：

A、对所述时间序列SAR单视复数影像中像元的时序幅度值进行排序，转换得到像元幅度值累积分布函数，其像元幅度值累积分布函数的无偏估计表示为:

$F_N\left(X\right)=\left\{\begin{array}{ccc}0,& if& X\<x_1\\ \frac{k}{N},& if& x_k\≤X\<x_{k+1}\\ 1,& if& X\≥x_N\end{array}\right.$

其中，F_N(X)为像元幅度值累积分布函数的无偏估计，X代表像元幅度值，N为配准后的SAR影像数量，x₁为时序幅度值排序中的第1个元素的幅度值，x_k为时序幅度值排序中的第k个元素的幅度值；x_k+1为时序幅度值排序中的第k+1个元素的幅度值；k取值范围为0…N-1。

B、再用KS检验或AD检验方法对相邻像元的无偏估计值进行判定，满足一定阈值的即为同质像元；

C、对于每个像元可获得其估计窗口内连通的所有同质像元，同质数量大于设定阈值的像元即为分布式散射体目标。