[发明专利]一种适用于大梯度地表沉降监测的时序差分雷达干涉方法

摘要

本申请公开了一种适用于大梯度地表沉降监测的时序差分雷达干涉方法，该方法通过短时间基线差分干涉图筛选、离散点相位解缠、基于短时间基线差分干涉图的形变分量建模和解算、形变分量可靠性检验、差分干涉图相位梯度修正、对上述过程进行迭代以确保形变分量解算正确，以及基于修正后差分干涉图的形变时间序列建模和解算等过程，并形成整体的技术方案，解决原有时序差分雷达干涉技术中由于形变和相位梯度较大导致形变相位模糊度解算困难或失败的问题，最终达到正确提取大梯度地表形变速率和形变时间序列的目的，并起到降低大梯度形变建模和解算所需合成孔径雷达影像数量的效果，节约时序差分雷达干涉的应用经济成本。

主权项

1.一种适用于大梯度地表沉降监测的时序差分雷达干涉方法，其特征在于，具体按照以下步骤实施：/n步骤1，对所有筛选后的SAR影像进行任意干涉组合并计算时间基线和空间基线；/n步骤2，设置时空基线阈值进行干涉对初始筛选，在保证干涉对数量的前提下限制时间失相干和空间失相干，按照干涉组合进行差分干涉处理得到差分干涉图集；/n步骤3，永久散射体PS和相干目标CT探测及点集合并，得到相干散射体CS点集，提取CS点集上的差分干涉相位，并以CS为节点构建不规则Delaunay三角网TIN；/n步骤4，基于Delaunay三角网和MCF方法的CS相位解缠；/n步骤5，线性形变速率和高程误差建模及解算；/n步骤6，基于数据模拟或地面测量数据的线性形变解算结果检验，若结果不可靠则执行步骤7，否则转向步骤8；/n步骤7，重新筛选参与计算的差分干涉图，并重新执行步骤5和步骤6，当达到计算终止条件时，执行步骤8；/n步骤8，从所有原始差分干涉图中减去线性形变相位分量，对残差相位进行重新解缠，然后再将线性形变相位分量加回重新解缠后的相位中，重新执行步骤5，得到更新后的线性形变速率和高程误差；/n步骤9，记录步骤8中所得到的线性形变速率，并从所有的原始差分干涉图中减去更新后的线性形变和高程误差相位分量，重新执行基于离散点的相位解缠，然后将步骤8中的线性形变分量重新加回新的解缠相位中，执行形变时间序列建模和解算过程，最终得到形变时间序列；/n步骤1具体实施方式为：排除受雨雪天气以及积雪影响的SAR影像后进行任意干涉组合配对，假设有N+1幅SAR影像，通过任意干涉组合形成N(N+1)/2个干涉对，每个干涉对由主、从两幅SAR影像构成，然后，根据每个干涉对主、从SAR影像的获取时间计算该干涉对的时间基线即SAR影像获取的时间差；根据主、从SAR影像的参数文件所记录的SAR传感器运行位置及相应的时间参数计算该干涉对的空间基线，主、从影像成像时SAR传感器的空间距离，取该空间距离在垂直于SAR传感器视线方向上的分量为空间基线，也即垂直基线。/n