1正电子发射断层扫描仪 PET 中的数据校正常用方法【摘要 】 介绍了正电子发射断层扫描（PET）中各种校正的意义及常见算法，这些校正包括归一化校正，衰变校正，散射与衰减校正，活度刻度等。对校正算法的最新进展和 PET 相关设备中的校正算法也做了些介绍。恰当的校正对提高 PET 的成象质量及定量分析的准确性非常重要。因此，校正算法是 PET 设备软件系统中所必不可少的组成部分。 【关键词】 PET； 正电子； 校正Abstract: Various corrections in clinical positron emission tomography (PET) are described, including detectors normalization, decay correction, scatter and attenuation correction, calibration, etc. The up to data developments in correction of PET and PET related devises are also introduced. Appropriate corrections are important for PET images with good quality and for quantitative analysis in clinics. PET correction algorithms consist of a necessary part in the software of PET.2Key words: PET; positron; correctionPET 是英文名称 Positron Emission Tomography 的缩写，即“正电子发射断层扫描仪” 。PET 是在分子水平上利用影像技术研究人体心脑代谢和受体功能的一种最先进的设备。PET 已成为肿瘤、心、脑疾病诊断的一种最有效的方法。在肿瘤学、心血管疾病学、神经系统疾病学、和新医药学开发等研究领域中已经显示出它的卓越性能。PET 是核医学领域中最先进的医疗诊断设备，是目前唯一可以在活体分子水平上完成生物学显示的影像技术，被誉为 20 世纪最伟大的十项发明之一。目前 PET 已经被广泛应用于核医学成像临床检查、疗效评价、药物开发和基础研究之中1 。引言PET 系统由机械系统，探测器环，电子线路，数据采集与处理系统和计算机系统，以及相应的图像重建和图像处理软件系统等构成2 ，见图 1。PET 利用其封闭环绕型探测器阵列对注入人体药物所引发的背对背的一对 光子进行符合测量，利用计算机对符合测量形成投影线 LOR（Line of Response）进行图像重建可求解出人体内部横断面放射性分布。但由于 PET 设备和测量过程存在误差，在日常使用中，采集到的原始数据在图象重建之前要根据仪器状况和病人状况进行校正。此外，PET 是复杂精密的大型医用诊断设备，3定期要做性能测试和（重新）标定，不断更新校正参数以确保系统处于良好的工作状态。本文概述了 PET 中常用的数据校正方法，主要是 2D（2 维采集模式）校正方法，包括：探测器灵敏度校正（归一化） ，同位素时间衰变校正，死时间校正，偶然符合校正，散射符合校正，衰减校正，几何校正及其他校正等。它们的顺序为图 2。1 各种校正的意义及常用方法*21.1 探测器效率的归一化（normalization） PET 中有成千上万个探测单元，受其各自几何位置和性能差异的影响，例如晶体条发光效率、晶体条与光电倍增管的耦合、晶体条对符合线的张角不同等，使其探测效率不尽一致。其校正方法是利用均匀分布的放射源，测量每个测量单元的计数响应 Di，i=1，M, M 是 PET 探测单元总数，算出归一化因子: NOPMi=Mi=1Di/MDi。这些因子以文件方式存于计算机，在对病人进行 PET 测量时，将测量值乘以相应归一化因子就实现了探测器效率校正。另外，由每天质控的空扫描检测探测器性能的漂移情况，必要时更新校正系数。1.2 同位素时间衰变校正正电子类核素的寿命都非常短（如 18F 为 110 分钟） ，放射性衰4变会使药物的强度随指数规律逐渐降低。特别是对于动态采集、全身扫描、门控采集和定量研究则必须考虑该项校正。根据指数衰变规律，注射时放射性强度为 A0、衰变系数为 的药物经过时间 t1采集到某一帧的时候，放射性强度下降到 A(t)=A0e-t1，据此，不难通过采集时刻的计数率求出注射时刻的药物强度。把 et1 作为刻度因子乘以该帧各个像素的计数值，就能将图象归一到注射时刻的情况。至于每一帧之间的差别，如果各帧的采集时间比药物的半衰期短，则可以忽略在每帧采集过程中放射性强度的变化。但在计算标准摄取值（SUV，standard uptake value）时，需根据帧采集周期的大小将计数率校正到药物注射时刻1 。1.3 死时间校正系统的死时间（dead time）是指系统处理每个事件所需的时间，它取决于探测器与电子学的时间特性以及数据处理器的速度、随机缓存器的性能等诸多因素。如果在后一个湮灭事件发生之前来不及处理完前一个事件，这两个事件就会丢失，这就是死时间损失。PET 出厂前都要进行死时间损失测量：根据测量结果画出计数率药物强度曲线，如图 3 所示。图 2 临床 PET 数据校正的顺序图图 3 计数率 药物强度曲线5在强度低的时候，计数率随药物强度正比增加，呈直线上升，当药物强度增加到某一限度后，曲线逐渐弯曲，它与直线的距离就是丢失的计数率，可以据此计算与记录校正参数以便进行死时间校正。死时间校正是有范围的，例如当上述曲线随药物强度呈下降趋势时，无法再进行校正。事实上有效评估 PET 计数特性的是噪声等效计数（NEC） 。NEC 定义为在无散射和偶然符合计数条件下达到同样的信噪比所需的真符合计数3 ，由于散射和偶然符合的存在，使NEC 先于计数率而饱和，因此要注意死时间校正的有效范围。1.4 偶然符合校正是指两个或两个以上没有关联的光子被同时探测到而造成的符合计数，也叫随机符合（random coincidence） ，见图 4，它与活度的平方成反比，它增加图象的噪声，影响图象的对比度。偶然符合校正硬件方法是使用延迟符合电路。只要延迟时间大于两倍的符合电路时间窗宽度，就能保证该符合电路输出中没有真的湮灭符合事件而只有偶然符合计数，然后再从总计数中减去。该方法简明有效、实时在线、速度快，易于实现，商用 PET 多采用这种方法。图 4 偶然符合与散射符合形成的示意图偶然符合校正软件方法：用公式 NAC=2N1N2 估计偶然符合6计数率，然后在总符合计数率扣除，其中 N1 和 N2 分别是参与符合的两个探测单元的单事例计数率， 是符合时间窗宽。1.5 散射符合校正散射符合：主要是指组织中正电子湮灭产生的两个光子在到达探测器之前其中之一或全部发生了康普顿散射而偏移了原来的运行轨迹，且无法用能量窗方法有效去除，造成错误的符合信息（如图 4所示） 。散射符合影响图象的对比度。散射校正有多种硬件与软件的校正方法，如双能量窗法，三能量窗法、卷积扣除法、人工神经网络法、MONTE CARLO.模拟法等。卷积扣除法（Convolution Subtraction）假设投影空间的散射符合分布可以通过真实符合分布的积分变换来近似表述。这种积分变换的核（kernel）函数一般以指数分布函数或者高斯分布函数的形式出现。如果以 T 表示真实符合，S 表示散射符合， R 表示实际测量的符合分布，则在投影空间 S=Th，因此真实符合就可以通过下式近似求取：T=R-SR-Rh ，或采用反卷积方法更精确地求解。核函数 h 的求取一般采用实验测量加函数拟合的方式。具体作法是用线源或点源放置一个模拟人体的水模中，在离中心轴线不同的距离上测量其符合投影值。而后对不同距离上的散射分布采取非线性最小二乘拟合的方法或尾端拟合方法求取核函数。通过对核函数积7分求出散射分数（散射所占比例） ，可对散射做进一步校正。双能窗法（Dual Energy Window）：所有的符合事例都可以在两个相邻的能量窗内获取。无论散射的、还是不散射的符合事例都收集在高窗（380850 keV）中，低窗（200380 keV）中只有散射事例。假设所有的散射符合均有相同的空间分布，则将高窗中的符合计数减去低窗中的符合计数，就可得到真实符合计数。而实际上，光子在低能部分对目标体的依赖性比在高能部分要大得多，因而该方法是近似的。商用 PET 使用的方法应该简洁、实用、易于实现，例如 ECAT EXCAT HR+型 PET 中，对 2D 采集采用卷积扣除法。1.6 衰减校正衰减校正是针对体内肌肉和骨骼等对光子的吸收衰减而进行的校正，从而得到真实的放射性药物分布图。软组织对 511 keV 的光子质量衰减系数约为 0.095 cm2/g，半衰减厚度约为 7.2 cm。对直径约 20 cm 的头部显像，超过 85%的光子被衰减，宽 40 cm 的躯干可将 95%以上的光子吸收掉，所以必须进行衰减校正，否则会造成PET 图象中外表组织影象过亮，内部组织影象过暗的现象，见图5。8图 5 衰减校正前后桶源重建图象在 PET 中，某一条符合线上的符合计数 A 表示如下：Ac=A0e-(x)dx=A0，该公式表明，某条符合线上的衰减因子 与源点的位置无关，即只要沿同一路径传播，不论湮灭点在哪里，测得的符合强度都相等。这样就可以用一条置于人体之外并与人体轴线平行的 68Ge 线（棒）源，来测量过源点的各条符合路径的衰减情况。线源绕人体一周完成透扫（transmission scan） ，就能取得沿所有符合线的衰减结果 Ac，与没有病人时的空扫结果 A0（blank scan）相比，就得到衰减校正因子 =A0/Ac=e-badl. 衰减量的测量精度取决于透扫的计数密度，它的统计噪声会进入校正后的发射图象，因此透射图象的计数一般要比被校正的发射图象的计数高十倍。对透扫数据进行平滑能够减少统计噪声，也能改善衰减校正的质量。图 5 显示了均匀桶源经衰减校正的效果。此外常用的衰减校正方法还有区域分割（segmentation ）法和使用 137Cs 源的单光子校正法等4 。1.7 几何校正PET 中的原始正弦（SINO）图是由探测器环上的探测器对通过事件符合，探测器编码、角度换算而得到的。由于探测器圆环结构，在某一角度下相邻符合线间的实际距离从中心到两边逐渐减小，空9间采样间距是不等间距的，也就是说直接所得的正弦图是错位的。因此应给予几何弧度校正，才能用以图象重建，否则重建的图象是畸形的。校正方法是通过线性插值计算或其他插值运算等分坐标位置上的计数值5 ，得到等物理间距的新的正弦图。迭代法图象重建可通过修正其系数矩阵而直接对原始正弦图进行重建，避免了线性插值计算，可提高重建精度。2 PET 中的刻度与定标原始 PET 数据经过上述校正与重建处理之后形成断层图象，其图象象素的绝对大小不反映相应组织中放射性核素浓度的大小。在定量分析时，要根据医师所选择的象素单位（如 Bq/ml、Ci/ml 等），将重建后的图象乘以刻度因子（scale factor）输出为计数率象素(countrate/pixel) 值，再乘以定标因子(calibration factor)，则输出比活度值。利用活度值和病人体重及用药量等可进一步计算SUV 值，这可用于肿瘤良恶的判断。PET 中的刻度因子与 SINO 图总测量计数、图象重建算法有关，是一种 SINO 图计数与图象象素值总和的除运算，每种 PET 具体计算方法不尽一样。而定标因子的计算如下：CF=刻度圆柱体源体积活度 刻度源分支比圆柱体重建图像感兴趣区的