CN110568440B - 一种基于DS-InSAR技术监测复杂山区形变的方法 - Google Patents

一种基于DS-InSAR技术监测复杂山区形变的方法 Download PDF

Info

Publication number
CN110568440B
CN110568440B CN201910854798.9A CN201910854798A CN110568440B CN 110568440 B CN110568440 B CN 110568440B CN 201910854798 A CN201910854798 A CN 201910854798A CN 110568440 B CN110568440 B CN 110568440B
Authority
CN
China
Prior art keywords
deformation
monitoring
area
research area
deformation rate
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910854798.9A
Other languages
English (en)
Other versions
CN110568440A (zh
Inventor
赵雪
薛廉
唐侨
潘国耀
王向东
吴森
陈超
齐得旭
谭超
陈渠波
王蕾
肖长波
袁伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Geological Engineering Exploration Institute Group Co ltd
Original Assignee
Sichuan Geological Engineering Exploration Institute Group Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan Geological Engineering Exploration Institute Group Co ltd filed Critical Sichuan Geological Engineering Exploration Institute Group Co ltd
Priority to CN201910854798.9A priority Critical patent/CN110568440B/zh
Publication of CN110568440A publication Critical patent/CN110568440A/zh
Application granted granted Critical
Publication of CN110568440B publication Critical patent/CN110568440B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904—SAR modes
    • G01S13/9058—Bistatic or multistatic SAR

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明公开了一种基于DS‑InSAR技术监测复杂山区形变的方法，包括：A.选取覆盖研究区域的多卫星平台多轨道的SAR影像数据；B.设置形变速率阈值，筛选出形变速率大于形变速率阈值的研究区域作为高形变速率区；C.在高形变速率区内选取典型监测点，根据典型监测点在每个监测时间点间的形变，绘制该监测点在整个监测时间段内的形变曲线；D.利用形变曲线计算研究区域的形变切线角；E.结合形变速率大小和形变切线角大小对研究区域的风险等级进行评价。本发明的方法可以在山区获得高密度的监测点，使得可靠数据量增大，极大的弥补了PS‑InSAR技术在复杂山区难以探测到PS点的问题，可以实现对复杂山区的灾害体进行初步识别和判断，提高排查效率。

Description

一种基于DS-InSAR技术监测复杂山区形变的方法
技术领域
本发明涉及地质灾害形变监测技术领域，特别涉及一种基于DS-InSAR技术监测复杂山区形变的方法。
背景技术
地质灾害，包括自然因素或者人为活动引发的危害人民生命和财产安全的山体崩塌、滑坡、泥石流、地面塌陷、地裂缝、地面沉降等与地质作用有关的灾害。地质灾害常常伴随着地表位移，常用的地表形变监测手段中，传统的地质灾害形变监测的方法为水准测量存在费时费力成本高、周期长、数据分布稀疏等问题，而且许多偏远地区难以到达，GPS测量虽然精度高，但是数据点稀疏，难以全面覆盖研究区域。
InSAR技术具有全天时、全天候、高精度的特点，能够快速获取面状、大范围、长时间序列的地表形变，为地质灾害体的判断和监测提供定量地表位移数据。现在常用的PS-InSAR技术可以对永久散射体(PS)进行监测，永久散射体指对雷达波的后向散射较强，在时序上较稳定的各种地物目标，如建筑物与构筑物的顶角、桥梁、栏杆、裸露的岩石以及人工布设的角反射器等目标。由于地质灾害多发生在植被茂密的山区，PS点较少，使用PS-InSAR技术难以获取大范围的形变数据。
发明内容
本发明的目的是克服上述背景技术中不足，提供一种基于DS-InSAR技术监测复杂山区形变的方法，采用DS-InSAR技术对地质灾害形变进行监测，从而获得高密度的野外监测点时序形变信息，通过使用多个平台和多个轨道的SAR影像数据可以获得更多维度的形变信息，使得形变信息更加全面，且减少野外工作量，降低监测成本，直接获取大范围、高密度、长时序的形变数据，最后根据监测点的形变曲线和形变变化趋势，对灾害体进行初步识别和判断，极大的弥补了PS-InSAR技术在复杂山区难以探测到PS点的问题。
为了达到上述的技术效果，本发明采取以下技术方案：
一种基于DS-InSAR技术监测复杂山区形变的方法，包括以下步骤：
A.选取覆盖研究区域的多卫星平台多轨道的SAR影像数据，根据监测区域将数据裁剪到监测的研究区域范围，并使用SqueeSAR算法获取研究区域范围内永久散射体PS点和分布式散射体DS点的形变速率以及每个监测时间之间的形变；
B.设置形变速率阈值，筛选出形变速率大于形变速率阈值的研究区域作为高形变速率区；其中，形变速率阈值的具体设定可以根据实际研究情况来设定；
C.在高形变速率区内选取典型监测点，根据典型监测点在每个监测时间点间的形变，绘制该监测点在整个监测时间段内的形变曲线；其中，典型监测点的选取具体是根据实际需要研究的点来确定，可直接选取需要研究的点作为典型监测点；
D.利用形变曲线计算研究区域的形变切线角；
E.结合形变速率大小和形变切线角大小对研究区域的风险等级进行评价，其中，形变速率越大，切线角越大，则研究区域地质风险性越高。
进一步地，所述步骤B中还包括将不同卫星平台和轨道监测到的高形变速率区作并集，获取到完整的高形变速率区分布情况。
进一步地，所述步骤D中所述形变切线角αi的计算公式如下：
Figure GDA0002533128800000021
其中，ti表示某一监测时刻，Δt表示单位时间段，S(i)表示ti时刻的形变量，ΔS表示单位时间段内S(i)的形变量。
进一步地，所述一个单位时间段为一个监测周期，如1周、1月等。
进一步地，当研究区域的形变切线角小于15°时则判定研究区域的风险等级较低，当研究区域的形变切线角大于等于85°时，则判定研究区域的风险等级较高。
本发明与现有技术相比，具有以下的有益效果：
本发明的基于DS-InSAR技术监测复杂山区形变的方法，使用DS-InSAR技术识别分布式散射体，分布式散射体一般为自然地面目标，利用DS-InSAR技术来监测复杂山区地表形变，可以在山区获得高密度的监测点，使得可靠数据量增大，极大的弥补了PS-InSAR技术在复杂山区难以探测到PS点的问题，可以实现对复杂山区的灾害体进行初步识别和判断，从而有导向性的对风险性高的区域进行重点排查，提高排查效率。
附图说明
图1是本发明的一个实施例的基于DS-InSAR技术监测复杂山区形变的方法的流程示意图。
图2是本发明的一个实施例中绘制的形变曲线示意图。
具体实施方式
下面结合本发明的实施例对本发明作进一步的阐述和说明。
实施例：
实施例一：
如图1所示，一种基于DS-InSAR技术监测复杂山区形变的方法，包括以下步骤：
步骤1.选取覆盖研究区域的多卫星平台(如Sentinel-1A、ALOS-2、TerraSAR-X等)多轨道(如包括升轨和降轨)的SAR影像数据，根据监测区域将数据裁剪到监测的研究区域范围，并使用SqueeSAR算法获取研究区域范围内永久散射体PS点和分布式散射体DS点的形变速率以及每个监测时间之间的形变。
具体的，SqueeSAR算法为现有技术中的一种算法，是在PS-InSAR算法基础上发展而来，PS-InSAR算法仅解算PS点的形变，SqueeSAR算法解算PS点和DS点的形变，SqueeSAR算法使用KS检验判定DS点，将判断出的DS点与PS点一起放入PS-InSAR处理流程中进行解算，得到研究区域的形变。
在影像上，DS点对应于具有相同分布特征的多个像素，这些像素点属于同一地物，地物的相干性为中等大小。SqueeSAR算法使用KS检验识别属于干涉数据堆栈的2个图像像素是否在统计上式均匀的，从而提取DS点。
具体的，KS检验判别像素P1、P2是否为统计均匀像素，基于2个数据集分布函数作差绝对值的最大值，分析DS点的随机变量在重新参数化下是稳定的。
其中，KS检验判别算式可表示为：
Figure GDA0002533128800000041
上式中：
Figure GDA0002533128800000042
和
Figure GDA0002533128800000043
分别为P1、P2的分布函数。DN的概率分布函数可用经验KS分布近似，可表示为：
Figure GDA0002533128800000044
KS检验通过判断DN是否大于设定阈值的概率来确定P1、P2是否符合同一分布。
步骤2.不同平台和不同轨道的数据监测到不同方向的形变和形变速率，设置形变速率阈值，筛选出形变速率大于形变速率阈值的研究区域作为高形变速率区；并将不同卫星平台和轨道监测到的高形变速率区作并集，获取到完整的高形变速率区分布情况。其中，形变速率阈值的具体设定可以根据实际研究情况来设定，本实施例中设定的形变速率阈值为20mm/year。
步骤3.在高形变速率区内选取典型监测点，根据典型监测点在每个监测时间点间的形变，绘制该监测点在整个监测时间段内的形变曲线；其中，典型监测点的选取具体是根据实际需要研究的点来确定，可直接选取需要研究的点作为典型监测点，如图2所示为本实施例中根据一个典型监测点在每个监测时间点间的形变，绘制的该监测点的形变曲线。
步骤4.利用形变曲线计算研究区域的形变切线角。
从图2中可看出，初始变形阶段，切线角总体表现为由大减小，逐渐平稳；当进入等速变形阶段，切线角基本稳定在45°；在加速变形阶段，切线角增长速率急剧扩大，直至90°。
具体的，图2中的形变曲线对应的形变切线角αi的计算公式如下：
Figure GDA0002533128800000051
其中，ti表示某一监测时刻，Δt表示单位时间段，S(i)表示ti时刻的形变量，ΔS表示单位时间段内S(i)的形变量。
其中，一个单位时间段为一个监测周期，如1周、1月等。
步骤5.结合形变速率大小和形变切线角大小对研究区域的风险等级进行评价，其中，形变速率越大，切线角越大，则研究区域地质风险性越高，为实地地质灾害勘察提供连续时间序列的定量地表形变数据。
本实施例中设定的在相同形变速率下，不同形变切线角对应的风险等级如下表所示：
Figure GDA0002533128800000052
本实施例中仅示意出了一种划分规则，实际中，也可根据具体情况设定各等级对应的角度值。
可以理解的是，以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式，然而本发明并不局限于此。对于本领域内的普通技术人员而言，在不脱离本发明的精神和实质的情况下，可以做出各种变型和改进，这些变型和改进也视为本发明的保护范围。

Claims (4)

1.一种基于DS-InSAR技术监测复杂山区形变的方法，其特征在于，包括以下步骤：
A.选取覆盖研究区域的多卫星平台多轨道的SAR影像数据，根据监测区域将数据裁剪到监测的研究区域范围，并使用SqueeSAR算法获取研究区域范围内永久散射体PS点和分布式散射体DS点的形变速率以及每个监测时间之间的形变；
B.设置形变速率阈值，筛选出形变速率大于形变速率阈值的研究区域作为高形变速率区，将不同卫星平台和轨道监测到的高形变速率区作并集，获取到完整的高形变速率区分布情况；
C.在高形变速率区内选取典型监测点，根据典型监测点在每个监测时间点间的形变，绘制该监测点在整个监测时间段内的形变曲线；
D.利用形变曲线计算研究区域的形变切线角；
E.结合形变速率大小和形变切线角大小对研究区域的风险等级进行评价，其中，形变速率越大，切线角越大，则研究区域地质风险性越高。
2.根据权利要求1所述的一种基于DS-InSAR技术监测复杂山区形变的方法，其特征在于，步骤D中所述形变切线角αi的计算公式如下：
Figure FDA0002533128790000011
其中，ti表示某一监测时刻，Δt表示单位时间段，S(i)表示ti时刻的形变量，ΔS表示单位时间段内S(i)的形变量。
3.根据权利要求2所述的一种基于DS-InSAR技术监测复杂山区形变的方法，其特征在于，一个单位时间段为一个监测周期。
4.根据权利要求1至3中任一所述的一种基于DS-InSAR技术监测复杂山区形变的方法，其特征在于，当研究区域的形变切线角小于15°时则判定研究区域的风险等级较低，当研究区域的形变切线角大于等于85°时，则判定研究区域的风险等级较高。
CN201910854798.9A 2019-09-10 2019-09-10 一种基于DS-InSAR技术监测复杂山区形变的方法 Active CN110568440B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910854798.9A CN110568440B (zh) 2019-09-10 2019-09-10 一种基于DS-InSAR技术监测复杂山区形变的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910854798.9A CN110568440B (zh) 2019-09-10 2019-09-10 一种基于DS-InSAR技术监测复杂山区形变的方法

Publications (2)

Publication Number Publication Date
CN110568440A CN110568440A (zh) 2019-12-13
CN110568440B true CN110568440B (zh) 2020-09-15

Family

ID=68779046

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910854798.9A Active CN110568440B (zh) 2019-09-10 2019-09-10 一种基于DS-InSAR技术监测复杂山区形变的方法

Country Status (1)

Country Link
CN (1) CN110568440B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112014841A (zh) * 2020-08-31 2020-12-01 中国矿业大学 一种基于DS-InSAR技术监测油田区地表形变分析方法
CN112269176B (zh) * 2020-10-14 2021-09-14 武汉工程大学 一种矿山地表沉陷早期识别监测方法
CN112526515A (zh) * 2020-11-05 2021-03-19 山西省交通科技研发有限公司 一种基于合成孔径雷达干涉测量技术的地表形变检测方法
CN113253271A (zh) * 2021-07-12 2021-08-13 中国测绘科学研究院 一种形变速率估计方法和系统
CN115585785B (zh) * 2022-10-26 2024-01-23 四川省公路规划勘察设计研究院有限公司 城市道路沉降InSAR评价方法
CN115638767B (zh) * 2022-11-07 2023-10-03 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) 地面沉降监测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105938193A (zh) * 2016-07-14 2016-09-14 中南大学 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法
CN106023157A (zh) * 2016-05-10 2016-10-12 电子科技大学 一种基于sar图像的山区地表微形变信息提取方法
CN106772377A (zh) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 一种基于InSAR的建筑物变形监测方法
CN109029344A (zh) * 2018-07-10 2018-12-18 湖南中科星图信息技术有限公司 一种基于高分影像和升降轨InSAR的堤坝沉降监测方法
CN109884635A (zh) * 2019-03-20 2019-06-14 中南大学 大范围高精度的InSAR形变监测数据处理方法
CN110058237A (zh) * 2019-05-22 2019-07-26 中南大学 面向高分辨率SAR影像的InSAR点云融合及三维形变监测方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102013150B (zh) * 2010-09-28 2013-03-27 浙江工业大学 基于雨强、斜坡土壤含水量和形变量的地质灾害预测系统
JP6349937B2 (ja) * 2014-05-09 2018-07-04 日本電気株式会社 変動検出装置、変動検出方法および変動検出用プログラム
CN104111457B (zh) * 2014-07-23 2016-10-12 中国国土资源航空物探遥感中心 一种升降轨PSInSAR地面沉降监测结果的互检验与时序融合方法
CN109031301A (zh) * 2018-09-26 2018-12-18 云南电网有限责任公司电力科学研究院 基于PSInSAR技术的山区地形形变提取方法
CN109238227B (zh) * 2018-10-31 2020-10-20 首都师范大学 一种表征地面沉降时空演变的方法
CN109991601A (zh) * 2018-11-29 2019-07-09 北京东方至远科技股份有限公司 一种基于PS-InSAR技术的房屋风险评估方法
CN109709550B (zh) * 2019-01-17 2020-10-30 武汉大学 一种基于InSAR影像数据的库岸边坡形变监测处理方法
CN109752715B (zh) * 2019-01-24 2023-05-12 深圳市数字城市工程研究中心 一种sar数据全散射体探测方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106023157A (zh) * 2016-05-10 2016-10-12 电子科技大学 一种基于sar图像的山区地表微形变信息提取方法
CN105938193A (zh) * 2016-07-14 2016-09-14 中南大学 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变的方法
CN106772377A (zh) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 一种基于InSAR的建筑物变形监测方法
CN109029344A (zh) * 2018-07-10 2018-12-18 湖南中科星图信息技术有限公司 一种基于高分影像和升降轨InSAR的堤坝沉降监测方法
CN109884635A (zh) * 2019-03-20 2019-06-14 中南大学 大范围高精度的InSAR形变监测数据处理方法
CN110058237A (zh) * 2019-05-22 2019-07-26 中南大学 面向高分辨率SAR影像的InSAR点云融合及三维形变监测方法

Also Published As

Publication number Publication date
CN110568440A (zh) 2019-12-13

Similar Documents

Publication Publication Date Title
CN110568440B (zh) 一种基于DS-InSAR技术监测复杂山区形变的方法
Massey et al. Landslides triggered by the 14 November 2016 Mw 7.8 Kaikōura earthquake, New Zealand
CN106023157B (zh) 一种基于sar图像的山区地表微形变信息提取方法
Li et al. Geological survey and unstable rock block movement monitoring of a post-earthquake high rock slope using terrestrial laser scanning
Mondini et al. Comparison of event landslide inventories: the Pogliaschina catchment test case, Italy
Trevisani et al. Surface texture analysis of a high-resolution DTM: Interpreting an alpine basin
CN102819023A (zh) 基于LiDAR的复杂地质背景区滑坡识别的方法及系统
Uprety et al. Damage detection using high-resolution SAR imagery in the 2009 L'Aquila, Italy, Earthquake
Wang et al. Automatic detection and update of landslide inventory before and after impoundments at the Lianghekou reservoir using Sentinel-1 InSAR
Song et al. Efficient water area classification using Radarsat-1 SAR imagery in a high relief mountainous environment
CN110988909A (zh) 基于tls进行高寒脆弱区沙地植被的植被盖度测定方法
Regard et al. Coastal chalk platform erosion modulated by step erosion and debris shielding: example from Normandy and Picardy (northern France)
CN115077489A (zh) 一种无人机倾斜摄影土方量计算的方法
CN109190593B (zh) 基于凹凸类别划分的山区道路沿线边坡稳定性初步判别方法
He et al. [Retracted] Integration and Fusion of Geologic Hazard Data under Deep Learning and Big Data Analysis Technology
CN116912068B (zh) 一种基于面域形变观测的滑坡预警方法
Oorthuis et al. Monitoring of rainfall and soil moisture at the Rebaixader catchment (Central Pyrenees)
Beck et al. Automated extraction of forest road network geometry from aerial LiDAR
CN112052274A (zh) 一种岩质地层地铁盾构施工地表沉降规律的数据挖掘方法
Zhang et al. Application of airborne LiDAR measurements to the topographic survey of the tidal flats of the Northern Jiangsu radial sand ridges in the Southern Yellow Sea
CN116721358A (zh) 一种InSAR地形产品建筑物提取方法
Gonçalves-Seco et al. Digital terrain model generation using airborne LiDAR in a forested area Galicia, Spain
Saepuloh et al. Detecting surface structures after large eruption of Mt. Merapi in 2010 using ALOS/PALSAR data
CN112907567A (zh) 基于空间推理方法的sar图像有序人造构筑物提取方法
Politz et al. Building change detection of airborne laser scanning and dense image matching point clouds using height and class information

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant