[发明专利]一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法

权利要求书

1.一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法，其特征在于，包括以下几个步骤：

步骤1：获取未覆盖整个SAR影像时序过程的时域离散InSAR干涉对；

将覆盖待监测矿区的所有SAR影像按照时间先后，利用差分合成孔径雷达干涉测量D-InSAR进行两两差分干涉，获得相干图组和解缠图组，从相干图组中找出无法干涉的InSAR干涉对，并剔除与之对应的相干图和解缠相位图，从而获得未覆盖时序过程的时域离散的InSAR干涉对及对应相干图组和解缠图组；

所述相干图，是评价两幅SAR影像相似程度的依据，在干涉处理中生成；

所述无法干涉，是指相干图的相干性小于设定的相干性阈值；

步骤2：获取时域离散InSAR干涉对中的高相干点；

依据设定的相干性阈值从未覆盖时序过程的时域离散的InSAR干涉对对应的相干图组中提取出每个时域离散InSAR干涉对中的高相干点；

所述高相干点是所有时域离散的相干图在该点的相干性均设定的相干性阈值；

步骤3：利用低通滤波削弱高相干点处的大气延迟和噪声，且忽略水平移动对雷达视线向形变的贡献，高相干点(i,j)的解缠相位δφ为：

δφ(i,j)=4πλcosθ(i,j)[W(tB,i,j)-W(tA,i,j)]+4πλB⊥Δh(i,j)rsinθ(i,j)]]>

其中，λ为雷达波长，θ(i,j)为高相干点的雷达入射角，r为雷达卫星距目标的距离，B_⊥为两幅SAR影像的垂直基线长度，t_B，t_A分别为两幅SAR影像的获取时间，上述参数均从对应的SAR影像的头文件中直接获取；

W为高相干点的地表下沉值，从时域离散InSAR干涉对的解缠相位图中获取，Δh为在高相干点的高程残差，为待求系数；

步骤4：选取矿区动态沉降模型为W(Δt)＝f(Δt,P)；

式中，Δt为相对于矿区初始沉降时刻t₀的间隔时间；f为动态沉降模型映射函数；P为模型待估参数，个数为numP；

步骤5：建立动态沉降模型待估参数P与解缠相位的关系方程；

δφ(i,j)=4πλcosθ(i,j)[f(ΔtB,P,i,j)-f(ΔtA,P,i,j)]+4πλB⊥Δh(i,j)rsinθ(i,j)]]>

其中，Δt_B＝t_B-t₀和Δt_A＝t_A-t₀为校正后的SAR影像时间；

步骤6：利用大于或等于numP+1个时域离散InSAR干涉对的解缠相位及对应雷达波长λ、入射角θ、斜距r和各干涉对的垂直基线长度B_⊥分别代入步骤5所述的公式中，将得到的多个方程联立方程组，计算出动态沉降模型的待估参数P和高相干点的高程残差Δh，所述时域离散InSAR干涉对的解缠相位从步骤1获得的解缠相位图中获得；将动态沉降模型的待估参数P代入步骤4选取的矿区动态沉降模型，计算出任意时刻的地表时序沉降，实现基于时域离散InSAR干涉对的矿区时序形变监测。