正电子断层扫描--PET的图像重建3.docVIP

中国医学物理学杂志 正电子断层扫描--PET的图像重建 刘力1 赵书俊2 （1中国科学技术大学近代物理系，安徽，合肥 230027；2协和医科大学阜外医院，北京100037） 摘要：随着计算机技术的迅速发展，PET图像重建方法得以很大的改进，其研究工作越来越受到人们的重视。特别是对能压抑噪声和提高空间分辨率的图像快速重建算法的研究是目前一个热门课题。本文简介了PET图像重建的基本原理与基本方法及其研究与应用现状。此外，对PET的三维重建也做了一些简介。 关键词：正电子断层扫描；图像重建；解析法；迭代法 中图分类号：R445.2, TP751 文章编号：1005202*（2000） 1. 引言 PET（positron emission tomography,正电子断层扫描仪）已广泛用于人体内部血流或代谢的三维显像，其空间分辨率可达到3~5mm。为简述PET图像重建的原理，考虑一个二维切面的成像，如图1所示。PET利用其封闭环绕型探测器阵列对正电子湮灭产生的背对背光子进行符合测量，符合测量中的投影线叫做线响应（line of response）。探测器阵列D在一个固定时间内采集沿 ( 方向的所有线响应，得该二维物体对应于 ( 方向的一个一维投影分布，测量值P（yr ,(）代表了该投影线上的被测放射性分布f(x, y)之和（线积分），yr 反映了入射光打在yr 轴上的横向位置。显然，某一个投影值给不出该投影线上各物体点的前后关系，为此需要围绕该二维物体的不同角度的投影分布测量值，即进行00到1800的旋转采样。取（yr ,(）为坐标，则由此表达的PET测量计数形式即为所谓的正弦图格式（sinogram）。利用这些投影数据可求解出该二维切片的待测放射性分布，这个过程就是图像重建。将重建出的诸断层顺次堆积在一起就得到了三维图像，或者说准三维图像。三维PET用三维空间的一次三维数据采集实现的三维空间图像重建，是真正的三维图像重建。 以上所述逐个断层分别进行重建的方法为传统的二维重建，其重建方法大致可分为两大类：解析（变换）法和代数（迭代）法，其中解析法中的滤波反投影法作为传统的重建方法，已广泛应用于各种断层扫描仪，而快速迭代算法作为可选项也在新型的核医学断层扫描仪中使用。三维PET能更有效地利用PET符合计数所采集到的数目巨大的投影数据（所谓过剩数据），可数倍地（3-5倍）提高探测灵敏度。 图1.一个断层的投影数据的采集示意图。（x, y）是二维物体空间的坐标， （xr, yr）是旋转采样的旋转坐标系。 2. 解析法图像重建 解析法是以中心切片定理（central slice theorem）为理论基础的求逆过程。根据具体计算过程不同，又可分为滤波反投影法（filtered back-projection, FBP），反投影滤波法（back-projection filtered, BPF）或 ( 滤波法（( -filtered layergram）,卷积反投影法(convolution-backprojection) 等算法。从严格的数学理论可以证明，在无穷取样的极限条件下，这些方法是完全相同的，都能给出断层图像的严格解。 其中FBP中的滤波是一维的，因运算速度快而为大多数断层扫描系统所采用[1]。 FBP 滤波反投影法首先在频率空间对 进行滤波，得滤波后的投影数据 （1） 其中P(,()是投影数据的一维傅立叶变换，变量 是与坐标yr 对应的频率坐标，式中因子 || 被称为斜坡函数（Ramp function, or Ramp filter）, 而W()是窗函数，用于控制噪声，窗函数的形状权衡着统计噪声和空间分辨。窗函数与斜坡函数的乘积形成滤波器。当W()取为常数1的情况下，即对高频信号不作抑制，空间分辨最好，但所重建的图像不平滑，易产生振荡。反之，过多压抑高频成分的低通窗函数会造成重建图像的模糊，故在变换法中低噪声和高分辨对滤波器的要求是矛盾的，需折衷选择。例如，典型的常用窗函数有General Hamming 窗（包含Hanning 窗与Hamming窗），Butterworth窗，以及可增强边缘和改善图像空间分辨率的Metz和Wiener窗等，见图2。各种滤波器都带有可调参数，可根据具体情况而取不同的值。滤波器的选取标准是既抑制低频本底和高频噪声，又尽可能地保留靶值（有用信息）以提高信噪比。临床中的视觉经验也是决定滤波参数的重要因素。 下一步对滤波后的投影 进行反投影，即对角度 ( 积分。所谓反投影（back-projection），就是对物体中任意一点（x，y），将所有经过该点的投影累加起来的运算，也就是将各投影值均匀（或加权）分配给投影线经过的每一个像素的