📚 课程目录第2章光学基础回顾

2、光学基础回顾：几何光学、波动光学、衍射极限、成像系统的基本参数（分辨率、景深、视场）

计算成像技术的核心目标，是突破传统光学系统在物理硬件层面所受到的制约。要理解这些突破是如何实现的，首先必须扎实地回顾经典光学的基础框架。本章并非简单复述大学物理，而是从计算成像的视角，重新审视几何光学、波动光学、衍射极限以及成像系统的基本参数。我们将重点厘清哪些限制是硬性的物理极限，哪些限制可以通过后期计算来缓解或规避，从而为后续章节中关于编码孔径、相位恢复、超分辨率重建等技术打下坚实的理论基础。

2.1 几何光学：成像系统的基石与近似边界

核心原理：光线模型与理想成像

几何光学将光视为沿直线传播的“光线”，其核心定律包括光的直线传播定律、独立传播定律以及反射和折射定律。在计算成像中，我们通常使用几何光学来建立成像系统的一阶近似模型，即忽略光的波动性，仅考虑光线的传播路径。薄透镜公式 \frac{1}{f} = \frac{1}{u} + \frac{1}{v} 是描述物距 u、像距 v 和焦距 f 关系的基础。然而，实际透镜并非理想，会产生球差、彗差、像散、场曲和畸变等单色像差，以及色差。这些像差在传统成像中是需要极力消除的，但在计算成像中，它们可以被建模为已知的“点扩散函数”空间变化，并通过算法进行部分补偿。

计算成像视角：几何光学提供了光线追迹的基础。在计算成像系统设计初期，我们利用光线追迹软件（如 Zemax、Code V）模拟光路，确定系统的视场角、F数、畸变量等参数。这些参数将作为后续算法建模的先验信息。例如，在计算摄影中，我们常利用已知的镜头畸变模型对图像进行校正，这本质上是对几何光学模型的逆运算。

判断方法：光线追迹与矩阵光学

对于复杂的多透镜系统，逐面应用折射定律进行光线追迹是标准方法。在计算成像中，更常用的是ABCD 矩阵（传输矩阵）方法，它将光线参数（高度 y 和角度 θ）通过矩阵乘法进行线性变换。这种方法不仅适用于共轴球面系统，也适用于分析谐振腔、光束传播等场景。判断一个系统是否可以用几何光学近似，主要依据是系统中光学元件的尺寸是否远大于光波长，且不考虑衍射效应。

# 薄透镜的ABCD矩阵示例（Python伪代码）
import numpy as np
# 自由空间传播距离d
def free_space(d):
    return np.array([[1, d], [0, 1]])
# 焦距为f的薄透镜
def thin_lens(f):
    return np.array([[1, 0], [-1/f, 1]])
# 组合系统：透镜+距离d+透镜
M = thin_lens(50) @ free_space(100) @ thin_lens(30)
print(M)

案例拆解：手机广角镜头的畸变校正

某手机摄像头采用 120° 超广角镜头，边缘畸变率高达 8%。在传统光学设计中，这被视为缺陷。但在计算成像流程中，我们首先通过标定板获取畸变场映射表，然后在 ISP（图像信号处理器）中利用几何校正算法（如双线性插值或双三次插值）对图像进行重映射。具体步骤包括：1）拍摄棋盘格标定板，提取角点坐标；2）计算理想坐标与实际坐标的映射关系；3）生成校正查找表（LUT）；4）对实时视频流应用 LUT 进行像素重采样。最终，虽然边缘分辨率因插值略有下降，但视觉上的“桶形畸变”被完全消除，且无需增加额外的光学镜片，这体现了计算成像对光学缺陷的“容忍-校正”设计哲学。

常见误区与总结

误区：认为几何光学是“过时”的理论。实际上，在计算成像中，几何光学是建立系统模型的基础。如果几何模型不准确，后续的逆问题求解将无从谈起。例如，在光场相机中，微透镜阵列的几何位置标定精度直接影响重聚焦质量。

复盘：几何光学提供了“光线从哪里来”的确定性描述，是计算成像中正向建模的基石。我们需要熟练掌握光线追迹、矩阵光学和像差理论，并理解其适用边界——当特征尺寸接近波长时，必须转向波动光学。

2.2 波动光学：从光线到复振幅的跨越

核心原理：光的电磁本质与标量衍射理论

波动光学将光视为电磁波，用复振幅 U(x,y) 描述光场的振幅和相位。这一视角是计算成像的核心，因为相位信息往往携带了物体深度、透明样本内部结构等关键信息，而传统强度探测器（CMOS/CCD）只能记录 |U|²。惠更斯-菲涅尔原理指出，波前上的每一点都可以看作新的球面波源，它们的相干叠加形成了后续的波前。基尔霍夫衍射公式和瑞利-索末菲衍射公式给出了严格的数学描述。在计算成像中，我们常用角谱传播或菲涅尔衍射近似来模拟光场的传播过程。

关键公式：菲涅尔衍射积分 U(x,y) = \frac{e^{jkz}}{j\lambda z} \iint U(\xi,\eta) \exp\left\{ \frac{jk}{2z} [(x-\xi)^2 + (y-\eta)^2] \right\} d\xi d\eta。该公式是计算全息、数字重聚焦等技术的理论基础。

判断方法：相干与非相干成像系统

成像系统可分为相干成像系统（如激光照明）和非相干成像系统（如自然光照明）。两者的传递函数性质截然不同：相干系统对复振幅线性，传递函数为相干传递函数（CTF），截止频率为 f_c = NA/\lambda；非相干系统对光强线性，传递函数为光学传递函数（OTF），其模为调制传递函数（MTF），截止频率为 2f_c。判断系统是相干还是非相干，决定了我们后续是处理复振幅还是处理强度，这直接影响重建算法的设计。

特性	相干成像	非相干成像
照明光源	激光（高相干性）	LED/卤素灯（低相干性）
线性对象	复振幅	光强
截止频率	NA/λ	2NA/λ
典型应用	全息显微、相干衍射成像	普通摄影、机器视觉

案例拆解：基于相位恢复的定量相位显微镜

在生物医学领域，观察未染色透明细胞（如活细胞）时，传统明场显微镜对比度极低。计算成像中的定量相位成像（如 Transport of Intensity Equation, TIE）方法，通过在不同焦面位置拍摄三张强度图像，利用光强传输方程求解出相位分布。具体流程：1）使用高精度电动载物台，在焦面 z0 处拍摄 I1，在 z0+Δz 处拍摄 I2，在 z0-Δz 处拍摄 I3（Δz 通常为 1-5 μm）；2）利用 TIE 方程 \nabla_\perp \cdot (I \nabla_\perp \phi) = -k \frac{\partial I}{\partial z} 进行求解；3）通过傅里叶变换法或多重网格法求解泊松方程，得到定量相位图。该相位图可以精确计算出细胞的光学厚度（干质量），无需染色即可实现无标记定量分析。此案例充分展示了波动光学中“相位”这一物理量在计算成像中的巨大价值。

常见误区与总结

误区：将波动光学简单地等同于“加一个衍射项”。实际上，波动光学引入了相干性这一核心概念。部分相干光场的建模（如使用互相干函数或维格纳分布）远比完全相干或完全非相干复杂，但更贴近实际工程场景。

复盘：波动光学是连接物理光学与信息处理的桥梁。我们必须习惯于用复振幅和传递函数来思考问题，理解衍射是光传播的固有属性，而非“缺陷”。计算成像的许多突破，如合成孔径、结构光照明，都是对波动光学规律的巧妙利用。

2.3 衍射极限与成像系统的基本参数

核心原理：瑞利判据与艾里斑

当一个点光源通过圆形孔径成像时，由于衍射效应，在像面上不会形成一个理想的点，而是形成艾里斑（Airy disk）。艾里斑的半径 r = 1.22 \lambda f / D = 0.61 \lambda / NA，其中 NA 为数值孔径。瑞利判据规定：当两个点光源的艾里斑中心间距等于第一暗环半径时，认为它们刚好可以被分辨。这就是经典的衍射极限。在传统成像中，这个极限由硬件参数（NA 和 λ）决定，无法通过改进探测器像素尺寸来突破。

重要推论：光学系统的截止频率 f_c = NA/\lambda（相干）或 2NA/\lambda（非相干）。这意味着物体通过光学系统后，高于此频率的空间信息将完全丢失。计算成像的目标之一，就是通过编码（如结构光、微扫描）或先验信息（如稀疏性）来恢复这些丢失的高频信息，即“超分辨率”。

判断方法：MTF曲线与系统级参数权衡

在实际工程中，我们通常使用MTF（调制传递函数）曲线来评价系统的分辨率性能。MTF 曲线描述了不同空间频率下对比度的衰减情况。一个成像系统的综合 MTF 是光学系统 MTF、探测器 MTF 和图像处理算法 MTF 的乘积。判断系统是否达到衍射极限，可以对比实际 MTF 曲线与理论衍射极限 MTF 曲线。此外，成像系统还有三个基本参数需要权衡：

分辨率（Resolution）：由光学衍射极限和探测器像素尺寸共同决定。在计算成像中，我们可以通过亚像素位移（超分辨率重建）来提升采样密度，但无法突破光学截止频率。
景深（Depth of Field, DoF）：指能在像面上清晰成像的物方深度范围。传统景深公式为 DoF ≈ 2Fc (s/f)²，其中 F 为 F数，c 为容许弥散圆直径。增大 F数可以增加景深，但会减少进光量并增加衍射影响。计算成像中的波前编码技术（如加入三次相位板）可以使 MTF 对离焦不敏感，从而通过算法恢复出大景深图像。
视场（Field of View, FoV）：指系统能够观察到的物方范围。视场与分辨率通常存在矛盾：视场越大，单个像素对应的物方尺寸越大，分辨率越低。计算成像中的拼接或扫描技术可以在保持高分辨率的同时扩大视场。

案例拆解：波前编码技术扩展景深

传统工业相机在检测高密度电路板时，由于焊点高度不一，常因景深不足导致边缘模糊。某视觉系统采用波前编码技术：在光学系统的孔径面插入一块非球面相位掩膜板（如三次相位板），使得离焦引起的 MTF 变化极小，但引入了固定的像差。具体实施：1）在镜头与传感器之间加入相位板，其相位函数为 φ(x,y) = α(x³ + y³)，α 为优化系数；2）拍摄得到的中间图像是模糊的，但模糊函数对离焦不敏感；3）通过维纳滤波或 Richardson-Lucy 算法进行去卷积，恢复出清晰图像。实验结果表明，在 F/2.0 的光圈下，传统系统景深约为 ±0.1mm，而波前编码系统通过算法重建后，有效景深扩展至 ±0.5mm，且分辨率损失控制在可接受范围内。这一案例展示了如何通过“光学编码+计算解码”来突破传统参数间的制约关系。

常见误区与总结

误区：认为像素数越高，分辨率越高。实际上，当探测器像素尺寸小于艾里斑半径的一半时，继续增加像素数对系统分辨率的提升微乎其微，此时系统处于“衍射受限”状态。计算成像中的超分辨率算法并不能创造超出光学截止频率的信息，只能恢复因采样不足而丢失的信息。

复盘：衍射极限是物理世界给传统成像画下的一条“红线”。但计算成像的思维是：既然高频信息在光学端已经丢失，我们是否可以通过改变光学编码方式（如用多个角度照明、改变孔径形状）来将高频信息“编码”到低频通道中，然后再通过算法解码？这正是傅里叶叠层显微术（FPM）和结构光照明显微术（SIM）的核心思想。理解衍射极限，不是为了接受它，而是为了找到绕过它的路径。

综上所述，几何光学、波动光学和衍射极限构成了计算成像的物理基础。几何光学告诉我们光线如何走，波动光学告诉我们光场如何演化，衍射极限则划定了传统系统的信息容量边界。计算成像的本质，是在这三者之间建立精确的数学模型，并通过计算手段进行逆问题求解，从而在分辨率、景深、视场等参数之间实现新的平衡，甚至突破传统极限。掌握这些基础概念，是深入理解后续编码策略与重构算法的必要前提。