4. 黑体辐射理论
辐射亮度测量，说白了就是搞清楚「有多少光，从哪个方向来，落在哪个面积上」。
我刚开始接触这行时，总觉得概念绕来绕去。后来亲手搭过一次测量系统，才真正明白——测量方程是骨架，立体角是眼睛，探测器响应是手脚。三者缺一，数据就废了。
4.1 测量方程：从光子到数字
辐射亮度测量的核心方程，其实就一行：
L = Φ / (Ω · A · cosθ)其中：
- L — 辐射亮度，单位 W/(sr·m²)
- Φ — 辐射通量，单位 W
- Ω — 立体角，单位 sr
- A — 源面积，单位 m²
- θ — 法线与观测方向的夹角
但实际测量时，探测器输出的是电压或数字值。所以工程上用的测量方程长这样：
DN = G · R(λ) · L · Ω · A_det · t_int参数说明：
| 符号 | 含义 | 常见坑点 |
|---|---|---|
| DN | 探测器输出数字值 | 注意暗电流扣除 |
| G | 系统增益 | 不同增益档位要分别标定 |
| R(λ) | 光谱响应度 | 波长相关，不能当常数 |
| L | 目标辐射亮度 | 我们要求的值 |
| Ω | 立体角 | 由光学系统决定 |
| A_det | 探测器像元面积 | datasheet 上写的是名义值 |
| t_int | 积分时间 | 线性区间内使用 |
关键认知：测量方程不是简单的公式代入。每个参数都有频率响应、温度漂移、非线性误差。我见过有人直接把 datasheet 上的 R(λ) 拿来用，结果反演出的 L 差了 30%。
4.2 立体角计算：别让几何坑了你
立体角是辐射亮度测量里最容易出错的地方。你想想看，一个锥形光束，到底覆盖了多少「角度面积」？
立体角的定义式：
Ω = A / r²其中 A 是球面上被锥体截出的面积，r 是球半径。单位是球面度（sr）。
实际工程中，我们常用小角度近似：
Ω ≈ π · (D/2f)²D 是入瞳直径，f 是焦距。这个近似在视场角小于 10° 时误差可忽略。
但要注意——一旦视场角超过 20°，必须用精确公式：
Ω = 2π · (1 - cos(θ/2))θ 是全锥角。
我的经验：做地物光谱测量时，镜头视场角通常 5°~15°，用近似公式没问题。但做室内 BRDF 测量，有时用到 30° 视场角，必须用精确公式。我曾经偷懒没换，结果反演出的辐射亮度偏高了 8%。
立体角计算的常见场景：
- 点光源测量 — 立体角由探测器孔径和距离决定
- 扩展源测量 — 立体角由光学系统视场决定
- 光纤耦合 — 立体角由数值孔径 NA 决定：Ω = π · NA²
4.3 探测器响应：从光子到电信号
探测器响应，就是把光子流变成可测量的电信号。这个过程不是线性的，也不是瞬时的。
响应度公式：
R(λ) = I_out / Φ_in单位是 A/W 或 V/W。但实际探测器有三个关键特性：
- 光谱响应 — 不同波长响应不同。硅探测器在 900nm 附近响应最高，到 1100nm 就几乎为零。
- 线性动态范围 — 信号太弱信噪比差，太强会饱和。我建议工作在线性区间的中间 60%。
- 响应时间 — 积分时间不能太短，否则信号建立不完全；也不能太长，否则暗电流累积过多。
避坑指南：我曾经用 InGaAs 探测器测 1550nm 激光，没注意温度对响应度的影响。室温 25°C 和 40°C 时，响应度差了 12%。后来我养成了习惯——每次测量前先做暗电流和参考源标定。
4.4 知识结构图：辐射亮度测量原理
下面这张图，把测量方程、立体角、探测器响应串在了一起。我建议你保存下来，做实验前看一眼。
4.5 实战中的测量流程
理论讲完了，咱们看看实际怎么做。我一般按这五步走：
- 确定几何关系 — 测量距离、视场角、立体角。这一步错了后面全白干。
- 标定探测器 — 用标准辐射源做响应度标定。我习惯用 NIST 可追溯的卤钨灯。
- 设置积分时间 — 先预估信号强度，让 DN 值落在满量程的 30%~80%。
- 采集暗信号 — 盖上镜头盖，采集暗电流。这一步不能省。
- 正式测量 — 采集目标信号，扣除暗信号，代入测量方程反演 L。
核心公式汇总：
辐射亮度：L = Φ / (Ω · A · cosθ)
探测器输出：DN = G · R(λ) · L · Ω · A_det · t_int
立体角（小角度）：Ω ≈ π · (D/2f)²
立体角（精确）：Ω = 2π · (1 - cos(θ/2))
嗯，到这里辐射亮度测量的核心原理就讲完了。记住一句话：测量方程是骨架，立体角是眼睛，探测器响应是手脚。下次你拿到一组辐射亮度数据，先问自己三个问题——几何关系对吗？探测器标定了吗？暗电流扣了吗？