4、杂散光关键参数：点源透过率(PST)、杂散光系数、辐射传输方程、BRDF/BSDF

杂散光分析的核心在于将“看不见的光”转化为“可量化的数”。在工程实践中，我们无法直接测量到达像面的每一束杂散光，但可以通过一系列标准化参数和数学模型，对杂散光的产生、传播和抑制效果进行精确评估。本章将围绕四个关键参数展开：点源透过率（PST）、杂散光系数、辐射传输方程以及双向反射分布函数（BRDF/BSDF）。这四个参数构成了杂散光分析从“表面属性”到“系统指标”的完整逻辑链：BRDF描述微观表面散射特性，辐射传输方程计算能量在系统内的传递路径，PST和杂散光系数则从系统层面评价杂散光的抑制能力。

上图展示了四个参数之间的逻辑关系：BRDF/BSDF是底层输入，辐射传输方程是计算工具，PST是设计验证手段，杂散光系数是最终验收指标。下面逐一深入解析。

4.1 点源透过率（PST）——系统杂散光抑制能力的“标尺”

核心原理

点源透过率（Point Source Transmittance, PST）定义为：离轴角为θ的点源，其辐射经过光学系统后在像面产生的辐照度，与轴上点源（θ=0）在像面产生的辐照度之比。数学表达式为：

PST(θ) = E_off-axis(θ) / E_on-axis(0)

其中E_off-axis(θ)是离轴角θ处的点源在像面产生的辐照度，E_on-axis(0)是轴上点源在像面产生的辐照度。PST是一个无量纲的比值，通常用对数坐标（dB）表示。PST值越小，说明系统对离轴杂散光的抑制能力越强。例如，一个高质量的空间相机在离轴角30°时的PST可能达到10⁻⁶量级（即-120dB）。

判断方法与操作流程

PST的测量和仿真需要遵循标准化流程：

1. 定义离轴角范围：根据系统视场和实际工作环境，确定需要评估的离轴角范围（通常为0°~90°）。
2. 建立仿真模型：在光学仿真软件（如TracePro、FRED、Zemax）中建立包含机械结构、光学元件表面属性（BRDF数据）的完整模型。
3. 设置点源：在无穷远处放置点光源，改变其离轴角，分别计算像面辐照度。
4. 归一化处理：将各离轴角的像面辐照度除以轴上点源的像面辐照度，得到PST曲线。
5. 实验验证：在暗室中使用平行光管和旋转平台，测量不同离轴角下的像面响应，与仿真结果对比。

案例拆解：星敏感器PST验证

某星敏感器要求在太阳离轴角30°时，PST≤1×10⁻⁶，以确保白天也能探测到暗弱恒星。初始设计采用遮光罩+内壁发黑处理，仿真PST为8.7×10⁻⁷，满足要求。但实测发现PST为9.2×10⁻⁷，接近指标边界。进一步分析发现，遮光罩内壁的螺纹结构在特定角度产生了镜面反射。解决方案：将遮光罩内壁改为阶梯状消光环，并增加一道挡光环。修改后实测PST降至2.1×10⁻⁷，余量充足。此案例说明，PST仿真必须考虑机械结构的二次反射，且实测验证不可或缺。

4.2 杂散光系数——工程验收的“及格线”

核心原理

杂散光系数（Veiling Glare Index, VGI）定义为：像面某点的杂散光照度E_stray与该点的总照度E_total（信号光+杂散光）之比，即VGI = E_stray / E_total。与PST不同，杂散光系数描述的是在正常工作视场内，杂散光对成像质量的污染程度。它直接反映了图像对比度的损失情况。例如，VGI=5%意味着图像最暗部分的亮度被抬高了5%，导致动态范围压缩。

判断方法与操作流程

杂散光系数的测量通常采用“黑斑法”：

1. 在均匀亮度背景中放置一个黑斑（吸收率≥99.9%）。
2. 用待测系统对黑斑成像，测量黑斑像中心的灰度值D_stray。
3. 移开黑斑，测量同一位置的灰度值D_total（此时为纯背景亮度）。
4. 计算VGI = D_stray / D_total。

对于仿真分析，杂散光系数可以通过在像面设置虚拟探测器，统计所有非设计光路到达像面的能量总和，除以总接收能量得到。

案例拆解：工业镜头杂散光优化

一款用于机器视觉的C接口镜头，在逆光拍摄时出现明显光晕，实测VGI达到8.5%。分析发现，镜筒内壁的螺纹反射和镜片边缘的未镀膜区域是主要贡献源。优化措施：①镜筒内壁改涂高吸收率消光漆（吸收率从85%提升至97%）；②镜片边缘涂黑漆；③在光阑片表面增加微结构。优化后VGI降至2.1%，图像对比度显著提升。该案例表明，杂散光系数优化往往需要多管齐下，单一措施难以奏效。

4.3 辐射传输方程——杂散光能量传递的“计算引擎”

核心原理

辐射传输方程（Radiative Transfer Equation, RTE）描述了辐射能在介质中传播时，由于吸收、散射、发射等过程导致的能量变化。在杂散光分析中，RTE被简化为几何光学框架下的能量追踪方程：

L(r, s) = L_emitted(r, s) + ∫ f_r(s_i, s) · L(r, s_i) · cos(θ_i) dω_i

其中L(r,s)是位置r处沿方向s的辐射亮度，f_r是BRDF函数，θ_i是入射方向与表面法线的夹角。该方程表明，任意表面出射的辐射亮度等于自身发射的辐射亮度加上所有入射方向散射贡献的积分。在杂散光分析中，我们通常忽略自身发射（除非是红外系统），重点计算散射项。

判断方法与操作流程

实际工程中，RTE的求解采用蒙特卡洛光线追迹法：

1. 表面属性赋值：为每个光学面和机械面赋予BRDF/BSDF数据（实测或模型库）。
2. 光线采样：从光源发射大量光线（通常10⁶~10⁸条），每条光线携带能量权重。
3. 追迹与分裂：光线遇到表面时，根据BRDF概率分布产生反射/透射/吸收，能量按概率分配。
4. 统计收敛：在像面统计到达光线的能量总和，当统计误差小于设定阈值（如1%）时停止追迹。

案例拆解：红外系统冷反射分析

某制冷型红外探测器系统，在低温环境下出现图像中心亮斑。通过RTE分析发现，冷反射路径为：冷探测器→透镜表面反射→回到探测器。该路径虽然能量极低（约10⁻⁵量级），但由于探测器自身温度低（80K），与背景温差大，导致冷反射信号明显。解决方案：在透镜表面镀高效增透膜（反射率从1%降至0.2%），并在镜筒内壁增加低温挡光环。优化后冷反射信号降低至原来的1/20。此案例说明，RTE不仅能计算杂散光总量，还能识别关键传播路径，为针对性优化提供依据。

4.4 BRDF/BSDF——表面散射特性的“指纹”

核心原理

双向反射分布函数（Bidirectional Reflectance Distribution Function, BRDF）定义为：出射方向辐射亮度与入射方向辐照度的比值，即f_r(θ_i, φ_i, θ_r, φ_r) = dL_r / dE_i。其单位是sr⁻¹。BRDF描述了表面在任意入射和出射方向组合下的散射特性。BSDF（Bidirectional Scattering Distribution Function）是BRDF和BTDF（透射）的统称。BRDF的测量通常在变角光度计上进行，通过改变入射角和探测角，获取完整的散射分布数据。

判断方法与操作流程

BRDF数据在杂散光分析中的应用流程：

1. 数据测量：使用变角光度计测量关键表面的BRDF，测量波长应与系统工作波长一致。
2. 数据拟合：将测量数据拟合为解析模型（如Harvey模型、ABg模型），便于仿真软件调用。
3. 模型验证：将拟合模型与原始数据对比，确保误差在可接受范围内（通常R²≥0.95）。
4. 库管理：建立表面属性库，记录不同材料、不同工艺（抛光、发黑、镀膜）的BRDF数据。

# Harvey模型示例
BRDF(θ_r) = b0 / (1 + (sin(θ_r) - sin(θ_i))² / l²)^(s/2)
# 参数：b0=0.01, l=0.1, s=2.5

案例拆解：遮光罩内壁材料选择

某空间望远镜的遮光罩内壁需要在宽光谱范围（0.4~2.5μm）内具有低散射特性。候选材料包括：黑色阳极氧化铝、碳纳米管涂层、聚酰亚胺黑膜。通过测量三种材料在632.8nm和1550nm波长的BRDF，发现碳纳米管涂层在10°散射角处的BRDF比黑色阳极氧化铝低两个数量级，但成本高出5倍。综合考虑性能与成本，最终选择在关键区域（直接可见路径）使用碳纳米管涂层，其他区域使用黑色阳极氧化铝。该案例表明，BRDF数据是材料选型的量化依据。

4.5 参数间的协同应用与工程决策

在实际工程中，四个参数并非孤立使用，而是形成一套完整的分析流程。首先，通过BRDF测量建立表面属性库；其次，利用辐射传输方程计算系统内各路径的能量贡献，识别主要杂散光路径；然后，通过PST曲线验证系统在不同离轴角下的抑制能力；最后，用杂散光系数作为最终验收指标。下表总结了各参数的应用场景和典型指标：

需要强调的是，这四个参数构成了一个闭环迭代过程：当杂散光系数超标时，通过PST定位问题离轴角范围，利用辐射传输方程分析该角度下的主要传播路径，再通过BRDF数据判断是哪个表面贡献最大，从而指导表面处理工艺改进或结构优化。这种系统化的分析方法，是解决复杂杂散光问题的关键。