CN112904337A - 一种基于Offset Tracking技术的边坡形变时序监测方法 - Google Patents
一种基于Offset Tracking技术的边坡形变时序监测方法 Download PDFInfo
- Publication number
- CN112904337A CN112904337A CN202110492614.6A CN202110492614A CN112904337A CN 112904337 A CN112904337 A CN 112904337A CN 202110492614 A CN202110492614 A CN 202110492614A CN 112904337 A CN112904337 A CN 112904337A
- Authority
- CN
- China
- Prior art keywords
- deformation
- vector
- target point
- axis
- slope
- 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
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
-
- 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
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种基于Offset Tracking技术的边坡形变时序监测方法,包括:步骤1:基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变;步骤2:结合地形数据和卫星轨道数据,反演目标区域的准三维形变;步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点;步骤4:综合多时间段的形变测量结果,计算目标点的时序形变。该发明能够提升Offset Tracking技术的适用性,为技术的工程化应用提供技术支持。
Description
技术领域
本发明涉及合成孔径雷达技术领域,特别涉及一种基于Offset Tracking技术的边坡形变时序监测方法。
背景技术
近年来,国内每年都会发生多起滑坡事故。一般来说,在边坡宏观失稳之前,其表面都会先发生微小位移。因此,为了更深入地研究边坡滑坡的触发机理,并在第一时间准确地实现临滑预报,我们需要对边坡进行全覆盖,长时间、连续、实时和高精度的形变监测。作为一种遥感技术手段,星载SAR技术具有监测覆盖范围大、全天时全天候观测等优势。基于星载SAR数据的InSAR技术和Offset Tracking技术都具备形变监测能力。
Offset Tracking技术的测量精度比InSAR技术低一个数量级,当图像分辨率为3m时,形变测量精度可达15cm。但是,与InSAR数据处理需要积累大量SAR图像的条件不同,Offset Tracking技术只需要两幅SAR图像就能开展形变监测,且不存在周期模糊现象,即使在干涉失相干的条件下,也能在方位向和距离向上获取几米至数十米的二维剧烈形变信息。
但是,在Offset Tracking技术的工程化应用过程中,还存在两个问题:(1)OffsetTracking技术获取的是距离向和方位向二维形变,而不是空间三维形变,难以较直观地展示空间形变特征;(2)Offset Tracking技术是采用计算窗口相关性的方式来测量形变,当窗口的选择不适当时,可能导致形变反演错误。
针对这两个问题,本发明首先基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变(距离向和方位向);然后,结合地形数据和卫星轨道数据,反演目标的准三维形变;接着,基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点;最后,综合多时间段的形变测量结果,计算目标点的时序形变。本发明能够提升OffsetTracking技术的适用性,为技术的工程化应用提供技术支持。
发明内容
本发明提供一种基于Offset Tracking技术的边坡形变时序监测方法,能够提升Offset Tracking技术的适用性,为技术的工程化应用提供技术支持。
本发明提供一种基于Offset Tracking技术的边坡形变时序监测方法,包括:
步骤1:基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变;
步骤2:结合地形数据和卫星轨道数据,反演目标区域的准三维形变;
步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点;
步骤4:综合多时间段的形变测量结果,计算目标点的时序形变。
在一种可能实现的方式中,步骤1中,基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变,包括:
基于SAR数据,并利用影像互相关算法,对形变前后的SAR图像进行粗配准和空间重采样;
选用滑动窗口,计算子像元间的相关性;
提取对应方位向和距离向的偏移量,同时拟合并去除由两幅图像成像时间和空间不同所带来的系统误差偏移值;
获得对应方位向和距离向的地表形变信息;
其中,二维形变指的是方位向和距离向的地表形变信息。
在一种可能实现的方式中,步骤2中,结合地形数据和卫星轨道数据,反演目标区域的准三维形变,包括:
以目标区域中心点为坐标原点,建立空间直角坐标系,其中,x轴正方向指向正东,y轴正方向指向正北,z轴正方向指向天空;
基于Offset Tracking技术的形变反演,获取目标区域中的目标点在距离向和方位向的二维形变,其中,距离向的形变向量dRa和方位向的形变向量dAz,在空间直角坐标系中分别表达为:
dRa=[xRa,yRa,zRa]=|dRa|[uxRa,uyRa,uzRa]=|dRa|uRa
dAz=[xAz,yAz,zAz]=|dAz|[uxAz,uyAz,uzAz]=|dAz|uAz
其中,uRa为距离向形变的单位向量,uAz为方位向形变的单位向量,|dRa|为距离向测得的形变值,|dAz|为方位向测得的形变值,uxRa为距离向形变单位向量uRa在x轴的分量,uyRa为距离向形变单位向量uRa在y轴的分量,uzRa为距离向形变单位向量uRa在z轴的分量;xRa为距离向形变向量dRa在x轴的分量,yRa为距离向形变向量dRa在y轴的分量,zRa为距离向形变向量dRa在z轴的分量;uxAa为方位向形变单位向量uAz在x轴的分量,uyAa为方位向形变单位向量uAz在y轴的分量,uzAa为方位向形变单位向量uAz在z轴的分量;xAa为方位向形变向量dAz在x轴的分量,yAa为方位向形变向量dAz在y轴的分量,zAa为方位向形变向量dAz在z轴的分量;
基于地形数据,计算每个目标点位置处的坡度角θ和坡向角φ,并在每个目标点的临近区域建立一个滑动平面,其中,所述滑动平面的方程表示为:
f(θ,φ)x+g(θ,φ)y+z=0
式中,f(θ,φ)和g(θ,φ)为滑动平面的方程系数,它们由坡度角θ和坡向角φ决定,根据空间几何知识,f(θ,φ)和g(θ,φ)通过以下代数方程求解:
基于上述理论,定义目标点的三维形变向量d=[x,y,z],基于计算的二维形变以及目标点形变位于滑动平面的基本理论假设,建立以下方程组:
通过解算上述方程组,可获取目标点的准三维形变向量d,进而反演得到目标区域的准三维形变。
在一种可能实现的方式中,步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点,包括:
基于三维形变反演结果进行第一校验;
基于SAR几何照射关系进行第二校验;
根据第一校验结果、第二校验结果,剔除测量误差点。
在一种可能实现的方式中,基于三维形变反演结果进行第一校验,包括:
获取目标点的坡度角θ和坡向角φ;
计算所述目标点沿坡向下的梯度矢量S,其数学表达式为:
S=[tanθcosφ,tanθsinφ,-1]
所述目标点的准三维形变量d与梯度矢量S的夹角φ由下式计算:
将所述夹角φ的绝对值与夹角阈值φTh进行比较,根据比较结果对测量误差点进行第一剔除。
在一种可能实现的方式中,基于SAR几何照射关系进行第二校验,包括:
基于地形数据和卫星轨道数据,获取每个目标点的对应的雷达方向向量aLOS和本地法向量aS;
确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态;
若是,对测量误差点进行第二剔除。
在一种可能实现的方式中,确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态,包括:
计算雷达LOS方向向量aLOS和本地法向量aS的夹角φ;
当夹角φ大于90°时,判定所述目标点在SAR图像中处于阴影状态;
当夹角φ的绝对值小于φTh时,判定所述目标点在SAR图像中处于叠掩状态。
在一种可能实现的方式中,步骤4中,综合多时间段的形变测量结果,计算目标点的时序形变,包括:
在目标区域内建立规则网格,并对网格上的点进行插值处理;
根据插值处理结果,获取形变序列结果;
形变序列结果,获取目标区域不同目标点的形变演化特征,并计算获取所述目标点的时序形变。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1基于OffsetTracking技术的边坡形变时序监测方法流程图;
图2由坡度和坡向角构成的目标点滑动平面示意图;
图3目标点阴影状态的判断方法;
图4目标点叠掩状态的判断方法;
图5基于OffsetTracking技术直接反演的三维形变结果;
图6误差点剔除后的三维形变结果;
图7某滑坡距离向空间累积形变量;
图8某滑坡方位向空间累积形变量;
图9某滑坡2018.11.10~2019.10.07准三维形变示意图;
图10为对应位置在距离向的空间累积形变历史曲线;
图11为对应位置在方位向的空间累积形变历史曲线。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
本发明提供一种基于Offset Tracking技术的边坡形变时序监测方法,如图1所示,包括:
步骤1:基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变;
步骤2:结合地形数据和卫星轨道数据,反演目标区域的准三维形变;
步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点;
步骤4:综合多时间段的形变测量结果,计算目标点的时序形变。
该实施例中,在进行步骤1-4的过程中,是基于DEM数据(地形高程数据)实现的。
该实施例中,通过对目标区域中的每个测量点进行校验,然后剔除其中的误差点,最后获取得到目标点。
上述技术方案的有益效果是:通过遥感测量技术手段,获取边坡区域的面状形变特征,剔除测量误差点,提升形变测量精度,特别是能适用于难以布设地面传感器的高陡边坡或危险边坡,最终为边坡形变监测和风险早期识别提供高效的技术支持。
在一种可能实现的方式中,步骤1中,基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变,包括:
基于SAR数据,并利用影像互相关算法,对形变前后的SAR图像进行粗配准和空间重采样;
选用滑动窗口,计算子像元间的相关性;
提取对应方位向和距离向的偏移量,同时拟合并去除由两幅图像成像时间和空间不同所带来的系统误差偏移值;
获得对应方位向和距离向的地表形变信息;
其中,二维形变指的是方位向和距离向的地表形变信息。
该实施例中,Offset Tracking技术利用影像互相关算法,通过追踪雷达幅度影像中特征目标的位置变化来监测地表形变。
在一种可能实现的方式中,如图2所示,步骤2中,结合地形数据和卫星轨道数据,反演目标区域的准三维形变,包括:
以目标区域中心点为坐标原点,建立空间直角坐标系,其中,x轴正方向指向正东,y轴正方向指向正北,z轴正方向指向天空;
基于Offset Tracking技术的形变反演,获取目标区域中的目标点在距离向和方位向的二维形变,其中,距离向的形变向量dRa和方位向的形变向量dAz,在空间直角坐标系中分别表达为:
dRa=[xRa,yRa,zRa]=|dRa|[uxRa,uyRa,uzRa]=|dRa|uRa
dAz=[xAz,yAz,zAz]=|dAz|[uxAz,uyAz,uzAz]=|dAz|uAz
其中,uRa为距离向形变的单位向量,uAz为方位向形变的单位向量,|dRa|为距离向测得的形变值,|dAz|为方位向测得的形变值,uxRa为距离向形变单位向量uRa在x轴的分量,uyRa为距离向形变单位向量uRa在y轴的分量,uzRa为距离向形变单位向量uRa在z轴的分量;xRa为距离向形变向量dRa在x轴的分量,为距离向形变向量dRa在y轴的分量,zRa为距离向形变向量dRa在z轴的分量;uxAa为方位向形变单位向量uAz在x轴的分量,uyAa为方位向形变单位向量uAz在y轴的分量,uzAa为方位向形变单位向量uAz在z轴的分量;xAa为方位向形变向量dAz在x轴的分量,yAa为方位向形变向量dAz在y轴的分量,zAa为方位向形变向量dAz在z轴的分量;
基于地形数据,计算每个目标点位置处的坡度角θ和坡向角φ,并在每个目标点的临近区域建立一个滑动平面,其中,所述滑动平面的方程表示为:
f(θ,φ)x+g(θ,φ)y+z=0
式中,f(θ,φ)和g(θ,φ)为滑动平面的方程系数,它们由坡度角θ和坡向角φ决定,根据空间几何知识,f(θ,φ)和g(θ,φ)通过以下代数方程求解:
基于上述理论,定义目标点的三维形变向量d=[x,y,z],基于计算的二维形变以及目标点形变位于滑动平面的基本理论假设,建立以下方程组:
通过解算上述方程组,可获取目标点的准三维形变向量d,进而反演得到目标区域的准三维形变。
该实施例中,结合滑坡的基本知识,目标点的形变向量是在该平面内的。
上述技术方案的有益效果是:便于获取准三维形变,为边坡形变监测和风险早期识别提供高效的基础。
在一种可能实现的方式中,步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点,包括:
基于三维形变反演结果进行第一校验;
基于SAR几何照射关系进行第二校验;
根据第一校验结果、第二校验结果,剔除测量误差点。
在一种可能实现的方式中,基于三维形变反演结果进行第一校验,包括:
获取目标点的坡度角θ和坡向角φ;
计算所述目标点沿坡向下的梯度矢量S,其数学表达式为:
S=[tanθcosφ,tanθsinφ,-1]
所述目标点的准三维形变量d与梯度矢量S的夹角φ由下式计算:
将所述夹角φ的绝对值与夹角阈值φTh进行比较,根据比较结果对测量误差点进行第一剔除。
在一种可能实现的方式中,基于SAR几何照射关系进行第二校验,包括:
基于地形数据和卫星轨道数据,获取每个目标点的对应的雷达方向向量aLOS和本地法向量aS;
确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态;
若是,对测量误差点进行第二剔除。
在一种可能实现的方式中,确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态,包括:
计算雷达LOS方向向量aLOS和本地法向量aS的夹角;
当夹角φ大于90°时,判定所述目标点在SAR图像中处于阴影状态;
当夹角φ的绝对值小于φTh时,判定所述目标点在SAR图像中处于叠掩状态。
该实施例中,在Offset Tracking的算法中,是基于主辅SAR图像窗口内数据的相关性来计算形变量。一般情况下,窗口的选取会影响算得的形变量。因此,Offset Tracking技术直接生成的形变量可能在局部区域包含错误的测量结果。
该实施例中,针对第一校验,一般来说,目标点主要的滑动方向是沿着滑动平面向下的。因此,每个目标点计算的准三维形变向量的方向与该点对应的地形梯度方向应基本相似,可通过这个条件来剔除形变测量明显出错的目标点。
通过设定合适的夹角阈值φTh,选出明显存在误差的形变测量结果(即|φ|>φTh),进而将其选出并剔除。
该实施例中,针对第二校验,考虑到处于阴影或叠掩状态下的目标信号质量较差,在后续分析过程中将剔除这些测量点。
如图3所示,在判断目标点阴影状态的过程中,可先计算雷达LOS方向向量aLOS和本地法向量aS的夹角φ。当夹角φ大于90°时,可认为该目标点在SAR图像中处于阴影状态。
如图4所示,在判断目标点叠掩状态的过程中,可先计算雷达LOS方向向量aLOS和本地法向量aS的夹角φ。当夹角φ的绝对值小于φTh时,可认为该目标点在SAR图像中处于叠掩状态,一般来说,φTh可设置为5°
上述技术方案的有益效果是:通过基于地形数据对形变监测结果进行校验,提升形变监测结果的可靠性。
在一种可能实现的方式中,步骤4中,综合多时间段的形变测量结果,计算目标点的时序形变,包括:
在目标区域内建立规则网格,并对网格上的点进行插值处理;
根据插值处理结果,获取形变序列结果;
形变序列结果,获取目标区域不同目标点的形变演化特征,并计算获取所述目标点的时序形变。
上述技术方案的有益效果是:将不同时期的测量结果进行时序处理,获取目标区域不同位置的形变演化特征,更好地研究目标区域在长时间范围内的形变特征。
针对上述技术方案,还具体通过下文对某边坡进行监测来进行具体说明:
SAR数据采用3m分辨率的COSMO-SkyMed数据,其具体参数信息如表1所示。
表1 COSMO-SkyMed数据参数
步骤一:
以20181110和20181122两期数据为例,首先采用OffsetTracking技术获得该时间间隔内的方位向和距离向的二维地表形变成果,基于此二维形变成果,结合准三维形变反演模型得到准三维形变反演结果,结果如图5所示。由图5可知,初步的准三维形变成果存在误差点。
步骤二:
基于本方法提出的模型,结合地形和卫星轨道数据可剔除掉形变测量误差点,进而能提取出目标区域准确的三维形变结果,结果如图6所示。时间上相邻的两期数据依次按照上述流程处理,得到目标区域相邻时间间隔内的准三维形变信息。
步骤三:
基于上述得到的各期二维和准三维形变成果,考虑到不同时间段的形变测量结果是由不同的SAR数据对生成,其测量点的位置可能存在差异,先在目标区域内建立规则网格,然后对网格上的点进行插值处理,最后获取形变序列结果。该滑坡在距离向和方位向的累积形变分别如图7和图8所示。合成的形变量及形变方向示意图如图9所示,其最大累积形变量达15m。在该滑坡中选择关键测量目标,其位置在图7中用蓝色五角星标识,该位置在距离向的空间累积形变历史曲线如图10所示,该位置在方位向的空间累积形变历史曲线如图11所示。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (8)
1.一种基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,包括:
步骤1:基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变;
步骤2:结合地形数据和卫星轨道数据,反演目标区域的准三维形变;
步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点;
步骤4:综合多时间段的形变测量结果,计算目标点的时序形变。
2.如权利要求1所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,步骤1中,基于SAR数据,采用Offset Tracking技术计算目标区域的二维形变,包括:
基于SAR数据,并利用影像互相关算法,对形变前后的SAR图像进行粗配准和空间重采样;
选用滑动窗口,计算子像元间的相关性;
提取对应方位向和距离向的偏移量,同时拟合并去除由两幅图像成像时间和空间不同所带来的系统误差偏移值;
获得对应方位向和距离向的地表形变信息;
其中,二维形变指的是方位向和距离向的地表形变信息。
3.如权利要求1所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,步骤2中,结合地形数据和卫星轨道数据,反演目标区域的准三维形变,包括:
以目标区域中心点为坐标原点,建立空间直角坐标系,其中,x轴正方向指向正东,y轴正方向指向正北,z轴正方向指向天空;
基于Offset Tracking技术的形变反演,获取目标区域中的目标点在距离向和方位向的二维形变,其中,距离向的形变向量dRa和方位向的形变向量dAz,在空间直角坐标系中分别表达为:
dRa=[xRa,yRa,zRa]=|dRa|[uxRa,uyRa,uzRa]=|dRa|uRa
dAz=[xAz,yAz,zAz]=|dAz|[uxAz,uyAz,uzAz]=|dAz|uAz
其中,uRa为距离向形变的单位向量,uAz为方位向形变的单位向量,|dRa|为距离向测得的形变值,|dAz|为方位向测得的形变值,uxRa为距离向形变单位向量uRa在x轴的分量,uyRa为距离向形变单位向量uRa在y轴的分量,uzRa为距离向形变单位向量uRa在z轴的分量;xRa为距离向形变向量dRa在x轴的分量,yRa为距离向形变向量dRa在y轴的分量,zRa为距离向形变向量dRa在z轴的分量;uxAa为方位向形变单位向量uAz在x轴的分量,uyAa为方位向形变单位向量uAz在y轴的分量,uzAa为方位向形变单位向量uAz在z轴的分量;xAa为方位向形变向量dAz在x轴的分量,yAa为方位向形变向量dAz在y轴的分量,zAa为方位向形变向量dAz在z轴的分量;
基于地形数据,计算每个目标点位置处的坡度角θ和坡向角φ,并在每个目标点的临近区域建立一个滑动平面,其中,所述滑动平面的方程表示为:
f(θ,φ)x+g(θ,φ)y+z=0
式中,f(θ,φ)和g(θ,φ)为滑动平面的方程系数,它们由坡度角θ和坡向角φ决定,根据空间几何知识,f(θ,φ)和g(θ,φ)通过以下代数方程求解:
基于上述理论,定义目标点的三维形变向量d=[x,y,z],基于计算的二维形变以及目标点形变位于滑动平面的基本理论假设,建立以下方程组:
通过解算上述方程组,可获取目标点的准三维形变向量d,进而反演得到目标区域的准三维形变。
4.如权利要求1所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,步骤3:基于地形数据,校验每个测量点的准三维形变测量结果,剔除测量误差点,包括:
基于三维形变反演结果进行第一校验;
基于SAR几何照射关系进行第二校验;
根据第一校验结果、第二校验结果,剔除测量误差点。
6.如权利要求4所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,基于SAR几何照射关系进行第二校验,包括:
基于地形数据和卫星轨道数据,获取每个目标点的对应的雷达方向向量aLOS和本地法向量aS;
确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态;
若是,对测量误差点进行第二剔除。
7.如权利要求6所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,确定雷达LOS方向向量aLOS和本地法向量aS的夹角关系,判断所述目标点是否处于阴影和叠掩状态,包括:
计算雷达LOS方向向量aLOS和本地法向量aS的夹角φ;
当夹角φ大于90°时,判定所述目标点在SAR图像中处于阴影状态;
当夹角φ的绝对值小于φTh时,判定所述目标点在SAR图像中处于叠掩状态。
8.如权利要求1所述的基于Offset Tracking技术的边坡形变时序监测方法,其特征在于,步骤4中,综合多时间段的形变测量结果,计算目标点的时序形变,包括:
在目标区域内建立规则网格,并对网格上的点进行插值处理;
根据插值处理结果,获取形变序列结果;
形变序列结果,获取目标区域不同目标点的形变演化特征,并计算获取所述目标点的时序形变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110492614.6A CN112904337B (zh) | 2021-05-07 | 2021-05-07 | 一种基于Offset Tracking技术的边坡形变时序监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110492614.6A CN112904337B (zh) | 2021-05-07 | 2021-05-07 | 一种基于Offset Tracking技术的边坡形变时序监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112904337A true CN112904337A (zh) | 2021-06-04 |
CN112904337B CN112904337B (zh) | 2021-11-05 |
Family
ID=76109006
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110492614.6A Active CN112904337B (zh) | 2021-05-07 | 2021-05-07 | 一种基于Offset Tracking技术的边坡形变时序监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112904337B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113848551A (zh) * | 2021-09-24 | 2021-12-28 | 成都理工大学 | 一种利用InSAR升降轨形变数据的滑坡深度反演方法 |
CN113960596A (zh) * | 2021-10-20 | 2022-01-21 | 苏州深蓝空间遥感技术有限公司 | 一种基于北斗和PS-InSAR的滑坡三维形变监测方法 |
CN114236541A (zh) * | 2021-12-08 | 2022-03-25 | 电子科技大学 | 基于Sentinel-1卫星SAR图像的大面积地表三维形变计算方法 |
CN117194928A (zh) * | 2023-11-07 | 2023-12-08 | 湖南中云图地理信息科技有限公司 | 基于gnss的地理形变监测系统 |
WO2024004215A1 (ja) * | 2022-07-01 | 2024-01-04 | 株式会社Synspective | 地盤変動解析装置および地盤変動解析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645476A (zh) * | 2013-12-18 | 2014-03-19 | 中国国土资源航空物探遥感中心 | 一种合成孔径雷达差分干涉图序列的时空同质滤波方法 |
CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
CN106066478A (zh) * | 2016-05-27 | 2016-11-02 | 中国矿业大学 | 融合像元偏移跟踪和短基线集的矿区地表形变解算方法 |
CN106203271A (zh) * | 2016-06-29 | 2016-12-07 | 南京大学 | 一种基于高分辨率sar图像相干性的高铁干线提取方法 |
CN109471104A (zh) * | 2018-11-15 | 2019-03-15 | 首都师范大学 | 一种从两平行轨道sar数据中获取地表三维移动量的方法 |
US20190131988A1 (en) * | 2017-11-02 | 2019-05-02 | Analog Devices, Inc. | Comparator error suppression |
-
2021
- 2021-05-07 CN CN202110492614.6A patent/CN112904337B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645476A (zh) * | 2013-12-18 | 2014-03-19 | 中国国土资源航空物探遥感中心 | 一种合成孔径雷达差分干涉图序列的时空同质滤波方法 |
CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
CN106066478A (zh) * | 2016-05-27 | 2016-11-02 | 中国矿业大学 | 融合像元偏移跟踪和短基线集的矿区地表形变解算方法 |
CN106203271A (zh) * | 2016-06-29 | 2016-12-07 | 南京大学 | 一种基于高分辨率sar图像相干性的高铁干线提取方法 |
US20190131988A1 (en) * | 2017-11-02 | 2019-05-02 | Analog Devices, Inc. | Comparator error suppression |
CN109471104A (zh) * | 2018-11-15 | 2019-03-15 | 首都师范大学 | 一种从两平行轨道sar数据中获取地表三维移动量的方法 |
Non-Patent Citations (1)
Title |
---|
姚丹丹: "基于偏移追踪法的D-InSAR矿区形变监测技术", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113848551A (zh) * | 2021-09-24 | 2021-12-28 | 成都理工大学 | 一种利用InSAR升降轨形变数据的滑坡深度反演方法 |
CN113960596A (zh) * | 2021-10-20 | 2022-01-21 | 苏州深蓝空间遥感技术有限公司 | 一种基于北斗和PS-InSAR的滑坡三维形变监测方法 |
CN114236541A (zh) * | 2021-12-08 | 2022-03-25 | 电子科技大学 | 基于Sentinel-1卫星SAR图像的大面积地表三维形变计算方法 |
CN114236541B (zh) * | 2021-12-08 | 2023-05-16 | 电子科技大学 | 基于Sentinel-1卫星SAR图像的大面积地表三维形变计算方法 |
WO2024004215A1 (ja) * | 2022-07-01 | 2024-01-04 | 株式会社Synspective | 地盤変動解析装置および地盤変動解析方法 |
CN117194928A (zh) * | 2023-11-07 | 2023-12-08 | 湖南中云图地理信息科技有限公司 | 基于gnss的地理形变监测系统 |
CN117194928B (zh) * | 2023-11-07 | 2024-01-26 | 湖南中云图地理信息科技有限公司 | 基于gnss的地理形变监测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112904337B (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112904337B (zh) | 一种基于Offset Tracking技术的边坡形变时序监测方法 | |
US20200103530A1 (en) | Method for extracting elevation control point with assistance of satellite laser altimetry data | |
Zhu et al. | Geodetic SAR tomography | |
Wang et al. | Improved SAR amplitude image offset measurements for deriving three-dimensional coseismic displacements | |
US8154435B2 (en) | Stability monitoring using synthetic aperture radar | |
CN106526593B (zh) | 基于sar严密成像模型的子像素级角反射器自动定位方法 | |
CN112986949B (zh) | 针对角反射器的sar高精度时序形变监测方法和装置 | |
CN106526590A (zh) | 一种融合多源sar影像工矿区三维地表形变监测及解算方法 | |
CN109116321B (zh) | 一种星载干涉成像高度计的相位滤波方法及高度测量方法 | |
EP3156818B1 (en) | Automatic three-dimensional geolocation of sar targets and simultaneous estimation of tropospheric propagation delays using two long-aperture sar images | |
CN110068817B (zh) | 一种基于激光测距和InSAR的地形测图方法、仪器和系统 | |
Ji et al. | An autofocus approach with applications to ionospheric scintillation compensation for spaceborne SAR images | |
CN112346027B (zh) | 用于确定合成孔径雷达图像的散射特性的方法和系统 | |
CN116908789B (zh) | 地基合成孔径雷达干涉测量建筑立面变形信息提取方法 | |
CN108646244B (zh) | 测量建筑五维形变的分析方法及系统 | |
Duro et al. | High resolution differential interferometry using time series of ERS and ENVISAT SAR data | |
CN116299466B (zh) | 一种输电通道地质形变监测方法和装置 | |
CN105403886A (zh) | 一种机载sar定标器图像位置自动提取方法 | |
CN117388851A (zh) | 一种InSAR变形监测精度可靠估计方法 | |
CN112505647A (zh) | 一种基于序贯子图像序列的动目标方位速度估计方法 | |
CN117073621A (zh) | 融合InSAR和北斗的城市沉降监测方法 | |
CN115201823B (zh) | 一种利用BDS-InSAR数据融合的地表形变监测方法 | |
Yan et al. | Determining suitable spaceborne SAR observations and ground control points for surface deformation study in rugged terrain with InSAR technique | |
CN115712095A (zh) | 基于单个角反射的sar卫星三维定位误差改正方法及系统 | |
Ge et al. | Integrating corner reflectors and PSInSAR technique to monitor regional land subsidence |
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 |