InSAR技术在形变灾害监测中的研究进展

  • 格式:pdf
  • 大小:235.85 KB
  • 文档页数:5

图2
InSAR 数据处理过程
2
InSAR 技术最新进展国内外学者源自过对大量 SAR 图像的研究, 发现
在城区和岩石地区, 在长时间范围内仍然存在相位 和幅度变化稳定的点。利用这些稳定点上的相位特 征 , 可以进行长时间尺度上的地表形变分析。最新
GPS 数据和水准数据进行了对比, 表明了 PS 测量的 [ 4] 准确性和有 效性。 W oraw attanam ateeku l等人 研 究了泰国 Bangkok 城区的形变 , Jun 等人研究了北京 地区的形变, 他们都取得了很精确的形变结果。 PS 技术的引入, 大大提高了 InSAR 技术监测地表形变 的精度, 使该方法的应用前景更加广阔。 2 3 小基线集 ( SBAS) 方法 Berard in o 和 L anari 等 人 提 出 了 小 基 线 集
用合成孔径雷达干涉测量技术 ( InS AR)监测地表 形变是目前地质灾害研究的热点。 InSAR 以 SAR 复数 据所承载的相位信息为信息源获取地表的三维信息和 变化信息, 它通过两副天线同时观测 (单轨模式 ), 或者 两次近平行观测 (重复轨道模式 ), 获取地面同一景观 的复图像对, 根据复图像上产生的相位差生成干涉纹 图, 干涉纹图中包含了斜距向上的点与两天线位置之 间的精确信息。因此, 利用传感器高度、 雷达波长、 波 束视向及天线基线距之间的几何关系, 可以精确地测 量出图像上每一点的三维位置和变化信息 。运用 InSAR 技术进行形变灾害监测, 不仅能够反映地表形 变的静态信息, 而且可以反映地表形变的动态变化规 律, 通过对形变区的相关测量, 可以探测地表形变的早 期信号, 有助于确定形变的范围和预测潜在事件发生 的可能性, 通过记录相关特征的变化史, 建立灾害的演 变模型, 有效地制定减灾措施。
3 InSAR 技术在形变灾害监测中的 应用现状
3 1 地面沉降监测 地面沉降是一个世界性的难题 , 均匀地地表下 沉不会对建筑物或城市基础设施产生大的影响 , 但 不均匀地地面下沉会严重威胁城市建筑和公路桥梁 的安全。由于城市周围一般地势比较平坦 , 植被相 对较少 , 雷达反射的相关性较好 , 容易获得质量较好 的干涉图像 , 因而是进行 InSAR 技术应用的理想领 [ 6] 域。国内外已有诸多实例。 B iegert等 ( 1997) 应 用美国加州 Be lridge和 L ost油田重复测量的合成孔 径雷达数据对该区的地而沉降进行了研究 , 结果显 示 70 d内沉降量达到 6 cm, 此结果与该区每年 30 cm 的地面沉降速率相吻合。 M arcovander( 2001) 对 该油田地面沉降的研究也证明了 InSAR 技术用于 地面沉降的可行性。李德仁等 ( 2000) 利用欧空局 ERS 1 和 ERS 2 相隔 1d 的重复轨 道 SAR 数据, 经 过差分处理对天津市地面沉降进行研究, 得到反应 地 面 沉 降 大 小 及 分 布 的 干 涉 条 纹 图。此 图 与
2 般通过曲线插值和拟合的方法来实现。
淮海工学院学报 ( 自然科学版 )
[ 2]
2009 年 12 月
( 2) 干涉图像形成 : 把两幅图像共扼相乘产生 干涉复数图像, 即干涉相位图, 其中包含目标点的高 度信息。为降低斑点嗓声 , 单视复数图像还要经过 多视平均 !处理得到多视复数图像 , 但这会降低空 间分辨率。 ( 3) 平地相位去除: I N SAR 干涉图中的相位由 两部分组成 , 一是地形的相对高度变化, 二是平地随 距离向位置的不同所引起的变化, 即平地在干涉条 纹图中所表现出来的距离向和方位向的周期变化。 为了简化相位展开处理 , 在相位解缠前先要把平地 效应去掉, 这个过程称为平地相位去除。 ( 4) 相位解缠: 由于干涉相位图的条纹是由相 位差形成的 , 它与地面位置直接相关的相位 以 2# 为模, 即只能测量出不足一个周期内的相位差 , 包含 了 2n#的模糊度。因此为了计算目标点的高程 , 必 须解出相位模糊度, 给相位侧量值加上 2n #的相位 周期。这种求解 2#模糊性问题的技术称为相位解 缠 , 它类似于 GPS 相位 测量中的整周模糊度问题 , 是一项重要 而有难度的工 作。 I N SAR 相位解缠方 法主要有分支切片法、 条纹检测法、 K al m an 滤波法 和最小二乘法等。 ( 5) 基线确定: 传统的基线确定方法是利用两 幅 SAR 图像之间的成像处理来补偿由卫星轨道参 数计算得到的基线。这种方法简单 , 但精度差。现 在采用地面标志点方法可以得到精确的基线估计。 ( 6) 数字高程模型 ( DEM ) 的建立 : 得到相位和 基线参数之后可计算出地形高度, 从而建立 DEM。 InSAR 数据处理流程如图 2 所示。
收稿日期 : 2009 10 11 ; 修订日期 : 2009 11 10 基金项目 : 江苏省测绘科研项目资助 ( JSCHKY 200902) , 江苏省高校自然科学研究项目资助 ( 09K JB420001) 作者简介 : 焦明连 ( 1963- ) , 男 , 河南商丘人 , 淮海工学院测绘工程学院教授 , 主要从事空间信息和测绘教育研究。
第 18 卷 专辑 2009 年 12 月
淮海工学院学报 ( 自然科学版 ) Journa l of Hua ihai Institute of T echno logy( N atural Scien ce E d it ion )
V o.l 18 S. I . D ec . 2008
DO I : 10 . 3969/ .j jssn. 1672 6685. 2009. S0 . 001
专辑
焦明连等 : InSAR 技术在形变灾害监测中的研究 进展
3
( SBAS)方法, 研究大尺度上的形变。 SBAS 方法将所 有获得的 SAR 数据组合成若干个集合, 原则是集合 内 SAR 图像基线距小, 集合间的 SAR 图像基线距大。 对每个小集合的地表形变时间序列可以很容易利用 LS 方法得到, 但是单个集合内时间采样不够。为了 解决这个问题, 利用奇异值分解 ( SVD) 方法将多个小 [ 5] 基线集联合起来求解 。与 LS 方法和 PS 方法相比, 小基线集方法限制了长基线导致的几何去相干, 而且 使更多的 SAR 图像参与到形变计算, 增加了时间上 的采样。因此, 它比前两种方法更先进。 Berardin o 等 人利用 SBAS 方法研究了意大利南部的 Cam pi F legrei 火山口和 Nables市区在空间低分辨率下 ( 约 100 m # 100 m ) 的时间序列形变。利用了 1992∀ 2000 年之间 的 44幅 ERS 数据, 组成了 3 个小基线子集 , 70 幅差 分干涉图, 集内垂直基线距小于 130 m。研究结果表 明, 在火山口区域 , 一直存在缓慢的沉降, 累积沉降达 到 20 cm。到 2000 年, 又开始呈隆起趋势, 这样的结 果与测量部门的实测数据很好的吻合。随着 Envisat 的升空, 他们还在同一地区探讨了利用 SBAS 方法将 ERS 和 Env isat数据集合起来求解形变的方法, 结果 表明 , Cam pi F legrei火山口区域在 2000∀ 2001年有轻 微隆起, 自 2002 年开始, 形变开始趋于稳定。实验结 果与水准测量结果相符。此外, M ora 等人结合了 PS 和 SBAS 方法的特点, 在少量的 SAR 图像的基础上, 进一步提出了进行形变分析的方法, 并在实验中取得 了很好的结果。
的研究方法主要有以下三种 。 2 1 最小二乘 ( LS)方法 1999 年 , U sai研究发现, 人造目标 ( 如桥梁、 道 路、 铁路和建筑物等 ) 在长时间序列的 SAR 图像上 仍然能够保持高相干性, 由此提出了一种利用多幅 干涉 SAR 图像反演地表形变的时间序列的方法, 该 方法将求单个 SAR 图像形变的问题转化到一个最 [ 2] 小二乘的问题。从原理 可知 , LS 方法将形变过程 考虑为线 性的 , 没有 考虑相 位解 缠中 的误差 h 和 DEM 误差 , 它给出的是整体意义上的最优解。对于 大气影响 , LS 方法并没有给出检测和分离的方法, 但是它在整体上类似于一个平均的过程 , 因此 U sa i 认为 , LS 方法对大气影响相位也是有一定消除作用 的。利用 LS 方法, U sa i研究了意大利 N aples 市的 Phleg rean 地区 在 1993∀ 1999 年期 间的形 变, 她用 20 幅 ERS 1 / 2 图像 , 生成了 43 幅干涉图。研究表 明整个地区一直处于下沉状态, 沉降中心的形变达 到 22 5 cm, 平均下沉速率为 0 3 cm /month 。这与 GPS 和水准测量的结果非常一致。 2 2 永久散射体 ( PS) 方法 永久散射体 ( PS)方法是通过选取一定时期内表现 出稳定干涉行为的孤立点, 克服许多阻碍传统雷达干 涉技术的分辨率、 空间及时间上基线限制等问题。在 城市或岩石裸露好的地区, 在雷达图像上可以识别出 许多不随时间改变其特征的后向散射点 (即永久散射 点 ), 这些点可以被用作估计地面累计变形的自然标 志。与传统方法相比, PS 技术的优点在于: ( 1) 非相关散射点中孤立的相关散射点可以被 识别出来 , 在像素对像素基础上, 地形位移可以被及 时追踪; ( 2) 临界基线和多普勒质心变化不再对 PS 有限制, 所有的 ERS 图像都可以被利用 ; ( 3) 由于 长的基线 , PS 的高程可以被精确确定 ; ( 4) 只要永 2 久散射点的密度足够大 ( > 5 个点 / k m ) , 应用适当 的处理方法, 可以将大气影响识别出来并从数据中 剔除。自从这种方法被提出以来 , 它己经取得了大 [ 3] 量的应用实例。 Ferretti等人 首先利用 PS 技术的 基本模型研究了意大利 Ancona 地区的滑坡以及美 国加州 Pom ona 地区的沉降 , 并将 PS 测量的结果与
[ 1]
图 1
InSAR 成像示意图