卷积与点扩散函数

任何线性平移不变成像系统都通过把物体与点扩散函数卷积来模糊,它在频域中即为传递函数

卷积描述了线性、平移不变系统如何变换输入。对一维信号 ff 与核 hh

(fh)(x)=f(x)h(xx)dx,(f * h)(x) = \int_{-\infty}^{\infty} f(x')\,h(x - x')\,dx',

并可自然推广到二维和三维。若系统是线性的(输出可叠加)且平移不变的(响应不随位置改变),则其全部行为都由它对单个点源的响应所刻画。这个响应就是点扩散函数 (PSF),而系统的输出即为物体与 PSF 的卷积。

逐项来读这个积分:xx 是我们要计算输出的位置;xx' 是扫过整个物体的求和变量;f(x)f(x') 是物体在 xx' 处的强度;h(xx)h(x-x') 是该点对位置 xx 的贡献权重,只取决于两者之差 xxx-x'——正是这一点体现了「平移不变」。核中的减号意味着 hh 在叠加前被翻转,这是卷积区别于相关的标志,对对称的 PSF 而言并无可见差别,但在一般情形下不可忽略。

直觉

把物体想象成一片由许多发光点组成的星空,而仪器无法把任何一个点成像为真正的点:每个点都被画成一小团固定形状的光斑,这团光斑就是 PSF。整幅图像,就是把每个点对应的、按其亮度缩放的光斑全部叠加起来。物体本身的信息没有丢失,只是被「这一团」反复盖章涂抹了。复原图像,就是反过来问:哪一片真实的星点,盖上这枚印章后,恰好得到我所看到的图像。

位置
输入(两个尖峰)高斯 PSF 核卷积输出(模糊后)

σ = 6.0 像素:两个尖峰被各自的 PSF 副本替换后相加。σ 越大模糊越强,相邻峰逐渐并为一团。由卷积定理,实空间卷积等于频域相乘,宽 PSF 的传递函数压低高频,细节随之损失。

这样一来,PSF 就像一枚理想化的模糊印章:物体上的每个点都被换成一份按比例缩放的 PSF 副本,再把所有副本叠加起来。窄的 PSF 保留精细细节,宽的则将其抹开。Cryo-ET 的成像在很好的近似下正是这样一个系统,故记录到的图像是投影势与仪器 PSF 的卷积。须留意,电镜的这枚印章并非温和的低通模糊:它的衬度传递函数随频率振荡并反复变号,在某些频带上将衬度整体翻转,因此其 PSF 带有振铃式的正负旁瓣,而非单调展宽的高斯钟形。

一个常用的解析模型是高斯 PSF,其宽度由标准差 σ\sigma 控制;σ\sigma 越大,两个相邻的尖锐特征越易被抹成一团而难以分辨。两个相距固定的点源经卷积后能否仍被分开,正是分辨率的直观体现。

把这件事算成数字会更清楚。设 PSF 为高斯 h(x)=exp(x2/2σ2)h(x)=\exp(-x^2/2\sigma^2),两个等强度的点源相距 dd。卷积后图像是两个中心相隔 dd 的高斯之和。当 dd 远大于 σ\sigma 时,两峰之间有明显的低谷,肉眼可分;当 dd 缩小到约 2σ2\sigma 时,中央的低谷被填平,两峰并成一个鼓包,无法再判定其为「两个」。于是 σ\sigma 就直接读作可分辨距离的尺度:σ\sigma22 个像素,约 44 像素以内的特征便混作一团。频域中这对应一个同样的高斯传递函数 H(k)=exp(2π2σ2k2)H(k)=\exp(-2\pi^2\sigma^2 k^2)——实空间越宽(σ\sigma 大),频域越窄,被压制的高频越多,正是细节流失的来由。

深入

由卷积定理,实空间中的卷积即频域中的相乘:

F{fh}(k)=F(k)H(k).\mathcal{F}\{f * h\}(k) = F(k)\,H(k).

这里 F(k)F(k)H(k)H(k) 分别是 ffhh 的傅里叶变换,kk 是空间频率。PSF 的变换 H(k)H(k) 是系统的传递函数,它说明每个空间频率被传递的强弱。一个有用的特例是 k=0k=0H(0)=h(x)dxH(0)=\int h(x)\,dx 是 PSF 的总「质量」,若它等于 11,卷积便不改变图像的整体亮度,只重新分配细节。对显微镜而言,这条传递函数就是衬度传递函数;它的振荡与零点,正是单凭一份实空间 PSF 难以直接读出的逐频率加权。在频域工作还使反卷积在概念上变得简单——除以 H(k)H(k)——但 HH 的零点使精确反演不可能:在 H(k)=0H(k)=0 的频率上,F(k)=Y(k)/H(k)F(k)=Y(k)/H(k) 成了 0/00/0,该频率的信息被仪器彻底丢弃,无从恢复;而在 H(k)H(k) 很小处,除法把那一带的噪声成倍放大,朴素的反演因此不切实际。

由此带来两点后果。复原图像,本质上就是要撤销一次卷积,这是一个反卷积问题,在传递函数很小或为零的地方并不适定。又因为卷积对应频域中的相乘,在频域里施加一个校正权重——锐化、平滑或 CTF 校正——回到实空间本身也就是一次卷积。这就把模糊、仪器响应,以及重构中处处用到的滤波,统一到了同一个线性系统的框架之下。

这条思路贯穿整个 Cryo-ET 的重构流程。前向成像被建模为「物体 * PSF」,于是重构就是反卷积;正因它在 CTF 零点附近不适定且对噪声敏感,单帧除以 H(k)H(k) 几乎不可用,实践中要么对多张焦距不同的图像合并以填补各自的零点,要么把 PSF 作为已知的前向算子嵌入迭代或生成式求解器,让先验信息去弥补仪器丢失的频段。换言之,理解卷积与 PSF,就是理解重构究竟在「撤销」什么、以及为什么必须借助先验才能撤销得动。

← 信号处理