地形照明校正要区分坡面受光差异与真实地表反射差异
说明坡度、坡向、太阳几何和地形阴影如何改变山地反射率,并比较经验与半经验校正的适用边界。
地形照明校正要区分坡面受光差异与真实地表反射差异
山地影像中,同一种地物会因坡度、坡向和太阳位置不同而呈现明显亮暗差异。朝向太阳的坡面接收更多直射辐射,背阴坡面则可能主要接收天空散射和邻近地形反射。地形照明校正试图减少这种几何效应,但不能把阴影中从未获得的信号完整恢复。
入射角描述局部受光
局部太阳入射角由太阳天顶角、太阳方位角、坡度和坡向共同决定。平地只需考虑太阳天顶角,山地则要把地表法向加入计算。坡度和坡向通常从 DEM 邻域估计,因此 DEM 分辨率、噪声和配准误差会直接进入校正。
影像与 DEM 必须在共同网格和坐标下对齐。山脊处一个像元的偏移就可能把错误坡向分配给高反差像元。先把低分辨率 DEM 任意上采样到影像网格,不会增加真实地形细节,反而可能产生过度精确的表面外观。
局部入射角只能描述直射几何,不包含大气散射、地形遮挡和邻坡反射。把所有辐射都按余弦关系缩放,常会使背阴坡过度增亮并放大噪声。
常见校正思路
余弦校正假设地表近似朗伯体,用太阳天顶角与局部入射角的余弦比调整反射率。它形式简单,但在低入射余弦和阴影区不稳定,也忽略地物方向性。
半经验方法通过影像样本拟合反射率与入射几何关系,增加常数或按地表类型调整强度,以减轻过校正。拟合系数依赖波段、地物、季节和太阳条件,不能从一个场景直接复制到另一个场景。
分层拟合可以按森林、裸地或其他覆盖分别估计,因为冠层结构和土壤的方向响应不同。但分类错误会把不合适的参数应用到像元。若分类本身由未校正影像产生,还要避免循环依赖。
阴影与不可恢复信息
地形自阴影发生在坡面背向太阳时,投射阴影则由远处地形遮挡形成。简单局部入射角不能完整计算投射阴影,需要视线或地形遮挡分析。云影是另一来源,应由云质量信息和时间条件处理。
深阴影中直射信号缺失,传感器记录可能接近噪声。把这些像元强行拉伸到阳坡亮度,会制造不可靠光谱。定量分析可把深阴影单独标记、限制校正系数或从特定任务排除。
雪、明亮岩石和水体也可能违反通用拟合关系。水面不应根据陆地坡面进行地形校正;积雪的高反射和方向性需要独立处理。掩膜规则应在拟合前执行。
与大气和 BRDF 的关系
大气校正处理辐射在大气中的吸收和散射,地形校正处理局部照明几何,两者不是同一步骤。先后顺序和算法假设需要一致。有些大气模型使用地形和邻近效应,重复校正可能过度补偿。
BRDF 描述地物反射随入射和观测方向变化。地形改变局部入射与观测几何,因此朗伯假设只是近似。多角度或高精度应用可能需要更完整的方向反射模型,但参数需求也更高。
地形校正目标也取决于任务。为了视觉镶嵌,可以追求坡向亮度均衡;为了反射率反演,应保持物理定义和不确定度。视觉上更平坦的影像不一定物理上更准确。
验证校正是否改善
常见诊断是比较同类地物校正前后反射率与入射角的相关性。相关降低说明地形依赖减弱,但并不能单独证明真实差异被保留。还应检查均值、方差、光谱关系和极值,防止所有坡面被过度压平。
可选择覆盖不同坡度、坡向和照明条件的稳定地物,使用独立场景或地面观测验证。分类任务则比较统一样本和空间划分下的结果,不能只看训练精度。时间序列还要检查校正是否减少季节性太阳角度造成的伪趋势。
输出应保留 DEM 来源、分辨率、坡度算法、太阳角度、校正公式、拟合样本、系数、阴影掩膜和限制条件。不同场景使用不同系数时,时间序列必须记录每期参数。
地形照明校正的合理目标是减少可解释的几何亮暗,同时承认阴影、方向反射和 DEM 误差留下的不确定性,而不是让山地影像看起来像平地。
处理顺序需要固定
山地流程中,正射校正、云雪筛选、大气校正、地形参数计算和照明校正彼此依赖。每次运行应固定处理顺序,并确认太阳角度对应场景获取时刻而不是处理时刻。镶嵌多景时,各景太阳几何不同,先统一拟合可能把景间差异混入地形系数;应先按场景处理和验证,再决定怎样拼接。
若产品已经执行过地形相关处理,二次校正前必须查阅处理说明并用未校正输入对照,避免重复补偿。
资料依据
- NASA ARSET:Fundamentals of Remote Sensing,https://appliedsciences.nasa.gov/sites/default/files/Fundamentals_RS_Session2_Land_Final.pdf
- Gonzalez、Woods:Digital Image Processing, Fourth Edition,https://www.pearson.com/en-us/subject-catalog/p/digital-image-processing/P200000003224/9780133356779