光谱辐射度、光度与色度计算原理
光谱辐射度、光度与色度计算原理
一、引言
在光学测量与照明工程领域，对"光"的量化描述存在三个相互关联却又各具独立意义的物理量体系：辐射度学（Radiometry）、光度学（Photometry） 与色度学（Colorimetry）。
辐射度学以纯粹的物理能量为尺度，描述电磁辐射在空间中传输、分布与接收的规律，其量值不依赖于任何观察者的感知特性。光度学则在辐射度学的基础上引入人眼视觉函数 $V(\lambda)$，将物理辐射能量加权转化为人眼主观感受的"明暗"度量，实现了从客观能量到主观感知的桥梁跨越。色度学更进一步，通过 CIE 标准色度观察者模型，将光谱功率分布映射为三维颜色空间中的坐标，使"颜色"这一主观体验得以被精确量化和传递。
三者的核心纽带是光谱功率分布 $\Phi(\lambda)$——它是辐射度量向光度量、色度量转换的共同起点。本文将沿着"辐射度量 → 光度量 → 色度量"的逻辑路径，系统阐述各物理量的定义、计算方法和内在联系，并结合实际光谱测量场景说明其工程应用价值。
二、辐射度学基础
辐射度学（Radiometry）研究波长范围涵盖紫外、可见与红外波段，其核心物理量均采用下标 e（energetic）标识，以区别于光度学中对应的下标 v（visual）。
2.1 辐射通量（Radiant Flux）
辐射通量 $\Phi_e$ 定义为单位时间内发射、传输或接收的辐射能量，单位为瓦特（W）。对于光谱测量而言，辐射通量通常表达为波长的函数——即光谱辐射通量 $\Phi_e(\lambda)$，单位为 $\mathrm{W/nm}$ 或 $\mathrm{W/}\mu\mathrm{m}$。
在有限波长区间内的积分辐射通量为：
$$ \Phi_e = \int_{\lambda_1}^{\lambda_2} \Phi_e(\lambda) \,\mathrm{d}\lambda $$
在实际光谱测量中，积分退化为离散求和：
$$ \Phi_e = \sum \Phi_e(\lambda) \,\Delta\lambda $$
2.2 辐射强度（Radiant Intensity）
辐射强度 $I_e$ 定义为点辐射源在单位立体角内发射的辐射通量，单位为 $\mathrm{W/sr}$（瓦特每球面度）：
$$ I_e = \frac{\mathrm{d}\Phi_e}{\mathrm{d}\Omega} $$
辐射强度描述了辐射源在空间中某一方向上的功率集中程度，是表征 LED、激光器等定向光源辐射特性的重要参数。
2.3 辐亮度（Radiance）
辐亮度 $L_e$ 是辐射度学中最核心的物理量之一，定义为面辐射源在单位投影面积、单位立体角内的辐射通量，单位为 $\mathrm{W/(sr \cdot m^2)}$：
$$ L_e = \frac{\mathrm{d}^2\Phi_e}{\mathrm{d}A \cdot \cos\theta \cdot \mathrm{d}\Omega} $$
其中 $\theta$ 为观察方向与表面法线的夹角，$\mathrm{d}A \cdot \cos\theta$ 为有效投影面积。
对于光谱测量中常见的面光源，辐亮度可由光谱辐射通量计算得到：
$$ L_e = \frac{1}{\pi S} \sum \Phi_e(\lambda) \,\Delta\lambda $$
其中 $S$ 为发光面积，系数 $1/\pi$ 来源于朗伯辐射体在半球空间内的积分关系。辐亮度的一个重要特性是：对于理想的朗伯辐射体，$L_e$ 与观察方向无关（即各方向亮度相等），这使得辐亮度成为表征显示屏、扩散面光源均匀性的首选参数。
2.4 辐照度（Irradiance）
辐照度 $E_e$ 定义为入射到单位接收面积上的辐射通量，单位为 $\mathrm{W/m^2}$：
$$ E_e = \frac{\mathrm{d}\Phi_e}{\mathrm{d}A} $$
辐照度表征的是"被照面"接收辐射能量的密度，与辐射强度（描述"源"）和辐亮度（描述"源的面亮度"）形成互补。在太阳辐射测量、紫外消毒剂量评估等场景中，辐照度是最直接相关的辐射度量。
2.5 辐射度学物理量小结
| 物理量 | 符号 | 定义 | 单位 |
|---|---|---|---|
| 辐射通量 | $\Phi_e$ | 单位时间的辐射能量 | $\mathrm{W}$ |
| 辐射强度 | $I_e$ | 单位立体角的辐射通量 | $\mathrm{W/sr}$ |
| 辐亮度 | $L_e$ | 单位投影面积·单位立体角的辐射通量 | $\mathrm{W/(sr \cdot m^2)}$ |
| 辐照度 | $E_e$ | 单位面积的入射辐射通量 | $\mathrm{W/m^2}$ |
三、光度学基础
光度学（Photometry）是辐射度学在可见光波段（约 380–780 nm）的"人眼加权"版本。它将物理辐射能量按照人眼对不同波长光的敏感程度进行加权积分，从而得到与人眼主观明暗感受相一致的量值。光度学物理量采用下标 v（visual）标识。
3.1 人眼视觉函数 $V(\lambda)$
CIE（国际照明委员会）于 1924 年定义了明视觉光谱光视效率函数 $V(\lambda)$（photopic luminous efficiency function），描述人眼在亮度较高（明视觉，主要是视锥细胞工作）条件下对不同波长单色光的相对敏感度。
$V(\lambda)$ 是一条钟形曲线，峰值位于 555 nm（黄绿色），峰值为 1。在偏离 555 nm 的波长处，$V(\lambda) < 1$，表示人眼对这些波长的光感受较"暗"。例如：
- 在 510 nm（青色），$V(\lambda) \approx 0.503$，即相同辐射功率下，510 nm 的光主观亮度约为 555 nm 的一半；
- 在 650 nm（红色），$V(\lambda) \approx 0.107$，即辐射功率需近 10 倍于 555 nm 才能产生相同亮暗感。
在暗视觉（低亮度，视杆细胞工作）条件下，则使用暗视觉光谱光视效率函数 $V'(\lambda)$，其峰值移至 507 nm。本文以下讨论均基于明视觉条件。
在实际工程计算中，常使用 $K(\lambda)$（或记为 $K_m V(\lambda)$）作为光谱光视效能函数，其峰值 $K_m = 683 \,\mathrm{lm/W}$。
3.2 光通量（Luminous Flux）
光通量 $\Phi_v$ 是辐射通量 $\Phi_e(\lambda)$ 经 $V(\lambda)$ 加权积分后的结果（乘上 $K_m$ 折算为流明），单位为流明（$\mathrm{lm}$）：
$$ \Phi_v = K_m \int_{380}^{780} \Phi_e(\lambda) \, V(\lambda) \,\mathrm{d}\lambda $$
离散形式为：
$$ \Phi_v = 683 \times \sum \Phi_e(\lambda) \, V(\lambda) \,\Delta\lambda $$
光通量是光度学中最基本的导出量，描述光源在可见光波段"对人眼产生的总明亮感觉"。
3.3 $K_m = 683 \,\mathrm{lm/W}$ 的由来
$K_m = 683 \,\mathrm{lm/W}$ 是明视觉条件下最大光谱光视效能（maximum spectral luminous efficacy）。它的定义可以追溯到 1979 年第 16 届 CGPM（国际计量大会）对坎德拉（candela）的重新定义：
坎德拉是发出频率为 $540 \times 10^{12}\,\mathrm{Hz}$（对应真空中波长约 555 nm）的单色辐射的光源在给定方向上的发光强度，该方向上的辐射强度为 $1/683 \,\mathrm{W/sr}$。
换言之，在 $V(\lambda)$ 峰值波长 555 nm 处，$1\,\mathrm{W}$ 的辐射通量恰好对应 $683\,\mathrm{lm}$ 的光通量。这一数值在 1979 年国际计量大会上被确认，并一直沿用至今。它是连接辐射度学（瓦特）与光度学（流明）两个量值体系的固定换算系数。
3.4 发光强度（Luminous Intensity）
发光强度 $I_v$ 是点光源在单位立体角内发出的光通量，单位为坎德拉（$\mathrm{cd}$），$\mathrm{1\,cd = 1\,lm/sr}$：
$$ I_v = \frac{\mathrm{d}\Phi_v}{\mathrm{d}\Omega} $$
坎德拉是国际单位制（SI）七个基本单位之一，也是光度学体系的基石。
3.5 亮度（Luminance）
亮度 $L_v$ 是面光源在单位投影面积、单位立体角内的光通量，单位为 $\mathrm{cd/m^2}$（有时称 nit）：
$$ L_v = \frac{\mathrm{d}^2\Phi_v}{\mathrm{d}A \cdot \cos\theta \cdot \mathrm{d}\Omega} $$
对于朗伯辐射体（单面发光），辐射度量与光度量之间存在简洁的转换关系：
$$ \Phi_v = 683 \times \Phi_e \times K \times 10^{-6} = \pi \times L_v \times S \times 10^{-6} $$
其中 $K$ 为光源的视觉函数加权平均系数，$S$ 为发光面积。该公式在 LED 显示屏、背光模组等朗伯型面光源的光通量评估中具有重要的工程实用价值。
3.6 照度（Illuminance）
照度 $E_v$ 是入射到单位面积上的光通量，单位为勒克斯（$\mathrm{lx}$），$\mathrm{1\,lx = 1\,lm/m^2}$：
$$ E_v = \frac{\mathrm{d}\Phi_v}{\mathrm{d}A} $$
照度是照明工程中最常用的设计指标，国家标准对不同环境和作业类型有明确的照度要求。
3.7 光度学物理量小结
| 物理量 | 符号 | 定义 | 单位 |
|---|---|---|---|
| 光通量 | $\Phi_v$ | 辐射通量的 $V(\lambda)$ 加权积分 | $\mathrm{lm}$ |
| 发光强度 | $I_v$ | 单位立体角的光通量 | $\mathrm{cd = lm/sr}$ |
| 亮度 | $L_v$ | 单位投影面积·单位立体角的光通量 | $\mathrm{cd/m^2}$ |
| 照度 | $E_v$ | 单位面积的入射光通量 | $\mathrm{lx = lm/m^2}$ |
四、辐射度量与光度量之间的转换
辐射度量与光度量之间的转换，本质上是将物理能量谱 $\Phi_e(\lambda)$ 通过人眼视觉函数 $V(\lambda)$ 进行加权积分，再乘以最大光谱光视效能 $K_m$。其通用关系式为：
$$ \Phi_v(\lambda) = 683 \times V(\lambda) \times \Phi_e(\lambda) $$
对整个可见光波段积分，得到累积光通量：
$$ \Phi_v = 683 \int_{380}^{780} \Phi_e(\lambda) \, V(\lambda) \,\mathrm{d}\lambda $$
离散形式为：
$$ \Phi_v = 683 \sum_{\lambda=380}^{780} \Phi_e(\lambda) \, V(\lambda) \,\Delta\lambda $$
4.1 转换流程示意
光谱辐射通量 Φ_e(λ) ──× V(λ)──▶ 加权光谱
│
▼
× 683 lm/W
│
▼
光谱光通量 Φ_v(λ)
│
Σ(积分/求和)
│
▼
总光通量 Φ_v (lm)
4.2 光谱光视效能 $K(\lambda)$
由 CIE 定义的光谱光视效能函数 $K(\lambda)$ 满足：
$$ K(\lambda) = 683 \times V(\lambda) \quad [\mathrm{lm/W}] $$
$K(\lambda)$ 直观地给出了每个波长处"1 W 辐射通量对应多少流明光通量"的转换系数。在 555 nm 处达到最大值 $683\,\mathrm{lm/W}$，在可见光波段边缘处趋近于零。
4.3 其他辐射度-光度对应量的转换
辐射度量与光度量之间存在一一对应的转换关系，如下表所示：
| 辐射度量（下标 e） | 光度量（下标 v） | 转换关系 |
|---|---|---|
| 辐射通量 $\Phi_e$ | 光通量 $\Phi_v$ | $\Phi_v = 683 \int \Phi_e(\lambda)V(\lambda)\mathrm{d}\lambda$ |
| 辐射强度 $I_e$ | 发光强度 $I_v$ | $I_v = 683 \int I_e(\lambda)V(\lambda)\mathrm{d}\lambda$ |
| 辐亮度 $L_e$ | 亮度 $L_v$ | $L_v = 683 \int L_e(\lambda)V(\lambda)\mathrm{d}\lambda$ |
| 辐照度 $E_e$ | 照度 $E_v$ | $E_v = 683 \int E_e(\lambda)V(\lambda)\mathrm{d}\lambda$ |
五、色度学基础
色度学（Colorimetry）将光谱功率分布映射为人类颜色感知空间的坐标。其理论基础是三色学说（Young–Helmholtz theory）——人眼视网膜中存在三种对不同波长敏感度各异的视锥细胞（L、M、S 型），任何颜色感知均可由三者响应强度的线性组合来描述。
5.1 CIE 1931 标准色度观察者
CIE 于 1931 年建立了 CIE 1931 标准色度观察者（CIE 1931 Standard Colorimetric Observer），定义了三条颜色匹配函数（Color Matching Functions, CMF）：$\bar{x}(\lambda)$、$\bar{y}(\lambda)$ 和 $\bar{z}(\lambda)$。
这三条函数是通过对大量观察者的颜色匹配实验进行统计平均得到的。值得注意的是，$\bar{y}(\lambda)$ 被刻意设计为与明视觉光谱光视效率函数 $V(\lambda)$ 完全一致——这一巧妙的设计使色度学中的 $Y$ 值（三刺激值之一）同时兼具光度学中光通量的意义，成为辐射度学、光度学与色度学交汇的关键节点。
CIE 1931 标准色度观察者基于 2° 视场角（小视场，对应视网膜中央凹区域的视锥细胞分布）。对于需要10° 视场角（大视场）的应用场景，应使用 CIE 1964 补充标准色度观察者，其颜色匹配函数有所不同，对短波长区域的权重更高。
5.2 三刺激值 $X$、$Y$、$Z$
三刺激值（Tristimulus Values）是光谱功率分布 $\Phi(\lambda)$ 与三条颜色匹配函数分别积分的结果：
$$ X = k \int_{380}^{780} \Phi(\lambda) \, \bar{x}(\lambda) \,\mathrm{d}\lambda $$
$$ Y = k \int_{380}^{780} \Phi(\lambda) \, \bar{y}(\lambda) \,\mathrm{d}\lambda $$
$$ Z = k \int_{380}^{780} \Phi(\lambda) \, \bar{z}(\lambda) \,\mathrm{d}\lambda $$
其中 $k$ 为归一化系数。对光源色（自发光体）而言，$k = 683 \,\mathrm{lm/W}$，此时 $Y$ 即为光通量；对物体色（反射或透射体）而言，$k$ 的选择使完全漫反射体的 $Y = 100$，此时 $\Phi(\lambda)$ 为相对光谱功率分布。
离散求和形式为：
$$ X = k \sum \Phi(\lambda) \, \bar{x}(\lambda) \,\Delta\lambda $$
$$ Y = k \sum \Phi(\lambda) \, \bar{y}(\lambda) \,\Delta\lambda $$
$$ Z = k \sum \Phi(\lambda) \, \bar{z}(\lambda) \,\Delta\lambda $$
5.3 色品坐标 $x$、$y$ 与 CIE 色度图
三刺激值 $X$、$Y$、$Z$ 包含了颜色的亮度和色度信息。将亮度分量分离后，可得到仅描述色彩品质（色相与饱和度）的色品坐标（Chromaticity Coordinates）：
$$ x = \frac{X}{X + Y + Z} $$
$$ y = \frac{Y}{X + Y + Z} $$
$$ z = \frac{Z}{X + Y + Z} = 1 - x - y $$
由定义可知，$x$ 与 $y$ 即可完全确定色品（$z$ 为冗余信息），所有物理可实现的颜色在 $(x, y)$ 坐标系中均落在一个马蹄形区域内——这就是著名的 CIE 1931 色度图（CIE Chromaticity Diagram）。
CIE 色度图的几个关键特征：
- 光谱轨迹（马蹄形边界）：由单色光（380–780 nm）的色品坐标连成的曲线，包含了最饱和的颜色；
- 紫红线：连接光谱轨迹短波端（380 nm）与长波端（780 nm）的直线，其上的颜色为两种极端波长的混合色，自然界中无对应的单色光；
- 等能白光点（E点）：$x = y = z = 1/3$，即 $x = 0.3333$、$y = 0.3333$；
- 色温轨迹（Planckian Locus）：黑体在不同温度下的色品坐标连成的曲线，白光源的色品通常以其最近的黑体色温（CCT，相关色温）来描述。
5.4 色品坐标计算示例
给定一组简化的光谱数据（以等间隔 $\Delta\lambda = 10\,\mathrm{nm}$ 采样），三刺激值的计算步骤为：
- 读取各波长的 $\bar{x}(\lambda)$、$\bar{y}(\lambda)$、$\bar{z}(\lambda)$ 值（CIE 标准数据表）；
- 计算加权积：$\Phi(\lambda) \cdot \bar{x}(\lambda)$、$\Phi(\lambda) \cdot \bar{y}(\lambda)$、$\Phi(\lambda) \cdot \bar{z}(\lambda)$；
- 求和并乘以 $k$ 与 $\Delta\lambda$ 得到 $X$、$Y$、$Z$；
- 计算 $x = X/(X+Y+Z)$、$y = Y/(X+Y+Z)$。
六、CIE 1976 $L^*a^*b^*$ 均匀色彩空间
CIE 1931 色度图虽然奠定了色度学基础，但其色品坐标 $(x, y)$ 在视觉上是不均匀的：在色度图的不同区域，相同几何距离对应的感知颜色差异并不相同。这给色差计算和色度公差设定带来了严重不便。
为了解决这一问题，CIE 于 1976 年推荐了 CIE 1976 $L^a^b^*$ 色彩空间（简称 CIELAB），它将三刺激值 $X$、$Y$、$Z$ 通过非线性变换映射到一个近似视觉均匀的三维空间。
6.1 坐标定义
CIELAB 空间的三个坐标定义为：
亮度：
$$ L^* = 116 \cdot f\left(\frac{Y}{Y_n}\right) - 16 $$
红-绿色度：
$$ a^* = 500 \left[ f\left(\frac{X}{X_n}\right) - f\left(\frac{Y}{Y_n}\right) \right] $$
黄-蓝色度：
$$ b^* = 200 \left[ f\left(\frac{Y}{Y_n}\right) - f\left(\frac{Z}{Z_n}\right) \right] $$
其中 $X_n$、$Y_n$、$Z_n$ 为参考白点（如 D65 标准照明体）的三刺激值，函数 $f(t)$ 为分段函数：
$$ f(t) = \begin{cases} t^{1/3}, & t > 0.008856 \[8pt] 7.787 \, t + \dfrac{16}{116}, & t \leq 0.008856 \end{cases} $$
6.2 坐标含义
- $L^*$（明度，Lightness）：取值范围 $0 \sim 100$，$L^=0$ 为纯黑，$L^=100$ 为参考白；
- $a^*$（红-绿轴）：正值偏红，负值偏绿；
- $b^*$（黄-蓝轴）：正值偏黄，负值偏蓝。
6.3 色差公式
CIELAB 空间的最大实用价值在于色差计算。两个颜色之间的色差 $\Delta E^*_{ab}$ 定义为欧氏距离：
$$ \Delta E^_{ab} = \sqrt{(L^_1 - L^_2)^2 + (a^_1 - a^_2)^2 + (b^_1 - b^*_2)^2} $$
通常情况下，$\Delta E^_{ab} < 1$ 表示色差人眼基本不可察觉；$\Delta E^{ab} \approx 2-3$ 为可察觉但较小的色差；$\Delta E^*{ab} > 5$ 则色差明显。该公式被广泛应用于 LED 分选（Binning）、显示器色准评估、印刷品色差控制等工业领域。
七、在实际光谱测量中的应用场景
7.1 LED 光源的综合评价
光谱辐射度、光度与色度计算在 LED 行业中的应用最为典型。通过光谱辐射计（Spectroradiometer）或积分球-光谱仪系统测量 LED 的绝对光谱功率分布 $\Phi_e(\lambda)$，可一次性导出全部关键参数：
- 辐射度量：辐射通量 $\Phi_e$（$\mathrm{W}$）；
- 光度量：光通量 $\Phi_v$（$\mathrm{lm}$）、发光效率（$\mathrm{lm/W}$）；
- 色度量：色品坐标 $(x, y)$、相关色温（CCT）、显色指数（CRI / $R_a$）；
- 色空间坐标：$L^a^b^$ 或 $L^u^v^$ 坐标，用于 LED 分 BIN。
7.2 显示屏与面光源测量
对于 OLED / LCD 显示屏、背光模组等面光源，辐亮度（亮度）测量是核心需求。通过成像亮度计或光谱辐射计配合望远镜头测量辐亮度 $L_e$（或亮度 $L_v$），再结合发光面积计算光通量。
利用朗伯辐射体的特性（$\Phi_v = \pi \cdot L_v \cdot S$），可在已知发光面积和亮度的情况下快速估算总光通量，这对于大面积光源的在线检测具有显著的工程便利性。
7.3 颜色品质控制与色差评估
在纺织印染、塑料着色、汽车漆面、印刷包装等行业，CIELAB 色差 $\Delta E^_{ab}$ 是颜色品质控制的通用语言。通过分光光度计或光谱辐射计测量样品的光谱反射率，计算 $X$、$Y$、$Z$ 和 $L^a^b^$，与目标颜色比较 $\Delta E^*_{ab}$，即可判断产品颜色是否在允许公差范围内。
7.4 太阳辐射与紫外测量
在光伏、气象、光生物安全等领域，辐照度 $E_e$（尤其是紫外波段的辐照度）是核心测量参数。通过光谱辐射计配合余弦校正器（漫射器）测量太阳或人工光源的绝对光谱辐照度，可计算 UV-A、UV-B、UV-C 各波段积分辐照度，用于光生物安全等级（IEC 62471）判定和光伏电池效率评估。
7.5 照明工程设计与评估
照明工程中，照度 $E_v$ 和亮度 $L_v$ 的测量与计算是设计的基础。通过光谱辐射计测量光源的光谱功率分布，可计算光通量和发光效率，进而结合配光曲线评估空间照度分布。同时，色品坐标和相关色温（CCT）用于评估照明光品质（色温适宜性、显色性），确保满足不同应用场景（办公室、商场、道路、体育场馆等）的照明标准。
八、总结
光谱辐射度、光度与色度计算共同构成了光学测量与照明工程的理论基石。三者以光谱功率分布为共同输入，沿不同路径输出不同维度的物理参数：
- 辐射度学提供与观察者无关的纯物理能量描述，是所有光辐射测量的原点；
- 光度学通过人眼视觉函数 $V(\lambda)$ 和转换系数 $683\,\mathrm{lm/W}$，将物理能量谱转化为人眼感知的"明暗"度量；
- 色度学通过 CIE 标准色度观察者的颜色匹配函数，将光谱映射为三维颜色空间坐标，并借助 CIELAB 均匀色彩空间实现色差的定量评估。
理解这三个体系各自的定义、相互间的转换逻辑以及离散积分的工程计算方法，是进行光谱测量系统设计、数据分析和结果判读的前提。随着光谱测量技术的不断进步（如阵列式光谱仪的普及、微型光谱传感器的小型化），辐射-光度-色度一体化测量已成为行业标准配置，对相关原理的深入掌握将直接提升测量方案设计的合理性和数据解释的可靠性。
本文由谱研互联整理编写。