Open Cryo-ET 试运行
中文 EN

缺失楔形

为什么 Cryo-ET 的断层重构总在竖直方向被拉长？

直觉

为了重建一个三维结构，你把样品在电镜里倾转，每个角度拍一张二维投影。但样品台不能转到 ±90°——倾角 θ\theta 越大，电子束穿过冰层平板的路径越长，按 1/cos⁡θ1/\cos\theta 增厚（60∘60^\circ 时已是两倍），信号被逐渐衰减（并非某个角度上的硬性截断），直到电子再也无法有效穿透。实际上一般只能转到 ±60° 左右。于是有一整片角度你从来没拍到过。

缺失楔形不是一个能调好的参数，而是几何决定的硬约束：它内部的傅里叶系数从未进入过任何一张投影，因此重构里关于它们的一切都是推断出来的，不是测出来的。这是 Cryo-ET 与普通显微成像最根本的区别。

这片”没拍到”的角度在傅里叶空间中的形状如下，倾转角可调：

倾转方向
傅里叶空间
已采样缺失楔形

倾转到 ±60°，傅里叶空间约有 33% 的角度从未被采样。这片缺口就是「缺失楔形」，它让重构在竖直方向上被拉长、模糊。

紫色是被采样到的区域，深色的两块楔形就是缺失楔形 (missing wedge)。

为什么是傅里叶空间里的一片楔形

深入

关键是中心切片定理 (central-slice theorem)：一张二维投影的傅里叶变换，等于三维物体傅里叶变换中、过原点且垂直于投影方向的那一个二维切片。

F2D[Pθ(x,y)]  =  f^(plane⊥θ)\mathcal{F}_{2D}\big[ P_\theta(x,y) \big] \;=\; \widehat{f}\big(\text{plane} \perp \theta\big)

逐项读这个式子：左边 Pθ(x,y)P_\theta(x,y) 是倾角 θ\theta 下拍到的那张二维投影，F2D\mathcal{F}_{2D} 是对它做二维傅里叶变换；右边 f^\widehat{f} 是三维物体 ff 的傅里叶变换，plane⊥θ\text{plane}\perp\theta 指过原点、法向沿投影方向 θ\theta 的那个二维平面。等号的意思是：拍一张投影，等于在三维频谱里抄下一个过原点的切片——而切片之外的频率，这张图什么都没告诉你。

每个倾转角 θ\theta 贡献一个过原点的切片。把 θ\theta 从 −θmax⁡-\theta_{\max} 扫到 +θmax⁡+\theta_{\max}，这些切片扫出的是一个双锥，留下的空缺是绕电子束轴的两块楔形。缺失楔形的半张角是

α=90∘−θmax⁡.\alpha = 90^\circ - \theta_{\max}.

这里 θmax⁡\theta_{\max} 是能达到的最大倾角，α\alpha 是单侧楔形从频率空间赤道面量起的张角。当 θmax⁡=60∘\theta_{\max}=60^\circ 时 α=30∘\alpha=30^\circ，约 1/3 的角度从未被采样；即便把样品台推到 θmax⁡=70∘\theta_{\max}=70^\circ，α\alpha 也只降到 20∘20^\circ，楔形仍在。要彻底消除它需要 θmax⁡=90∘\theta_{\max}=90^\circ，而这在平板冰样品上做不到。

直觉

为什么是”楔形”而不是别的形状？因为缺的不是某些孤立的频率点，而是一整段方向——所有近乎竖直（沿电子束轴）的频率分量都没被采样。频率的方向对应实空间里结构的取向：竖直方向的频率编码的是水平界面（如平躺的膜）上下变化的细节。缺了这段方向，凡是”上下”方向的边界都会变虚。换句话说，缺失楔形不是均匀地降低分辨率，而是只在某些取向上降低——这正是它伪影如此有特征性的原因。

后果

傅里叶空间竖直方向缺信息，等价于实空间里点扩散函数沿 Z 轴被拉长。可以这样理解：重构等于真值与一个”模糊核”卷积，而这个模糊核的形状正是缺失楔形的傅里叶逆变换——它在 Z 方向被抹开成一条尾巴。直接的可见后果：

一个常用的缓解手段是双轴断层 (dual-axis tomography)：采完一条倾转序列后把样品在平面内转 90∘90^\circ，再采一条。两条序列的缺失楔形方向正交，叠加后空缺缩小为一个缺失锥 (missing pyramid)，伪影减轻但并未消除，且代价是双倍剂量。

这些都不是噪声，而是确定性的信息缺失：缺失楔形里的频率从未被测量，任何线性滤波都无法把它们造出来。这正是 CryoGEN 这类方法试图用学到的先验和自监督学习去”填补”的东西——把缺失楔形当成一个待补全的逆问题，靠对真实生物结构的统计先验去推断那些从未拍到的频率。在管线中，缺失楔形破坏的恰是 重构 阶段输出的那个观测体 yy，而 倾转序列采集 的几何上限正是它的根源。

← 电子断层成像