一种稳健的基于雷达分布式目标的地面沉降监测方法与流程

技术特征：

1.一种稳健的基于雷达分布式目标的地面沉降监测方法，其特征在于：具体包括以下步骤：

s1：sar时间序列影像的获得：利用合成孔径雷达全天时、全天候对地实施观测，按照时间序列获得多组单视复数影像(slc)，选取一组单视复数影像(slc)为主影像，剩余的各组单视复数影像(slc)为从影像，将各组从影像按照同名像素点的方式进行匹配，实现配准过程，获得sar时间序列影像；

s2：ds候选点的识别：

s201：定义目标函数：假设ap和aq分别为sar时间序列影像的p和q两个样本的幅度值，两者的表达式为：

①和②式中，f(g)是累积分布函数，当有n组sar时间序列影像数据，ap的经验累积密度分布函数(cdf)表达式为：

③式中，是一个指标函数，当ap,i≤a时为1，反之则为0；

s202：ks检验：定义p和q两个样本的经验分布函数之间的最大距离，两个样本ap和aq及其对应的经验累积密度分布函数(cdf)和ks检验统计量定义式为：

④式中，是距离集合的最大值，当n>8时检验统计量tks分布收敛于柯尔莫哥洛夫分布，给定显著性水平α，

当满足⑤式时，ap假设被拒，⑤式中kα是柯尔莫哥洛夫分布的α分位数，kα通过柯尔莫哥洛夫cdf计算，计算式为：

s203：ds候选点：当不满足⑤式时，ap和aq的假设满足要求，sar时间序列影像的p和q两个样本符合相同的概率分布函数，即p和q两个样本均为后向散射特性相同的sar时间序列影像像元点，识别为ds候选点；

s3：协方差矩阵的构建：步骤s203中，sar时间序列影像的某一像元点协方差矩阵的数学表达式为：

⑦式中，e[]表示为期望，y为像元时间序列上的标准化后的复数观测向量，y＝[y1，y2，y3…yn]，且e[|yj|²]＝1，在⑦的基础上，对步骤s203获得的ds候选点的相位信息进行自适应滤波，构建协方差矩阵，协方差矩阵数学表达式为：

⑧式中，非对角线上对应的元素的绝对值为空间自适应滤波后的相干性值，相位为空间自适应滤波后的干涉相位，获得更加精确的干涉相位和相干性；

s4：估算协方差矩阵最优相位值：对步骤s3中获得的协方差矩阵⑧进行特征分解，特征分解公式为：

⑨式中，λi为特征值，按降序排序λ1≥λ2≥,...,≥λn,ui为对应的特征向量，假设u1是最大特征分量，约束条件为则可表述为最大值约束问题，u1表述公式为：

特征分解的相位估算值表达式为：

式中，θm,n＝θm-θn，rm,n和ζm,n分别是协方差矩阵t和矩阵|u1||u1|^t中第m行n列的元素，符号o代表哈达玛积，

按照上述分解方式，从协方差矩阵中提取与最大特征值对应的特征向量的相位来估计最优相位；

s5：确定最终的ds点目标：假设估算后的最优相位为θ＝[θ1，θ2，...，θn]，为第n景和第k景slc的对应的原始干涉相位，则

式中，r为给定的阈值，大于给定阈值的ds候选点被保留并视为最终的ds点目标；

s6：获得ps点目标：对永久散射体合成孔径雷达干涉测量，获得ps点目标；

s7：最终形变速率和时序形变序列的确定：对步骤s5中获得的ds点目标和ps点目标联合进行时序insar处理，经误差改正、基线精化和相位解缠得到最终形变速率和时序形变序列。

2.根据权利要求1所述的一种稳健的基于雷达分布式目标的地面沉降监测方法，其特征在于：所述步骤s7中，误差改正包括地形改正和大气改正，地形改正的方式为：时间序列干涉相位与影像空间基线进行线性回归模型分析，其斜率即是地形改正，大气改正的方式：为通过对干涉相位在时间域和空间域进行不同的带通滤波消除大气相位，进行大气改正。

3.根据权利要求1-2任一项所述的一种稳健的基于雷达分布式目标的地面沉降监测方法，其特征在于：所述步骤s7中，基线精化的方式为：结合地面控制点，通过最小二乘等算法对空间基线进行精化估计，得到更为准确的空间基线。

4.根据权利要求1-3任一项所述的一种稳健的基于雷达分布式目标的地面沉降监测方法，其特征在于：所述步骤s7中，相位解缠的方式为：获得的干涉相位位于[-π，+π]之间，与地面真实相位之间相差nπ个整周模糊度，恢复干涉相位的整周模糊度，并生成一副连续的解缠相位图，实现相位解缠。