4. 目标辐射特性与背景环境建模

光电载荷成像链路的第一性输入,是进入入瞳的辐射亮度。要判断一幅图像里目标是否“亮得起来”、对比度是否够用,必须先把目标、地表背景和大气三者各自的辐射贡献算清楚,再按几何关系叠加成场景辐射亮度。本节围绕这条主线,给出可落地的建模方法。

一、辐射传输的基本分解

在可见—短波红外(0.4~2.5 μm)与中波—长波红外(3~5 μm、8~12 μm)两个区间,入瞳辐射亮度的构成不同,但都可写成“路径辐射 + 背景/目标反射或自身发射经大气衰减”的形式。对反射波段,传感器入瞳处光谱辐射亮度可写为:

L(λ) = L_path(λ) + τ(λ)·[ ρ(λ)·E_g(λ)/π + L_adj(λ) ]

其中 L_path 为大气程辐射,τ 为大气透过率,ρ 为地表或目标反射率,E_g 为到达地面的总辐照度(直射 + 漫射),L_adj 为邻近像元贡献。热红外波段则把 ρ·E_g/π 替换为 ε·B(λ,T),即目标自身发射项,ε 为发射率,B 为普朗克函数。理解这一分解,是后续所有参数取值的依据。

二、地表反射特性建模

地表反射不是单一常数,而是随波长、入射角、观测角变化的函数,即双向反射分布函数 BRDF。工程上常用简化模型:朗伯体假设 ρ 与角度无关,适合植被、粗糙土壤的粗略估算;对水体、金属屋面、结冰路面等镜向性强的目标,必须引入方向反射,否则会严重低估或高估辐亮度。

实际建模时,先按地物类型给出典型光谱反射率曲线。例如:健康植被在 0.55 μm 附近有绿峰(ρ≈0.10~0.15),在 0.68 μm 有叶绿素吸收谷(ρ≈0.05),在 0.75~1.1 μm 近红外平台迅速抬升(ρ≈0.4~0.6);干沙在可见—短波红外缓慢上升,ρ 从 0.3 增至 0.5;沥青路面整体偏低,ρ≈0.05~0.12。这些数值属于教学示例量级,实际任务应查光谱库或现场实测。

若目标尺寸小于像元,还需考虑面积加权:像元反射率等于目标与背景反射率按投影面积占比的线性混合,这是亚像元目标对比度下降的直接原因。

三、大气背景建模

大气对成像的影响有三项:衰减、程辐射、背景辐射。衰减由气溶胶和分子散射、气体吸收共同决定,用透过率 τ 描述;程辐射 L_path 是大气自身散射进入视线的辐射,在短波蓝光段最显著,是图像“发灰”的主因;热红外段还要叠加大气上行、下行辐射。

工程上常用 MODTRAN、6S 等工具,输入大气模式(如中纬度夏季)、气溶胶类型(乡村/城市/海洋)、能见度、观测几何(太阳天顶角、方位角、传感器高度与天底角),输出各波段 τ 与 L_path。以 0.55 μm、能见度 23 km、太阳天顶角 30° 为例,τ 约 0.8~0.9,L_path 明显抬升暗目标亮度,使对比度下降。若不做大气校正,直接用表观辐亮度反演反射率,误差可达数十个百分点。

四、目标光谱特征提取

目标光谱特征建模的关键,是区分“本征特征”与“场景特征”。本征特征来自材料发射率/反射率,如高温尾焰在 3~5 μm 呈强辐射峰,伪装网在近红外反射率被压低;场景特征则叠加了大气与背景。建模步骤为:先获取目标材料光谱,再按观测几何与大气条件传播到入瞳,最后与背景辐亮度相减得到目标—背景对比度。

对比度常用两种度量:辐射对比度 C=(L_t−L_b)/(L_t+L_b),以及信噪比意义上的对比度噪声比 CNR=(L_t−L_b)/NEDL,NEDL 为载荷噪声等效辐亮度。前者反映场景本身,后者决定能否被探测。

五、场景辐射亮度构建流程

把上述要素组织成可执行的构建流程,是本节的核心交付。下面用一张关系图表达各输入如何汇聚为场景辐亮度立方体。

场景辐射亮度构建关系图 地表反射率 ρ(λ) BRDF / 光谱库 / 实测 大气参数 τ, L_path MODTRAN / 6S 输出 目标光谱 ε/ρ(λ) 材料本征 + 温度 观测几何 太阳/传感器角度 辐射传输计算 反射波段:L = L_path + τ·(ρ·E_g/π + L_adj) | 热红外:L = L_path + τ·ε·B(λ,T) 按波段逐 λ 计算,输出入瞳光谱辐亮度 目标辐亮度 L_t(λ) 含大气衰减与程辐射 背景辐亮度 L_b(λ) 地表 + 大气背景 对比度评估 C = (L_t − L_b)/(L_t + L_b) | CNR = (L_t − L_b)/NEDL 关键误差源 忽略 L_path / 朗伯假设 亚像元混合未加权 输出成果 场景辐亮度立方体 (x, y, λ) 逐像元

六、一个可复算的教学示例

设某反射波段中心波长 0.66 μm,太阳天顶角 30°,地表为沥青路面 ρ=0.08,植被背景 ρ=0.35,大气透过率 τ=0.85,程辐射 L_path=8 W/(m²·sr·μm),到达地面总辐照度 E_g=600 W/m²。按朗伯近似:

目标:L_t = 8 + 0.85×(0.08×600/π) ≈ 8 + 13.0 = 21.0
背景:L_b = 8 + 0.85×(0.35×600/π) ≈ 8 + 56.8 = 64.8
对比度:C = (21.0 − 64.8)/(21.0 + 64.8) ≈ −0.51

负号表示目标比背景暗,这是沥青路面的典型情形。若忽略程辐射,C 会变成 −0.63,明显高估对比度;若把沥青误当朗伯体而实际存在镜向反射,结果还会进一步偏离。可见程辐射与 BRDF 两项是误差的主要来源。

七、建模检查要点

  • 波段划分是否与载荷光谱通道一致,避免用宽波段平均值代替逐 λ 计算。
  • 大气模式、气溶胶类型、能见度是否与任务季节、地域匹配。
  • 目标与背景反射率是否来自同一光谱库或同一批次实测,避免量纲与基准不一致。
  • 亚像元目标是否按面积占比做线性混合,混合后再算对比度。
  • 热红外目标是否同时给出温度与发射率,二者缺一不可。
  • 最终是否输出逐像元、逐波段的辐亮度立方体,而非单一平均值。
建模时先固定观测几何与大气条件,再替换目标/背景反射率,可快速得到同一场景下不同地物的对比度排序,用于载荷波段选择与曝光设计。

完成上述步骤后,场景辐亮度立方体即可作为成像链路仿真的输入,用于评估入瞳辐亮度动态范围、量化位数需求以及目标可探测性。下一环节的载荷响应建模,将直接建立在这一辐亮度基准之上。