TS-InSAR原理与应用(概括)

1. InSAR/D-InSAR技术原理

与传统遥感技术通过幅度信息来获取地表信息不同,合成孔径雷达干涉测量(InSAR)通过对信号的相位信息进行干涉来得到地表的高精度DEM或形变信息。

InSAR利用对相同区域拍摄的两景影像,经过干涉后获取该区域的干涉图,干涉条纹中带有地表的高程信息和两景影像拍摄期间的地表形变信息。如果两景影像拍摄的时间点很接近,则可以忽略掉干涉条纹中的地表形变信息,进而得到地表的高程。

相反,如果已有地表的高程信息,则可以将地表起伏引起的干涉条纹从原始条纹中减去,得到只包含地表形变信息的干涉条纹,即为差分InSAR(D-InSAR)技术。理论上,D-InSAR获取的地表形变精度可达到毫米(mm)级。

2. TS-InSAR技术原理

虽然理论上通过D-InSAR获取的地表形变精度可达mm级,但其只能获得离散的地表形变。并且,由于D-InSAR没有多余观测,得到的结果误差较大。

时序InSAR(TS-InSAR)技术首先获取覆盖同一区域的多景影像,将所有影像配准到同一幅主影像后,根据多个影像干涉对之间的多余观测来获取地表的位移时间序列,并减弱结果的误差。目前较为常用的TS-InSAR技术有持久散射体技术(PS-InSAR)、小极限集技术(SBAS-InSAR)、分布式散射体技术(DS-InSAR)。

2.1 PS-InSAR

PS-InSAR的思想是选取在整个观测时间内始终保持高相干性的点(如人工建筑物、道路、岩石等)构建观测网,两个PS点之间的连接线成为弧段。在对相邻PS点进行差分后,以弧段作为基础观测值,进行地表位移序列的反演。目前较为常用的PS-InSAR软件为StaMPS、GAMMA软件的IPTA模块。

2.2 SBAS-InSAR

在植被茂盛的区域,可使用的PS点较少,PS-InSAR的结果精度受限。SBAS-InSAR首先根据预先设定的时间和空间阈值构建小基线集合,之后通过最小二乘方法获取地表形变序列。

3.3 DS-InSAR

DS-InSAR认为分布式散射体(如裸露的地表、灌木丛等)有着相同的散射特性,因此对分布式散射体进行自适应滤波可以提高信号的信噪比,进而提高结果的精度。目前,DS-InSAR是主流的TS-InSAR技术。

3 TS-InSAR应用

TS-InSAR的应用范围非常广泛,如:滑坡、地震、采矿、地面沉降、冰川、工程建设等。在此只举一些例子,在以后的文章中会展开描述。

第一次发表文章,有问题随时可以私信我,常驻!

  • 1
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
D-InSAR(Differential Synthetic Aperture Radar Interferometry)是一种通过SAR(Synthetic Aperture Radar)图像以及多次采集的SAR图像的相位信息来测量地形变的技术。其基本原理是通过比较两张SAR图像的相位差异,从而计算出地表的形变情况。 D-InSAR主要分为以下几个步骤: 1. 预处理 预处理是D-InSAR的第一步,主要目的是去除大气影响、地球表面的影响、以及减小地球电离层的影响。预处理主要包括以下步骤: - 平场校正(Flat Earth Removal),消除由于地球曲率引起的相位误差; - 相位解缠(Phase Unwrapping),将相位展开,使得相位差在$[-\pi,\pi]$之间; - 干涉滤波(Interferometric Filtering),去除非地面信号,如大气、杂波等; - 水平位移校正(Horizontal Offset Correction),校正由于不同的SAR成像几何条件产生的影响,如地形和地物的不同造成的影响。 2. 生成干涉图 生成干涉图的主要目的是计算两张SAR图像之间的差异,从而得到形变的信息。干涉图的公式为: $$I(x,y)=\frac{S_1(x,y)S_2^*(x,y)}{|S_1(x,y)S_2^*(x,y)|}$$ 其中,$S_1(x,y)$和$S_2(x,y)$分别表示两张SAR图像在$(x,y)$处的复振幅,$^*$表示共轭复数运算。 3. 相位差计算 相位差是指两张SAR图像在同一地点的相位差异。相位差的公式为: $$\Delta \phi(x,y)=\phi_2(x,y)-\phi_1(x,y)$$ 其中,$\phi_1(x,y)$和$\phi_2(x,y)$分别表示两张SAR图像在$(x,y)$处的相位。 4. 三维变形反演 通过对相位差进行分析,可以得到地表的形变情况。三维变形反演的公式为: $$\Delta h(x,y)=\frac{\Delta \phi(x,y)\lambda}{4\pi \sin\theta}$$ 其中,$\Delta h(x,y)$表示在$(x,y)$处的高程变化量,$\lambda$表示雷达波长,$\theta$表示入射角。 以上就是D-InSAR的基本原理及各步骤公式。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值