在山区、丘陵等地形复杂区域,地形起伏会引发几何畸变和辐射畸变两类核心问题。几何畸变表现为透视收缩、叠掩和阴影,导致地物位置偏移、形状扭曲甚至信息丢失;辐射畸变则源于地形引起的局部入射角变化、照射面积差异和遮挡效应,使得相同地物在不同地形条件下呈现出不同的后向散射系数,严重影响地物分类、变化检测等定量应用的精度。本文系统阐述了SAR地形影响的产生机制,详细梳理了几何校正和辐射校正两大核心技术体系,对比分析了不同方法的适用场景与精度差异,为
SAR数据采集的高精度处理与应用提供技术参考。
一、SAR地形影响的产生机制
1. 几何畸变的产生
SAR采用距离-多普勒成像原理,目标的位置由斜距(雷达天线到目标的直线距离)和多普勒频率决定。当地形存在起伏时,同一距离向的不同高程目标会被投影到同一斜距单元,同一方位向的不同高程目标会因多普勒频移差异产生位置偏移,从而引发三种典型的几何畸变:
(1)透视收缩:当雷达波束照射到面向雷达的山坡时,山坡上不同点的斜距差小于其实际水平距离差,导致山坡在SAR图像上沿距离向被压缩。坡度越陡,透视收缩越严重;当坡度等于雷达入射角时,透视收缩达到最大值。
(2)叠掩:当面向雷达的山坡坡度大于雷达入射角时,山顶的斜距小于山脚的斜距,导致山顶在图像上先于山脚成像,形成山顶与山脚影像重叠的现象。叠掩是SAR特有的严重几何畸变,无法通过简单的几何变换消除,只能通过DEM进行校正。
(3)阴影:当雷达波束被地形遮挡时,背向雷达的山坡或山谷无法接收到雷达信号,形成无回波的阴影区域。阴影区域的长度与地形坡度和雷达入射角成正比,坡度越陡、入射角越小,阴影越长。
2. 辐射畸变的产生
地形起伏对SAR辐射特性的影响主要体现在三个方面:
(1)局部入射角变化:雷达波束的入射角是指雷达视线与地面法线之间的夹角。在平坦地区,局部入射角等于雷达入射角;但在地形起伏地区,局部入射角会随坡度和坡向的变化而显著变化。后向散射系数与局部入射角密切相关,一般来说,入射角越大,后向散射系数越小。
(2)照射面积差异:SAR图像的像素值代表单位斜距面积内的后向散射能量。当地形倾斜时,单位斜距面积对应的实际地面面积会发生变化。面向雷达的山坡,单位斜距面积对应的实际地面面积较小,像素值偏高;背向雷达的山坡,单位斜距面积对应的实际地面面积较大,像素值偏低。
(3)遮挡与阴影效应:地形遮挡不仅会导致几何上的阴影,还会使阴影区域的后向散射能量为零,同时会影响相邻区域的二次散射和多次散射,进一步加剧辐射畸变。
二、SAR地形几何校正方法
地形几何校正的目标是将SAR图像从斜距投影转换为地理坐标投影(如UTM投影),消除地形起伏引起的几何畸变,得到正射校正后的SAR图像。根据校正原理的不同,主要可分为多项式校正法、距离-多普勒模型校正法和基于DEM的正射校正法三类。
1. 多项式校正法
多项式校正法是一种基于地面控制点(GCP)的经验校正方法,其基本原理是通过建立SAR图像坐标与地理坐标之间的多项式映射关系,实现图像的几何校正。常用的多项式模型有一次、二次和三次多项式,模型阶数越高,校正精度越高,但需要的控制点数量也越多,且容易出现过拟合现象。
多项式校正法的优点是原理简单、易于实现,不需要卫星轨道参数和DEM数据;缺点是校正精度较低,且无法消除叠掩和阴影等SAR特有的几何畸变。该方法仅适用于地形平坦、精度要求不高的
SAR数据采集处理。
2. 距离-多普勒模型校正法
距离-多普勒(R-D)模型是SAR成像的基本物理模型,它精确描述了目标的斜距、多普勒频率与地理坐标之间的关系。R-D模型由三个基本方程组成:
(1)距离方程:描述目标到雷达天线的斜距与时间的关系;
(2)多普勒方程:描述目标的多普勒频率与雷达和目标相对运动的关系;
(3)地球椭球方程:描述目标在地球表面的位置关系。
通过联立求解这三个方程,可以得到目标的地理坐标。R-D模型校正法的优点是物理意义明确,校正精度较高,不需要大量的地面控制点;缺点是需要精确的卫星轨道参数和姿态参数,且在地形起伏较大的地区,由于忽略了地形高程的影响,校正精度会显著下降。
3. 基于DEM的正射校正法
基于DEM的正射校正法是目前最常用、精度最高的SAR几何校正方法。它在R-D模型的基础上,引入DEM数据提供的地形高程信息,通过迭代求解R-D方程组,得到每个DEM格网点对应的SAR图像坐标,然后通过重采样将SAR图像转换为地理坐标投影。
该方法的核心步骤包括:
(1)轨道参数精化:利用少量地面控制点对卫星轨道参数进行精化,提高轨道精度;
(2)正射变换:对于每个DEM格网点,根据R-D模型计算其在SAR图像中的坐标;
(3)重采样:采用最近邻法、双线性插值法或三次卷积法等重采样方法,将SAR图像的像素值赋值给对应的地理坐标格网;
(4)叠掩与阴影处理:根据DEM和雷达成像几何,检测并标记叠掩和阴影区域。
基于DEM的正射校正法能够有效消除地形起伏引起的所有几何畸变,包括叠掩和阴影,校正精度主要取决于DEM的分辨率和精度、卫星轨道参数的精度以及地面控制点的数量和分布。目前,该方法已被广泛应用于ENVI、SNAP、SARscape等主流遥感图像处理软件中。
三、SAR地形辐射校正方法
地形辐射校正(TRC)的目标是消除地形起伏引起的辐射畸变,恢复地物真实的后向散射系数。根据校正原理的不同,主要可分为经验统计法、半经验模型法和物理模型法三类。
1. 经验统计法
经验统计法通过分析同一地物在不同地形条件下的后向散射系数与地形参数(如坡度、坡向、局部入射角)之间的统计关系,建立校正模型。常用的经验统计模型有线性回归模型、多项式回归模型等。
经验统计法的优点是原理简单、计算量小,不需要复杂的物理模型;缺点是需要大量的训练样本,且校正模型具有较强的区域依赖性,泛化能力较差。该方法仅适用于地形相对简单、地物类型单一的区域。
2. 半经验模型法
半经验模型法结合了物理模型的理论基础和经验统计的灵活性,是目前应用最广泛的地形辐射校正方法。其中,最具代表性的是余弦校正法和SCS+G校正法。
(1)余弦校正法:余弦校正法是最简单的半经验模型,其基本假设是后向散射系数与局部入射角的余弦成正比。校正公式为:
σ⁰_corr = σ⁰_raw * (cosθ / cosθ_local)
其中,σ⁰_raw 是原始后向散射系数,σ⁰_corr 是校正后的后向散射系数,θ 是雷达入射角,θ_local 是局部入射角。
余弦校正法的优点是计算简单、效率高;缺点是忽略了地形遮挡和照射面积差异的影响,在地形起伏较大的地区校正效果较差。
(2)SCS+G 校正法:SCS校正法最初是为光学遥感数据设计的,后来被引入到 SAR 数据的地形辐射校正中。SCS+G 校正法在 SCS 模型的基础上,增加了增益因子 G,考虑了地形坡度和坡向对后向散射系数的综合影响。校正公式为:
σ⁰_corr = σ⁰_raw * (cosθ * cosα / (cosθ_local * cosα_ref)) + G
其中,α 是地形坡度,α_ref 是参考坡度,G 是增益因子。
SCS+G 校正法能够有效消除地形起伏引起的辐射畸变,校正效果优于余弦校正法,是目前应用最广泛的地形辐射校正方法之一。
3. 物理模型法
物理模型法基于电磁波与地表相互作用的物理原理,建立精确的辐射传输模型,模拟地形起伏对后向散射系数的影响。常用的物理模型有基尔霍夫模型、小扰动模型、积分方程模型等。
物理模型法的优点是物理意义明确,校正精度高,泛化能力强;缺点是模型复杂,计算量大,需要大量的地表参数(如地表粗糙度、介电常数、植被覆盖度等),且在复杂地形和植被覆盖区域,模型的精度会受到影响。该方法主要用于科学研究和高精度定量应用。
四、不同校正方法的对比与选择
不同的地形校正方法具有不同的原理、精度、计算复杂度和适用场景,在实际应用中,需要根据数据特点、研究区域地形条件和应用需求选择合适的校正方法。表1对比了主要地形校正方法的优缺点和适用场景。
表1 SAR地形校正方法对比
|
校正类型
|
方法名称
|
优点
|
缺点
|
适用场景
|
|
几何校正
|
多项式校正法
|
原理简单,易于实现,无需轨道和 DEM 数据
|
精度低,无法消除叠掩和阴影
|
地形平坦,精度要求不高
|
|
|
距离 - 多普勒模型法
|
物理意义明确,精度较高,无需大量控制点
|
忽略地形高程影响,山区精度低
|
地形起伏较小,有精确轨道参数
|
|
|
基于 DEM 的正射校正法
|
精度高,可消除所有几何畸变
|
需要高精度 DEM 和轨道参数
|
地形复杂,精度要求高
|
|
辐射校正
|
经验统计法
|
原理简单,计算量小
|
区域依赖性强,泛化能力差
|
地形简单,地物类型单一
|
|
|
余弦校正法
|
计算简单,效率高
|
忽略遮挡和面积差异影响
|
地形起伏较小
|
|
|
SCS+G 校正法
|
校正效果好,应用广泛
|
仍存在一定的经验性
|
大多数地形条件
|
|
|
物理模型法
|
精度高,泛化能力强
|
模型复杂,计算量大,参数多
|
科学研究,高精度定量应用
|
一般来说,对于地形复杂的山区,必须采用基于DEM的正射校正法进行几何校正,并结合SCS+G校正法或物理模型法进行辐射校正;对于地形平坦的平原地区,可以采用多项式校正法或距离-多普勒模型法进行几何校正,辐射校正可以简化或省略。
地形影响校正是
SAR数据采集预处理中不可或缺的关键环节,直接关系到SAR数据的解译精度和应用效果。本文系统阐述了SAR地形影响的产生机制,详细梳理了几何校正和辐射校正两大核心技术体系,对比分析了不同方法的优缺点和适用场景,并探讨了当前的研究热点与未来发展方向。
MiniSAR聚焦于微型合成孔径雷达(SAR)制造研发,为用户提供定制化机载SAR、轻型MiniSAR、无人机载MiniSAR、SAR数据采集服务、SAR飞行服务等。如您有相关业务需求,欢迎联系!