本页目录

探测 III · 成像质量与计算成像

层次:本科核心 + 业界标准语言 + 前沿方向。 "这个镜头好不好"需要一个客观答案。MTF 提供了它——把成像系统当作空间频率的滤波器(承第四页 4f 系统)。本页先讲清 MTF 这套语言,再讲计算成像如何模糊了"光学"与"算法"的边界,以及它的真实边界在哪。

学习层:清晰度不是一个旋钮,而是一条频率账

1. 具体谜题:更密的像素一定让细节更清楚吗?

一枚像素间距 \(p=5\,\mu\mathrm{m}\) 的传感器观察 \(80\,\mathrm{lp/mm}\) 的黑白条纹。镜头 MTF 约为 \(0.67\),像素孔径 MTF 约为 \(0.76\)。如果把镜头换成 MTF 更高的版本,系统一定还能分辨更细的条纹吗?再把频率提高到 \(120\,\mathrm{lp/mm}\),你看到的是“更弱的真实细节”,还是已经进入混叠区的伪纹理?

2. 先预测:先写下三个判断

在打开实验结果前,先判断:

  1. \(p=5\,\mu\mathrm{m}\) 时奈奎斯特频率是多少?\(80\,\mathrm{lp/mm}\) 是否已经混叠?
  2. 若镜头、像素、离焦三个 MTF 分别为 \(0.67,0.76,0.87\),系统 MTF 会不会高于其中任何一个?
  3. 把 \(120\,\mathrm{lp/mm}\) 的输入换成一个低频条纹,是否就能“恢复”原来没有采到的高频相位?

3. 最小模型:卷积把每一环变成乘法

点扩散函数的傅里叶变换是 OTF,线性、平移不变的级联系统满足

\[ \mathrm{OTF}_{sys}(f)=\mathrm{OTF}_{lens}(f)\,\mathrm{OTF}_{pixel}(f)\,\mathrm{OTF}_{blur}(f),\qquad \mathrm{MTF}=|\mathrm{OTF}|. \]

这里的 OTF/MTF 级联采用非相干照明下的线性强度成像模型。若讨论相干成像,线性传递的是复振幅,瞳函数的场频谱支持上限是 \(\mathrm{NA}/\lambda\);探测到的强度还要取 \(|E|^2\),不能把它直接当成同一个线性强度 OTF。还要分清三条不同的界线:非相干强度 OTF 的衍射截止 \(f_c=2\mathrm{NA}/\lambda\)、有限像素孔径的首个零点 \(1/p\),以及采样奈奎斯特频率 \(f_N=1/(2p)\)。

有限宽度像素在这个强度模型中对光强做空间积分,因此

\[ \mathrm{MTF}_{pixel}(f)=\left|\operatorname{sinc}(\pi f p)\right|,\qquad f_N=\frac1{2p}. \]

这给出两个不变量:级联的 MTF 不会凭空超过每一环;\(f>f_N\) 时采样会把频率折回低频,而不是把原始波形无损搬进数字图像。实验中的 lensQuality 和离焦尺度是高斯预采样 MTF 代理,不输入 \(\mathrm{NA}\) 与 \(\lambda\),所以它不模拟硬的衍射截止;超过奈奎斯特的 120 lp/mm 例子只是在说明采样混叠,不是在声称高于光学截止仍有可靠信息。

4. 实验:在乘法账本上先押注,再看曲线

先选择你对当前条纹频率的判断,核对后再拖动像素间距、镜头质量和频率。SVG 同时显示镜头、像素孔径、离焦和系统 MTF,并把奈奎斯特线与混叠频率标出来。

无 JavaScript 时的静态读法:默认 \(p=5\,\mu\mathrm{m}\)、\(f=80\,\mathrm{lp/mm}\);示例模型给出镜头 MTF \(0.67\)、像素 MTF \(0.76\)、离焦 MTF \(0.87\),所以系统预采样 MTF 约为 \(0.44\),奈奎斯特频率为 \(100\,\mathrm{lp/mm}\)。

输入频率 镜头 MTF 像素 MTF 系统预采样 MTF 采样判断
\(40\,\mathrm{lp/mm}\) 0.91 0.94 0.82 低于奈奎斯特,尚未混叠
\(80\,\mathrm{lp/mm}\) 0.67 0.76 0.44 低于奈奎斯特,但已被级联衰减
\(120\,\mathrm{lp/mm}\) 0.41 0.50 0.15 超过奈奎斯特,会折叠成伪低频

表中 120 lp/mm 的非零值是预采样响应,不代表采样后能可靠恢复该频率;若要计算具体光学截止,还必须给出与波长同单位的 NA 和 \(\lambda\)。

5. 失败边界:MTF 高不等于信息可恢复

  • MTF 只记录幅度传递;PTF 的相位反转仍可能把条纹位置和形状改掉。
  • 反卷积只能放大已经进入系统且有足够 SNR 的频率;超过截止频率或被噪声淹没的成分没有可逆账本。
  • AI 超分辨可用先验生成“像真的”高频,但不能把生成细节冒充为测量证据。医学、司法和科学测量必须区分测得的频率与猜出的纹理。

6. 迁移任务:为一条真实规格写出验收条件

给定镜头 MTF 曲线、像素间距、曝光抖动和目标条纹频率,写出一份三行验收表:系统 MTF、奈奎斯特频率、以及是否需要 OLPF 或更高光子数。然后说明在“检测缺陷”与“生成好看照片”两种任务中,为什么同一个 MTF 结论不能直接替代 SNR 和误报率。

一、MTF:成像系统的频率响应

承第四页:成像是线性滤波,\(\tilde{I} = \tilde{O}\cdot\mathrm{OTF}\)。

读 MTF 曲线的方法:

MTF:衍射极限、像差与探测器

图 11-1 系统 MTF 的构成。衍射极限(虚线)给出理论上限,在截止频率 \(f_c=2\mathrm{NA}/\lambda\) 处降为零;实际镜头(蓝)因像差而更低;探测器因像素有限尺寸而有自身 MTF(其零点在 \(1/p\));系统总 MTF 是各环节相乘(紫)。注意奈奎斯特频率 \(1/2p\)——高于它的信息会混叠成伪信号。系统性能由最差的环节主导,这就是"镜头与传感器必须匹配"的定量含义。

级联规则:

\[\mathrm{MTF}_{\text{系统}} = \mathrm{MTF}_{\text{镜头}}\times\mathrm{MTF}_{\text{探测器}}\times\mathrm{MTF}_{\text{离焦/抖动}}\times\cdots\]

这个乘法关系给出一条重要的工程判断:任何一环很差,整体就好不了。 给一个平庸的镜头配一个高像素传感器,或反之,都是浪费——匹配比单项优秀更重要(承本课导论 §3)。

探测器的 MTF:像素是有限大小的方形,对光强做空间积分 → 相当于一个 sinc 型低通,首个零点在 \(1/p\)(\(p\) 为像素间距)。

混叠:奈奎斯特频率 \(f_N = 1/(2p)\)。高于它的空间频率会折叠成低频伪信号(摩尔纹、伪色)。 对策:光学低通滤波器(OLPF,故意轻微模糊)——但现代高像素传感器多已取消它,因为像素已足够密,且宁可保留锐度、用算法处理偶发摩尔纹。这是一个典型的"从光学解决转向算法解决"的迁移。

二、其他评价指标

三、计算成像:重新划分光学与算法的边界

传统思路:光学系统尽力做到完美 → 探测器忠实记录 → 图像即最终产品。

计算成像的思路:光学与算法联合设计——光学系统不必产生"好看"的图像,只需编码足够的信息,由算法解码出所需结果。

代表性技术:

技术 思想
相位编码/波前编码 故意引入特定像差,使 PSF 对离焦不敏感 → 算法反卷积后获得极大景深
光场相机 微透镜阵列记录光线的方向信息 → 可事后重聚焦(代价:空间分辨率大幅下降)
压缩感知成像 利用信号稀疏性,用远少于奈奎斯特的测量重建(单像素相机)
编码孔径 用特定图案孔径,可做深度估计或 X 射线/伽马成像(无法用透镜时)
多帧/多曝光合成 HDR、超分辨、去噪(承第十页)
无透镜成像 传感器紧贴掩模,全靠计算重建 → 极薄,用于生物芯片

成功的关键在于"联合设计":单独看,编码光学产生的原始图像往往很糟;但光学 + 算法作为整体,性能可超越任何一方单独优化的结果。

这与 🔗 微电子站第九页"把模拟难题转成数字问题"、第二页"畸变交给软件"是同一条主线——把约束从难的域转移到便宜的域。

四、诚实的边界:算法能做什么、不能做什么

这是本页最需要清醒的部分(承第四页):

能做:

不能做:

关于"AI 超分辨":生成模型能产生看起来清晰的图像,但那些细节是根据训练先验推测的,不是测量到的。

这在不同场合的性质截然不同:

判读此类工作的标准问题:这些高频信息来自测量还是来自先验?有没有在"真值已知"的独立数据上验证过? 这与 🔗 生物站、语言站反复强调的判读纪律完全一致。

五、系统工程视角

设计一个成像系统时,需要联合考虑:

\[\text{光源/照明} \to \text{光学} \to \text{探测器} \to \text{读出电子学} \to \text{算法}\]

常见的失衡(都是实际项目中反复出现的):

综合指标:单看 MTF 或单看 SNR 都不够,可用噪声等效调制(NEM)或可探测量子效率(DQE)等把两者结合(医学影像领域尤其常用 DQE)。

六、要点


线三结束。下一线转向光的另一大用途——传输信息。光纤是人类建造过的最重要的信息基础设施。