2. 遥感载荷辐射定标基础

辐射定标是把遥感载荷输出的原始数字量(DN)与真实物理辐射量建立定量映射的过程。没有这层映射,影像上的灰度值只是无量纲的计数,无法用于反演地表反射率、亮温或辐射收支。要理解定标,先要理清三个层次:用什么物理量描述辐射、辐射从目标到载荷经历了什么、以及用什么数学关系把两端连起来。

一、辐射量定义:从辐射源到探测器的物理量链

描述辐射场最基础的量是辐射亮度 L,单位为 W·m⁻²·sr⁻¹·μm⁻¹,它表示单位面积、单位立体角、单位波段宽度内通过的辐射功率,是遥感中最常用的量,因为它与观测方向无关地刻画了视场内某点的辐射强度。与之相关的还有辐照度 E(W·m⁻²·μm⁻¹),表示入射到单位面积上的辐射功率,常用于描述太阳光到达地表或入瞳处的能量;两者关系为 E = ∫L·cosθ·dΩ,即辐照度是辐射亮度在半球立体角上的积分。

当辐射进入载荷后,探测器把光信号转换为电信号,再经量化得到数字量 DN。DN 是无量纲整数,其位数由量化器决定,例如 12 bit 对应 0~4095,14 bit 对应 0~16383。DN 与入瞳辐射亮度之间通常近似线性,但存在暗电流偏置和增益差异,这正是定标要标定的对象。此外,光谱辐射亮度强调单位波长间隔,而波段积分辐射亮度则考虑载荷光谱响应函数 RSR(λ) 的加权,实际定标系数往往对应后者。

二、辐射传输链路:从太阳到 DN 的完整路径

以太阳反射波段为例,辐射从太阳出发,经大气衰减后到达地表,被地物反射后再次穿过大气进入载荷入瞳,最终被探测器接收。这条链路可拆解为:太阳辐照度 → 大气上行透过率 → 地表反射(含方向性) → 大气上行透过率与程辐射 → 入瞳辐射亮度 → 光学系统透过率 → 探测器响应 → 量化输出 DN。其中大气程辐射是叠加项,它不携带地表信息,却是入瞳辐射的组成部分,因此定标时必须区分“入瞳辐射亮度”与“地表反射率”两个层级。

热红外波段的链路不同:地表自身发射辐射,经大气上行衰减与大气自身发射叠加后进入入瞳,定标目标通常是入瞳辐亮度或亮温。可见近红外与热红外的定标方程形式相似,但物理含义和参考源不同,这是后续绝对定标与相对定标讨论的前提。

太阳 E_sun(λ) 大气层顶 TOA 大气层:吸收 / 散射 / 程辐射 地表(反射率 ρ 或发射率 ε) 太阳下行辐照度 地表反射辐射 大气程辐射 L_path 入瞳辐射亮度 L_TOA(λ) 光学系统 τ_opt 探测器响应 R(λ) 数字量 DN 12/14 bit 量化 定标方程 L = Gain·DN + Offset Gain:辐射 响应斜率 Offset:暗电流 偏置 L 为入瞳辐射 亮度 辐射传输链路:太阳 → 大气下行 → 地表反射/发射 → 大气上行+程辐射 → 入瞳 → 光学系统 → 探测器 → DN

三、定标方程:把 DN 翻译成辐射亮度

在探测器线性响应区间内,入瞳辐射亮度 L 与 DN 的关系可写成:

L(λ) = Gain(λ) · DN + Offset(λ)

其中 Gain 的单位是 W·m⁻²·sr⁻¹·μm⁻¹/DN,Offset 的单位与 L 相同。若载荷响应存在非线性,则需用二次多项式或分段线性拟合:L = a·DN² + b·DN + c。定标的任务就是通过已知辐射源确定 Gain 与 Offset(或 a、b、c)。

实际工程中,定标方程常以波段积分形式给出,即把光谱响应函数 RSR(λ) 纳入:

L_band = ∫ L(λ)·RSR(λ) dλ / ∫ RSR(λ) dλ
DN = (L_band - Offset) / Gain

这意味着同一 DN 值在不同波段对应的辐射亮度不同,定标系数必须逐波段给出。以某 12 bit 可见近红外载荷为例,假设实验室积分球标定得到某波段 Gain = 0.025 W·m⁻²·sr⁻¹·μm⁻¹/DN,Offset = -3.2 W·m⁻²·sr⁻¹·μm⁻¹,当 DN = 2000 时,L = 0.025×2000 - 3.2 = 46.8 W·m⁻²·sr⁻¹·μm⁻¹。该数值属于示例假设,用于说明计算过程。

四、绝对辐射定标与相对辐射定标

绝对辐射定标解决的是“DN 对应多少物理辐射量”的问题,目标是建立 DN 与入瞳辐射亮度(或亮温)之间的绝对映射,使不同时间、不同载荷的数据可以统一到同一物理量纲下比较。它依赖已知辐射基准,如实验室积分球、标准黑体、太阳漫反射板或场地同步实测。绝对定标系数通常包含 Gain 和 Offset,精度用不确定度衡量,典型要求为 3%~5%(反射波段)或 0.5~1 K(热红外)。

相对辐射定标解决的是“同一载荷不同探元之间响应是否一致”的问题,目标是消除探测器阵列各像元之间的响应差异,使均匀辐射场景下输出均匀的 DN。它不追求与绝对物理量对齐,而是追求内部一致性。相对定标通常用均匀光源(积分球、大面积均匀场地)获取各探元的响应均值,再计算归一化系数,使所有探元向同一参考水平看齐。

对比维度 绝对辐射定标 相对辐射定标
目标 DN → 入瞳辐射亮度/亮温 消除探元间响应差异
参考源 积分球、黑体、太阳板、场地实测 均匀光源、均匀场景
输出 Gain、Offset(逐波段) 归一化系数(逐探元)
精度指标 绝对不确定度 3%~5% 探元间非均匀性 <1%
典型应用 定量反演、跨载荷比较 影像去条纹、拼接匀色
相对定标是绝对定标的前置步骤:先让所有探元“步调一致”,再给整机赋予绝对物理量纲。若相对定标未做好,绝对定标系数会掩盖探元间差异,导致影像出现条纹。

五、从概念到操作:一次定标的数据流

以实验室积分球绝对定标为例,操作流程为:设置积分球输出 5~8 个辐射亮度等级,每个等级采集 100 帧影像取平均以抑制噪声;对每帧扣除本底(快门或暗场);按探元提取 DN 均值;用最小二乘拟合 L-DN 直线,得到 Gain 和 Offset;最后用独立等级验证残差,若残差超过 2% 则检查积分球稳定性或探测器非线性。相对定标则在同一组均匀数据上,计算各探元 DN 与全阵列均值的比值,得到归一化系数矩阵。

理解这两个概念的差异,关键在于区分“量纲对齐”与“内部一致”:绝对定标回答“是多少”,相对定标回答“齐不齐”。二者共同构成遥感载荷定量化应用的基础。