合成孔径雷达(Synthetic Aperture Radar,SAR)的概念始于20世纪50年代,是正在发展中的极具潜力的微波遥感技术。SAR具有全天时、全天候的工作能力,能够穿透云层,对某些地物具有一定的穿透能力,并能够在不同的微波频段、不同极化状态下得到地面目标的高分辨率图像,因此很快引起各国的热切关注,得到飞速发展。

合成孔径雷达干涉技术(Interferometric Synthetic Aperture Radar,InSAR)出现于20世纪60年代末,它是SAR与射电天文学干涉测量技术结合的产物。当扫过地面同一目标区域时,利用成像几何关系,通过成像、一些特殊的数据处理和几何转换,即可提取地表目标区域的高程信息和形变信息。由于InSAR技术有效利用了SAR的回波相位信息,测高精度为米级甚至亚米级,而一般雷达立体测量方法只利用灰度信息来实现三维制图,测高精度仅能达到数十米,因此该技术迅速引起了地学界及相关领域科研工作者的极大兴趣,现已成为微波遥感领域的研究热点

D-InSAR(Different InSAR,差分干涉)技术是在InSAR的基础上发展起来的,它以合成孔径雷达复数据提供的相位信息为信息源,可从包含日标区域地形和形变等信息的一幅或多幅干涉纹图中提取地面目标的微小形变信息。D-InSAR具有高形变敏感度、高空间分辨率、几乎不受云雨天气制约和空中遥感等突出的技术优势,因而有人认为它是独特的基于面观测的空间大地测量新技术,可补充已有的基于点观测的低空间分辨率大地测量技术如全球定位系统(GPS)、甚长基线干涉(VLBI)和精密水准等,从而可以揭示出更多的地球物理现象,最终为地球物理学提供一种全新的动态研究途径。

2 InSAR、D-InSAR基本原理

InSAR测量模式主要有两种:一种是双天线单轨(Single Pass)模式,主要用来生成数字高程模型,一般用于机载SAR;另一种是双轨(Two Pass)模式,主要用于获取地表变形,一般用于星载SAR。下面以重复轨道干涉测量为例,简要介绍InSAR技术的基本原理(见图1)。假设卫星以一定的时间间隔和轨道偏离(通常为几十米到一公里左右)重复对某一区域成像,并在两次飞行过程中处于不同的空间位置S1和S2,则空间干涉基线向量为B,长度用B表示;基线向量B与水平方向的夹角为基线倾角α;S1和S2至地面点P的斜距分别为R和R+△R;将基线沿视线方向分解,得到平行于和垂直于视线向的分量Bp、B;H为S1到参考面的高度;从S1发射波长为λ的信号经目标点P反射后被S1接收,得到测量相位Φ1

同样,另一空间位置S2上测量到相位Φ2

式中arg{u1}和arg{u2}表示不同散射特性造成的随机相位。假设两幅图中随机相位的贡献相同,则S1和S2关于目标P点的相位差

Φ也称为干涉相位(Interferometric Phase),可由经过配准的两幅SAR SLC(Single Look Complex,单视复数)图共扼相乘得到。根据图1中的几何关系并利用余弦定理可得:

(4)、(5)两式即为InSAR确定高程的原理性公式

监测形变的D-InSAR技术的基本思想如下:首先选取重复时间间隔尽量短的SLC图像对,用来生成反映数字高程模型(Digital Elevation Model,DEM)的干涉相位图;然后选择跨越同一目标区域形变发生时间段的SLC图像对生成包含地形和形变信息的相位图;将前后两幅相位图相减,即得到反映视线向形变的相位图;再根据相应的成像几何关系计算垂直、水平方向上的形变。

理论上,由相位解缠(相位模糊度解算)、成像几何和轨道参数就可以重建DEM,但实际的处理过程相当复杂,一般包括SAR信号预处理、图像配准、生成干涉图、去除平地效应、相位解缠、基线估计、地理编码、DEM建立等过程。长期以来,成像处理、复数影像的精配准、相位解缠、大气效应的改正等都是InSAR数据处理的重点和难点。从InSAR理论提出以来,大部分研究工作都集中在这些方面。但是,至今相位解缠以及大气效应的改正仍未得到很好解决。

3 InSAR、D-InSAR在地面沉降监测中的应用

作为一种新兴的地面形变研究方法,InSAR技术在地面沉降监测方面发挥了愈来愈明显的作用,国内外已有诸多实例。Biegert等(1997)应用不同卫星在美国加利福尼亚州Belridge和Lost山油田重复测量的合成孔径雷达数据对该区的地面沉降进行了研究,结果显示70天内沉降量达到6厘米,此结果与该区每年30厘米的地面沉降速率相吻合。Marco van der (2001)对该油田地面沉降的研究也证明了InSAR技术用于地面沉降的可行性。李德仁等(2000)利用欧空局ERS-1和ERS-2相隔1天的重复轨道SAR数据,经过差分处理对天津市地面沉降进行研究,得到反映地面沉降大小及分布的干涉条纹图。此图与1995~1997年重复水准测量求得的地面沉降等值线图比较,具有明显的一致性和相似性。

刘国祥等(2001)用卫星雷达差分干涉技术研究在近海回填地基上建立香港赤腊角机场的稳定性,获得该机场在近1年内的非均匀沉降场,地面分辨率为20米x20米,在填海区域内下沉量呈0~50米的空间分布,与离散水准测量结果吻合较好(相关系数0.89)。证实了ERS-2干涉系统对微小地面沉降敏感度高,精度小于1厘米。

4发展前景

GPS是20世纪70年代由美国陆海空三军联合研制的新一代空间卫星导航定位系统。随着GPS接收机硬件性能的提高和软件处理技术的进步,GPS相对定位的精度从以前的10-7提高到10-9量级,其在地面沉降应用中的精度已达到亚毫米级。

GPS和InSAR起源于大地测量和遥感两个不同的领域,但它们技术上的互补性使得它们可以共同应用于在获得空间三维信息和地表变形监测等领域。目前世界上许多国家建立了CGPS网(Continuous GPS networks),用于连续和动态地监测本地区的地面沉降情况。大部分CGPS网的时间分辨率为30秒,基本满足实时动态监测的要求,而空间分辨率却只有几十公里,不能满足小区域地面沉降监测的要求。而InSAR技术虽然具有很高的空间分辨率,星载SAR已达到约20米,但是雷达卫星因其固有的运行周期,所以时间分辨率还满足不了要求。另外大气参数的变化,卫星轨道误差等问题单纯依靠InSAR数据本身难以解决,必须加入其它的辅助数据和必要的技术手段来加以改善。

由GPS网络得到的数据可以计算出对流层的水蒸气含量和电离层干扰,这些结果可以修正InSAR中的大气模型。由于GPS的坐标系统与大地坐标系统是联系在一起的,所以在这种意义上,GPS获得的坐标系可以认为是“绝对的”。相反,InSAR结果可认为是“相对”测量值。另外,由于InSAR影像的高空间分辨率性,它在空间意义上可以用来加密GPS结果。所以显然两种技术是互补的。

对于地面沉降监测来说,具体的做法如下:第一步,根据GPS数据获得大气模型改正,给出水蒸气可降水量的预测值,从而得出电离层延迟。第二步,以GPS定位结果作为约束条件来减少或消除SAR卫星轨道误差。第三步,进行内插,首先在空间领域利用InSAR结果对GPS网内插,然后在时间领域利用动态模型对加密的网格点进行内插,可以使用自适应滤波的方法建立动态模型。第四步,在双内插结果的基础上,利用卡尔曼滤波往前滤波方法对格网中所有点进行估计,从而得到某一时刻的变形值,这就达到了时间域和空间域的双估计。由于入射角的关系,InSAR技术对高程信息特别敏感,理论上探测精度可达亚厘米级,上已述及,实际中受各种因素的影响,即使利用常规观测结果进行评价其探测精度只能达到厘米级或毫米级。提出用GPS观测数据修正InSAR结果中的诸如对流层、电离层、卫星轨道等误差,必然会提高InSAR技术的实际探测精度,使其接近或达到理论精度

InSAR在地面沉降监测中的应用及发展前景的更多相关文章

  1. Linux在Hyper-V中实现与Windows的宽带共享上网

    相信不少读者都会做(或者曾经做过)这么一项活动——把一个Linux系统安装到Windows系统下的虚拟机软件中,然后在Windows这个大环境中对Linux进行学习或者一些实验操作.在进行这么一项活动 ...

  2. Linux中open函数以及退出进程的函数

    open函数的flag详解1 读写权限:O_RDONLY O_WRONLY O_RDWR (1)linux中文件有读写权限,我们在open打开文件时也可以附带一定的权限说明 (譬如O_RDONLY就表 ...

  3. zabbix web监测

    web monitoring(监测)属于业务监控,用来监控Web站点多方面的可用性,可以监控Web站点的下载速度.返回码和响应时间.Zabbix能够检测HTML中包含的预先定义的字符串,也可以模拟登录 ...

  4. KAFA 监测| Kafka监测的方法和工具

    1.目标 在我们上一篇Kafka教程中,我们讨论了Kafka Tools.今天,我们将看到Kafka Monitoring.在此,我们将学习如何监控Apache Kafka的概念.此外,我们将涵盖在故 ...

  5. javascript中的BOM

    浏览器对象模型BOM,提供了访问浏览器的接口.这些功能大多和网页内容无关,多年来,由于缺乏规范导致BOM中的不同方法在不同浏览器中的实现有所差异,直到html5,才将BOM的主要方面纳入规范. BOM ...

  6. intouch与PLC之间通讯状态监测和设置

    intouch与PLC进行通讯状态监测中,一般做法需要PLC来实施主动脉冲计数,或者bool变化来实现.本文通过上位机自带参数设置,实现对intouch通讯状态监视,将画面恢复初始状态,并及时弹窗报警 ...

  7. Java 中的监控与管理原理概述

    点赞再看,动力无限.Hello world : ) 微信搜「程序猿阿朗 」. 本文 Github.com/niumoo/JavaNotes 和 程序猿阿朗博客 已经收录,有很多知识点和系列文章. 当前 ...

  8. Echarts(一)

    echarts3.61.<!-- 为ECharts准备一个具备大小(宽高)的Dom --> <div id="barMain" style="heigh ...

  9. 【遥感专题系列】微波遥感(三、SAR图像特征)

    SAR是主动式侧视雷达系统,且成像几何属于斜距投影类型.因此SAR图像与光学图像在成像机理.几何特征.辐射特征等方面都有较大的区别.在进行SAR图像处理和应用前,需要了解SAR图像的基本特征. 本文主 ...

随机推荐

  1. OWIN katana注册中间件的几种写法

    首先特别说明下在startup中注册完中间件的两个注意事项,看到有人写的东西有误导人的作用.关于startup启动发现类的内容,参照这里 http://www.asp.net/aspnet/overv ...

  2. MySQL索引的查看创建和删除

    1.索引作用 在索引列上,除了上面提到的有序查找之外,数据库利用各种各样的快速定位技术,能够大大提高查询效率.特别是当数据量非常大,查询涉及多个表时,使用索引往往能使查询速度加快成千上万倍. 例如,有 ...

  3. Entity Framework4.0 (七) EF4的存储过程

    前面了解了EF4的CRUD的操作,你会发现EF4使用起来比较简单的.呵呵,之前我们使用数据库的时候,有时会使用存储过程代替在代码中直接使用SQL语句. 使用存储过程的好处: 提高效率:因为存储过程是经 ...

  4. POJ3616 Milking Time 简单DP

    注意0,1,.....,N是时间点,i~i+1是时间段 然后就是思路:dp[i]代表到时间点 i 获得的最大价值, 1:dp[i]=max(dp[i],dp[s-r]+e),表示有以s为开头,i为结尾 ...

  5. centos ssh 免密码登录

    最近在学习的过程中遇到这个问题: A主机和B主机: A 免密码登录B: 首先在A的 ~/.ssh 目录中 执行 ssh-keygen -t rsa 一路回车 最后生成连个文件: 将id_rsa.pub ...

  6. zabbix监控域名带宽

    代码地址:https://github.com/Ma-Jing/python/blob/master/ngxv2_traffic_daemon.py READ.md里有使用说明! #!/usr/bin ...

  7. POJ 1113&&HDU 1348

    题意:凸包周长+一个完整的圆周长.因为走一圈,经过拐点时,所形成的扇形的内角和是360度,故一个完整的圆. 模板题,之前写的Graham模板不对,WR了很多发....POJ上的AC代码 #includ ...

  8. uvalive 4992 Jungle Outpost

    题意:一个凸边型,目标在凸边型内且最优.问最多删除几个点使目标暴露在新凸边型外面. 思路:二分+半平面相交. #include<cstdio> #include<cmath> ...

  9. Bzoj 2006: [NOI2010]超级钢琴 堆,ST表

    2006: [NOI2010]超级钢琴 Time Limit: 20 Sec  Memory Limit: 552 MBSubmit: 2222  Solved: 1082[Submit][Statu ...

  10. mysql 清空表的两种方法

    一.Delete DELETE FROM `table`; 二.Truncate TRUNCATE `table`; 第一种方法其实就是去掉where条件,没有了条件,也就是删除掉表里面的所有记录了: ...