1.一种SAR偏移量二维形变时序计算方法,其特征在于,包括以下步骤:步骤1,获取待监测地表区域在不同时间点的SLC影像,每两个不同时间点对应的SLC影像形成一组SLC影像对,得到多组SLC影像对;求取多组SLC影像对的二维形变数据L1;
步骤2,对待监测地表区域的每个监测点的形变建立复线性模型:A1X1=L1 (1)式中,A1表示设计矩阵;X1表示监测点的二维形变时间序列;M表示SLC影像对的数量;N表示SLC影像的个数;
(1) T ‑1 T
步骤3,根据X =(A1P1A1) A1 P1L1和 求解监测点的二维形变时间序(1)
列X1的估计值X 及X1的协因数矩阵
T
其中,A1表示设计矩阵;P1表示L1的先验信息;A1表示A1的转置矩阵;
步骤4,动态更新待监测地表区域的二维形变时间序列:步骤4.1,当每增加一景新的SLC影像,将该新的SLC影像分别与步骤1中的SLC影像进行组队,形成多组SLC影像;求取多组SLC影像对的二维形变数据L2;
步骤4.2,对待监测地表区域的每个监测点的形变建立如下复线性模型:式中,A2和B均表示设计矩阵;N′表示更新后SLC影像的个数;M′表示更新后SLC影像对的(2)数量;X 表示更新后的二维形变时间序列;Y表示增加的一景新的SLC影像对应的时间点累计形变;
(2)
步骤4.3,利用公式(4)求解得到X 和Y,(2) T ‑1
式中,Jx是增益矩阵,QJ是X 的协因数矩阵;B表示B的转置矩阵;QJ 表示QJ的逆矩阵;
P2表示L2的先验信息。
2.如权利要求1所述的SAR偏移量二维形变时序计算方法,其特征在于,所述的步骤1中,采用强度互相关系数最大化方法求取多组SLC影像对的二维形变数据L1。
3.如权利要求1所述的SAR偏移量二维形变时序计算方法,其特征在于,所述的步骤4.1中,采用强度互相关系数最大化方法求取多组SLC影像对的二维形变数据L2。
4.一种SAR偏移量二维形变时序计算系统,其特征在于,包括:SLC影像配准模块,用于获取待监测地表区域在不同时间点的SLC影像,每两个不同时间点对应的SLC影像形成一组SLC影像对,得到多组SLC影像对;二维形变数据获取模块,用于求取多组SLC影像对的二维形变数据L1;
库存数据复线性模型建立模块,用于对待监测地表区域的每个监测点的形变建立复线性模型:A1X1=L1 (1)式中,A1表示设计矩阵;X1表示监测点的二维形变时间序列;M表示SLC影像对的数量;N表示SLC影像的个数;
(1) T ‑1 T
库存数据二维形变时间序列获取模块,用于根据X =(A1 P1A1) A1 P1L1和(1)求解监测点的二维形变时间序列X1的估计值X 及X1的协因数矩阵T
其中,A1表示设计矩阵;P1表示L1的先验信息;A1表示A1的转置矩阵;
二维形变时间序列更新模块,包括:
更新后的SLC影像配准模块,用于每增加一景新的SLC影像,将该新的SLC影像分别与SLC影像配准模块中的SLC影像进行组队,形成多组SLC影像对;更新后的二维形变数据获取模块,用于求取多组SLC影像对的二维形变数据L2;
更新数据复线性模型建立模块,用于对待监测地表区域的每个监测点的形变建立如下复线性模型:式中,A2和B均表示设计矩阵;N′表示更新后SLC影像的个数;M′表示更新后SLC影像对的(2)数量;X 表示更新后的二维形变时间序列;Y表示增加的一景新的SLC影像对应的时间点累计形变;
(2)
更新数据二维形变时间序列获取模块,用于利用公式(4)求解得到X 和Y,(2) T ‑1式中,Jx是增益矩阵,QJ是X 的协因数矩阵;B表示B的转置矩阵;QJ 表示QJ的逆矩阵;
P2表示L2的先验信息。
5.如权利要求4所述的SAR偏移量二维形变时序计算系统,其特征在于,所述的二维形变数据获取模块中,采用强度互相关系数最大化方法求取多组SLC影像对的二维形变数据L1。
6.如权利要求4所述的SAR偏移量二维形变时序计算系统,其特征在于,所述的更新后的SLC影像配准模块中,采用强度互相关系数最大化方法求取多组SLC影像对的二维形变数据L2。