中心切片定理
投影的傅里叶变换等于物体频谱中过原点的一个切片,把层析倾转角与频域覆盖联系起来
中心切片定理,又称傅里叶切片定理,把物体的投影与其傅里叶变换联系起来。它指出:二维物体平行投影的一维傅里叶变换,等于物体二维变换中过原点、且与投影方向垂直的一条直线。三维版本把直线换成平面:一张投影图像的二维变换,是物体三维频谱中过原点的一个中心切片。这条定理把实空间里的一次求和,等价成频域中沿一条直线或一片平面的采样,从而为层析采集与重构提供了统一的几何图景。
这件事为什么值得专门立一条定理?因为它把一个看似无望的反问题变得可解。从一堆投影还原内部三维结构,在实空间里没有直接的下手处;但只要承认每张投影都对应频域里一片确定的切片,重构就变成一道”拼图”:把各角度的切片摆进同一个三维频谱,填满之后做一次逆变换即可。难点也随之换了位置 —— 不再是”怎么反演”,而是”切片够不够、摆得准不准”。
每张投影在二维频谱中沿一条过原点的直线采样(中心切片定理)。倾角范围为 ±60° 时,竖直轴两侧各留下半角 30° 的楔形从未被测量;当最大倾角趋近 90° 时楔形闭合。
每个倾转角所对应的切片都过频谱原点,因此低频信息被所有投影反复采样,而高频细节只沿各自的切片方向被记录。这条规律有两个直接后果:低频(物体的整体形状、大尺度衬度)总是冗余而稳健,几乎不依赖采集多少角度;高频(边缘、细小颗粒、原子级特征)则脆弱得多,只有恰好落在某条切片方向上的那部分才被测到,旋转 90° 看同一个特征,得到的信息可能完全不同。
形式上,设 为 沿方向 的投影,则
这里 是待求的物体(密度), 是在倾角 下沿一条直线把物体加起来得到的一维投影, 是这条投影上的坐标。等号左边 是这条投影的一维傅里叶变换, 是沿投影方向的空间频率。等号右边 是物体的二维傅里叶变换, 是频域里一条过原点、与投影方向成角 的直线上的点。换句话说:把一维投影做变换,得到的恰好是物体二维频谱沿某条过原点直线上的取值,二者一一对应、互不损失。因此每个投影沿一条过原点的直线(二维)或平面(三维)采样物体的频谱。在许多角度上采集投影即可填满频域,随后一次精确逆变换便能还原物体。这正是从倾转序列做层析重构的原理。
投影沿电子束方向对物体求和,从而丢弃了该方向上的一切变化。用频率的语言说,沿某轴求和只保留了与该轴正交的那片频谱切片。因此每个倾转角贡献一片新切片,角度越多,频域被填得越满。
可以这样记:求和就是只取频率为零的那一档。沿某个方向求和,等于把那个方向上所有非零频率清零,只剩下与该方向垂直的频率分量幸存下来 —— 这些幸存者拼起来,正好是过原点、与求和方向垂直的那片切片。
定理为什么成立,一行就能看清。投影沿 方向求和:。对它做一维傅里叶变换,
关键的一步是把对 的内层求和等同于在二维变换里令 :因为 ,沿 积分恰好就是二维变换在 上的取值。于是沿 投影对应频谱中 这条线,即与投影方向垂直、过原点的切片。把坐标系旋转角 ,这条线就转成 ,得到上面的一般式。傅里叶变换的旋转不变性(旋转物体等于同样旋转其频谱)保证了这一步合法。
由此能读出几个量化事实。其一,每片切片的角宽随频率反比收窄:在半径 处,相邻两个倾角 之间留下的弧长约为 ,所以高频处切片之间的空隙更大,采样更稀 —— 这就是为什么仅靠加密角度也难以恢复最高频细节。其二,离散采集下重构常用滤波反投影:逆变换里的极坐标雅可比带来一个 因子,等价于对每张投影先乘一个斜坡滤波器(高频增益正比于 )再反投影,用以补偿低频被过采样、抵消其重复计数。其三,相位远比幅度重要:物体的边缘位置由各频率的相位对齐决定,切片采样误差里相位的错配比幅度的衰减更致命。
这条定理给 Cryo-ET 带来一个直接而无法回避的后果。样品只能在有限范围内倾转,通常约 到 ,再往上,有效厚度和载网几何就让更高的倾角无法使用了。可用角度贡献的切片,在频域中留下一块楔形区域始终没被采样到 —— 这就是缺失楔形。一个具体的数字感受:若最大倾角为 ,则在垂直于电子束的平面内,两侧各有 、合计约 的扇区从未落进任何一片切片,约占该平面频率内容的三分之一完全缺测。这块楔形里的信息从未被测到,于是分辨率变得各向异性:沿电子束轴被拉伸,特征也出现种种典型的畸变(球状物拉成椭球、平行于缺口的边缘被抹淡)。由于这一缺口定义在频域中,单凭实空间里的任何操作都填不回来 —— 你能内插、能加权,却无法凭空造出从未采到的频率。如何应对它,正是 CryoGEN 等方法所要面对的核心问题:它们不去”填洞”,而是用一个学到的先验,去推断哪一种完整结构在缺测频率之外、又与所有已测切片一致。