光学与光子学
📐 进阶推导:照度 (E) 与亮度 (L) 的数学本质与转化之谜(含 π 的由来)
文章目录
在光学设计中,我们经常遇到这样的场景: 已知探测器接收到的照度是 \(E\) (Lux),想推算出光源看起来有多亮 \(L\) (Nits)。
很多工程师直接套用公式 \(E = \pi L\)。 但你有没有想过:立体角明明是 \(2\pi\)(半球),为什么系数却是 \(\pi\)?
本文将从微分几何的角度,硬核推导两者的转化关系。
01. 两个定义的再回顾
照度 (Illuminance, \(E\))
亮度 (Luminance, \(L\))
02. 关键连接点:朗伯体 (Lambertian Surface)
要建立 \(E\) 和 \(L\) 的简单关系,必须引入一个理想模型:朗伯体(如理想的漫反射白板、纸张、墙壁)。
定义:一个表面向各个方向反射(或发射)的光的亮度 \(L\) 是常数,不随观察角度 \(\theta\) 变化。
注意:虽然亮度 \(L\) 不变,但发光强度 \(I\) 却遵循余弦定律:
\[I(\theta) = I_0 \cos\theta\](这就是为什么你侧着看纸张,觉得亮度没变,但实际上侧面发出的总能量变少了,只是因为你的投影面积也按 \(\cos\theta\) 变小了,两者抵消。)
[Image of Lambertian cosine law diagram]
03. 核心推导:\(\pi\) 是怎么来的?
假设我们有一个朗伯发光面 \(dA\),亮度为 \(L\)。我们要计算它向半球空间辐射的总光出射度(也可以理解为对外的照度贡献)\(M\)。
根据定义,光通量是亮度在立体角上的积分:
我们要对半球空间 (Hemisphere) 进行积分。 在球坐标系下,立体角微元 \(d\Omega\) 的表达式为:
代入积分公式:
因为是朗伯体,\(L\) 是常数,可以提到积分号外面:
第一步:积方位角 \(\phi\)
第二步:积天顶角 \(\theta\) 利用三角恒等式 \(\sin\theta \cos\theta = \frac{1}{2}\sin(2\theta)\):
第三步:合并结果
结论: 对于朗伯体,光出射度(或等效照度)与亮度的关系是:
这就是为什么系数是 \(\pi\) 而不是 \(2\pi\)。那个丢失的 \(1/2\) 来自于 \(\cos\theta\) (投影面积修正) 在积分过程中的衰减作用。
04. 现实世界的复杂性:BRDF
如果表面不是朗伯体(比如镜子、金属、有光泽的塑料),\(L\) 就不是常数,而是角度 \(\theta\) 的函数。
这时候,简单的 \(\pi\) 就不管用了,我们需要引入 BRDF (双向反射分布函数):
- 反射亮度 = 入射照度 \(\times\) BRDF
\[L_{out} = \int_{\Omega} f_r \cdot L_{in} \cos\theta_i d\omega_i\]
在 Zemax 或 TracePro 等光学软件中,当你定义一个表面属性时,实际上就是在定义这个 \(f_r\) 函数(是 Gaussian 散射,还是 Lambertian 散射,还是 ABg 散射模型)。
05. 总结
-
公式记忆:
- 一般情况:\(L\) 是 \(E\) 在角度上的微分。
- 朗伯体特例:\(E = \pi L\)。
-
单位换算:
- 如果你看到单位 Lambert (L),它等于 \(1/\pi \ cd/cm^2\)。
- 如果你看到单位 Foot-Lambert (fL),它等于 \(1/\pi \ cd/ft^2\)。
- 这些古老的单位里包含 \(\pi\),就是为了让 \(E\) 和 \(L\) 的数值在朗伯体下相等(为了偷懒)。但在现代 SI 单位(Lux 和 Nit)中,你必须手动乘除这个 \(\pi\)。
-
应用场景:
- 屏幕校色:色度仪测的是 \(L\) (nits),但计算屏幕总辐射功率时,如果假设屏幕是朗伯体,可以乘以 \(\pi\) 估算总光通量密度。
- 环境光传感器:手机上的光感测的是 \(E\) (lux),它无法分辨光是集中来自一个灯泡,还是来自均匀的天花板。
CONTINUE EXPLORING
DISCUSSION
留言与讨论
分享你的理解,或留下一个值得继续讨论的问题。
还没有留言,欢迎成为第一个参与讨论的人。