利用DEM进行SAR图像模拟及地形辐射校正

薛东剑，李增元，郑 洁，李成绕，李婉秋

(1. 中国林业科学研究院资源信息研究所，北京100091； 2. 成都理工大学地球科学学院，四川 成都 610059； 3. 地学空间信息技术国土资源部重点实验室，四川 成都 610059)



利用DEM进行SAR图像模拟及地形辐射校正

薛东剑1,2,3，李增元1，郑 洁2,3，李成绕2,3，李婉秋2,3

(1. 中国林业科学研究院资源信息研究所，北京100091； 2. 成都理工大学地球科学学院，四川 成都 610059； 3. 地学空间信息技术国土资源部重点实验室，四川 成都 610059)

在山区获取地面控制点比较困难，运用模拟SAR进行配准、校正，具有较大的优势。本文在分析SAR成像几何结构及ALOS PALSAR卫星轨道参数特征的基础上，运用RD定位模型对DEM的每个网格点进行雷达成像点的位置计算，模拟SAR图像,并提取当地入射角、投影角及规则化因子等；模拟出的SAR图像与真实SAR图像纹理吻合，有利于控制点的自动配准。在此基础上对ALOS PALSAR进行编码，构建基于规则化因子及入射角的地形辐射校正模型，消除面积效应及地形起伏造成的畸变问题，从结果中分析，校正后的图像明暗差异明显减少，这对雷达定量反演研究具有一定的现实意义。

PALSAR；SAR模拟；数字高程；地形辐射校正

合成孔径雷达(SAR)能够全天候、全天时工作，对地表具有一定的穿透力，相比于光学图像具有较大的优势与应用前景，但SAR特有的成像几何结构对其处理具有一定的难度。SAR数据处理的好坏直接关系到SAR的应用，而校正是SAR应用的前提。目前SAR图像校正主要分为3种：椭球表面校正地理编码(GET)、地形校正地理编码(GTC)及地形辐射校正。GET是没有进行地形校正的，GTC是采用DEM进行了地形校正的。由于SAR图像会存在透视收缩、叠掩、阴影等几何与辐射畸变，在运用SAR进行分类及定量参数反演前(如土壤湿度、森林生物量)需进行地形辐射校正，以消除由于地形起伏造成的畸变问题。1989年Freeman A[1]提出了基于当地入射角的地形辐射校正方法，随后国内外大量学者进行了地形辐射校正方法研究[2-5]。但是对于叠掩较为严重的区域，这些方法还一直处于改进与发展中[6-7]。本文在分析总结以往地形辐射校正、SAR图像模拟的基础上，采用一种基于规则化因子、入射角等相结合的方法，对ALOS PALSAR数据进行地形辐射校正，对运用SAR图像进行定量参数反演等具有一定的现实意义。

1 SAR系统参数及几何定位模型

1.1 ALOS PALSAR系统参数介绍

ALOS卫星是日本的对地观测卫星，轨道高度为691.65 km，周期为98.7 min，重复周期46 d，带有3个传感器，其中PALSAR有3种观测模式：精细、扫描、极化。精细又分为精细波束单极化模式(FBS)、精细波束双极化模式(FBD)。数据选取2010年9月11日的ALOS PALSAR数据(见表1)，地理坐标为：103°35′47″E—104°22′40″E，31°23′48″N—32°00′13N，位于四川安县与茂县的接合部位；区内地形复杂多样，主要为平坝、丘陵(台地)、中低山3种类型。即以大光包斜冲断层和北川冲断层为界，西北部属四川西部地槽区的前龙门山褶断带，系龙门山脉，地势较高，高川乡境内的大光包海拔3047 m，为区内最高点，东南部地势相对较低，最低点位于界牌镇石安村与绵阳市交界的安昌河，海拔为490 m。

表1 ALOS PALSAR影像主要参数

1.2 SAR图像几何定位模型

SAR作为一种主动遥感成像方式，提供了传感器到目标的距离和多普勒历史信息。运用这些信息可以精确地将卫星和地面坐标相联系，从而构建SAR定位模型，通过解算定位模型就可以得到每个像元的地理坐标。R-D模型是由Brown首先提出，Curlander等发展了该模型，并给出了作为分析问题出发点的3个基本方程。R-D定位原理是利用等距离线、等多普勒线在地球等高面上的交点确定像元位置，如图1所示。

图1 SAR成像几何结构

(1)

R2=(X-Xs)2+(Y-YS)2+(Z-ZS)2

(2)

(3)

2 SAR图像模拟

在对SAR图像进行辐射校正前需先进行几何校正，模型参数迭代优化过程需要利用控制点。对于山区，人工寻找控制点需要消耗大量的人力物力，有时很难实现。采用模拟影像与真实SAR影像可以实现快速、自动的高精度配准。目前模拟SAR图像主要分3类：相干系统模拟、非相干的图像模拟和基于图像的模拟[8-11]。在此采用DEM建立地面坐标系与斜距坐标系的映射关系，进行SAR模拟。首先根据PALSAR成像参数,采用RD定位模型,计算DEM每个像元的SAR影像坐标(i,j),并将结果保存到查询列表(LUTA)；根据模拟出的SAR影像，计算每个DEM像元的当地入射角(如图2所示)和投影角；根据纹理及几何特征对模拟SAR影像和待校正SAR进行自动配准，建立多项式方程，计算模拟SAR影像对应的真实SAR影像坐标，并更新、精化查询列表(LUTA),求算每个DEM像元对应的真实SAR影像坐标(i,j),将其像元值重采样到DEM坐标空间,得到与DEM大小相同的模拟SAR影像(如图3所示)。模拟的图像纹理特征与实际的SAR图像非常逼近，反映了成像所独有的透视收缩、叠掩、阴影等基本特征。由于模拟影像的后向散射强度是通过计算像元和局部入射角的经验函数相乘得到的，在模拟时没考虑地表后向散射的独立性，对于平原区大部分为常数值，而研究区大部分为山区，并不影响整个图像的自动配准。利用模拟影像与真实影像的匹配可以弥补轨道数据的不准确，从而实现精确的几何校正。此时校正只是根据地形和定位模型将SLC影像的像元值经重采样正射到其地面实际位置,与原始图像相比发生了明显的变化，尤其是叠掩区域，明显拉伸，并没有对像元值进行地形辐射校正。

图2 当地入射角

图3 DEM模拟的SAR

3 地形辐射校正

在地形复杂区，雷达成像的辐射特性受地形的影响较大，几何校正后还需对SAR图像进行地形辐射校正，校正的好坏直接影响定量化反演的精度。图1中n为地面单元的法线，m为成像面的法线，θ为当地入射角(n与入射向量的夹角)，φ为投影角(向量n与m之间的夹角)，δ为参考入射角，则

(4)

(5)

投影角φ又可表示为

(6)

(7)

(8)

式中,RAO和n可以通过已知的DEM求算,但要求算对应目标矢量RAO的SAR影像坐标(i,j)及卫星位置矢量RSA,需要通过SAR卫星成像参数及建立R-D定位模型来确定。设地距图像中距离向与方位向间隔为：ρg、ρa则在雷达视线投影向上的面积A为

(9)

式中，a为最大坡度角；a1为距离向坡度角；a2为方位向坡度角。推导后其关系为：tan2a=tan2a1+tan2a2。为了消除面积效应，对灰度图像乘以规则化散射因子，其因子为

(10)

结合规则化因子及入射角，建立校正模型

(11)

式中，β0为成像面散射单元对应的后向散射系数。

为了评价校正的效果一般常采用影像校正前后方差减少的百分比来衡量，即

(12)

校正前影像方差为0.15,校正后为0.12,减少了20%，校正效果较明显。从图4分析，校正后朝向雷达入射方向与背向入射方向的明暗差异明显减少、纹理特征更加均匀，地形引起的后向散射单位面积变化对散射系数的影响得到了明显改善。另外经统计分析，校正前散射系数与当地入射角具有明显的线性相关性，校正后相关系数显著下降，由此验证了本文运用的方法是可行的。

图4 地形辐射校正前后SAR图像

4 结 语

由于雷达特有的成像方式受地形影响较大，对SAR图像进行地形辐射校正具有重要的意义，尤其对定量反演研究。文中以ALOS PALSAR数据为例，在对SAR成像几何结构分析的基础上，重点研究了基于DEM的SAR图像模拟的地形辐射校正。模拟后的SAR图像在山区纹理特征与实际的SAR图像非常逼近，有利于运用自动控制点进行配准。通过模拟的SAR图像可以精化查询列表，结合RD定位模型对原始SAR图像作了GTC校正，将SLC影像的像元值经重采样正射到其地面实际位置。结合当地入射角及投影角，在计算规则化散射因子的基础上进行了地形辐射校正，经过校正，朝向雷达入射方向与背向入射方向的明暗差异明显减少、纹理特征更加均匀，散射系数与当地入射角的线性相关性显著降低。

[1] FREEMAN A,CURLANDER J.Radiometric Correction and Calibration of SAR Images[J].Photogrammetric Engineering and Remote Sensing,1989,55:1295-1301.

[2] LOEW A,MAUSER W.Generation of Geometrically and Radiometrically Terrain Corrected SAR Image Products[J].Remote Sensing of Environment,2007,106(3):337-349.

[3] 张永红,张继贤,林宗坚,等.地形引起的雷达辐射畸变及其校正[J].测绘科学,2002,27(4):23-26.

[4] 陈尔学,李增元,田昕,等.星载SAR 地形辐射校正模型及其效果评价[J].武汉大学学报(信息科学版),2010,35(3):322-327.

[5] SMALL D.Flattening Gamma:Radiometric Terrain Correction for SAR Imagery[J].IEEE Transaction on Geoscience and Remote Sensing,2011,49(8):3081-3093.

[6] SHIMADA M.Ortho-rectification and Slope Correction of SAR Data Using DEM and Its Accuracy Evaluation[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2010,3(4):657-671.

[7] 刘文祥,张继贤,黄国满,等.地形辐射校正在极化SAR影像分类中的应用[J].测绘科学,2014,39(3): 98-102.

[8] 吴涛,王超,张红.雷达图像模拟研究综述[J].遥感信息,2005(4):237-240.

[9] FRANCESCHETTI G，MIGLIACCIO M，RICCIO D． SARAS: A Synthetic Aperture Radar (SAR) Raw Signal Simulator[J].IEEE Transactions on Geoscience and Remote Sensing，1992，30(1):110-123.

[10] 刘智,张永贺,王冬红,等.基于DEM 的SAR 图像模拟[J].测绘科学技术学报, 2008,25(5): 317-320.

[11] 秦志远,张宝印,靳国旺,等.基于DEM 和地物分类信息的SAR 图像模拟[J].系统仿真学报,2009，21(8):2437-2440.

[12] XUE Dongjian,ZHENG Jie,LI Chengrao,et al.Landslides Monitoring with SAR/InSAR Technolog[C]∥The Third Imaging Radar Senior Academic Seminar of Earth Observation.[S.l.]:[s.n.],2016:542-545.

[13] 王庆,曾琪明,焦健,等. 基于DEM 和查找表的高分辨率机载SAR 图像正射校正[J].测绘通报,2013(11):10-13.

[14] SMALL D F.Lattening Gamma:Radiometric Terrain Correction for SAR Imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2011,49(8):3081-3093.

Simulating SAR Images and Correcting Terrain Radiation Using DEM

XUE Dongjian1,2,3,LI Zengyuan1,ZHENG Jie2,3,LI Chengrao2,3,LI Wanqiu2,3

(1. Research Institute of Forest Resources Information Techniques,Chinese Academy of Forestry,Beijing 100091,China； 2. College of Earth Sciences ,Chengdu University of Technology,Chengdu 610059,China； 3. Key Laboratory of Geoscience Spatial Information Technology,Ministry of Land and Resources of the P.R.China,Chengdu 610059，China)

The simulated SAR image is beneficial to the automatic registration of the control points and image correction,especially in the mountains, where it is very difficult to get the ground control point. Based on the analysis of SAR imaging geometry and characteristics of ALOS PALSAR satellite orbit parameters, the range doppler (RD)geo-location model is used to calculate the position of the radar imaging point for each grid point of the DEM. The SAR image is simulated, the local incidence angle, projection angle，regularization factor are extracted. The simulated SAR image coincides with the real SAR image texture, which is beneficial to the automatic registration of the control points. And the terrain radiation correction model based on the regularization factor and angle of incidence is constructed to eliminate the distortion caused by the area effect and terrain undulation. From the analysis results, the difference of the brightness and darkness of the image after correction is obviously reduced. It has certain practical significance to radar quantitative inversion research.

PALSAR; SAR simulation; digital elevation model; terrain radiation correction

薛东剑，李增元，郑洁，等.利用DEM进行SAR图像模拟及地形辐射校正[J].测绘通报，2017(7)：14-17.

10.13474/j.cnki.11-2246.2017.0215.

2016-12-13

四川省教育厅重点项目(16ZA0100)；国家973计划(2013CB733404)；国土资源部地学空间信息技术重点实验室开放基金(KLGSIT2013-06)

薛东剑(1977—)，男，博士生，副教授，主要研究方向为SAR图像处理及干涉测量。E-mail:xdj101@sina.com

P23

A

0494-0911(2017)07-0014-04