4.2 辐射定标与大气校正

4.2 辐射定标与大气校正：从 DN 值到地表反射率

本节摘要：辐射定标用厂商的增益与偏置参数把 DN 值换算成辐亮度，再除以太阳辐照度得到表观反射率；大气校正进一步剥离程辐射与透过率损耗，得到可与地面实测对话的地表反射率。本节完成整套数值演算，并给出云检测、去薄雾、严格校正三级方法的适用边界。

从计数值到物理量

第 2 章的辐射传输方程在此兑现。传感器输出的 DN 值只是"带编号的计数值"——不同传感器的 DN 值没有可比性，同一传感器不同时期也可能因增益调整而系统性漂移。定标就是用随产品发布的参数把 DN 值翻译成物理量，分两步。第一步定标到辐亮度：线性换算，辐亮度等于增益乘 DN 值加偏置，参数随每期产品一起发布，传感器老化时参数会更新，用错时期的参数是"系统性偏移但不报错"的典型暗伤。第二步除以该波段的大气顶太阳辐照度（按波长加权），把辐亮度归一化为无量纲的表观反射率，它修正了太阳照度随季节的距离变化与波段差异，让不同日期的影像有了初步可比性。

# 辐射定标完整演算：DN -> 辐亮度 -> 表观反射率
import math
dn, gain, offset = 512, 0.0125, -0.2       # 假设的某波段计数值与定标参数
L = gain * dn + offset                      # 辐亮度 W/(m2·sr·um)
print(f"辐亮度 L = {L:.4f}")

esun, d, theta_s = 1550.0, 1.0163, 32.0     # 波段辐照度、日地距离(天文单位)、太阳天顶角
rho = L * math.pi * d * d / (esun * math.cos(math.radians(theta_s)))
print(f"表观反射率 rho = {rho:.4f}")
# 表观反射率已消除太阳照度与日地距离差异，跨日期粗可比；
# 但它仍含大气项：程辐射抬高暗端、透过率压低亮端
​

表观反射率与地表反射率之间隔着大气的两笔账：程辐射把暗目标抬亮（水体在蓝光波段读数明显偏高），透过率与散射把亮目标压暗。把这两笔账退掉，才是第 6 章反演真正想要的地表反射率。

校正的三级阶梯

大气校正方法按代价与精度分三级，任务日志里按需取用。零级替代：云与薄雾的保守处理。 云检测算法用亮温与可见光反射率的组合把云像元标记出来，后续分析按掩膜剔除；轻度薄雾可用"暗目标减除"快速去雾——把深净水体的读数视为程辐射本底，整景减去这个底数。这类处理不追求物理严格，但把"云污染"从数据里隔离出来，是所有业务流程的第一道闸。一级工程：暗目标法与经验线性法。 暗目标法假设影像内存在反射率近零的浓密植被或水体，其读数全部记为程辐射，逐波段扣除后再按标准大气近似估计透过率；经验线性法用同步实测的亮暗两块地面按两点法定标，工程上简单可靠。二级严格：辐射传输模型法。 输入气溶胶类型与光学厚度、水汽含量、观测几何，数值求解完整辐射传输方程，输出物理意义上严格的地表反射率；定量产品线与跨传感器联合分析的标准配置，代价是需要大气参数输入与可观算力。

# 暗目标法去程辐射：整景逐波段减除暗像元读数
import statistics
band_nir = [12, 11, 10, 38, 45, 40, 9, 12, 11]   # 某景近红外读数示例
dark = statistics.quantiles(band_nir, n=100)[1]   # 约取 2% 分位作为暗本底
print("暗本底 =", dark)
clean = [max(0, v - dark) for v in band_nir]
print("去程辐射后:", clean)
# 注意：减除假设暗目标真实反射率为零，城市区与土壤主导区并不满足
# 该法适合快速业务去雾，定量结论仍需模型法背书
​

图 4-2：从 DN 值到地表反射率的换算链

定标参数的时效与交叉验证

定标系数不是一劳永逸的数字。探测器在真空与辐射环境里缓慢老化，增益逐年漂移，运营方靠星上定标器与地面校验场定期更新参数，发布新版本系数表。工程上有两条纪律。其一，核对参数版本与成像时间的对应关系：下载产品时把定标文件与影像一起归档，注明版本号；多年期分析若跨参数版本，要么统一用反演后的表观反射率产品，要么分版本换算后再拼接。其二，做一次便宜的交叉验证：取一景已知地表反射率产品覆盖重叠区，抽样对比自己换算的数值，偏差应在几个百分点以内——这一步十分钟，能拦住"参数用错版本"这类静默事故。

常见问答

问：表观反射率与地表反射率，日常该用哪个？ 记住分界线：做比值类指数、快速摸底、单期目视分析，表观反射率够用且省事；做跨期比较、参数反演、多源拼接，必须用地表反射率。拿不准时往严里选——返工的成本永远高于多算一步。

问：产品说明里的增益偏置每年都变，用哪一版？ 用影像成像日期对应版本。定标参数随传感器老化滚动更新，版本与日期对不上就引入系统性偏移。稳妥做法是把定标文件与影像一起归档，处理日志里记下版本号——这些元数据将来就是审计依据。

问：有没有"一把尺子"检验大气校正做没做好？ 有两个便宜手段：一是过水体检验——严格校正后的深净水体在各可见光波段反射率应接近零，蓝光仍明显偏高就是程辐射没退干净；二是跨期一致性检验——同一稳定地物（大型广场、机场跑道）在多期产品上的反射率曲线应当重合，发散就是有期数据没洗干净。

参数速查：三级校正的选择表

这张表值得贴在工位上——"该用哪一级"是预处理车间每天被问的问题，答案永远从任务精度预算倒推，而不是从工具菜单里挑。

要点回顾：DN 值是计数值不是物理量，定标的线性换算把计数翻译成辐亮度，归一化成表观反射率后跨日期粗可比；程辐射抬高暗端、透过率压低亮端，退掉这两笔账才得到地表反射率；云掩膜与去薄雾是业务第一道闸，暗目标法与经验线性法是工程中档，模型法是定量标准配置；定标参数随时间漂移，版本核对与重叠区抽检是防静默事故的便宜保险。下一节完成几何账的收尾——配准、镶嵌与质量验收。

附：辐射账本速记卡

速记四条。一、三步换算：增益偏置到辐亮度、除以太阳照度成表观反射率、退掉大气两项到地表反射率——每一步解决一个可比性问题。二、参数对日期：定标版本必须与成像时间对应，归档时文件同行。三、指数可偷懒、反演不能偷：比值指数可在表观反射率上算，反演与跨期检测必须走到地表反射率。四、水体是试纸：深净水体的可见光反射率应近零，蓝光不近零就是程辐射没退净——一行检查，胜过十页报告。另记一个业务细节：云掩膜要区分云与云影，影子的反射率低、会把暗目标检验带偏，掩膜产品没分影时，太阳几何加几何关系的自动判别值得一做。

延伸：跨传感器一致性的一半之功

多传感器联合分析日益普遍，辐射一致性成了日常问题，给你一份"一半之功"清单——做完这几步，多数跨源偏差能消掉大半：一、统一到地表反射率，且大气校正用同一套模型与参数体系；二、做光谱响应匹配：两传感器同名义波段的光谱响应函数不同，用响应函数卷积或线性回归做波段等效；三、统一空间网格与聚合方式（面积权重聚合，不用最近邻粗暴重采样）；四、留出永久地物检验区做残差评估，残差超限再建线性校正。做完这四步仍剩下的差异，才轮到"训练时加传感器编码"这类模型手段兜底。顺序别颠倒——用模型硬吃本可物理消除的偏差，是当前多源项目里最常见的偷懒与翻车点。