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。其中大气程辐射是叠加项,它不携带地表信息,却是入瞳辐射的组成部分,因此定标时必须区分“入瞳辐射亮度”与“地表反射率”两个层级。
热红外波段的链路不同:地表自身发射辐射,经大气上行衰减与大气自身发射叠加后进入入瞳,定标目标通常是入瞳辐亮度或亮温。可见近红外与热红外的定标方程形式相似,但物理含义和参考源不同,这是后续绝对定标与相对定标讨论的前提。
三、定标方程:把 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 与全阵列均值的比值,得到归一化系数矩阵。
理解这两个概念的差异,关键在于区分“量纲对齐”与“内部一致”:绝对定标回答“是多少”,相对定标回答“齐不齐”。二者共同构成遥感载荷定量化应用的基础。