CN115469308A - 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质 - Google Patents

多轨道InSAR震间形变速率场拼接方法、装置、设备及介质 Download PDF

Info

Publication number
CN115469308A
CN115469308A CN202211025070.3A CN202211025070A CN115469308A CN 115469308 A CN115469308 A CN 115469308A CN 202211025070 A CN202211025070 A CN 202211025070A CN 115469308 A CN115469308 A CN 115469308A
Authority
CN
China
Prior art keywords
deformation rate
target
track
deformation
insar
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
CN202211025070.3A
Other languages
English (en)
Other versions
CN115469308B (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.)
INSTITUTE OF GEOLOGY CHINA EARTHQUAKE ADMINISTRATION
Original Assignee
INSTITUTE OF GEOLOGY CHINA EARTHQUAKE ADMINISTRATION
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 INSTITUTE OF GEOLOGY CHINA EARTHQUAKE ADMINISTRATION filed Critical INSTITUTE OF GEOLOGY CHINA EARTHQUAKE ADMINISTRATION
Priority to CN202211025070.3A priority Critical patent/CN115469308B/zh
Publication of CN115469308A publication Critical patent/CN115469308A/zh
Application granted granted Critical
Publication of CN115469308B publication Critical patent/CN115469308B/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/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/9094—Theoretical aspects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
    • 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
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/35—Constructional details or hardware or software details of the signal processing chain
    • G01S19/37—Hardware or software details of the signal processing chain

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本公开实施例涉及一种多轨道InSAR震间形变速率场拼接方法、装置、设备及介质，方法包括：获取目标区域的GNSS的三维形变速率和SAR数据，将三维形变速率投影至目标向得到第一目标向(目标向为LOS向或地距向)形变速率，将SAR数据处理得到每个轨道的第二目标向形变速率，基于第一目标向形变速率对第二目标向形变速率进行参考基准统一得到每个轨道的第三目标向形变速率后进行入射角校正得到每个轨道的最终目标向形变速率，在地理坐标系下将每个轨道的最终目标向形变速率进行拼接。使用GNSS数据对InSAR形变速率场进行基准统一，在目标向进行大范围单一轨道形变速率场拼接时进行入射角校正，获取高精度、高空间分辨率、大范围目标向形变速率场。

Description

多轨道InSAR震间形变速率场拼接方法、装置、设备及介质
技术领域
本公开涉及数据处理技术领域，尤其涉及一种多轨道InSAR(InterferometricSynthetic Aperture Radar，合成孔径雷达干涉测量)震间形变速率场拼接方法、装置、设备及介质。
背景技术
InSAR是一种高精度、高效率的地表形变测量方法。随着合成孔径雷达干涉测量技术的飞速发展，海量高质量的干涉图使得大面积的地表形变监测成为可能。
然而，星载SAR的标准景幅宽有限，常常无法覆盖整个形变区域，需要将多轨InSAR数据进行拼接来开展大范围地表形变的监测，以青藏高原地区为例，其内部分布许多东西走向的大型走滑断裂，走向长度往往达到上百公里，需要多个同一轨道(升轨或者降轨)SAR数据条带才能完全覆盖。
相关技术中，通过邻轨PS(永久散射体，Persistent Scatterer)-InSAR地面沉降结果拼接处理的方式是假设地表只发生了垂直运动，但是实际情况中，地表的运动方向是三维的，导致拼接结果不够准确。
发明内容
为了解决上述技术问题或者至少部分地解决上述技术问题，本公开提供了一种多轨道InSAR震间形变速率场拼接方法、装置、设备及介质。
本公开实施例提供了一种多轨道InSAR震间形变速率场拼接方法，所述方法包括：
获取目标区域的全球导航卫星系统GNSS(Global Navigation SatelliteSystem，全球导航卫星系统)的三维形变速率和多轨道合成孔径雷达SAR数据；
将所述三维形变速率投影至目标向，得到第一目标向形变速率；
将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；
基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；
对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接。
本公开实施例还提供了一种多轨道InSAR震间形变速率场拼接装置，所述装置包括：
获取模块，用于获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径SARSAR数据；
投影模块，用于将所述三维形变速率投影至目标向，得到第一目标向形变速率；
处理模块，用于将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；
基准统一模块，用于基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；
校正模块，用于对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率；
拼接模块，用于在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接。
本公开实施例还提供了一种电子设备，所述电子设备包括：处理器；用于存储所述处理器可执行指令的存储器；所述处理器，用于从所述存储器中读取所述可执行指令，并执行所述指令以实现如本公开实施例提供的多轨道InSAR震间形变速率场拼接方法。
本公开实施例还提供了一种计算机可读存储介质，所述存储介质存储有计算机程序，所述计算机程序用于执行如本公开实施例提供的多轨道InSAR震间形变速率场拼接方法。
本公开实施例提供的技术方案与现有技术相比具有如下优点：本公开实施例提供的多轨道InSAR震间形变速率场拼接方案，获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径雷达SAR数据，将三维形变速率投影至目标向，得到第一目标向形变速率，将SAR数据进行时序处理，得到每个轨道的第二目标向形变速率，基于第一目标向形变速率对第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率，对第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将每个轨道的最终目标向形变速率进行拼接。采用上述技术方案，使用GNSS数据对InSAR形变速率场进行基准统一，在目标向进行大范围单一轨道形变速率场拼接时进行入射角校正，获取高精度、高空间分辨率、大范围目标向比如LOS向或地距向的形变速率场。
附图说明
结合附图并参考以下具体实施方式，本公开各实施例的上述和其他特征、优点及方面将变得更加明显。贯穿附图中，相同或相似的附图标记表示相同或相似的元素。应当理解附图是示意性的，原件和元素不一定按照比例绘制。
图1为本公开实施例提供的一种多轨道InSAR震间形变速率场拼接方法的流程示意图；
图2为本公开实施例提供的一种多轨道InSAR震间形变速率场拼接装置的结构示意图。
具体实施方式
下面将参照附图更详细地描述本公开的实施例。虽然附图中显示了本公开的某些实施例，然而应当理解的是，本公开可以通过各种形式来实现，而且不应该被解释为限于这里阐述的实施例，相反提供这些实施例是为了更加透彻和完整地理解本公开。应当理解的是，本公开的附图及实施例仅用于示例性作用，并非用于限制本公开的保护范围。
应当理解，本公开的方法实施方式中记载的各个步骤可以按照不同的顺序执行，和/或并行执行。此外，方法实施方式可以包括附加的步骤和/或省略执行示出的步骤。本公开的范围在此方面不受限制。
本文使用的术语“包括”及其变形是开放性包括，即“包括但不限于”。术语“基于”是“至少部分地基于”。术语“一个实施例”表示“至少一个实施例”；术语“另一实施例”表示“至少一个另外的实施例”；术语“一些实施例”表示“至少一些实施例”。其他术语的相关定义将在下文描述中给出。
需要注意，本公开中提及的“第一”、“第二”等概念仅用于对不同的装置、模块或单元进行区分，并非用于限定这些装置、模块或单元所执行的功能的顺序或者相互依存关系。
需要注意，本公开中提及的“一个”、“多个”的修饰是示意性而非限制性的，本领域技术人员应当理解，除非在上下文另有明确指出，否则应该理解为“一个或多个”。
本公开实施方式中的多个装置之间所交互的消息或者信息的名称仅用于说明性的目的，而并不是用于对这些消息或信息的范围进行限制。
为了实现多轨道InSAR震间形变速率场拼接，将垂直形变监测领域扩展到更为复杂的构造形变监测领域中，既要统一多轨道InSAR形变速率场的参考基准，又要克服雷达卫星入射角的影响，本公开提出一种多轨道InSAR震间形变速率场拼接方法，通过使用GNSS数据对InSAR形变速率场进行基准统一，在目标向进行大范围单一轨道形变速率场拼接时进行入射角校正，获取高精度、高空间分辨率、大范围目标向比如LOS向或地距向的形变速率场。
具体地，图1为本公开实施例提供的一种多轨道InSAR震间形变速率场拼接方法的流程示意图，该方法可以由多轨道InSAR震间形变速率场拼接装置执行，其中该装置可以采用软件和/或硬件实现，一般可集成在电子设备中。如图1所示，该方法包括：
步骤101、获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径雷达SAR数据。
其中，目标区域可以根据应用场景选择设置，比如青藏高原地区等。
本公开实施例中，获取GNSS的三维形变速率的方式有很多种，在一些实施方式中，收集目标区域对应的SAR影像覆盖的GNSS三维形变速率或者收集目标区域原始GNSS数据，经过差分处理，计算出对于站点(GNSS站)的三维形变速率，最终得到GNSS站的三维形变速率。
具体地，使用GAMIT/GLOBK软件处理原始GNSS数据，采用双差电离层组合观测值，比如采用时间间隔30s，时间段24h，高度截止角定为10，对流层每1h估计天顶对流层延迟及每天估计1个NS(北南)和EW(东西)梯度，固体潮汐模型采用IERS2003模型，极潮模型采用FES2004模型，通过上述模型与方法，可以进行单日松弛解处理，得到一系列待估参数及其方差-协方差阵。利用GLOBK进行多时段综合解算，获得网平差结果。在进行GLOBK解算时，把基准站单日松弛解与IGS(International GNSS Service，国际GNSS服务)站单日松弛解联合平差，选取全球均匀分布的IGS站为基准点，最后得到GNSS的三维形变速率。
在本公开实施例中，获取SAR数据的方式有很多种，在一些实施方式中，收集目标区域的升轨或者降轨SAR数据。
步骤102、将三维形变速率投影至目标向，得到第一目标向形变速率。
其中，目标向可以是LOS向、地距向等，具体根据应用场景需要选择计算。第一目标向形变速率指的是将GNSS的三维形变速率投影至目标向后的形变速率，比如第一LOS向形变速率、第一地距向形变速率。
在本公开实施例中，将三维形变速率投影至目标向，得到第一目标向形变速率的方式有很多种，在一些实施方式中，当SAR卫星包括水平运动和垂直运动，获取三维形变速率对应的垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角，基于第一计算公式对垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角进行计算，得到LOS视线向形变速率作为第一目标向形变速率。
在另一些实施方式中，SAR卫星包括水平运动，获取三维形变速率对应的南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角，基于第二计算公式对南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角进行计算，得到第一地距向形变速率第一目标向形变速率。以上两种方式仅为将形变速率投影至目标向，得到第一目标向形变速率的示例，本公开实施例不对将形变速率投影至目标向，得到第一目标向形变速率的具体方式进行限定。
步骤103、将SAR数据进行时序处理，得到每个轨道的第二目标向形变速率。
其中，对SAR数据进行时序处理，获取对应的LOS向形变速率或者是将LOS向形变速率投影到地距向，得到地距向形变速率。比如可以根据D-InSAR技术处理获取对应LOS向三维形变速率，然后通过相位叠加(即将同一轨道相位叠加除以总的时间差)获取每个轨道的LOS向形变速率。第二目标向形变速率指的是SAR数据对应的LOS向形变速率、地距向形变速率。
其中，根据D-InSAR技术处理获取对应LOS向形变速率包括：基于数字高程模型(Digital Elevation Model，DEM)、扫描行校正器(Scan Line Corrector，SLC)主影像模拟SAR影像后进行影像配准，将DEM转换到雷达坐标系后获取地形相位，同时SLC主影像和SLC辅影像进行影像配准后辅影像进行重采样，得到复共轭相乘生成干涉图后进行基线估算，得到去平地相位，基于地形相位和去平地相位得到差分干涉图后进行干涉图滤波、相位解缠和地理编码得到LOS向形变速率。
步骤104，基于第一目标向形变速率对第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率。
其中，第三目标向形变速率指的是对SAR数据对应的LOS向形变速率、地距向形变速率按照GNSS的三维形变速率投影至目标向后的形变速率进行参考基准统一后的目标向形变速率。
在一些实施例中，使用GNSS数据投影的目标向形变速率，对SAR数据的目标向形变速率比如LOS向三维形变速率和地距向三维形变速率分别进行参考基准统一(即GNSS数据投影的LOS向三维形变速率和时序技术得到的LOS向三维形变速率存在一个差值)。
在一个具体实施例中，对第一目标向形变速率和LOS向形变速率进行差值处理，得到第一目标差值，使用二阶二次多项式以第一目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并基于最小二乘求解第一待估参数，基于第一待估参数确定第一形变速率比值，基于第一形变速率比值和LOS向形变速率进行计算，得到每个轨道的第三目标向形变速率。
步骤105，对第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将每个轨道的最终目标向形变速率进行拼接。
其中，最终目标向形变速率指的是对第三目标向形变速率进行入射角校正后的目标向形变速率。
在本公开实施例中，对参考基准统一后的目标向比如LOS向形变速率进行入射角校正(入射角指的是SAR卫星的入射角，即基于多项式改正的方法实现入射角校正)。
在一些实施方式中，基于SAR卫星的方位角，将各轨道的第三目标向形变速率转换到SAR坐标系下，并获取相邻轨道第三目标向形变速率之间的第二目标差值，使用二阶二次多项式以第二目标差值作为因变量、SAR卫星视线向作为自变量进行拟合，并采用最小二乘求解第二待估参数，基于第二待估参数确定第二形变速率比值，基于第二形变速率比值和第三目标向形变速率进行计算，得到每个轨道的最终目标向形变速率。
在另一些实施方式中，将第二目标向形变速率投影到地距向，得到第二地距向形变速率，计算第一地距向形变速率和第二地距向形变速率之间的第三目标差值，使用二阶二次多项式以第三目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并采用最小二乘求解第三待估参数，基于第三待估参数确定第三形变速率比值，基于第三形变速率比值和第二地距向形变速率进行计算，得到每个轨道的最终目标向形变速率。
在一些实施方式中，在目标向比如分别LOS向和地距向进行大范围形变速率场拼接(即在地理坐标系下将各个轨道的最终目标向形变速率直接进行拼接，其中，地面上的点坐标必然在一个参考坐标系下，地理坐标系是指以经纬度来表示地面点位置)。
本公开实施例提供的多轨道InSAR震间形变速率场拼接方案，获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径雷达SAR数据，将三维形变速率投影至目标向，得到第一目标向形变速率，将SAR数据进行时序处理，得到每个轨道的第二目标向形变速率，基于第一目标向形变速率对第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率，对第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将每个轨道的最终目标向形变速率进行拼接。采用上述技术方案，使用GNSS数据对InSAR形变速率场进行基准统一，在目标向进行大范围单一轨道形变速率场拼接时进行入射角校正，获取高精度、高空间分辨率、大范围目标向比如LOS向或地距向的形变速率场。
在一些实施例中，SAR卫星包括水平运动和垂直运动，将三维形变速率投影至目标向，得到第一目标向形变速率，包括：获取三维形变速率，即获取对应的垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角；基于第一计算公式对垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角进行计算，得到LOS视线向形变速率作为第一目标向形变速率。
具体地，将GNSS三维形变速率(东西南北垂直形变速率)投影到LOS向，比如第一计算公式为：
Figure BDA0003815332210000091
其中，dlos表示InSAR的LOS向形变速率，du、dn和de分别表示地表垂直向、南北向和东西向的三维形变速率，θinc和αh分别表示雷达卫星相对采集位置的入射角和方位角。
在一些实施例中，基于第二计算公式对南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角进行计算，得到第一地距向形变速率作为第一目标向形变速率。
具体地，将GNSS三维形变速率(东西南北垂直形变速率)投影到地距向，比如第二计算公式为：
dald＝-[sin(θinc)DN(range_line)+sin(θinc)DE(range_line)]
其中，dald表示地距向的三维形变速率；dn和de分别表示南北向和东西向的形变速率，θinc和αh分别表示雷达卫星相对采集位置的入射角和方位角；
Figure BDA0003815332210000101
Figure BDA0003815332210000102
在一些实施例中，对第一目标向形变速率和LOS向形变速率进行差值处理，得到第一目标差值，使用二阶二次多项式以所述第一目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并基于最小二乘求解第一待估参数，基于第一待估参数确定第一形变速率比值，基于第一形变速率比值和LOS向形变速率进行计算，得到每个轨道的第三目标向形变速率。
具体地，使用GAMMA软件对SAR数据进行了主辅影像配准、地形去除、平地相位去除、干涉图生成和滤波、相位解缠、stacking集成、地理编码等，获取了地表LOS三维形变速率。后续步骤都是在地理坐标系下进行的。
具体地，基于SAR卫星的入射角和方位角，将GNSS的三维形变速率投影到LOS向，并在各个轨道InSAR的形变速率中，基于地理坐标搜索与GNSS站对应的InSAR数据LOS向形变速率场，然后GNSS投影的对LOS向形变速率和搜索到的LOS向形变速率进行差值处理，并以差值作为因变量，地理坐标(即经纬度)作为自变量，使用二阶二次多项式进行拟合，并基于最小二乘求解对应的待估参数，最后使用该二次多项式，并加上InSAR形变速率场，实现对InSAR形变速率进行基准的绝对统一。
其中，使用该二次多项式，并加上InSAR形变速率场的公式为：
Figure BDA0003815332210000103
其中，
Figure BDA0003815332210000111
表示入射角修正后的LOS向形变速率；dlos(x,y)表示入射角修正前的LOS向形变速率。
其中，
Figure BDA0003815332210000112
其中，
Figure BDA0003815332210000113
表示相邻轨道的重叠区域点位上形变速率差相对于主轨道的形变信号差异的比值；(x,y)表示重叠区域像素点的东西和南北向距离，a，b，c是第一待估参数。
在一些实施例中，基于SAR卫星的方位角，将各轨道的第三目标向形变速率转换到SAR坐标系下，并获取相邻轨道第三目标向形变速率之间的第二目标差值，使用二阶二次多项式以第二目标差值作为因变量、SAR卫星视线向作为自变量进行拟合，并采用最小二乘求解第二待估参数，基于第二待估参数确定第二形变速率比值，基于第二形变速率比值和第三目标向形变速率进行计算，得到每个轨道的最终目标向形变速率。
具体地，在基于GNSS数据对InSAR数据LOS向形变速率场进行基准统一后，相邻轨道重叠区域仍然存在入射角的差异，且不可忽略。首先基于SAR卫星的方位角，将各轨道InSAR形变速率场转换到雷达坐标系下，并获取相邻轨道InSAR形变速率场的差值，以差值作为因变量，距离向(即SAR卫星视线向)作为自变量，使用二阶多项式进行拟合，并采用最小二乘求待估参数。以该二阶多项式消除或减弱入射角的影响。最后，对入射角修正后的InSAR三维形变速率场重新进行地理编码实现在地理坐标系下的InSAR大范围形变速率场的拼接。
在一些实施例中，对第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，包括：将第二目标向形变速率投影到地距向，得到第二地距向形变速率；计算第一地距向形变速率和第二地距向形变速率之间的第三目标差值；使用二阶二次多项式以第三目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并采用最小二乘求解第三待估参数；基于第三待估参数确定第三形变速率比值；基于第三形变速率比值和第二地距向三维形变速率进行计算，得到每个轨道的最终目标向形变速率。
具体地，若目标区域以水平运动为主，即可以忽略垂直形变的影响，可以实现在地距向大范围InSAR形变速率的拼接，从而避免入射角的影响。即基于SAR卫星的入射角和方位角，将GNSS形变速率场投影到地距向，InSAR数据LOS向形变速率场投影到地距向，然后搜索GNSS站对应的InSAR数据LOS向三维形变速率场。
进一步对地距向的GNSS数据和搜索到的地距向三维形变速率数据进行差值处理，并以差值作为因变量，地理坐标作为自变量，同样可以使用二阶二次多项式进行拟合，并基于最小二乘求解对应的待估参数，最后使用该二次多项式，并加上InSAR形变速率，实现对InSAR形变速率进行基准的绝对统一，最后在地理坐标系下，实现大范围在地距向大范围InSAR形变速率的拼接。
由此，如果仅仅以水平运动为主，可以考虑只针对地距向投影计算，如果水平和垂直都要考虑的话是通过LOS方向投影计算。
由此，通过使用GNSS数据对InSAR形变速率进行了基准统一，在LOS向进行大范围单一轨道形变速率拼接时，进行入射角校正，成功获取高精度、高空间分辨率、大范围LOS向形变速率数据。考虑到在构造形变领域，地表水平运动一般较垂直运动更为强烈，通过以GNSS数据为参考数据，对InSAR形变速率在地距向进行基准统一，实现了在地距向进行大范围单一轨道形变速率拼接。
综上，本公开提出的技术方案考虑了影响大范围单一轨道InSAR形变速率场的主要因素，即参考点和入射角，在保证各个轨道SAR数据时间覆盖基本一致的条件下，解决参考基准和入射角的差异后,在地距向和LOS向拼接的图像结果很好，成功实现了高精度、大范围InSAR三维形变速率重建，对研究大范围、高精度构造形变的空间变化具有重要意义，方式简单且满足应用需求。
图2为本公开实施例提供的一种多轨道InSAR震间形变速率场拼接装置的结构示意图，该装置可由软件和/或硬件实现，一般可集成在电子设备中。如图2所示，该装置包括：
获取模块201，用于获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径SARSAR数据；
投影模块202，用于将所述三维形变速率投影至目标向，得到第一目标向形变速率；
处理模块203，用于将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；
基准统一模块204，用于基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；
校正模块205，用于对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率；
拼接模块206，用于在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接。
可选的，SAR卫星包括水平运动和垂直运动，所述投影模块202具体用于：
获取所述三维形变速率对应的垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角；
基于第一计算公式对所述垂直向形变速率、所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到LOS视线向形变速率作为所述第一目标向形变速率。
可选的，所述SAR卫星仅包括水平运动，所述投影模块202具体还用于：
基于第二计算公式对所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到第一地距向形变速率作为所述第一目标向形变速率。
可选的，所述基准统一模块204具体用于：
对所述第一目标向形变速率和所述LOS向形变速率进行差值处理，得到第一目标差值；
使用二阶二次多项式以所述第一目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并基于最小二乘求解第一待估参数；
基于所述第一待估参数确定第一形变速率比值；
基于所述第一形变速率比值和所述LOS向形变速率进行计算，得到所述每个轨道的第三目标向形变速率。
可选的，所述校正模块205具体用于：
基于SAR卫星的方位角，将各轨道的所述第三目标向形变速率转换到SAR坐标系下，并获取相邻轨道所述第三目标向形变速率之间的第二目标差值；
使用二阶二次多项式以所述第二目标差值作为因变量、SAR卫星视线向作为自变量进行拟合，并采用最小二乘求解第二待估参数；
基于所述第二待估参数确定第二形变速率比值；
基于所述第二形变速率比值和所述第三目标向形变速率进行计算，得到所述每个轨道的最终目标向形变速率。
可选的，所述校正模块205具体用于：
将所述第二目标向形变速率投影到地距向，得到第二地距向形变速率；
计算所述第一地距向形变速率和所述第二地距向形变速率之间的第三目标差值；
使用二阶二次多项式以所述第三目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并采用最小二乘求解第三待估参数；
基于所述第三待估参数确定第三形变速率比值；
基于所述第三形变速率比值和所述第二地距向形变速率进行计算，得到所述每个轨道的最终目标向形变速率。
可选的，所述拼接模块206具体用于：
基于所述任意相邻两个轨道在所述地理坐标系下的最终目标向形变速率进行叠加，得到拼接结果。
本公开实施例所提供的多轨道InSAR震间形变速率场拼接装置可执行本公开任意实施例所提供的多轨道InSAR震间形变速率场拼接方法，具备执行方法相应的功能模块和有益效果。
本公开实施例还提供了一种计算机程序产品，包括计算机程序/指令，该计算机程序/指令被处理器执行时实现本公开任意实施例所提供的多轨道InSAR震间形变速率场拼接方法。
根据本公开的一个或多个实施例，本公开提供了一种电子设备，包括：
处理器；
用于存储所述处理器可执行指令的存储器；
所述处理器，用于从所述存储器中读取所述可执行指令，并执行所述指令以实现如本公开提供的任一所述的多轨道InSAR震间形变速率场拼接方法。
根据本公开的一个或多个实施例，本公开提供了一种计算机可读存储介质，所述存储介质存储有计算机程序，所述计算机程序用于执行如本公开提供的任一所述的多轨道InSAR震间形变速率场拼接方法。
以上描述仅为本公开的较佳实施例以及对所运用技术原理的说明。本领域技术人员应当理解，本公开中所涉及的公开范围，并不限于上述技术特征的特定组合而成的技术方案，同时也应涵盖在不脱离上述公开构思的情况下，由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本公开中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案。
此外，虽然采用特定次序描绘了各操作，但是这不应当理解为要求这些操作以所示出的特定次序或以顺序次序执行来执行。在一定环境下，多任务和并行处理可能是有利的。同样地，虽然在上面论述中包含了若干具体实现细节，但是这些不应当被解释为对本公开的范围的限制。在单独的实施例的上下文中描述的某些特征还可以组合地实现在单个实施例中。相反地，在单个实施例的上下文中描述的各种特征也可以单独地或以任何合适的子组合的方式实现在多个实施例中。
尽管已经采用特定于结构特征和/或方法逻辑动作的语言描述了本主题，但是应当理解所附权利要求书中所限定的主题未必局限于上面描述的特定特征或动作。相反，上面所描述的特定特征和动作仅仅是实现权利要求书的示例形式。

Claims (10)

1.一种多轨道InSAR震间形变速率场拼接方法，其特征在于，包括：
获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径雷达SAR数据；
将所述三维形变速率投影至目标向，得到第一目标向形变速率；
将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；
基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；
对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接。
2.根据权利要求1所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，基于时序InSAR技术获取的SAR卫星的形变速率包括水平运动和垂直运动，所述将所述GNSS的三维形变速率投影至目标向，得到第一目标向形变速率，包括：
获取所述三维形变速率对应的垂直向形变速率、南北向形变速率、东西向形变速率和SAR卫星的入射角和方位角；
基于第一计算公式对所述垂直向形变速率、所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到LOS视线向三维形变速率作为所述第一目标向形变速率。
3.根据权利要求2所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述SAR卫星仅包括水平运动，还包括：
基于第二计算公式对所述南北向形变速率、所述东西向形变速率和所述SAR卫星的入射角和方位角进行计算，得到第一地距向形变速率作为所述第一目标向形变速率。
4.根据权利要求2所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率，包括：
对所述第一目标向形变速率和所述LOS向形变速率进行差值处理，得到第一目标差值；
使用二阶二次多项式以所述第一目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并基于最小二乘求解第一待估参数；
基于所述第一待估参数确定第一形变速率比值；
基于所述第一形变速率比值和所述LOS向形变速率进行计算，得到所述每个轨道的第三目标向形变速率。
5.根据权利要求2所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，包括：
基于SAR卫星的方位角，将各轨道的所述第三目标向形变速率转换到SAR坐标系下，并获取相邻轨道所述第三目标向形变速率之间的第二目标差值；
使用二阶二次多项式以所述第二目标差值作为因变量、SAR卫星视线向作为自变量进行拟合，并采用最小二乘求解第二待估参数；
基于所述第二待估参数确定第二形变速率比值；
基于所述第二形变速率比值和所述第三目标向形变速率进行计算，得到所述每个轨道的最终目标向形变速率。
6.根据权利要求3所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率，包括：
将所述第二目标向形变速率投影到地距向，得到第二地距向形变速率；
计算所述第一地距向形变速率和所述第二地距向形变速率之间的第三目标差值；
使用二阶二次多项式以所述第三目标差值作为因变量、地理经纬度坐标作为自变量进行拟合，并采用最小二乘求解第三待估参数；
基于所述第三待估参数确定第三形变速率比值；
基于所述第三形变速率比值和所述第二地距向形变速率进行计算，得到所述每个轨道的最终目标向形变速率。
7.根据权利要求1所述的多轨道InSAR震间形变速率场拼接方法，其特征在于，所述在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接，包括：
基于任意相邻两个轨道在所述地理坐标系下的最终目标向形变速率进行叠加，得到拼接结果。
8.一种多轨道InSAR震间形变速率场拼接装置，其特征在于，包括：
获取模块，用于获取目标区域的全球导航卫星系统GNSS的三维形变速率和多轨道合成孔径SARSAR数据；
投影模块，用于将所述形变速率投影至目标向，得到第一目标向形变速率；
处理模块，用于将所述SAR数据进行时序处理，得到每个轨道的第二目标向形变速率；
基准统一模块，用于基于所述第一目标向形变速率对所述第二目标向形变速率进行参考基准统一，得到每个轨道的第三目标向形变速率；
校正模块，用于对所述第三目标向形变速率进行入射角校正，得到每个轨道的最终目标向形变速率；
拼接模块，用于在地理坐标系下将所述每个轨道的最终目标向形变速率进行拼接。
9.一种电子设备，其特征在于，所述电子设备包括：
处理器；
用于存储所述处理器可执行指令的存储器；
所述处理器，用于从所述存储器中读取所述可执行指令，并执行所述指令以实现上述权利要求1-7中任一所述的多轨道InSAR震间形变速率场拼接方法。
10.一种计算机可读存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序用于执行上述权利要求1-7中任一所述的多轨道InSAR震间形变速率场拼接方法。
CN202211025070.3A 2022-08-25 2022-08-25 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质 Active CN115469308B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211025070.3A CN115469308B (zh) 2022-08-25 2022-08-25 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211025070.3A CN115469308B (zh) 2022-08-25 2022-08-25 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质

Publications (2)

Publication Number Publication Date
CN115469308A true CN115469308A (zh) 2022-12-13
CN115469308B CN115469308B (zh) 2023-07-04

Family

ID=84368712

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211025070.3A Active CN115469308B (zh) 2022-08-25 2022-08-25 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质

Country Status (1)

Country Link
CN (1) CN115469308B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116465302A (zh) * 2023-03-31 2023-07-21 中国地震局地质研究所 一种断层运动的监测方法、装置、设备及存储介质
CN117109426A (zh) * 2023-08-28 2023-11-24 兰州交通大学 一种融合GNSS/InSAR观测资料的三维形变场建模方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594224A (zh) * 2018-03-30 2018-09-28 中国电力工程顾问集团中南电力设计院有限公司 融合不同平台和轨道sar数据的三维时序形变监测方法
CN109029344A (zh) * 2018-07-10 2018-12-18 湖南中科星图信息技术有限公司 一种基于高分影像和升降轨InSAR的堤坝沉降监测方法
CN112284332A (zh) * 2020-08-31 2021-01-29 北京四象爱数科技有限公司 基于高分辨率insar的高层建筑沉降监测结果三维定位方法
CN112540369A (zh) * 2020-11-27 2021-03-23 武汉大学 融合GNSS与升降轨时序InSAR的滑坡三维形变解算方法及系统
CN113405447A (zh) * 2020-05-19 2021-09-17 湖南北斗微芯产业发展有限公司 一种综合InSAR和GNSS的轨道交通形变监测方法、装置及设备
CN114660602A (zh) * 2022-03-24 2022-06-24 中南大学 一种广域InSAR形变速率自适应拼接融合方法、装置、设备及介质

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108594224A (zh) * 2018-03-30 2018-09-28 中国电力工程顾问集团中南电力设计院有限公司 融合不同平台和轨道sar数据的三维时序形变监测方法
CN109029344A (zh) * 2018-07-10 2018-12-18 湖南中科星图信息技术有限公司 一种基于高分影像和升降轨InSAR的堤坝沉降监测方法
CN113405447A (zh) * 2020-05-19 2021-09-17 湖南北斗微芯产业发展有限公司 一种综合InSAR和GNSS的轨道交通形变监测方法、装置及设备
CN112284332A (zh) * 2020-08-31 2021-01-29 北京四象爱数科技有限公司 基于高分辨率insar的高层建筑沉降监测结果三维定位方法
CN112540369A (zh) * 2020-11-27 2021-03-23 武汉大学 融合GNSS与升降轨时序InSAR的滑坡三维形变解算方法及系统
CN114660602A (zh) * 2022-03-24 2022-06-24 中南大学 一种广域InSAR形变速率自适应拼接融合方法、装置、设备及介质

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116465302A (zh) * 2023-03-31 2023-07-21 中国地震局地质研究所 一种断层运动的监测方法、装置、设备及存储介质
CN116465302B (zh) * 2023-03-31 2023-11-10 中国地震局地质研究所 一种断层运动的监测方法、装置、设备及存储介质
CN117109426A (zh) * 2023-08-28 2023-11-24 兰州交通大学 一种融合GNSS/InSAR观测资料的三维形变场建模方法
CN117109426B (zh) * 2023-08-28 2024-03-22 兰州交通大学 一种融合GNSS/InSAR观测资料的三维形变场建模方法

Also Published As

Publication number Publication date
CN115469308B (zh) 2023-07-04

Similar Documents

Publication Publication Date Title
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
Xu et al. High-rate precise point positioning (PPP) to measure seismic wave motions: an experimental comparison of GPS PPP with inertial measurement units
CN115469308B (zh) 多轨道InSAR震间形变速率场拼接方法、装置、设备及介质
Catalão et al. Merging GPS and atmospherically corrected InSAR data to map 3-D terrain displacement velocity
CN104597471B (zh) 面向时钟同步多天线gnss接收机的定向测姿方法
CN102968631B (zh) 山区多光谱遥感卫星影像的自动几何纠正与正射校正方法
CN110058236A (zh) 一种面向三维地表形变估计的InSAR和GNSS定权方法
CN101241011B (zh) 激光雷达平台上高精度定位、定姿的装置和方法
Fielding et al. Surface deformation related to the 2019 Mw 7.1 and 6.4 Ridgecrest earthquakes in California from GPS, SAR interferometry, and SAR pixel offsets
CN105627991A (zh) 一种无人机影像实时全景拼接方法及系统
US20160259044A1 (en) Three-dimensional positioning method
CN108983232B (zh) 一种基于邻轨数据的InSAR二维地表形变监测方法
CN111650579B (zh) 一种岩移参数自适应获取的InSAR矿区三维形变估计方法、装置及介质
CN111398980B (zh) 一种机载LiDAR数据处理的方法及装置
CN101876701A (zh) 一种侧视雷达遥感影像定位方法
CN109471104B (zh) 一种从两平行轨道sar数据中获取地表三维移动量的方法
Cao et al. Bundle adjustment of satellite images based on an equivalent geometric sensor model with digital elevation model
Tang et al. Geometric accuracy analysis model of the ZiYuan-3 satellite without GCPs
CN105444778A (zh) 一种基于成像几何反演的星敏感器在轨定姿误差获取方法
Ma Very long baseline interferometry applied to polar motion, relativity and geodesy
KR20120009186A (ko) Ｓａｒ데이터를 이용하여 수치고도모형을 제작하기 위한 방법
Cramer et al. On the use of GPS/inertial exterior orientation parameters in airborne photogrammetry
Ye et al. Initial orbit determination of BDS-3 satellites based on new code signals
CN113761808A (zh) 一种基于gps和经验潮汐模型的地表潮汐位移获取方法及应用方法和系统
CN115201825A (zh) 一种InSAR震间形变监测中的大气延迟校正方法

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