一种融合像元偏移跟踪和短基线集的矿区地表形变解算方法。适用于地表矿区地表形变研究使用。其步骤为：对矿区地表形变区域进行预估；选取时间序列上相干图中的高相干点；利用时序偏移量跟踪算法恢复大形变区域高相干点的相位整周数N；大形变区域和小形变区域的高相干点联合建立短基线集解算模型；矿区地表高程误差及形变速率的解算。其监测精度高、范围大，实现操作过程简单，费用低，克服了传统时序InSAR方法无法正确获取矿区大形变梯度下的地表沉降问题，也解决了基于SAR幅度信息的时序像元偏移跟踪算法和基于差分相位的短基线集技术难以联合解算地表形变速率和高程误差的问题，具有广泛的实用性。
【技术实现步骤摘要】
本专利技术涉及一种矿区地表形变解算方法，尤其适用于一种融合像元偏移跟踪和短基线集的矿区地表形变解算方法。技术背景合成孔径雷达干涉测量技术(InSAR，Interferometric Synthetic Aperture Radar)具有全天候、覆盖范围大、精度高、具有穿透力等优点，在地形测绘、形变监测、地表参数反演等领域的应用越来越广。DInSAR(Differential InSAR)主要用于地表形变监测，在火山监测、地面沉降、地震形变场获取、滑坡监测等领域具有极大的应用潜能。但该技术受到时间基线、空间基线、大气延迟等因素的影响，仅对部分高相干点或永久散射点进行处理分析的时序InSAR方法应运而生。该类方法主要包括：PS-InSAR(Permanent Scatterer-InSAR)、SBAS(Small Baseline Subsets)、IPTA(Interferometric Point Target Analysis)等。时序InSAR处理技术中，因主影像不唯一，时空基线较短，因此相干性较高的短基线集技术(SBAS)应用较广。但该方法在解算时首先要对差分相位进行解缠，而矿区开采沉陷速度快、形变量大，相位解缠只能解算相邻点形变梯度在(-π,π]间的沉降值，大形变区域的地表下沉难以通过相位解缠方法正确得到，导致这些区域监测量级较小，难以满足矿区大形变监测要求。同时，像元偏移跟踪算法能够解算地表大形变，目前主要应用于冰川移动、地震形变等，但利用SAR幅度信息得到的地表形变精度较低，并且非常耗时，解算结果中的地形误差与基于差分相位建立的短基线集技术中的地形误差并不等价，无法将其带入短基线集模型进行直接解算。因此，目前没有将像元偏移跟踪和短基线集融合解算地表形变和高程误差的方法。
技术实现思路
针对上述技术问题，提供一种克服了时序InSAR方法缺点，计算量小，精度高的融合像元偏移跟踪和短基线集的矿区地表形变解算方法。为实现上述技术目的，本专利技术的融合像元偏移跟踪和短基线集的矿区地表形变解算方法步骤如下：步骤1：预估矿区地表形变区域：选取被测因煤炭开采导致的矿区，利用卫星雷达获取选取矿区的SAR影像，利用DInSAR或开采沉陷预计得到的地表沉降结果，将SAR影像区域中相邻像元形变差大于SAR影像波长一半的下沉盆地初步划分成大形变区域，将像元形变量小于10mm设定为无形变区域，其余区域设定为小形变区域；步骤2：选取时间序列上相干图中的高相干点：设定SAR影像的时间和空间基线阈值，选取小于阈值的SAR影像进行干涉处理，再依据相干性获取整个或大部分时间序列上具有高相干性的点位；步骤3：利用时序偏移量跟踪算法恢复大形变区域高相干点的相位整周数N：在划分出的无形变区域内选取相干性大于0.3的高相干点计算整体偏移量，利用最小二乘原理建立偏移量的双线性多项式内插函数，将形变区域内相干点的影像坐标带入内插函数计算形变区域内相干点的整体偏移量；利用小窗口最大频谱配准对大形变区域的高相干点进行精确配准，获取局部区域偏移量；从局部精确偏移量中去除整体偏移量，并去除未下沉区域因地形起伏引起的残差，从而获取最终的视线向像元偏移量offset；利用视线向像元偏移量offset得到大形变区域高相干点的相位整周数N，利用相位整周数N和大形变区域内高相干点的原差分缠绕相位联合恢复大形变区域高相干点的真实相位差δΦb；步骤4：利用大形变区域和小形变区域的高相干点联合建立短基线集解算模型：在随机时间点t1、随机时间点t2分别采集了两幅被测矿区的SAR影像信息，判断步骤2选取的相干点位于SAR影像的小形变和无形变区域时，通过传统最小费用流方法解缠差分干涉图中高相干点(l,m)的真实相位差δΦs；判断高相干点(l,m)位于大形变区域时，不能采用传统最小费用流方法解缠差分干涉图中高相干点(l,m)真实相位差，此时则由步骤3获取高相干点(l,m)的真实相位差δΦb；最终，利用δΦs、δΦb联合构建短基线集的解算模型；步骤5，解算矿区地表高程误差及形变速率：利用短基线集处理模型，采用最小二乘解算矿区高相干点处地表高程误差及形变速率，在获取矿区地表形变后，通过反演矿区概率积分法预计参数，最终获取矿区地表沉陷预计及分析模型，为矿区沉陷工程应用提供基础数据。所述未下沉区域因地形起伏引起的残差利用公式：计算得到；式中，offsettopo为地形起伏引起的偏移量误差；a0，a1，a2为最小二乘的拟合系数；H为对应位置的高程；为研究区域的平均高程；利用公式：N＝ceil(2offset·Δ/λ)，计算得到大形变区域高相干点的相位整周数N，式中，ceil()为向下取整运算，△为斜距向像元尺寸，λ为SAR影像波长；从而克服像元偏移跟踪算法求解的偏移量仅用到SAR影像幅度信息，不能同传统短基线技术得到的差分相位进行最小二乘建模的问题。利用公式：δΦb＝2Nπ+δφdiff，计算大形变区域高相干点的真实相位差，式中，δφdiff为缠绕的差分相位；利用真实相位差δΦb建立短基线集处理模型A方法为：首先利用真实相位差δΦb构建解缠相位矩阵δΦ＝[δΦsδΦb]，再利用公式：建立短基线集处理模型；式中：Δhl,m为相干点(l,m)的高程误差；V为相干点的下沉速率；W为残余相位；λ为卫星雷达波长；B⊥,l,m为两幅SAR影像相干点处的垂直基线；θ(l,m)为相干点处的雷达入射角；rl,m为卫星距离点目标(l,m)的斜距；局部区域偏移量还可以通过最小二乘匹配方法对大形变区域的高相干点进行精确配准获取；所述矿区高相干点处地表高程误差及形变速率利用公式：计算，式中P为权阵，可取单位阵。有益效果：本专利技术将采用幅度获取地表大形变的时序偏移量跟踪算法和利用相位获取地表微小形变的时序短基线算法有机融合，将两者纳入到一个模型中统一计算，利用两种方法的互补性，通过利用真实相位差建立短基线集处理模型解决了两者无法统一建模计算和传统时序InSAR方法无法正确获取矿区大形变梯度下的地表沉降问题，改变了原有短基线集算法中只能解算地表微小形变的方式，现有时序偏移量方法不能有效解算地形误差的问题也得到了的解决。方法简单，监测精度高、范围大，费用低，对指导矿区生产、预警地质灾害、治理生态环境、维护矿群关系等具有重要的实际意义和应用价值。附图说明图1为本专利技术实施例一种融合像元偏移跟踪和短基线集的矿区地表形变解算方法流程图图2为本专利技术实施例采用相干性阈值选取的高相干点分布图图3为本专利技术实施例解算的相干点处地表年形变速率图图4为本专利技术实施例估计的相干点处DEM误差分布图具体实施方式下面将结合图和具体实施过程对本专利技术做进一步详细说明：如图1所示，本专利技术的融合像元偏移跟踪和短基线集的矿区地表形变解算方法，包括如下步骤：步骤1：矿区地表形变区域预估：本实例选取被测因煤炭开采导致的矿区，采用的是裁切后的利用卫星获取选取矿区的11景1m分辨率的Terrasar-X卫星影像，范围为2000×1600像元。利用DInSAR或开采沉陷预计得到的地表沉降结果，将SAR影像区域中相邻像元形变差大于SAR影像波长一半的下沉盆地初步划分成大形变区域，将像元形变量小于10mm设定为无形变区域，其余区域设定为小形变区域；步骤2：本文档来自技高网...
【技术保护点】
一种融合像元偏移跟踪和短基线集的矿区地表形变解算方法，其特征在于步骤如下：步骤1：预估矿区地表形变区域：选取被测因煤炭开采导致的矿区，利用卫星雷达获取选取矿区的SAR影像，利用DInSAR或开采沉陷预计得到的地表沉降结果，将SAR影像区域中相邻像元形变差大于SAR影像波长一半的下沉盆地初步划分成大形变区域，将像元形变量小于10mm设定为无形变区域，其余区域设定为小形变区域；步骤2：选取时间序列上相干图中的高相干点：设定SAR影像的时间和空间基线阈值，选取小于阈值的SAR影像进行干涉处理，再依据相干性获取整个或大部分时间序列上具有高相干性的点位；步骤3：利用时序偏移量跟踪算法恢复大形变区域高相干点的相位整周数N：在划分出的无形变区域内选取相干性大于0.3的高相干点计算整体偏移量，利用最小二乘原理建立偏移量的双线性多项式内插函数，将形变区域内相干点的影像坐标带入内插函数计算形变区域内相干点的整体偏移量；利用小窗口最大频谱配准对大形变区域的高相干点进行精确配准，获取局部区域偏移量；从局部精确偏移量中去除整体偏移量，并去除未下沉区域因地形起伏引起的残差，从而获取最终的视线向像元偏移量offset；利用视线向像元偏移量offset得到大形变区域高相干点的相位整周数N，利用相位整周数N和大形变区域内高相干点的原差分缠绕相位联合恢复大形变区域高相干点的真实相位差δΦb；步骤4：利用大形变区域和小形变区域的高相干点联合建立短基线集解算模型：在随机时间点t1、随机时间点t2分别采集了两幅被测矿区的SAR影像信息，判断步骤2选取的相干点位于SAR影像的小形变和无形变区域时，通过传统最小费用流方法解缠差分干涉图中高相干点(l,m)的真实相位差δΦs；判断高相干点(l,m)位于大形变区域时，不能采用传统最小费用流方法解缠差分干涉图中高相干点(l,m)真实相位差，此时则由步骤3获取高相干点(l,m)的真实相位差δΦb；最终，利用δΦs、δΦb联合构建短基线集的解算模型；步骤5，解算矿区地表高程误差及形变速率：利用短基线集处理模型，采用最小二乘解算矿区高相干点处地表高程误差及形变速率，在获取矿区地表形变后，通过反演矿区概率积分法预计参数，最终获取矿区地表沉陷预计及分析模型，为矿区沉陷工程应用提供基础数据。...
【技术特征摘要】
1.一种融合像元偏移跟踪和短基线集的矿区地表形变解算方法，其特征在于步骤如下：步骤1：预估矿区地表形变区域：选取被测因煤炭开采导致的矿区，利用卫星雷达获取选取矿区的SAR影像，利用DInSAR或开采沉陷预计得到的地表沉降结果，将SAR影像区域中相邻像元形变差大于SAR影像波长一半的下沉盆地初步划分成大形变区域，将像元形变量小于10mm设定为无形变区域，其余区域设定为小形变区域；步骤2：选取时间序列上相干图中的高相干点：设定SAR影像的时间和空间基线阈值，选取小于阈值的SAR影像进行干涉处理，再依据相干性获取整个或大部分时间序列上具有高相干性的点位；步骤3：利用时序偏移量跟踪算法恢复大形变区域高相干点的相位整周数N：在划分出的无形变区域内选取相干性大于0.3的高相干点计算整体偏移量，利用最小二乘原理建立偏移量的双线性多项式内插函数，将形变区域内相干点的影像坐标带入内插函数计算形变区域内相干点的整体偏移量；利用小窗口最大频谱配准对大形变区域的高相干点进行精确配准，获取局部区域偏移量；从局部精确偏移量中去除整体偏移量，并去除未下沉区域因地形起伏引起的残差，从而获取最终的视线向像元偏移量offset；利用视线向像元偏移量offset得到大形变区域高相干点的相位整周数N，利用相位整周数N和大形变区域内高相干点的原差分缠绕相位联合恢复大形变区域高相干点的真实相位差δΦb；步骤4：利用大形变区域和小形变区域的高相干点联合建立短基线集解算模型：在随机时间点t1、随机时间点t2分别采集了两幅被测矿区的SAR影像信息，判断步骤2选取的相干点位于SAR影像的小形变和无形变区域时，通过传统最小费用流方法解缠差分干涉图中高相干点(l,m)的真实相位差δΦs；判断高相干点(l,m)位于大形变区域时，不能采用传统最小费用流方法解缠差分干涉图中高相干点(l,m)真实相位差，此时则由步骤3获取高相干点(l,m)的真实相位差δΦb；最终，利用δΦs、δΦb联合构建短基线集的解算模型；步骤5，解算矿区地表高程误差及形变速率：利用短基线集处理模型，采用最小二乘解算矿区高相干点处地表高程误差...
【专利技术属性】
技术研发人员：范洪冬，杜森，黄继磊，闫世勇，邓喀中，汪云甲，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。