CN114185048A - 地基InSAR提取滑坡位移向量的方法、系统及存储介质 - Google Patents

地基InSAR提取滑坡位移向量的方法、系统及存储介质 Download PDF

Info

Publication number
CN114185048A
CN114185048A CN202210135030.8A CN202210135030A CN114185048A CN 114185048 A CN114185048 A CN 114185048A CN 202210135030 A CN202210135030 A CN 202210135030A CN 114185048 A CN114185048 A CN 114185048A
Authority
CN
China
Prior art keywords
slope
landslide
insar
gradient
displacement
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.)
Granted
Application number
CN202210135030.8A
Other languages
English (en)
Other versions
CN114185048B (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.)
Hunan Jihe Information Technology Co ltd
Original Assignee
Hunan Jihe Information Technology 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 Hunan Jihe Information Technology Co ltd filed Critical Hunan Jihe Information Technology Co ltd
Priority to CN202210135030.8A priority Critical patent/CN114185048B/zh
Publication of CN114185048A publication Critical patent/CN114185048A/zh
Application granted granted Critical
Publication of CN114185048B publication Critical patent/CN114185048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)

Abstract

本发明提供一种地基InSAR提取滑坡位移向量的方法、系统及存储介质，其中方法包括如下步骤：首先根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 889307DEST_PATH_IMAGE001
和坡度
Figure 241790DEST_PATH_IMAGE002
；再基于单视图像的地基InSAR，测量得到雷达到目标位置的视线（LOS）方向位移
Figure 122022DEST_PATH_IMAGE003
；最后将步骤S2获得的视线（LOS）方向位移投影到步骤S1中坡向
Figure 559956DEST_PATH_IMAGE001
和坡度

Description

地基InSAR提取滑坡位移向量的方法、系统及存储介质
技术领域
本发明涉及合成孔径雷达技术领域，尤其涉及一种基于单视图像的地基InSAR沿最大坡度方向提取滑坡位移向量的方法、系统及存储介质。
背景技术
地基干涉合成孔径雷达（Interferometric Synthetic Aperture Radar，InSAR）是一种有效的边坡稳定监测遥感技术手段，相比星载、机载干涉合成孔雷达和三维激光、摄影测量等手段，具有远距离非接触监测、形变测量精度高、数据更新率高、不受气候光照影响等优点，成为灾害隐患现场持续监测和早期预警的重要装备。
但地基InSAR由于其原理性的限制，即测量得到的相位干涉差仅能获取视线（LOS）方向的位移量。而实际滑坡发生的位移方向往往与雷达LOS方向不一致，此时雷达测量到的方向实际上是实际滑坡位移向量在LOS方向的投影。因为这一局限性，地基InSAR往往用于趋势性边坡形变测量，难以对边坡形变进行定量分析。因此，地基InSAR如何获取实际边坡位移大小和方向（或位移向量）一直是研究的热点之一。
传统的获取实际边坡位移向量的方法，在条件许可下，一般可以通过多视图像，也即多个SAR轨道从不同视角对同一场景生成的图像，对同一位置处提取不同视角的位移值，再通过三角运算获得二维或三维位移值，然后合成实际位移向量。多轨图像一般通过同一平台按不同轨道运行、或不同平台按不同轨道运行获得。在地基InSAR应用中，比较典型的包括《一种基于多基地MIMO-SAR的三维形变监测系统》（CN 104849712 A）和《Multi-staticMIMO-SAR Three Dimensional Deformation Measurement System》（10.1109/APSAR.2015.7306212）提出利用三个空间上分布的MIMO雷达对同一场景进行监测，从而通过三角投影关系解算出该场景的三维形变场的方法；但是利用地基InSAR进行边坡监测，特别是对地质灾害进行应急监测的应用中，从成本和部署难度考虑，往往难以布置多个设备或一个设备在多个位置对同一场景进行监测。
在星载InSAR应用中，可以采用多孔径干涉（MAI）技术通过方位向滤波器将SAR孔径分成若干子孔径，并利用子孔径形成方位向的多视图像估计方位和距离二维位移值。但地基SAR孔径很短，即使分为若干子孔径，子孔径的视角变化非常小，实际上方位向估计精度很差；
也有研究表明，采用偏移量跟踪（OT）方法通过时序图像中相邻图像的互相关获取场景在方位和距离二维上亚分辨率级位移，一般位移测量精度可以达到分辨率的1/30到1/10。但由于地基SAR孔径很短，方位向分辨率往往不高（如5mrad方位角分辨率在1公里处对应方位向分辨率5米），此方法位移测量精度无法满足边坡稳定监测毫米级的需求；
在地下采矿引起的地表沉降、冰川移动等特殊情况下，也可以通过基于裂缝发育等先验知识对三维位移进行估计，但该方法对于大多数边坡并不适用。
因此，本领域存在极大不足，考虑到地基InSAR多视图像获取工程上的困难性，亟需一种基于单视图像的地基InSAR提取实际滑坡位移向量的方法、系统，既能够降低成本和部署难度，又能保证测量精度。
发明内容
鉴于上述现有技术的不足之处，本发明提供一种基于单视图像的地基InSAR提取实际滑坡位移向量的方法、系统及存储介质，能够通过LOS方向位移结合三维地形对实际滑坡大小和方向进行估计，且设备成本低和部署难度低，还能保证测量精度。
为了达到上述目的，本发明采取了以下技术方案：
地基InSAR提取实际滑坡位移向量的方法，包括步骤：
S1、根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 296269DEST_PATH_IMAGE001
和坡度
Figure 776929DEST_PATH_IMAGE002
；
S2、基于单视图像的地基InSAR，测量得到雷达到目标位置的视线（LOS）方向位移
Figure 272633DEST_PATH_IMAGE003
；
S3、将步骤S2获得的视线（LOS）方向位移
Figure 864151DEST_PATH_IMAGE003
投影到步骤S1中坡向
Figure 994656DEST_PATH_IMAGE001
和坡度
Figure 911796DEST_PATH_IMAGE002
所确定的最大坡度方向上，从而获得实际滑坡位移向量。
作为上述方案进一步的改进，在所述步骤S1中，提取滑坡场景的坡向和坡度时，三维地形的坡度图和坡向图利用GDAL（https://gdal.org）的gdaldem slope/aspect命令生成，或ArcGIS、QGIS相关工具软件获得。
作为上述方案进一步的改进，对获得的滑坡场景的坡度图和坡向图进行平滑滤波，如合适窗大小的均值、高斯等滤波器，以抑制小尺度坡度和坡向带来的噪声。
作为上述方案进一步的改进，在步骤S3中，
Figure 629216DEST_PATH_IMAGE004
坡向
Figure 24426DEST_PATH_IMAGE001
和坡度
Figure 838798DEST_PATH_IMAGE002
所确定的最大坡度方向的方法如下，雷达监测的目标位置的单位向量的三维坐标
Figure 598944DEST_PATH_IMAGE005
表示如下：
Figure 803660DEST_PATH_IMAGE006
则坡向
Figure 2560DEST_PATH_IMAGE001
和坡度
Figure 131095DEST_PATH_IMAGE002
所确定的最大坡度方向即为指向单位向量的三维坐标
Figure 124458DEST_PATH_IMAGE005
的方向。
作为上述方案进一步的改进，通过地理角度定义或几何角度定义来获取坡向
Figure 816471DEST_PATH_IMAGE001
；
坡向以地理角度定义顺时针增加，以正北方向为0°角度；或者坡向以几何角度定义逆时针增加，以东方向为0°角度；且地理角度和几何角度存在以下关系：
Figure 553483DEST_PATH_IMAGE007
，其中
Figure 14551DEST_PATH_IMAGE008
为目标位置的坡向地理角度，
Figure 444395DEST_PATH_IMAGE009
为目标位置的坡向几何角度。
作为上述方案进一步的改进，在步骤S2中，视线（LOS）方向的获取步骤如下：
S21：假设雷达位置的经纬高坐标为
Figure 623704DEST_PATH_IMAGE010
，其中
Figure 898828DEST_PATH_IMAGE011
分别为雷达的经度、纬度和高度对应的地理坐标；雷达监测的目标位置的经纬高坐标为
Figure 978517DEST_PATH_IMAGE012
，其中
Figure 313683DEST_PATH_IMAGE013
分别为目标位置的经度、纬度和高度对应的地理坐标；
S22：以目标位置为直角坐标系为零点
Figure 980288DEST_PATH_IMAGE014
，重建三维直角坐标系，则雷达的直角坐标
Figure 59102DEST_PATH_IMAGE015
为：
Figure 494763DEST_PATH_IMAGE016
其中
Figure 266410DEST_PATH_IMAGE017
为地球周长，则视线（LOS）方向即为指向
Figure 154731DEST_PATH_IMAGE018
的方向。
作为上述方案进一步的改进，视线（LOS）方向和最大坡度方向的夹角
Figure 37237DEST_PATH_IMAGE019
为：
Figure 825939DEST_PATH_IMAGE020
。
作为上述方案进一步的改进，沿最大坡度方向的实际位移量为
Figure 768487DEST_PATH_IMAGE021
，则
Figure 144105DEST_PATH_IMAGE022
其中
Figure 564722DEST_PATH_IMAGE023
为单视图像的地基InSAR测量得到的视线（LOS）方向位移量；由此，对于滑坡场景中任何一点地基InSAR获得的视线（LOS）方向位移，都可得到其沿最大坡度方向的实际位移。
本发明还提供一种地基InSAR提取实际滑坡位移向量的系统，包括：
坡向和坡度获取模块，用于获取目标位置的坡向和坡度；
最大坡度方向生成模块，与所述坡向和坡度获取模块连接，用于根据坡向和坡度得到目标位置的单位向量的三维坐标
Figure 974975DEST_PATH_IMAGE024
；
视线（LOS）方向获取模块，用于获取雷达到目标位置的视线（LOS）方向；
视线（LOS）方向和最大坡度方向的夹角
Figure 88424DEST_PATH_IMAGE025
计算模块，分别与所述最大坡度方向生成模块和视线（LOS）方向获取模块连接，用于计算得到视线（LOS）方向和最大坡度方向的夹角
Figure 951338DEST_PATH_IMAGE025
；
实际位移量获取模块，与所述视线（LOS）方向和最大坡度方向的夹角
Figure 175646DEST_PATH_IMAGE026
计算模块连接，用于计算实际滑坡位移向量。
本发明还提供一种存储介质，所述存储介质上存储有地基InSAR提取滑坡位移向量的方法，所述地基InSAR提取滑坡位移向量的程序被处理器执行时实现上述地基InSAR提取实际滑坡位移向量的方法的步骤。
需要说明的是，考虑到目标位置的单位向量的三维坐标
Figure 673361DEST_PATH_IMAGE027
仅仅与三维地形及雷达位置有关，当监测开始后不再变化，因此可以根据数字高程信息，在监测前对滑坡场景中每一个点从经纬度
Figure 957712DEST_PATH_IMAGE028
提取高度
Figure 307922DEST_PATH_IMAGE029
，并根据雷达位置
Figure 335921DEST_PATH_IMAGE030
生成对应的投影系数
Figure 517503DEST_PATH_IMAGE031
，在监测时直接调用即可。
由于本发明采用了以上技术方案，使本申请具备的有益效果在于 ：
相较于现有技术，本发明提供的地基InSAR提取滑坡位移向量的方法，首先根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 644859DEST_PATH_IMAGE032
和坡度
Figure 482365DEST_PATH_IMAGE033
；再基于单视图像的地基InSAR，测量得到雷达到目标位置的视线（LOS）方向位移
Figure 48476DEST_PATH_IMAGE003
；最后将步骤S2获得的视线（LOS）方向位移投影到步骤S1中坡向
Figure 350144DEST_PATH_IMAGE032
和坡度
Figure 146937DEST_PATH_IMAGE033
所确定的最大坡度方向上，从而获得滑坡位移向量。与传统的多轨图像测量雷达LOS方向相比，本发明提供的实际滑坡位移向量提取方法，仅需在现场部署一个雷达，大大减低了现场的部署难度以及整体设备成本，尤其适合地质灾害进行应急监测；与传统的偏移量跟踪（OT）方法相比，地基SAR孔径短对测试结果没有影响，从而避免了由于地基SAR孔径短造成的位移测量精度不高的问题，也即，提高了实际滑坡位移向量的测量精度。
具体实施方式
基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
另外，本发明各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
以下对本发明作进一步描述:
实施例1：
大部分滑坡是在重力作用下发生和发展的，特别滑坡的初始变形和等速变形阶段中重力往往是主要影响因素，而这两个阶段是地基InSAR边坡稳定监测的重要对象，因为地基InSAR的亚毫米精度和分钟级数据更新率非常适用于识别早期稳定性隐患。而此类滑坡的一个重要特点是，位移向量往往与边坡的最大坡度方向一致，即滑坡沿最大坡度方向进行，基于此，在重力性滑坡的早期阶段，本发明提供一种地基InSAR提取实际滑坡位移向量的方法，包括步骤：
S1、根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 534056DEST_PATH_IMAGE034
和坡度
Figure 841540DEST_PATH_IMAGE035
，具体的，提取滑坡场景的坡向和坡度时，三维地形的坡度图和坡向图利用GDAL（https://gdal.org）的gdaldem slope/aspect命令生成，或ArcGIS、QGIS相关工具软件获得，同时对获得的滑坡场景的坡度图和坡向图进行平滑滤波，如合适窗大小的均值、高斯等滤波器对坡度图和坡向图进行平滑滤波，以抑制小尺度坡度和坡向带来的噪声；
S2、基于单视图像的地基InSAR，测量得到雷达到目标位置的视线（LOS）方向位移
Figure 732136DEST_PATH_IMAGE003
；
S3、将步骤S2获得的雷达LOS方向位移投影到步骤S1中坡向
Figure 466874DEST_PATH_IMAGE034
和坡度
Figure 75709DEST_PATH_IMAGE035
所确定的最大坡度方向上，从而获得滑坡位移向量。
本发明提供的地基InSAR提取滑坡位移向量的方法，首先根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 249202DEST_PATH_IMAGE034
和坡度
Figure 197566DEST_PATH_IMAGE035
；再基于单视图像的地基InSAR，测量得到雷达到目标位置的视线（LOS）方向位移
Figure 165522DEST_PATH_IMAGE003
；最后将步骤S2获得的视线（LOS）方向位移投影到步骤S1中坡向
Figure 697872DEST_PATH_IMAGE034
和坡度
Figure 409477DEST_PATH_IMAGE035
所确定的最大坡度方向上，从而获得实际滑坡位移向量。与传统的多轨图像测量雷达LOS方向相比，本发明提供的实际滑坡位移向量提取方法，仅需在现场部署一个雷达，大大减低了现场的部署难度以及整体设备成本，尤其适合地质灾害进行应急监测；与传统的偏移量跟踪（OT）方法相比，地基SAR孔径短对测试结果没有影响，从而避免了由于地基SAR孔径短造成的位移测量精度不高的问题，也即，提高了实际滑坡位移向量的测量精度。
作为优选的实施例，在步骤S3中，
Figure 274664DEST_PATH_IMAGE004
坡向
Figure 351205DEST_PATH_IMAGE034
和坡度
Figure 934633DEST_PATH_IMAGE035
所确定的最大坡度方向的方法如下，雷达监测的目标位置的单位向量的三维坐标
Figure 122032DEST_PATH_IMAGE036
表示如下：
Figure 107305DEST_PATH_IMAGE037
则坡向
Figure 354747DEST_PATH_IMAGE034
和坡度
Figure 861689DEST_PATH_IMAGE035
所确定的最大坡度方向即为指向单位向量的三维坐标
Figure 852779DEST_PATH_IMAGE038
的方向。
作为优选的实施例，通过地理角度定义或几何角度定义来获取坡向
Figure 692559DEST_PATH_IMAGE034
；坡向以地理角度定义顺时针增加，以正北方向为0°角度；或者坡向以几何角度定义逆时针增加，以东方向为0°角度；且地理角度和几何角度存在以下关系：
Figure 110902DEST_PATH_IMAGE039
，其中
Figure 403343DEST_PATH_IMAGE040
为目标位置的坡向地理角度，
Figure 198124DEST_PATH_IMAGE041
为目标位置的坡向几何角度。
作为优选的实施例，在步骤S2中，所述视线（LOS）方向的获取步骤如下：
S21：假设雷达位置的经纬高坐标为
Figure 892410DEST_PATH_IMAGE042
，其中
Figure 980190DEST_PATH_IMAGE043
分别为雷达的经度、纬度和高度对应的地理坐标；雷达监测的目标位置的经纬高坐标为
Figure 25506DEST_PATH_IMAGE044
，其中
Figure 358399DEST_PATH_IMAGE045
分别为目标位置的经度、纬度和高度对应的地理坐标；
S22：以目标位置为直角坐标系为零点
Figure 907192DEST_PATH_IMAGE014
，重建三维直角坐标系，则雷达的直角坐标
Figure 995233DEST_PATH_IMAGE046
为：
Figure 199950DEST_PATH_IMAGE047
其中
Figure 336533DEST_PATH_IMAGE048
为地球周长，则视线（LOS）方向即为指向
Figure 739833DEST_PATH_IMAGE049
的方向。
作为优选的实施例，视线（LOS）方向和最大坡度方向的夹角
Figure 998776DEST_PATH_IMAGE050
为：
Figure 189323DEST_PATH_IMAGE051
。
作为优选的实施例，沿最大坡度方向的实际位移量为
Figure 926335DEST_PATH_IMAGE021
，则
Figure 387403DEST_PATH_IMAGE052
其中
Figure 817248DEST_PATH_IMAGE053
为单视图像的地基InSAR测量得到的视线（LOS）方向位移量；由此，对于滑坡场景中任何一点地基InSAR获得的视线（LOS）方向位移，都可以得到其沿最大坡度方向的实际位移。
需要说明的是，考虑到目标位置的单位向量的三维坐标
Figure 730977DEST_PATH_IMAGE054
仅仅与三维地形及雷达位置有关，当监测开始后不再变化，因此可以根据数字高程信息，在监测前对滑坡场景中每一个点从经纬度
Figure 271680DEST_PATH_IMAGE055
提取高度
Figure 852834DEST_PATH_IMAGE056
，并根据雷达位置
Figure 188000DEST_PATH_IMAGE057
生成对应的投影系数
Figure 353140DEST_PATH_IMAGE058
，在监测时直接调用即可。
实施例2：
本发明还提供一种地基InSAR提取滑坡位移向量的系统，包括：
坡向和坡度获取模块，用于获取目标位置的坡向和坡度；
最大坡度方向生成模块，与所述坡向和坡度获取模块连接，用于根据坡向和坡度得到目标位置的单位向量的三维坐标
Figure 431955DEST_PATH_IMAGE059
；
视线（LOS）方向获取模块，用于获取雷达到目标位置的视线（LOS）方向；
视线（LOS）方向和最大坡度方向的夹角
Figure 867615DEST_PATH_IMAGE060
计算模块，分别与所述最大坡度方向生成模块和视线（LOS）方向获取模块连接，用于计算得到视线（LOS）方向和最大坡度方向的夹角
Figure 311366DEST_PATH_IMAGE060
；
实际位移量获取模块，与所述视线（LOS）方向和最大坡度方向的夹角
Figure 465267DEST_PATH_IMAGE060
计算模块连接，用于计算实际滑坡位移向量。
实施例3：
本发明还提供一种存储介质，所述存储介质上存储有地基InSAR提取滑坡位移向量的程序，所述地基InSAR提取滑坡位移向量程序被处理器执行时实现上述地基InSAR提取实际滑坡位移向量的方法的步骤。
上述存储介质可能包含高速RAM存储器，也可能还包括非易失性存储器(non-volatilememory)，例如至少一个磁盘存储器。可以理解的是，存储介质可以为随机存储器(Random Access Memory，RAM)、磁碟、硬盘、固态硬盘(solid state disk，SSD)或者非易失性存储器等各种可以存储程序代码的机器可读介质。
本领域内的技术人员应明白，本发明实施例可提供为方法或存储介质。因此，本发明实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上是本发明的详细的介绍，本文中应用了具体个例对本发明的原理以及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法以及核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (10)

1.一种地基InSAR提取滑坡位移向量的方法，其特征在于，包括步骤：
S1、根据数字高程信息提取滑坡场景中目标位置的坡向
Figure 289922DEST_PATH_IMAGE001
和坡度
Figure 665540DEST_PATH_IMAGE002
；
S2、基于单视图像的地基InSAR，测量得到雷达到目标位置的视线方向位移
Figure 289419DEST_PATH_IMAGE003
；
S3、将步骤S2获得的视线方向位移
Figure 496409DEST_PATH_IMAGE003
投影到步骤S1中坡向
Figure 46077DEST_PATH_IMAGE001
和坡度
Figure 971308DEST_PATH_IMAGE002
所确定的最大坡度方向上，从而获得滑坡位移向量。
2.根据权利要求1所述的地基InSAR提取滑坡位移向量的方法，其特征在于，在步骤S3中，根据坡向
Figure 195616DEST_PATH_IMAGE001
和坡度
Figure 460375DEST_PATH_IMAGE002
所确定的最大坡度方向的方法如下，雷达监测的目标位置的单位向量的三维坐标
Figure 479147DEST_PATH_IMAGE004
表示如下：
Figure 891673DEST_PATH_IMAGE005
则坡向
Figure 60618DEST_PATH_IMAGE006
和坡度
Figure 242200DEST_PATH_IMAGE007
所确定的最大坡度方向即为指向单位向量的三维坐标
Figure 431873DEST_PATH_IMAGE008
的方向。
3.根据权利要求2所述的地基InSAR提取滑坡位移向量的方法，其特征在于，通过地理角度定义或几何角度定义来获取坡向
Figure 66117DEST_PATH_IMAGE009
；
坡向以地理角度定义顺时针增加，以正北方向为0°角度；或者坡向以几何角度定义逆时针增加，以东方向为0°角度；
且地理角度和几何角度存在以下关系：
Figure 68446DEST_PATH_IMAGE010
，其中
Figure 370114DEST_PATH_IMAGE011
为目标位置的坡向地理角度，
Figure 933951DEST_PATH_IMAGE012
为目标位置的坡向几何角度。
4.根据权利要求1-3任意一项所述的地基InSAR提取滑坡位移向量的方法，其特征在于，在步骤S2中，视线方向的获取步骤如下：
S21：假设雷达位置的经纬高坐标为
Figure 55490DEST_PATH_IMAGE013
，其中
Figure 425292DEST_PATH_IMAGE014
分别为雷达的经度、纬度和高度对应的地理坐标；雷达监测的目标位置的经纬高坐标为
Figure 519150DEST_PATH_IMAGE015
，其中
Figure 316204DEST_PATH_IMAGE016
分别为目标位置的经度、纬度和高度对应的地理坐标；
S22：以目标位置为直角坐标系为零点
Figure 862723DEST_PATH_IMAGE017
，重建三维直角坐标系，则雷达的直角坐标
Figure 770637DEST_PATH_IMAGE018
为：
Figure 951957DEST_PATH_IMAGE019
其中
Figure 326438DEST_PATH_IMAGE020
为地球周长，视线方向即为指向
Figure 953728DEST_PATH_IMAGE021
的方向。
5.根据权利要求4所述的地基InSAR提取滑坡位移向量的方法，其特征在于，视线方向和最大坡度方向的夹角
Figure 337436DEST_PATH_IMAGE022
为：
Figure 202624DEST_PATH_IMAGE023
。
6.根据权利要求5所述的地基InSAR提取滑坡位移向量的方法，其特征在于，沿最大坡度方向的实际位移量为
Figure 279164DEST_PATH_IMAGE024
，则
Figure 862592DEST_PATH_IMAGE025
其中
Figure 814106DEST_PATH_IMAGE026
为单视图像的地基InSAR测量得到的视线方向位移量；由此，对于滑坡场景中任何一点地基InSAR获得的视线方向位移，都可得到其沿最大坡度方向的实际位移。
7.根据权利要求1-3任意一项所述的地基InSAR提取滑坡位移向量的方法，其特征在于，所述步骤S1中，提取滑坡场景的坡向和坡度时，三维地形的坡度图和坡向图利用GDAL的gdaldem slope/aspect命令生成，或ArcGIS、QGIS相关工具软件获得。
8.根据权利要求7所述的地基InSAR提取滑坡位移向量的方法，其特征在于，对获得的滑坡场景的坡度图和坡向图进行平滑滤波，以抑制小尺度坡度和坡向带来的噪声。
9.一种地基InSAR提取滑坡位移向量的系统，其特征在于，包括：
坡向和坡度获取模块，用于获取目标位置的坡向和坡度；
最大坡度方向生成模块，与所述坡向和坡度获取模块连接，用于根据坡向和坡度得到目标位置的单位向量的三维坐标
Figure 533800DEST_PATH_IMAGE027
；
视线方向获取模块，用于获取雷达到目标位置的视线方向；
视线方向和最大坡度方向的夹角
Figure 781242DEST_PATH_IMAGE028
计算模块，分别与所述最大坡度方向生成模块和视线方向获取模块连接，用于计算得到视线方向和最大坡度方向的夹角
Figure 851966DEST_PATH_IMAGE028
；
实际位移量获取模块，与所述视线方向和最大坡度方向的夹角
Figure 843055DEST_PATH_IMAGE028
计算模块连接，用于计算实际滑坡位移向量。
10.一种存储介质，其特征在于，所述存储介质上存储有如权利要求1-8任意一项所述的地基InSAR提取滑坡位移向量的方法，所述地基InSAR提取滑坡位移向量程序被处理器执行时实现上述地基InSAR提取实际滑坡位移向量的方法的步骤。
CN202210135030.8A 2022-02-15 2022-02-15 地基InSAR提取滑坡位移向量的方法、系统及存储介质 Active CN114185048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210135030.8A CN114185048B (zh) 2022-02-15 2022-02-15 地基InSAR提取滑坡位移向量的方法、系统及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210135030.8A CN114185048B (zh) 2022-02-15 2022-02-15 地基InSAR提取滑坡位移向量的方法、系统及存储介质

Publications (2)

Publication Number Publication Date
CN114185048A true CN114185048A (zh) 2022-03-15
CN114185048B CN114185048B (zh) 2022-05-31

Family

ID=80545911

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210135030.8A Active CN114185048B (zh) 2022-02-15 2022-02-15 地基InSAR提取滑坡位移向量的方法、系统及存储介质

Country Status (1)

Country Link
CN (1) CN114185048B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966601A (zh) * 2022-08-01 2022-08-30 南京隼眼电子科技有限公司 基于毫米波雷达的山体滑坡预测方法和电子设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110048925A (ko) * 2009-11-03 2011-05-12 연세대학교 산학협력단 3차원 지표 변위 추출 방법 및 그 장치
CN102645650A (zh) * 2012-03-06 2012-08-22 北京北科安地科技发展有限公司 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术
CN102967851A (zh) * 2012-11-14 2013-03-13 电子科技大学 一种双基地sar的空间同步方法
CN108387899A (zh) * 2018-04-17 2018-08-10 南京师范大学 合成孔径雷达干涉测量中地面控制点自动选取方法
US20200225360A1 (en) * 2017-09-29 2020-07-16 Saab Ab A method for determining the base line for a synthetic aperture of a sar using gnss
CN113534154A (zh) * 2021-09-16 2021-10-22 成都理工大学 一种sar视线向变形与坡度坡向敏感度计算方法
CN113848551A (zh) * 2021-09-24 2021-12-28 成都理工大学 一种利用InSAR升降轨形变数据的滑坡深度反演方法
CN113866764A (zh) * 2021-08-13 2021-12-31 成都理工大学 基于InSAR和LR-IOE模型的滑坡易发性改进评估方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20110048925A (ko) * 2009-11-03 2011-05-12 연세대학교 산학협력단 3차원 지표 변위 추출 방법 및 그 장치
CN102645650A (zh) * 2012-03-06 2012-08-22 北京北科安地科技发展有限公司 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术
CN102967851A (zh) * 2012-11-14 2013-03-13 电子科技大学 一种双基地sar的空间同步方法
US20200225360A1 (en) * 2017-09-29 2020-07-16 Saab Ab A method for determining the base line for a synthetic aperture of a sar using gnss
CN108387899A (zh) * 2018-04-17 2018-08-10 南京师范大学 合成孔径雷达干涉测量中地面控制点自动选取方法
CN113866764A (zh) * 2021-08-13 2021-12-31 成都理工大学 基于InSAR和LR-IOE模型的滑坡易发性改进评估方法
CN113534154A (zh) * 2021-09-16 2021-10-22 成都理工大学 一种sar视线向变形与坡度坡向敏感度计算方法
CN113848551A (zh) * 2021-09-24 2021-12-28 成都理工大学 一种利用InSAR升降轨形变数据的滑坡深度反演方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
XIE HU ET AL.: "Internal kinematics of the Slumgullion landslide (USA) from high-resolution UAVSAR InSAR data", 《REMOTE SENSING OF ENVIRONMENT》 *
徐克科 等: "联合GPS、InSAR建立同震地表三维位移场", 《大地测量与地球动力学》 *
马娟: "一种大地坐标变换近似算法", 《安徽建筑大学学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966601A (zh) * 2022-08-01 2022-08-30 南京隼眼电子科技有限公司 基于毫米波雷达的山体滑坡预测方法和电子设备
CN114966601B (zh) * 2022-08-01 2022-10-21 南京隼眼电子科技有限公司 基于毫米波雷达的山体滑坡预测方法和电子设备

Also Published As

Publication number Publication date
CN114185048B (zh) 2022-05-31

Similar Documents

Publication Publication Date Title
CN113091599B (zh) 融合无人机dom和星载sar影像的地表三维形变提取方法
Lissak et al. Remote sensing for assessing landslides and associated hazards
Caduff et al. A review of terrestrial radar interferometry for measuring surface change in the geosciences
Rosenau et al. A processing system to monitor Greenland outlet glacier velocity variations at decadal and seasonal time scales utilizing the Landsat imagery
Kääb et al. Coseismic displacements of the 14 November 2016 M w 7.8 Kaikoura, New Zealand, earthquake using the Planet optical cubesat constellation
Zhang et al. Satellite SAR geocoding with refined RPC model
Crosetto et al. Urban subsidence monitoring using radar interferometry
Vassilakis et al. Post-event surface deformation of Amyntaio slide (Greece) by complementary analysis of Remotely Piloted Airborne System imagery and SAR interferometry
Jones et al. Pre-event and post-formation ground movement associated with the Bayou Corne sinkhole
CN114185048B (zh) 地基InSAR提取滑坡位移向量的方法、系统及存储介质
Ma et al. Decision-making fusion of InSAR technology and offset tracking to study the deformation of large gradients in mining areas-Xuemiaotan mine as an example
Ittycheria et al. Time series analysis of surface deformation of Bengaluru city using Sentinel-1 images
Dai et al. Characterization of the 2008 phreatomagmatic eruption of okmok from ArcticDEM and InSAR: deposition, erosion, and deformation
Sefercik et al. Country-scale discontinuity analysis of AW3D30 and SRTM Global DEMS: case study in Turkey
Atanasova et al. Application of interferometry in landslide area analysis
Walter et al. Fully automatic geo-coding of CARABAS-II VHF SAR images
Pawluszek-Filipiak et al. Comparison of PSI and DInSAR approach for the subsidence monitoring caused by coal mining exploitation
Hlaváčová et al. Sentinel-1 INSAR processing of corner reflector information in the northern-bohemian coal basin
Aziz et al. Vertical Accuracy Assessment of Improvised Global Digital Elevation Models (MERIT, NASADEM, EarthEnv) Using GNSS and Airborne IFSAR DEM
Bolkas On the estimation of uncertainties and fusion of multi-platform geodetic data
Yu et al. Urban monitoring using persistent scatterer InSAR and photogrammetry
Haugen Qualitative and quantitative comparative analyses of 3D lidar landslide displacement field measurements
Palamà et al. Analysis of mining-induced terrain deformation using multitemporal distributed scatterer SAR interferometry
Osmanoglu Applications and development of new algorithms for displacement analysis using InSAR time series
Weissgerber et al. LabSAR, a one-GCP coregistration tool for SAR–InSAR local analysis in high-mountain regions

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