D - InSAR原理及其数据处理流程-

1、第38卷增刊2008年11月吉林大学学报(地球科学版Journal of J ilin University (Earth Science Editi on Vol 138Sup 1Nov 12008D 2I nSAR 原理及其数据处理流程吕乐婷,陈圣波,车大为,刘银萍吉林大学地球探测科学与技术学院,长春130026摘要:合成孔径雷达是近20年发展起来的一种空间对地观测技术,目前已被广泛应用于获取地表信息及地表形变。介绍了D 2I nS AR 基本原理,发展以及在其地学方面的应用。并对三轨D 2I nS AR 数据处理的算法及流程作了简要介绍。关键词:D 2I nS AR;地表形变;相位解缠中

2、图分类号:TP75文献标识码:A作者简介:吕乐婷(1984,女,辽宁大连人,硕士研究生,主要从事遥感与地理信息系统研究,E 2mail:ayier1220 。The Pr i n c i ple and Da t a Processi n g Flow of I n terferom etrySyn theti c Aperture RadarL Le 2ting ,CHEN Sheng 2bo,CHE Da 2wei,L I U Yin 2p ingCollege of GeoExploration Science and Technology,J ilin U niversity,Chan

3、gchun 130026,ChinaAbstract:Synthetic aperture radar interfer ometry is a ne w devel oped and potential technol ogy for observati on of earth fr om s pace,and now has been widely used in obtaining the surface infor mati on of the earth .I n this paper,the devel op ing and the app licati on of D 2I nS

4、 AR is intr oduced,and als o the algorithm and data p r ocessing fl ow .Key words:D 2I nS AR;surface def or mati on;phase un wrapp ing0引言合成孔径雷达S AR (synthetic aperture radar 是近20年发展起来的一种空间对地观测技术,目前已被广泛应用于获取地表起伏的信息。与传统光学遥感相比,S AR 具有全天时全天候,穿透力强,分辨率高等诸多优点。用S AR 图像通过干涉测量技术I n 2S AR (interfer ometry synt

5、hetic aperture radar 和差分干涉测量技术D 2I nS AR (differential interfer ometry synthetic aperture radar 可以进一步得到地面高程信息和地表形变信息。在国外,合成孔径雷达雷达干涉测量已经从早期的论证研究阶段发展到了比较成熟的应用研究阶段,积累了大量经验和研究成果,并且已取得了厘米级的应用精度,而我国在这方面的应用正在起步。虽然在具体应用中还有许多关键技术和限制有待突破与克服,S AR 技术的绝对优势在垂直变形监测中已毋庸置疑。1合成孔径雷达的发展与地学应用111I nSAR 的发展合成孔径雷达S AR 是20世

6、纪50年代末研制成功的一种微波传感器,早期的雷达遥感大多基于单张S AR 图像的灰度信息来进行地质调查、极地冰川、土地利用、植被和生态环境监测等。70年代,射电天文领域发展成熟的干涉技术被引入,将覆盖同一地区的两张S AR 图像联合处理并提取对应像素的相位差(干涉相位信息,以此恢复地物信息,从而诞生了I nS AR 技术,即合成孔径雷达干涉测量。差分干涉测量技术D 2I nS AR 是在雷达干涉测量的基础上发展起来的,它利用遥感卫星多时相的复雷达图像相干信息进行地表的垂直形变量提取,其精度已达到毫米级,且具有前所未有的连续空间覆盖、高度自动化和高精度监测地表变形的能力,为地表变形的自动化监测提

7、供了全新的方法。112I nSAR 及D 2I nSAR 的地学应用1978年美国国家航空航天局发射的海洋卫星SE AS AT 上首次装载了合成孔径雷达,由此,星载S AR 技术由实验室成功的进入应用领域。发展至今,I nS AR 测量技术已成为获取地形、地表形变和地物信息的最具潜力的新技术,他的应用领域也得以不断扩展,包括:数字高程模型DE M 的获取;地图测绘;生成城市建筑三维透视图以及地球动力学应用。D 2I nS AR 技术通过消除I nS AR 获取的相位图中地形相位部分,可以获得由于地表运动而造成的地表目标在雷达视线方向上位移所对应的相位贡献。因此广泛用于:(1地震现象的研究:包括

8、同震、震后和抗震构造的机理研究;(2火山的下陷与抬升研究:通过对火山的运动规律分析,进行火山爆发的预测研究;(3冰川研究:通过I nS AR 技术获取完整的,高分辨率的,高精度的地形数据,以及测量冰流和其他变化;(4细微地形变化:包括滑坡,地面沉降等地表形变。2D 2I nSAR 基本原理及数据处理流程211D 2I nSAR 基本原理假定卫星两次对同一地区成像,两副天线接收所信号的路径分别为和+r ,这时生成干涉图纹两次成像的相位差为=4,因此我们可在基线距B 、两次观测的信号路径长度、视角、基线与水平方向夹角sin (-几个参数之间建立几何关系,并忽略(r 2,得:r =B sin (-=

9、B p ,=4B p 。(1其中B p 代表基线距在辅图像视线方向上分量的大小。若在同一地区保持主图像相同(即和保持不变,选择该地区另一幅S AR 影像生成新的干涉纹图,则可以将先后生成的两副干涉纹图进行比较和计算。应有=4B p,所以,两幅干涉纹图之间相位的比值应该等于各基线距辅图像视线分量的比值,而与地形起伏无关,即=B pB r。如果在第3次观测期间,地表发生了形变,并且这种形变与雷达的分辨单元相比很小,则可认为这三幅S AR 影像的雷达信号仍然相关。这时其干涉纹图所提供的相位信息除了地形信息之外,还包含有雷达视线方向的地表形变分量。那么第2幅干涉纹图的相位就应为=4(B p +,进一步

10、可推出雷达视线方向上的地形变量表达式:-B p B p=4。(2式中左边各量可以根据卫星的轨道参数和干涉相位推算得到,借此进一步获得S AR 图像上各个位置处雷达视线方向上的地形变量。212D 2I nSAR 数据处理流程可以看出,合成孔径雷达干涉测量的几何原理并不复杂,但是由于对于同一场景获得两幅相干性高的S AR 图像本身具有复杂性,并且从复图像无法直接获得干涉相位的真实值,使得I nS AR 测量成为一个相当复杂的过程。三轨法差分干涉测量是利用三景图像生成两幅干涉纹图,一幅反映地形信息,一幅反映地表形变信息,进行平地效应消除后,分别进行相位解缠,最后利用差分干涉原理计算得到地表形变信息。

11、三轨法的主要优点是无需地面信息,数据间的配准较易实现。缺点是相位解缠的好坏将影响最终结果。其基本步骤为:(1两幅干涉图像对生成干涉图1,进行平地消除,干涉图l 中包含地形信息和形变信息。(2主图像与第三幅图像形成干涉图像对,生成干涉图2,进行平地效应消除,干涉图2仅包含地形信息。(3对去平地效应的干涉图1和2进行相位解缠。(4利用公式(2进行形变计算。(5将形变图投影至地理坐标系中。3结语I nS AR 是一个交叉性很强的新领域,涉及对地观测,电磁波传播,信号处理,影像处理与模式识别,空间大地测量和数字摄影测量等多个领域。用D 2I nS AR 技术监测地表变形的优点更是不言而喻的,但要实际用

12、来取代水准测量监测微小的地面垂直变形还需要克服许多技术难点,其中最重要的是克服相干性的限制和突破相位解缠的障碍。经过了20多年的研究,随着多颗遥感卫星的成功发射并投入使用,I nS AR 技术本身及其应用都得到了快速发展,应用领域不断扩大。随着电子技术的不断改进以及旨在提高图像清晰度的信号处理算法的不断成熟,高分辨率的S AR 系统将必成为今后发展的趋势。332增刊吕乐婷,等:D 2I nS AR 原理及其数据处理流程图1星载合成孔径雷达差分干涉测量数据处理流程图F i g11the flow chart of D2I nSARs da t a processi n g参考文献:1M ass

13、onnet D,RossiM,Car mona C,et al.The dis p lace2ment field of the Landers earthquake mapped by rader in2terfer ometryJ.Nature,1993,364:1382142.2Zebker H A,Rosen P A,Goldstein R M,et al.On thederivati on of coseis m ic dis p lace ment fields using differen2tial radar interfer ometry:the landers earthquakeJ.Geo23王超,张红,刘智.星载合成孔径雷达干涉测量M.北京:科学出版社,2002.4廖明生,林晖.雷达干涉测量原理与信号处理M.北京:测绘