第三章 光学基础：傅里叶光学、成像系统模型与传递函数
压缩感知光学系统的核心任务，是在物理域完成对场景信息的编码与采样。要设计好这样的系统，必须深刻理解光场在成像链路中的传播规律。本章从傅里叶光学的标量衍射理论出发，建立成像系统的线性系统模型，并重点剖析点扩散函数（PSF）与调制传递函数（MTF）这两个贯穿光学设计、仿真与实验评价的核心工具。只有掌握了这些基础，才能理解后续的编码孔径、散射介质成像以及计算重建算法为何有效，以及它们各自的能力边界在哪里。
3.1 傅里叶光学核心原理：从衍射到空间频谱
傅里叶光学将光场视为复振幅的二维分布，并利用傅里叶变换将空间域信号分解为不同方向传播的平面波分量。在近轴近似下，光场传播满足菲涅尔衍射积分，而夫琅禾费远场衍射则直接对应光场的傅里叶变换。这一性质使得透镜天然成为光学傅里叶变换器：在透镜后焦面上，光场分布正比于入射光场的二维傅里叶变换。
对于压缩感知光学系统，这一原理具有双重意义。第一，它提供了设计编码掩膜的理论工具——通过在频谱面放置振幅或相位调制器件，可以直接对场景的空间频率进行加权或随机化。第二，它解释了成像系统分辨率极限的物理根源：有限孔径相当于在频域施加了一个低通滤波，丢失了高于截止频率的空间信息，而这些信息无法通过任何计算手段恢复。
关键公式：相干成像系统中，像面复振幅 U_i 与物面复振幅 U_o 满足：
U_i(x,y) = U_o(x,y) ⊗ h(x,y)
H(u,v) = ℱ{h(x,y)} = P(-λd_i u, -λd_i v)其中 h 为相干点扩散函数，H 为相干传递函数（CTF），P 为光瞳函数，λ 为波长，d_i 为像距。非相干成像系统则对光强线性叠加，传递函数为 OTF = CTF ⊗ CTF*。
3.2 成像系统模型与点扩散函数（PSF）
一个完整的成像系统模型包含物平面、光学系统、像平面和探测器采样四个环节。在空间域，非相干成像系统是光强线性系统，像面照度等于物面照度与PSF的卷积。PSF描述了理想点光源经过系统后在像面产生的光斑分布，它综合了衍射、像差、离焦以及编码掩膜的影响。
在压缩感知成像中，PSF的设计直接决定了测量矩阵的性质。传统成像追求PSF尽可能接近δ函数以获得清晰图像；而压缩感知成像（如编码孔径相机）则刻意设计宽域、随机或伪随机的PSF，使不同场景点的光在探测器上混叠，再通过稀疏重建算法解混叠。此时PSF的形态决定了测量矩阵的相干性、条件数以及重建质量上限。
| PSF类型 | 频域特征 | 压缩感知适用性 | 典型应用 |
|---|---|---|---|
| 艾里斑（衍射极限） | 低通，截止频率明确 | 差（信息丢失不可逆） | 传统显微/摄影 |
| 编码孔径PSF | 频谱宽且起伏大 | 优（满足RIP或低相干） | 单像素相机、编码孔径成像 |
| 散射介质PSF | 近似随机相位分布 | 优（高随机性） | 透过散射介质成像 |
| 离焦PSF | 频谱有零点 | 中（需多帧或编码） | 深度估计、扩展景深 |
实际测量PSF时，通常使用亚像素尺寸的点光源（如针孔或荧光微球），在像面记录其强度分布。对于压缩感知系统，还需要测量不同视场位置的PSF，因为大视场下PSF往往具有空间移变性。
3.3 调制传递函数（MTF）与系统评价
MTF定义为OTF的模，即不同空间频率下系统对比度的传递能力。MTF曲线从低频的1（或归一化值）逐渐下降，在截止频率处降为0。MTF不仅用于评价传统光学系统的成像质量，在压缩感知光学中同样重要：它决定了系统对场景高频信息的物理采集能力，即使算法再强大，也无法恢复MTF为零的频率分量。
MTF的测量方法包括：点光源法（对PSF做傅里叶变换取模）、狭缝法（对线扩散函数做一维傅里叶变换）、刀口法（对边缘响应求导得到LSF再变换）以及正弦光栅法。在压缩感知系统中，由于PSF可能不是中心对称的，需要测量二维MTF或不同方向的一维MTF剖面。
工程经验：在搭建压缩感知成像系统时，建议在系统集成后立即测量实际MTF，而不是依赖设计仿真值。常见偏差来源包括：透镜装配应力导致的像散、编码掩膜与探测器之间的间距误差（引入额外衍射）、以及环境温度变化引起的焦面漂移。实测MTF应作为重建算法中前向模型校正的依据。
3.4 实战案例：编码孔径成像系统的PSF标定与MTF分析
背景与目标：某实验室搭建了一套可见光编码孔径成像系统，目标是对1m距离处的场景进行压缩感知重建，压缩比为10:1。系统由主透镜（焦距50mm，F/#=2.8）、编码掩膜（随机二值图案，特征尺寸50μm）和CMOS探测器（像元尺寸5.5μm，2048×2048）组成。
初始现象：系统集成后，直接拍摄点光源（10μm针孔照明）得到的PSF与仿真结果差异显著。实测PSF主瓣宽度约为仿真值的1.8倍，且旁瓣分布不对称，导致重建图像出现明显的条纹伪影。
分析判断过程：首先检查掩膜与探测器间距，设计值为2mm，实测为2.15mm（装配误差）。利用菲涅尔衍射公式计算该间距误差对PSF的影响，发现主瓣展宽约30%。其次，检查掩膜图案的加工精度，发现最小特征尺寸实际为58μm而非设计值50μm，这进一步改变了衍射分布。最后，通过测量不同视场位置的PSF，确认存在轻微的场曲，边缘视场PSF比中心展宽约20%。
实际操作步骤：
- 使用精密垫片将掩膜-探测器间距调整至2.00±0.05mm，并用激光干涉仪验证平行度。
- 重新加工掩膜，确保特征尺寸公差控制在±2μm以内。
- 在3×3共9个视场位置分别采集PSF，建立PSF库。
- 对每个PSF做傅里叶变换得到OTF，取模得到MTF，并拟合为解析模型（如高斯-洛伦兹混合模型）供重建算法调用。
- 将校正后的前向模型代入压缩感知重建框架（TVAL3算法），重新处理测试图像。
结果对比：校正前，重建图像的峰值信噪比（PSNR）为21.3dB，结构相似度（SSIM）为0.72；校正后PSNR提升至27.8dB，SSIM达到0.91。条纹伪影基本消失，边缘清晰度接近衍射极限。
风险与代价：PSF标定过程耗时约4小时，且需要高精度位移台和点光源。此外，PSF库的存储和调用增加了重建算法的内存开销（9个PSF，每个512×512浮点数，约9.4MB）。
复盘经验：压缩感知光学系统的性能对PSF精度极其敏感，任何毫米级的机械误差都会导致重建质量显著下降。建议在系统设计阶段预留PSF在线标定接口（如内置LED点光源），并建立定期校准制度。
3.5 常见误区与总结复盘
误区一：认为MTF越高越好。在压缩感知系统中，MTF的“形状”比“高度”更重要。一个在宽频带内保持中等且非零的MTF（如编码孔径）往往优于一个在低频很高但高频迅速截止的MTF（如普通透镜），因为前者提供了更丰富的频率混叠信息。
误区二：忽略PSF的空间移变性。许多初学者假设PSF在全视场恒定，直接使用中心视场PSF进行重建。对于大视场或大孔径系统，这种近似会导致边缘区域重建失败。正确做法是分区标定或使用视场相关的PSF模型。
误区三：混淆相干与非相干传递函数。激光照明系统必须使用CTF（振幅线性），而自然光或LED照明系统使用OTF（强度线性）。错误使用会导致重建算法中的前向模型错误，产生无法收敛或伪影严重的结果。
误区四：忽视噪声对MTF的影响。实测MTF在接近截止频率处往往被噪声抬高，因为噪声功率谱是平坦的。建议在MTF测量时进行多次采集平均，并扣除暗场噪声，或使用正则化方法平滑MTF曲线。
总结：傅里叶光学为压缩感知成像提供了物理可实现性的理论边界。PSF和MTF不仅是评价工具，更是设计工具——它们直接决定了测量矩阵的性质和重建算法的性能上限。在实际工程中，务必以实测数据为准，将光学标定作为系统集成的必要环节，而非可选项。