分析浅地表空间利用差异模式下地面沉降演变规律的方法技术

本发明专利技术公开一种分析浅地表空间利用差异模式下地面沉降演变规律的方法，对设定期间内反映地表形变信息的图像序列进行差分干涉处理，并对处理得到的多幅差分干涉相位图像进行光谱滤波得到SBAS数据集和PS数据集；采用空间相关组分、空间非相关视角误差评估方法，从得到的所述SBAS数据集和所述PS数据集中选择误差小于设定值的一组SBAS点和PS点；对所选的一组SBAS点和PS点进行数据集融合处理，并对融合后的数据集进行相位解缠，通过奇异值分解及范数约束提取像元形变信息，得到相对于参考区域的主图像的时间序列形变信息。

【技术实现步骤摘要】

本专利技术涉及InSAR(Synthetic Aperture Radar Interferometry，合成孔径雷达 干涉)技术下的地面沉降监测及时序演化规律领域，具体而言，涉及一种分析浅地表空间利 用差异模式下地面沉降演变规律的方法。
技术介绍
相对于水准测量、分层标测量等常规监测方法，InSAR技术可以实时动态获取大范 围、高精度的地表形变细节信息。永久散射体干涉测量(Permanent Scatterer for SAR Interferometry，PS_InSAR)技术、小基线集雷达干涉测量(Small Baseline-InSAR，SBAS) 是目前具有代表性的时序InSAR监测技术，克服了传统D-InSAR技术空间、时间去相干影响， 减弱大气延迟带来的误差组分，提高形变监测信息的时、空分辨率及数据处理的精度，在城 区或植被覆盖稀少地区可以达到毫米级的地表形变监测精度。国内外众多专家学者不断改 进该两种方法，并取得了成功的应用。H 〇〇per(2007)提出一种新的融合小基线（SBAS)和PS 技术监测地表形变的算法，证明了小基线（SBAS)技术和PS方法虽然用不同的离散模型，但 二者存在优势互补之处，融合两者数据集不仅能够增加相干像元数量，而且提高了像元的 信噪比，大幅度减小了解缠空间的偏差，这一优势是单一小基线技术或PS技术都无法比拟 的。 在时序InSAR技术获取地面沉降监测信息基础上，借助地理信息系统强大的空间 数据管理和分析能力，可以更全面地揭示地面沉降时空演化过程。因此，时序InSAR与GIS技 术的集成应用，是地面沉降监测和演化研究的重要发展方向。
技术实现思路
本专利技术提供一种，用以 研究浅地表空间的地面沉降演变规律。 为达到上述目的，本专利技术提供了一种分析浅地表空间利用差异模式下地面沉降演 变规律的方法包括以下步骤： 对设定期间内反映地表形变信息的图像序列进行差分干涉处理，并对处理得到的 多幅差分干涉相位图像进行光谱滤波得到SBAS数据集和PS数据集； 采用空间相关组分、空间非相关视角误差评估方法，从得到的所述SBAS数据集和 所述PS数据集中选择误差小于设定值的一组SBAS点和PSA; 对所选的一组SBAS点和PS点进行数据集融合处理，并对融合后的数据集进行相位 解缠，通过奇异值分解及范数约束提取像元形变信息，得到相对于参考区域的主图像的时 间序列形变信息。 进一步地，上述还包括 以下步骤： 采用设定大小的移动窗口从反映时间序列形变信息的主图像上选择多个典型区 域，对于所选的每个典型区域，采用GIS空间分析方法，根据该典型区域的时间序列形变信 息获取该典型区域的时间序列的不均匀沉降的演化规律。 进一步地，上述还包括 以下步骤： 采用GIS空间分析方法，根据多个典型区域中每个典型区域的时间序列形变信息 获取多个典型区域内的时间序列的不均匀沉降的演化规律。 进一步地，所述多个典型区域为5个。进一步地，所述移动窗口为6km2大小的正方形。 进一步地，该典型区域的时间序列的不均匀沉降的演化规律包括年际地面沉降 量。 进一步地，多个典型区域内的时间序列的不均匀沉降的演化规律包括地面沉降速 率极值和年平均地面沉降速率。 本专利技术采用的永久散射体干涉测量（Permanent Scatterer for SAR Interferometry，PS_InSAR)技术、小基线集雷达干涉测量(Small Baseline-InSAR,SBAS) 是目前具有代表性的时序InSAR监测技术，克服了传统D-InSAR技术空间、时间去相干影响， 减弱大气延迟带来的误差组分，提高形变监测信息的时、空分辨率及数据处理的精度，在城 区或植被覆盖稀少地区可以达到毫米级的地表形变监测精度。同时，结合GIS空间分析方 法，在区域浅表层空间（地铁、城市密集建筑群(CBD)、立体交通网络设施)不同的变异模式 下，选取多个典型地面沉降区域，结合土地利用分类信息、可见光遥感影像、地质资料等，分 析时间序列的不均匀沉降的演化规律。【附图说明】 为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。 图1为本专利技术一个实施例的融合PS和小基线InSAR技术数据处理流程图； 图2为本专利技术一个实施例的5个典型区域不均匀地面沉降演变过程分析图； 图3为5个典型的移动窗口的位置分布及沉降速率情况图； 图4为典型区域1的沉降分布(左)与年际形变演化过程(右)图； 图5为典型区域2的沉降分布(左)与年际形变演化过程(右)图； 图6为典型区域3的沉降分布(左)与年际形变演化过程(右)图； 图7为典型区域4的沉降分布(左)与年际形变演化过程(右)图； 图8为典型区域5的沉降分布(左)与年际形变演化过程(右)图； 图9为5个典型区域PS点和小基线点不同沉降速率比例分布图； 图10为5个典型区域相干点沉降速率极值比较分布图。【具体实施方式】下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于 本专利技术中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他 实施例，都属于本专利技术保护的范围。 图1为本专利技术一个实施例的融合PS和小基线InSAR技术数据处理流程图；图2为本 专利技术一个实施例的5个典型区域不均匀地面沉降演变过程分析图；如图所示，分析浅地表空 间利用差异模式下地面沉降演变规律的方法包括以下步骤： 对设定期间内反映地表形变信息的图像序列进行差分干涉处理，并对处理得到的 多幅差分干涉相位图像进行光谱滤波得到SBAS数据集和PS数据集； 采用空间相关组分、空间非相关视角误差评估方法，从得到的所述SBAS数据集和 所述PS数据集中选择误差小于设定值的一组SBAS点和PS点；从得到的所述SBAS数据集和所 述PS数据集选择设定数量的SBAS点和PS点，并对所选的SBAS点和PS点进行空间相关组分、 空间非相关视角误差评估； 对所选的一组SBAS点和PS点进行数据集融合处理，并对融合后的数据集进行相位 解缠，通过奇异值分解及范数约束提取像元形变信息，得到相对于参考区域的主图像的时 间序列形变信息。 进一步地，上述还包括 以下步骤：采用设定大小的移动窗口（如6km2大小的正方形）从反映时间序列形变信息的主 图像上选择多个(如5个)典型区域，对于所选的每个典型区域，采用GIS空间分析方法，根据 该典型区域的时间序列形变信息获取该典型区域的时间序列的不均匀沉降的演化规律(如 包括年际地面沉降量）。 进一步地，上述还包括 以下步骤： 采用GIS空间分析方法，根据多个典型区域中每个典型区域的时间序列形变信息 获取多个典型区域内的时间序列的不均匀沉降的演化规律(如包括地面沉降速率极值和年 平均地面沉降速率）。 其中，图1、图2实施例中是选取2003~2009年的覆盖研究区的29景降轨Envisat ASAR数据，以2005年1本文档来自技高网...

【技术保护点】
一种分析浅地表空间利用差异模式下地面沉降演变规律的方法，其特征在于，包括以下步骤：对设定期间内反映地表形变信息的图像序列进行差分干涉处理，并对处理得到的多幅差分干涉相位图像进行光谱滤波得到SBAS数据集和PS数据集；采用空间相关组分、空间非相关视角误差评估方法，从得到的所述SBAS数据集和所述PS数据集中选择误差小于设定值的一组SBAS点和PS点；对所选的一组SBAS点和PS点进行数据集融合处理，并对融合后的数据集进行相位解缠，通过奇异值分解及范数约束提取像元形变信息，得到相对于参考区域的主图像的时间序列形变信息。

【技术特征摘要】

【专利技术属性】
技术研发人员：陈蓓蓓，宫辉力，李小娟，梁悦，雷坤超，赵文吉，
申请(专利权)人：首都师范大学，
类型：发明
国别省市：北京;11