衬度传递函数 (CTF)

电镜不是忠实成像 —— 它把每个频率按一条振荡曲线"打折",还会反转衬度

直觉

电镜并不忠实成像。受离焦和球差影响,显微镜对不同空间频率的传递强度并不一致 —— 有的频率被增强,有的被压到零(彻底丢失),还有的衬度被整个反转(本该黑的地方变白)。描述这条传递曲线的,就是衬度传递函数 (CTF)

换个说法:把一张图拆成各种粗细的条纹(空间频率),电镜对每一档条纹都给一个介于 1-1+1+1 之间的增益。+1+1 表示这一档原样传过来,00 表示这一档被抹掉,1-1 表示这一档黑白互换地传过来。CTF 就是”增益随条纹粗细变化”的那条曲线,而它不是平的、也不是单调的,而是来回振荡、反复变号。

为什么显微镜会这样?冰中的生物分子几乎不吸收电子,主要是给穿过的电子波叠加相位(相位物体)。但探测器只能记录强度,记不下相位。要把看不见的相位差变成看得见的明暗,靠的是故意把镜头离焦:离焦让不同频率的散射波在到达探测器时积累不同的相位延迟,这些延迟有的让波相长(亮)、有的相消(暗),从而把相位信息”翻译”成衬度。代价是这个翻译过程对每个频率的效果都不同,于是就有了那条振荡曲线。

CTF 随空间频率的振荡如下,欠焦量与振幅衬度可调:

第一个零点空间频率 k (1/Å)+1−1
CTF(k)衬度反转区 (CTF < 0)

CTF 在零点处把对应频率完全丢失,在负值区把衬度反转。欠焦越大,振荡越快、第一个零点越靠低频——这就是为什么单一欠焦无法覆盖所有频率,需要跨欠焦合并。

紫色阴影区是 CTF < 0、衬度被反转的频段;黄色虚线是第一个零点 —— 在它之后曲线开始反复穿过零,对应频率被反复丢失。试着把欠焦调大:第一个零点会向低频移动,曲线整体变密 —— 这意味着大欠焦能给低频(粗大特征)很强的衬度,却更早开始丢高频(精细细节)。把欠焦调小则相反:低频衬度弱,但第一个零点推到更高频,能保住更多细节。这正是采图时要在”看得见”和”看得清”之间做的取舍。

数学形式

深入

CTF 由波像差函数 χ(k)\chi(k) 决定。对空间频率 kk、欠焦 Δf\Delta f、球差 CsC_s、电子波长 λ\lambda

χ(k)=πλΔfk2    π2Csλ3k4\chi(k) = \pi\,\lambda\,\Delta f\,k^2 \;-\; \frac{\pi}{2}\,C_s\,\lambda^3\,k^4

逐项读这个相位:χ(k)\chi(k) 是频率为 kk 的散射波相对未散射波多走的相位(弧度)。kk 是空间频率(单位 1/A˚1/\text{Å},越大代表越精细的特征);Δf\Delta f 是欠焦量(习惯上欠焦取正);CsC_s 是球差系数,刻画镜头对大角度散射波的额外相位误差;λ\lambda 是电子波长,由加速电压决定。第一项 πλΔfk2\pi\lambda\,\Delta f\,k^2 来自离焦,随 k2k^2 增长,在中低频就主导整条曲线;第二项 π2Csλ3k4\tfrac{\pi}{2}C_s\lambda^3 k^4 来自球差,随 k4k^4 增长,只有到高频才显著。

加入振幅衬度 ww 后:

CTF(k)=[1w2sinχ(k)+wcosχ(k)]\mathrm{CTF}(k) = -\Big[\sqrt{1-w^2}\,\sin\chi(k) + w\,\cos\chi(k)\Big]

这里 sinχ\sin\chi 项是相位衬度:相位 χ\chi 每增加 π\pisinχ\sin\chi 就翻一次号,这就是曲线反复变号的来源。ww(典型约 0.070.10.07\text{–}0.1)是振幅衬度占比,对应真正被样品吸收/散射掉、不靠离焦也能成像的那一小部分;它通过 cosχ\cos\chi 项在零频附近给一点非零衬度,所以 k0k\to 0CTFw\mathrm{CTF}\to -w 而非 00。前面的负号是符号约定,让欠焦下低频衬度为负、与常见显示一致。

上面的演示用 300 kV(λ0.0197\lambda \approx 0.0197 Å)、Cs=2.7C_s = 2.7 mm。k2k^2 项由欠焦主导、随频率二次增长,k4k^4 项由球差主导、在高频才显现 —— 两者叠加造就了那条越往高频越密的振荡。

一个有用的具体数字:相位每累积 π\pi 就出现一个新零点。由 sinχ=0\sin\chi=0χ(kn)=nπ\chi(k_n)=n\pi,在只保留离焦项时第一个零点约在 k11/(λΔf)k_1\approx\sqrt{1/(\lambda\,\Delta f)}。代入 Δf=1 μm=104 A˚\Delta f=1\ \mu\text{m}=10^4\ \text{Å}λ0.0197 A˚\lambda\approx 0.0197\ \text{Å},得 k11/14 A˚1k_1\approx 1/14\ \text{Å}^{-1},即首个零点落在约 14 A˚14\ \text{Å} 的特征尺度上。把欠焦增到 2 μm2\ \mu\text{m}k1k_1 缩小 2\sqrt 2 倍,第一个零点退到约 20 A˚20\ \text{Å} —— 这就量化了”大欠焦更早丢细节”。

为了把相位翻译成衬度,离焦其实带来了双重代价。其一是上面的振荡:从某个频率往后,sinχ\sin\chi 一会儿正一会儿负,于是同一张图里粗特征和细特征的明暗关系会对不上。其二是探测器只测得到强度 2|\,\cdot\,|^2,丢掉了波的相位,CTF 正是这次”相位→强度”翻译留下的逐频率指纹。理解这两点,就能明白为什么后面非做校正不可。

为什么要在意

还有一个实际细节:真实样品往往有厚度、还可能轻微像散,于是同一张图里不同位置、不同方向的有效欠焦并不一致。处理流程要先估计每张图(乃至每个区域)的欠焦与像散参数,把那条 CTF 曲线拟合出来,才能做校正 —— 这一步通常用图像功率谱里的同心环(Thon 环)来标定,环的疏密和椭圆度直接读出欠焦和像散。

CTF 和缺失楔形是 Cryo-ET 成像的两块基石:一个在频率方向”打折”,一个在角度方向”挖空”。它和采样与奈奎斯特共同决定了一张图里哪些频率”可信”:奈奎斯特设上限、CTF 在上限以内逐频率加权。CTF 的逐频率衰减也是图像信噪比随频率下降的一大来源。重构前的 CTF 校正属于滤波的范畴 —— 在频域施加一个逐频率权重。而 CryoGEN 这类生成式方法不把校正当成孤立的预处理,而是把 CTF(频率方向的退化)和缺失楔形(角度方向的退化)一起写进成像模型,让复原与重构在同一个框架里同时应对这两种退化。

← 信号处理