中心切片定理

投影的傅里叶变换等于物体频谱中过原点的一个切片,把层析倾转角与频域覆盖联系起来

中心切片定理,又称傅里叶切片定理,把物体的投影与其傅里叶变换联系起来。它指出:二维物体平行投影的一维傅里叶变换,等于物体二维变换中过原点、且与投影方向垂直的一条直线。三维版本把直线换成平面:一张投影图像的二维变换,是物体三维频谱中过原点的一个中心切片。这条定理把实空间里的一次求和,等价成频域中沿一条直线或一片平面的采样,从而为层析采集与重构提供了统一的几何图景。

这件事为什么值得专门立一条定理?因为它把一个看似无望的反问题变得可解。从一堆投影还原内部三维结构,在实空间里没有直接的下手处;但只要承认每张投影都对应频域里一片确定的切片,重构就变成一道”拼图”:把各角度的切片摆进同一个三维频谱,填满之后做一次逆变换即可。难点也随之换了位置 —— 不再是”怎么反演”,而是”切片够不够、摆得准不准”。

已采样切片缺失楔形 (半角 30°)

每张投影在二维频谱中沿一条过原点的直线采样(中心切片定理)。倾角范围为 ±60° 时,竖直轴两侧各留下半角 30° 的楔形从未被测量;当最大倾角趋近 90° 时楔形闭合。

每个倾转角所对应的切片都过频谱原点,因此低频信息被所有投影反复采样,而高频细节只沿各自的切片方向被记录。这条规律有两个直接后果:低频(物体的整体形状、大尺度衬度)总是冗余而稳健,几乎不依赖采集多少角度;高频(边缘、细小颗粒、原子级特征)则脆弱得多,只有恰好落在某条切片方向上的那部分才被测到,旋转 90° 看同一个特征,得到的信息可能完全不同。

实空间:物体与一次投影投影1D 傅里叶变换傅里叶空间:过原点的切片中心切片

形式上,设 pθ(s)p_\theta(s)f(x)f(\mathbf{x}) 沿方向 θ\theta 的投影,则

pθ^(k)=F(kcosθ,  ksinθ),\widehat{p_\theta}(k) = F\big(k\cos\theta,\; k\sin\theta\big),

这里 f(x)f(\mathbf{x}) 是待求的物体(密度),pθ(s)p_\theta(s) 是在倾角 θ\theta 下沿一条直线把物体加起来得到的一维投影,ss 是这条投影上的坐标。等号左边 pθ^(k)\widehat{p_\theta}(k) 是这条投影的一维傅里叶变换,kk 是沿投影方向的空间频率。等号右边 FF 是物体的二维傅里叶变换,(kcosθ,  ksinθ)\big(k\cos\theta,\;k\sin\theta\big) 是频域里一条过原点、与投影方向成角 θ\theta 的直线上的点。换句话说:把一维投影做变换,得到的恰好是物体二维频谱沿某条过原点直线上的取值,二者一一对应、互不损失。因此每个投影沿一条过原点的直线(二维)或平面(三维)采样物体的频谱。在许多角度上采集投影即可填满频域,随后一次精确逆变换便能还原物体。这正是从倾转序列做层析重构的原理。

直觉

投影沿电子束方向对物体求和,从而丢弃了该方向上的一切变化。用频率的语言说,沿某轴求和只保留了与该轴正交的那片频谱切片。因此每个倾转角贡献一片新切片,角度越多,频域被填得越满。

可以这样记:求和就是只取频率为零的那一档。沿某个方向求和,等于把那个方向上所有非零频率清零,只剩下与该方向垂直的频率分量幸存下来 —— 这些幸存者拼起来,正好是过原点、与求和方向垂直的那片切片。

深入

定理为什么成立,一行就能看清。投影沿 yy 方向求和:p(x)=f(x,y)dyp(x)=\int f(x,y)\,dy。对它做一维傅里叶变换,

p^(kx)= ⁣ ⁣f(x,y)e2πikxxdydx= ⁣ ⁣f(x,y)e2πi(kxx+0y)dxdy=F(kx,0).\widehat{p}(k_x)=\int\!\!\int f(x,y)\,e^{-2\pi i k_x x}\,dy\,dx = \int\!\!\int f(x,y)\,e^{-2\pi i (k_x x + 0\cdot y)}\,dx\,dy = F(k_x, 0).

关键的一步是把对 yy 的内层求和等同于在二维变换里令 ky=0k_y=0:因为 e2πi0y=1e^{-2\pi i\cdot 0\cdot y}=1,沿 yy 积分恰好就是二维变换在 ky=0k_y=0 上的取值。于是沿 yy 投影对应频谱中 ky=0k_y=0 这条线,即与投影方向垂直、过原点的切片。把坐标系旋转角 θ\theta,这条线就转成 (kcosθ,  ksinθ)\big(k\cos\theta,\;k\sin\theta\big),得到上面的一般式。傅里叶变换的旋转不变性(旋转物体等于同样旋转其频谱)保证了这一步合法。

由此能读出几个量化事实。其一,每片切片的角宽随频率反比收窄:在半径 kk 处,相邻两个倾角 Δθ\Delta\theta 之间留下的弧长约为 kΔθk\,\Delta\theta,所以高频处切片之间的空隙更大,采样更稀 —— 这就是为什么仅靠加密角度也难以恢复最高频细节。其二,离散采集下重构常用滤波反投影:逆变换里的极坐标雅可比带来一个 k|k| 因子,等价于对每张投影先乘一个斜坡滤波器(高频增益正比于 k|k|)再反投影,用以补偿低频被过采样、抵消其重复计数。其三,相位远比幅度重要:物体的边缘位置由各频率的相位对齐决定,切片采样误差里相位的错配比幅度的衰减更致命。

这条定理给 Cryo-ET 带来一个直接而无法回避的后果。样品只能在有限范围内倾转,通常约 ±60\pm 60^\circ±70\pm 70^\circ,再往上,有效厚度和载网几何就让更高的倾角无法使用了。可用角度贡献的切片,在频域中留下一块楔形区域始终没被采样到 —— 这就是缺失楔形。一个具体的数字感受:若最大倾角为 ±60\pm 60^\circ,则在垂直于电子束的平面内,两侧各有 3030^\circ、合计约 6060^\circ 的扇区从未落进任何一片切片,约占该平面频率内容的三分之一完全缺测。这块楔形里的信息从未被测到,于是分辨率变得各向异性:沿电子束轴被拉伸,特征也出现种种典型的畸变(球状物拉成椭球、平行于缺口的边缘被抹淡)。由于这一缺口定义在频域中,单凭实空间里的任何操作都填不回来 —— 你能内插、能加权,却无法凭空造出从未采到的频率。如何应对它,正是 CryoGEN 等方法所要面对的核心问题:它们不去”填洞”,而是用一个学到的先验,去推断哪一种完整结构在缺测频率之外、又与所有已测切片一致。

← 信号处理