基于SAR的遥感图像几何校正研究

摘要：高分辨率遥感卫星和合成孔径雷达(Synthetic Aperture Radar，SAR)的空前进步，能够为地形图测绘提供重要的数据源。SAR作为一种主动式微波遥感，其影像的高几何分辨率对制图非常有利，可用于测绘地形图、制作正射影像图、编绘各种专题图。本文讲述了雷达遥感原理及图像的几何校正方法。

关键词：合成孔径雷达:几何校正；数字高程模型；

Abstract: With the country's economic development, there is a growing demand for the topographic maps . High-resolution remote sensing satellites and SAR have had an unprecedented progress and it becomes an important data source for the topographic mapping. SAR as an positive microwave remote sensing, high geometric resolution of its images are very beneficial to mapping. It can be used for mapping topographic maps, produced orthophoto maps, compilation of various thematic maps. This paper describes the principles of radar remote sensing and image geometric correction method.

Keyword: Synthetic Aperture Radar; Geometric Correction; Digital elevation model

中图分类号：TP7文献标识码：A 文章编号：

一、引言

高分辨率遥感卫星以及影像处理系统的相继出现使得困难地区的地形图测绘和快速更新大比例尺地形图成为可能。合成孔径雷达作为一种特殊的微波遥感器，其影像分辨率在不断提高，且利用SAR测图，具有仅用少量控制点、测图自动化程度高、工作效率高等优势。因此，深刻理解SAR成像原理，探索如何对SAR图像进行几何校正具有很强的理论意义和现实意义。

二、SAR的成像原理

SAR，是用多普勒频移这一物理现象来改善雷达成像的方位向分辨率的，它利用一个小天线作为单个辐射单元，将此单元沿直线不断移动，在移动过程中选择若干位置，在每一个位置上发射一个信号，接收相应发射位置的回波信号存储记录下来，同时保存接收信号的相位和幅度。如图1，假设一个长度为L的真实孔径雷达天线从点a移动到点b再到点c，被成像点D的雷达斜距则由大变小再变大，雷达接收到从地面点D反射回来的雷达脉冲频率也会产生变化，即频率漂移由大变小。通过精确测定接收脉冲的雷达相位延迟并跟踪频率漂移，最后可以合成一个脉冲，使方位向的目标被锐化，即提高了方位向分辨率。

图1 SAR成像几何原理

三、SAR图像的特点

雷达侧视斜距投影受到地形起伏的影响，使得SAR图像存在几何畸变，主要畸变特征有：斜距显示的近距离压缩、雷达图像透视收缩、雷达迭掩、雷达阴影、影像位移等。

(1)斜距显示的近距离压缩。在斜距显示的图像上，地面上等间距的地物目标间距离均被缩短了，但近距端(即雷达波束照射在距雷达近的一端)要比远距端缩短得更多，使图像产生几何畸变，这种现象称为图像沿斜距向的近距离压缩。

(2)SAR图像的透视收缩。SAR图像上斜坡的长度按比例尺换算后总有比实际长度短的现象，称为透视收缩。如图2所示，斜坡AB在SAR图像上的构像A1B1，显著的缩小了，而BC线段的构像BlC缩小的比较少。

(3)雷达叠掩。SAR成像时，地距大的地物目标的斜距小于或等于地距小的地物目标的斜距，在SAR图像上表现为斜距小的地物目标先于倾斜大的地物目标成像(如图3，B点和C点所成的像b和a)，或者表现为一个以上的地物目标点成像为一个像点(如图3中所示，A点和C点所成的像a)，这种图像变形称为雷达叠掩。

图2 SAR 图像的透视收缩图3 雷达叠掩

（4）雷达阴影。雷达波束在山区除了会造成透视收缩和雷达叠掩以外，还会形成阴影，即雷达阴影。在山的后坡雷达波束不能到达，因而就不可能有回波信号，在图像上的相应位置出现暗区，没有信息。当侧视角与地面坡度α之和大于90º时，在斜坡的背部形成雷达盲区，即有阴影形成。阴影的长度L与地物高度H和侧视角有关。

（5）地形起伏引起的影像位移。由于地形起伏或高大建筑物等具有相对高程，其顶部的雷达回波先于底部被天线接收，故产生影像向底点方向移位的现象，在此称为影像位移。雷达图像上地形起伏引起的像点位移与中心投影产生的像点位移相反。

四、SAR图像的几何校正方法

目前对SAR图像进行几何校正主要有基于地面控制点的校正方法和基于DEM来模拟SAR影像的校正方法。基于地面控制点的校正方法根据校正变换模型的不同又可以分为多项式校正法、共线方程校正法以及基于SAR成像原理的距离多普勒模型校正法。

1、基于地面控制点的多项式校正原理

这类校正方法的主要思路是通过在待校正图像选择地面控制点，并获取其相应的地理坐标，从而在图像空间与地理坐标空间之间建立一种变换关系模型，实现图像坐标空间向地理坐标空间的变换。

（1）多项式校正

多项式校正法的基本思想就是回避成像的空间几何过程，而直接对图像变形的本身进行数学模拟。它认为雷达图像的总体变形可以看作是平移、缩放、旋转、仿射、偏扭、弯曲以及更高次的基本变形的综合作用结果，因而校正前后图像相应点之间的坐标关系可以用一个适当的多项式来表达。

一般多项式校正变换公式可表达为:

（式1）

其中：x,y为某像素的原始图像坐标；

X,Y为同名像素的地理坐标；

（i=0,1…9）,（i=0,1…9）为多项式的待定系数。

多项式的待定系数可用最小二乘法原理求解。先根据(式1)确定所需要的最少控制点数目(N，不小于待定系数的一半)，再按照最小二乘原理求解系数。

这种方法适合于地形比较平坦的地区，但由于简单，因而利用率最高。

（2）基于地面控制点的共线方程校正原理

传感器的共线方程本身就是共线法的校正公式。1988年的第16届国际摄影测量与遥感学会上，国际摄影测量学者G.Konency利用类似的共线方程式构造SAR图像点与地面点之间的关系，称之为G.Konency公式

（式2）

其中：

为姿态参数 的方向余弦。

类似式1中对系数的求解过程，利用8个以上的控制点用最小二乘法就可以求出各个待定系数。

共线方程校正法比多项示校正法在理论上较为严密，因为它是建立在恢复实际成像条件的基础之上的。在地形起伏的影响比较大的情况下，因为在共线方程中可考虑高程的影响，所以相对多项示校正有其优越性。