insar技术提取dem数据的应用

《insar技术提取dem数据的应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、InSAR技术提取DEM数据的应用ok3etricSyntheticApertureRadar,InSAR)技术是极具潜力的微波遥感新技术。本文介绍了InSAR的基本原理、数据处理流程,通过对卫星提供的真实SAR数据进行干涉处理。结果表明,利用InSAR技术可以有效获取大范围的地面数字高程模型。【关键字】InSAR数据处理流程干涉图像数字高程模型1引言合成孔径雷达干涉测量(InterferometricSyntheticApertureRadar,InSAR)技术是极具潜力的微波遥感新技术,它是合成孔径雷达遥感成像与电磁波干涉两大技术融合的结果。其基本思路是通过两副天线同时观

2、测(单轨双天线模式)或两次平行观测(重复轨道模式),获得同一地区的重复观测数据(复数影像对),提取同一目标对应的两个回波信号之间的相位差,结合观测飞行平台的轨道数据提取高程信息,可以获取高精度、高分辨率的地面高程信息。在InSAR基础上进一步发展了差分干涉技术(DifferentialInterferometricSyntheticAertureRadar,D_InSAR),利用D_InSAR技术可以精密测定地表的微小形变,精度可达厘米级甚至毫米级[1]。InSAR技术在地形测绘、城市目标探测和城市形态分析、海洋表面及其流速场探测、极地冰川监测、农业和资源调查等领域有着极其广

3、阔的应用前景。而D_InSAR技术在地震、地面沉降、火山活动、冰川漂移、山体滑坡等大范围地面形变的监测和研究方面表现出极大的潜力[2]。2InSAR基本原理InSAR技术主要利用SAR影像的后向散射强度进行两幅SAR图像的相关性计算和图像配准,而利用其相位信息,根据雷达波长和周期以及雷达天线位置推算出地面散射单元的位置,即利用SAR干涉影像的相位成分来算出地面散射单元的位置从而形成有关地球表面的地形图。如图1所示,设A1和A2是卫星两天线位置,天线之间的距离即基线距用B表示,其平行和垂直于视线向分量分别为和。基线与水平方向的夹角为α,H表示天线飞行高度,地面一点O到天线A1,

4、A2的路径分别为R1,R2(R2=R1+ΔR),θ是雷达波入射角,O点高程为z。设S1是在位置A1被拍摄到的相位为Φ1的影像,S2是在位置A2被拍摄到的相位为Φ2的影像。根据合成孔径雷达成像原理有:图1InSAR成像示意图在由天线A1,A2和地面点O组成的三角形中,利用余弦定理,经计算整理有:由于在星载系统中,R>>B,所以有ΔR≈,则(2)式可近似的写成:于是,干涉测量中的雷达波入射角计算公式应该为:在合成孔径雷达干涉测量处理中,雷达天线飞行高度H、基线B和基线倾角α都可以从雷达系统参数辅以卫星星历等信息获得,ΔΦ和斜距R可以从干涉图上求得,通过公式(5)可求得

5、精确的θ值,由公式(6)可计算得到目标的精确高程。其精度在理论上可达到波长的量级(cm级)。以上给出的解释模型是针对两部天线分别独立的发射与接收回波的系统推导而来的,对于一部天线发射、两部天线同时接收回波的系统,只要把公式中的4π代之以2π就可以了。3InSAR数据处理过程经过多年的研究,InSAR数据的主要处理过程已基本形成[3],主要包括以下几个部分:SAR数据的读取、基线估计、图像配准和重采样、生成干涉图、相位解缠、地理编码。从预处理后的两幅复数图像出发,整个过程大致如图2所示:[1][2][3]下一页ok3建模前,均需要计算基线的长度和倾角,所以一般先进行基线估计。基

6、线值估计结果的准确性受轨道数据精度影响。目前基线估计的方法有轨道参数几何法、基于地面控制点(GCP)精确估算法和干涉条纹图频谱估算法等基线估计方法。3.2复图像配准和重采样SAR图像配准通常有两类方法:一、基于精密卫星星历和相干系数的配准;二、基于地面控制点的配准。方法一不需要外部信息,直接通过卫星的基本参数就可以达到精度在1/10像元的配准要求,但是这需要两幅图像有很好的相关性,卫星的轨道信息很精确等较强的约束条件;方法二需要在图像上选取控制点,根据控制点的坐标来配准两幅影像,通常这些控制点是SAR图像上的一些容易辨认的特征点或者布设的角反射器。配准精度通常限制在误差小于0

7、.2个像元。复图像配准之后需要对从图像进行重采样,使从图像上的每一个像元值精确的对应于主图像上的每一个像元。重采样的方法通常有最近临法,线性内插法,(4点、6点)立方体卷积插值法,(6点、8点或16点)有限长sinc函数法。3.3干涉图和相关图的生成对于地面上任意一点在两幅已配准的SAR图像上的像元值可用复数表示,两复数共轭相乘即可得到该像元干涉后的复数表示形式。对所有像元进行共轭相乘即可生成干涉图。由于SAR图像在距离向和方位向上的分辨率不同,在干涉后一般需将方位向像元做多视处理,使得干涉图符合实际的