📚 课程目录第4章光学遥感物理基础

4. 光学遥感物理基础

光学遥感载荷设计的第一步，不是选探测器，而是把“目标—大气—太阳—载荷”这条能量链路算清楚。载荷最终能收到多少能量，决定了口径、焦距、积分时间、探测器灵敏度和信噪比指标的分配。本节围绕太阳反射式被动成像这一主流体制，梳理从太阳辐照到入瞳辐亮度的完整物理链条，并给出可落地的建模方法。

一、辐射传输的基本量与链路

描述这条链路需要几个基本辐射量。辐亮度 L（W·m⁻²·sr⁻¹·μm⁻¹）是单位投影面积、单位立体角、单位波段的辐射功率，它是贯穿整个链路的核心量，因为它在无损耗自由空间中沿程守恒。辐照度 E（W·m⁻²·μm⁻¹）是到达某表面的辐射通量密度。反射率 ρ 描述目标对入射能量的反射能力，分为方向—半球反射率与二向反射分布函数 BRDF。

太阳反射链路可概括为：太阳以约 1361 W/m² 的太阳常数（大气层外，日地平均距离）照射地球，经大气上行衰减后到达地表；地表目标按自身反射特性将部分能量散射到上半球；目标反射能量再经大气上行衰减进入载荷入瞳；入瞳辐亮度经光学系统成像到焦面探测器。整条链路中，大气在“下行”和“上行”两段各衰减一次，同时叠加自身的程辐射。

二、太阳反射与大气窗口

太阳辐射近似 5800 K 黑体，其能量集中在 0.3～3 μm，峰值约在 0.5 μm，这正是可见光—近红外遥感的工作区间。太阳常数在日地平均距离处约 1361 W/m²，实际到达地表的辐照度还要乘以日地距离修正因子和太阳天顶角余弦。

大气对太阳辐射的作用包括吸收、散射和发射。吸收主要由水汽、二氧化碳、臭氧、氧气等分子造成，形成一系列强吸收带；散射包括瑞利散射（分子，与 λ⁻⁴ 成正比，蓝光散射强）和米氏散射（气溶胶，与波长关系较弱）。这些作用把光谱切成若干“窗口”——透过率高、可有效观测地面的波段。工程上常用的窗口有：可见光—近红外（0.4～0.9 μm）、短波红外（1.0～1.75 μm、2.0～2.4 μm）、中波红外（3～5 μm）和长波红外（8～14 μm）。设计波段时必须避开 1.4 μm、1.9 μm 等强水汽吸收带，否则信号会被大气“吃掉”。

三、目标反射特性

目标反射特性用反射率光谱和 BRDF 描述。不同地物的光谱曲线差异是遥感识别的物理基础：植被在可见光因叶绿素吸收而反射率低，在 0.7 μm 附近出现“红边”陡升，近红外反射率高；水体在近红外几乎全吸收，反射率极低；干土、雪、云则各有特征。BRDF 描述反射的方向性，分为镜面反射、朗伯漫反射和方向性反射。多数自然地物近似朗伯体，可用半球反射率简化；但水面、金属屋顶、冰面等存在明显镜面或热点效应，若仍按朗伯处理会带来显著误差。

建模时通常把目标反射率写成 ρ(λ, θ_s, θ_v, φ)，其中 θ_s 为太阳天顶角，θ_v 为观测天顶角，φ 为相对方位角。对朗伯假设，反射辐亮度为 L_t = ρ·E_d·cosθ_s/π，这个 1/π 来自朗伯面把入射辐照度转换为辐亮度的积分关系，是初学者最容易漏掉的因子。

四、入瞳辐亮度建模与示例

把上述环节串起来，入瞳辐亮度可写为：

L_TOA(λ) = [ ρ(λ)·E_d(λ)·cosθ_s / π ] · τ↑(λ) + L_p(λ)

其中：
  E_d(λ) = E_sun(λ)·cosθ_s·τ↓(λ) + E_sky(λ)   （直射 + 天空漫射）
  E_sun(λ) 为大气层外太阳光谱辐照度
  τ↓、τ↑ 分别为下行、上行大气透过率
  L_p(λ) 为大气程辐射（路径辐射）

下面给出一个可见光波段的示例计算（数值为教学假设，非实测）。设 0.66 μm 处大气层外太阳光谱辐照度 E_sun = 1550 W·m⁻²·μm⁻¹，太阳天顶角 θ_s = 30°（cosθ_s ≈ 0.866），下行透过率 τ↓ = 0.80，天空漫射 E_sky = 80 W·m⁻²·μm⁻¹，目标反射率 ρ = 0.25，上行透过率 τ↑ = 0.85，程辐射 L_p = 8 W·m⁻²·sr⁻¹·μm⁻¹。

先算到达地表的直射辐照度：1550 × 0.866 × 0.80 ≈ 1074 W·m⁻²·μm⁻¹；加上天空漫射 80，得 E_d ≈ 1154 W·m⁻²·μm⁻¹。目标反射辐亮度 L_t = 0.25 × 1154 × 0.866 / π ≈ 79.5 W·m⁻²·sr⁻¹·μm⁻¹。经上行衰减并叠加程辐射：L_TOA = 79.5 × 0.85 + 8 ≈ 75.6 W·m⁻²·sr⁻¹·μm⁻¹。可以看到，程辐射占了约 10%，在低反射率目标（如水体）观测中占比会更高，必须扣除或校正。

建模时先固定参照约定：太阳天顶角 θ_s 从当地天顶量起，观测天顶角 θ_v 从目标天顶量起，方位角 φ 以太阳方位为基准。全文与图示保持一致，避免正负号与角度定义混乱。

五、工程落地要点

波段选择：优先落在高透过率窗口内，避开 1.4 μm、1.9 μm 水汽吸收带和 0.76 μm 氧气吸收带。
反射率输入：用实测或标准地物光谱库（如植被、水体、土壤、雪）作为 ρ(λ) 输入，不要用单一常数代替全波段。
大气参数：τ↓、τ↑、L_p 随气溶胶光学厚度、水汽含量、观测几何变化，可用 MODTRAN 或 6S 按实际大气模式计算。
程辐射校正：低反射率目标场景中 L_p 占比大，需在反演中显式扣除，否则反演反射率偏高。
与载荷指标衔接：把 L_TOA 作为入瞳辐亮度输入，结合光学系统 F 数、透过率、探测器量子效率与积分时间，推算到达焦面的电子数，进而评估 SNR 是否满足任务要求。

把这条链路算通，后续的口径、焦距、积分时间、探测器选型才有物理依据。入瞳辐亮度建模不是一次性的公式套用，而是随波段、观测几何、大气条件和目标类型反复迭代的过程，也是载荷指标分配中最需要与用户任务反复确认的环节。