4. 辐射传输基础：朗伯余弦定律、立体角计算、点光源与面光源模型
辐射传输，说白了就是研究光能量怎么从A点跑到B点的。你想想看，一个光源发出的光，经过一段空间，最后落到探测器上——这中间发生了什么？能量怎么分配？方向怎么变化？这就是本章要聊的核心。
我个人习惯把辐射传输拆成三块来理解：光源怎么发光（朗伯定律）、光在空间怎么分布（立体角）、光源怎么建模（点/面光源）。这三块搞清楚了，后面那些复杂的辐射度计算才有根基。
4.1 朗伯余弦定律
先讲一个我在项目中踩过的坑。有一次做均匀光源标定，我拿一个LED面板当朗伯体用，结果测出来的角度分布完全不对。后来才发现——不是所有漫射面都是朗伯体。
朗伯余弦定律说的是：一个理想漫射面，其辐射亮度在各个方向上恒定，但辐射强度随观察角度的余弦变化。公式很简单：
I(θ) = I₀ · cos(θ)其中 I₀ 是法线方向的辐射强度，θ 是观察方向与法线的夹角。
核心要点：
- 朗伯面的辐射亮度 L 与方向无关（常数）
- 辐射强度 I 随 cos(θ) 衰减
- 法线方向辐射最强，切向方向为零
为什么会这样？你可以想象一个微面元 dA，从侧面看过去，它的投影面积是 dA·cos(θ)。单位立体角内发出的功率自然就按这个比例衰减。嗯，这里要注意：朗伯定律是理想模型，真实材料只能近似满足。
工程经验：我在做积分球涂层选型时，发现硫酸钡涂层在 0°~60° 范围内近似朗伯体，超过 60° 偏差就开始明显了。如果你需要高精度测量，记得做角度修正。
4.2 立体角计算
立体角是辐射度学里绕不开的概念。说白了，它就是二维平面角在三维空间的推广。平面角用弧度衡量，立体角用球面度（sr）衡量。
定义式：
Ω = A / r²其中 A 是球面上被锥面截取的表面积，r 是球半径。整个球面的立体角是 4π sr。
4.2.1 常见立体角计算
| 几何形状 | 立体角公式 | 说明 |
|---|---|---|
| 整个球面 | 4π sr | 全空间 |
| 半球面 | 2π sr | 上半空间 |
| 圆锥半角 α | 2π(1 - cosα) sr | 锥顶在球心 |
| 矩形锥（θ₁, θ₂, φ₁, φ₂） | ∫∫ sinθ dθ dφ | 需数值积分 |
我记得有一次做LED配光曲线测试，需要计算一个窄光束的立体角。当时直接用圆锥近似，结果误差很大。后来发现——LED的光束形状往往不是标准圆锥，得用数值积分。
避坑指南：我曾经在计算小立体角时直接用 sinθ ≈ θ 近似，结果在 θ=30° 时误差已经超过 5%。小角度近似只适用于 θ < 10° 的情况。
4.2.2 立体角在辐射传输中的应用
立体角的核心作用是把辐射强度和辐射照度联系起来：
E = ∫ L · cosθ · dΩ这个积分式看着简单，实际计算时经常让人头疼。我建议你记住几个常用场景：
- 点光源对平面：E = I · cosθ / r²
- 朗伯面元对平面：E = L · π · sin²α（α 为面元半张角）
- 均匀天空对水平面：E = π · L
4.3 点光源模型
点光源是辐射度学里最基础的模型。它假设光源尺寸远小于观测距离，可以忽略光源的几何形状。你想想看，一颗星星离我们几光年远，它就是个完美的点光源。
点光源的特性：
- 辐射强度 I 是方向函数 I(θ, φ)
- 辐射照度按距离平方反比衰减：E = I / r²
- 没有空间扩展，没有自遮挡
工程判断准则：当光源最大尺寸 d 与观测距离 r 满足 d/r < 0.1 时，可以用点光源近似。误差通常在 1% 以内。
我在做汽车前照灯仿真时，一开始把LED当成点光源，结果近场照度分布完全不对。后来才意识到——近场必须用面光源模型，点光源只适用于远场。
4.4 面光源模型
当面光源的尺寸不可忽略时，就得用面光源模型了。面光源的辐射传输计算比点光源复杂得多，因为每个面元都对接收点有贡献。
4.4.1 面光源的辐射传输方程
一个面光源 dAₛ 对接收面 dAᵣ 的辐射照度贡献为：
dE = Lₛ · cosθₛ · cosθᵣ / r² · dAₛ其中：
- Lₛ 是光源面元的辐射亮度
- θₛ 是光源面元法线与连线方向的夹角
- θᵣ 是接收面元法线与连线方向的夹角
- r 是两个面元之间的距离
整个面光源对接收点的总照度就是对这个方程做面积分：
E = ∫∫ Lₛ · cosθₛ · cosθᵣ / r² · dAₛ4.4.2 常见面光源模型对比
| 模型类型 | 辐射亮度分布 | 典型应用 | 计算复杂度 |
|---|---|---|---|
| 均匀朗伯面 | L = 常数 | 积分球、标准光源 | 低（解析解） |
| 非均匀朗伯面 | L = L(x,y) | LED阵列、背光模组 | 中（数值积分） |
| 各向异性面 | L = L(θ, φ) | 反射面、散射面 | 高（蒙特卡洛） |
实用技巧：我处理大面积均匀朗伯光源时，常用“立体角投影法”快速估算。把面光源对接收点的立体角投影算出来，乘以 L 和 π，就能得到近似照度。这个方法在 d/r < 2 时精度还不错。
4.5 本章知识体系
下面这张图把本章的核心逻辑串起来了。你顺着箭头看，就能明白辐射传输的脉络：
从这张图你能看到，辐射传输基础其实就三条线：光源怎么发光（朗伯定律决定方向分布）、光怎么在空间铺开（立体角决定能量密度）、光源怎么简化（点/面模型决定计算方式）。三条线最终汇到工程应用上。
重要提醒：我见过不少工程师在仿真时直接套用点光源模型，结果近场误差大到 30% 以上。记住一个原则：先判断 d/r 比值，再选模型。d/r < 0.1 用点光源，0.1 < d/r < 1 用面光源数值积分，d/r > 1 必须用完整面光源模型。
好了，这一章的内容就这些。朗伯定律、立体角、光源模型——这三样东西是辐射传输的基石。你把它搞透了，后面那些复杂的辐射度方程、BRDF、蒙特卡洛仿真，学起来就会轻松很多。