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