2. 遥感载荷辐射传输链路解析
辐射定标的核心任务，是把传感器记录的数字量化值（DN）还原为具有物理意义的入瞳辐亮度或表观反射率。要做到这一点，必须先理清从太阳到传感器的完整辐射传输链路，明确每一个环节引入的物理量、作用机制与可测量参数。本节沿着“太阳—大气—地表—大气—传感器”这条路径，逐段建立辐射定标所需的物理量链条与测量方程。
一、链路全景：能量从哪里来，到哪里去
对于工作在可见光—短波红外（VIS-SWIR）波段的被动光学载荷，太阳是唯一的初始辐射源。太阳以约5778 K的等效黑体温度向外辐射，在日地平均距离处，大气层顶（TOA）垂直于太阳光方向的太阳辐照度记为 E₀(λ)，即太阳常数光谱分布。太阳辐射以天顶角 θₛ 斜射进入大气，穿过大气时受到分子散射（瑞利散射）、气溶胶散射（米氏散射）和气体吸收（水汽、臭氧、氧气、二氧化碳等）的共同作用，一部分被衰减，一部分被散射到其他方向。
到达地表的太阳辐射分为两部分：一是直接透射的太阳直射光，二是被大气散射后从天空各方向到达地表的漫射光。两者共同构成地表总下行辐照度。地表按自身反射特性（方向反射率或双向反射分布函数 BRDF）将入射能量反射回上行方向，反射能量再次穿过大气，经历上行路径的散射与吸收，最终到达传感器入瞳。此外，大气自身的散射（程辐射）也会直接进入传感器视场，形成路径辐射项。因此，传感器入瞳处的辐亮度是“地表反射贡献”与“大气程辐射贡献”的叠加。
二、逐段物理量链条与测量方程
把上述链路写成可计算的方程，是辐射定标建模的关键。在朗伯体地表假设下，传感器入瞳辐亮度 L(λ) 可表示为：
L(λ) = [ E₀(λ) · cosθₛ · T↓(λ) · ρ(λ) · T↑(λ) ] / π + Lp(λ)式中各项含义如下：E₀(λ) 为大气层顶太阳光谱辐照度，单位 W/(m²·μm)；θₛ 为太阳天顶角；T↓(λ) 为太阳—地表下行路径的大气透过率；ρ(λ) 为地表半球反射率（朗伯假设下等于方向反射率）；T↑(λ) 为地表—传感器上行路径的大气透过率；Lp(λ) 为大气程辐射，即未到达地表而直接散射进入传感器视场的辐亮度。除以 π 是因为朗伯面反射辐亮度与反射辐照度之间满足 L = ρE/π 的余弦关系。
若把地表反射贡献与程辐射合并，可写成更紧凑的形式：
L(λ) = L₀(λ) + [ E₀(λ)·cosθₛ·T↓(λ)·ρ(λ)·T↑(λ) ] / π其中 L₀(λ) 即程辐射项。这一形式在辐射定标中非常实用：L₀(λ) 与地表无关，只取决于大气状态和观测几何，可通过暗像元法或大气辐射传输模型（如 6S、MODTRAN）计算；而含 ρ 的项则与地表耦合。定标时，若已知 L(λ) 与 ρ(λ)，即可反推 T↓、T↑、Lp 等大气参数；反之，若大气参数已知，则可由 L(λ) 反演地表反射率。
三、从入瞳辐亮度到 DN：传感器端的转换
入瞳辐亮度 L(λ) 进入传感器后，经光学系统聚焦、分光、探测器光电转换、放大与量化，最终输出 DN 值。在线性响应区间内，这一过程可表示为：
DN = G · L(λ) + DN₀其中 G 为绝对辐射定标增益（DN/(W·m⁻²·sr⁻¹·μm⁻¹)），DN₀ 为暗电流与偏置对应的零点输出。定标的任务就是通过实验室积分球或星上定标器，确定 G 与 DN₀，使 DN 能够被还原为 L(λ)。把传感器方程与链路方程联立，就得到完整的“太阳—大气—地表—传感器”测量方程：
DN = G · { [ E₀(λ)·cosθₛ·T↓(λ)·ρ(λ)·T↑(λ) ] / π + Lp(λ) } + DN₀这个方程是遥感定量反演与在轨定标评估的共同基础。它清楚地表明：DN 不仅取决于载荷自身的增益和偏置，还受到太阳角度、大气透过率、程辐射和地表反射率的共同影响。因此，在轨定标时必须把载荷响应与大气、地表贡献分离，否则会把大气变化误判为载荷衰减。
四、各环节的测量与获取方式
要让上述方程可解，需要为每个物理量提供实际获取途径。E₀(λ) 通常采用标准太阳光谱（如 ASTM E490 或 WRC 谱），按日地距离进行修正；θₛ 由观测时间与经纬度计算。T↓、T↑、Lp 依赖大气状态，可通过同步的大气参数测量（太阳光度计、微波辐射计）结合辐射传输模型计算，也可用暗像元、对比度法等图像方法估算。ρ(λ) 则通过地面同步实测（地物光谱仪）或已知反射率靶标获得。
在星上定标场景中，常用太阳漫反射板作为已知反射率目标：太阳以固定角度照射漫反射板，板面反射光直接进入传感器，此时 ρ 已知、路径大气影响可忽略或可建模，从而直接标定 G 与 DN₀。这种“太阳—漫反射板—传感器”的短链路，正是对上述长链路方程的简化与验证。
五、链路解析对定标实践的意义
理解这条链路，能帮助我们在定标中做出正确判断。例如，当星上定标器显示载荷响应下降时，需要区分是探测器真实衰减，还是漫反射板污染导致反射率变化；当同一地物在不同日期 DN 差异明显时，应检查太阳角度、气溶胶光学厚度和水汽含量是否发生变化。只有把链路中每个环节的贡献量化清楚，才能把“载荷性能变化”从“观测条件变化”中剥离出来，这正是辐射定标与性能评估的核心逻辑。
实际工作中，可先固定观测几何与大气条件，用辐射传输模型生成一组 L(λ) 与 ρ(λ) 的对应关系，再与传感器 DN 做线性回归，得到 G 与 DN₀ 的初值；随后用独立样本验证残差是否随 ρ 或波长呈现系统性偏差，以此判断大气校正或定标系数是否需要迭代修正。这一过程把链路方程从理论公式转化为可执行、可检验的定标流程。