正电子发射断层扫描仪PET中的数据校正常用方法.docx

正电子发射断层扫描仪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.Key words: PET; positron; correctionPET是英文名称Positron Emission Tomography的缩写，即"正电 子发射断层扫描仪”。

PET是在分子水平上利用影像技术研究人体心 脑代谢和受体功能的一种最先进的设备PET已成为肿瘤、心、脑疾 病诊断的一种最有效的方法在肿瘤学、心血管疾病学、神经系统疾 病学、和新医药学开发等研究领域中已经显示出它的卓越性能PET 是核医学领域中最先进的医疗诊断设备，是目前唯一可以在活体分子 水平上完成生物学显示的影像技术，被誉为20世纪最伟大的十项发 明之一目前PET已经被广泛应用于核医学成像临床检查、疗效评价、 药物开发和基础研究之中［1］引言PET系统由机械系统，探测器环，电子线路，数据采集与处理系统 和计算机系统，以及相应的图像重建和图像处理软件系统等构成［2］， 见图1PET利用其封闭环绕型探测器阵列对注入人体药物所引发的 背对背的一对Y光子进行符合测量，利用计算机对符合测量形成投影 线LOR (Line of Response)进行图像重建可求解出人体内部横断面 放射性分布但由于PET设备和测量过程存在误差，在日常使用中， 采集到的原始数据在图象重建之前要根据仪器状况和病人状况进行 校正此外，PET是复杂精密的大型医用诊断设备，定期要做性能测 试和(重新)标定，不断更新校正参数以确保系统处于良好的工作状 态。

本文概述了 PET中常用的数据校正方法，主要是2D （2维采集模式） 校正方法，包括：探测器灵敏度校正（归一化），同位素时间衰变校 正，死时间校正，偶然符合校正，散射符合校正，衰减校正，几何校 正及其他校正等它们的顺序为图21各种校正的意义及常用方法〖*2〗1.1探测器效率的归一化（normalization） PET中有成千上万个探测单元，受其各自几何位 置和性能差异的影响，例如晶体条发光效率、晶体条与光电倍增管的 耦合、晶体条对符合线的张角不同等，使其探测效率不尽一致其校 正方法是利用均匀分布的放射源，测量每个测量单元的计数响应Di， i=1，…，M, M是PET探测单元总数，算出归一化因子： NOPMi二苍Mi=1Di/MDi这些因子以文件方式存于计算机，在对病人进 行PET测量时，将测量值乘以相应归一化因子就实现了探测器效率校 正另外，由每天质控的空扫描检测探测器性能的漂移情况，必要时 更新校正系数1.2同位素时间衰变校正正电子类核素的寿命都非常短（如18F为110分钟），放射性衰变会 使药物的强度随指数规律逐渐降低特别是对于动态采集、全身扫描、 门控采集和定量研究则必须考虑该项校正。

根据指数衰变规律，注射 时放射性强度为A0、衰变系数为入的药物经过时间t1采集到某一 帧的时候，放射性强度下降到A（t）二A0e-入t1，据此，不难通过采集 时刻的计数率求出注射时刻的药物强度把e入t1作为刻度因子乘以 该帧各个像素的计数值，就能将图象归一到注射时刻的情况至于每 一帧之间的差别，如果各帧的采集时间比药物的半衰期短，则可以忽 略在每帧采集过程中放射性强度的变化但在计算标准摄取值(SUV， standard uptake value)时，需根据帧采集周期的大小将计数率校 正到药物注射时刻［1］1.3死时间校正系统的死时间(deadtime )是指系统处理每个事件所需的时间，它 取决于探测器与电子学的时间特性以及数据处理器的速度、随机缓存 器的性能等诸多因素如果在后一个湮灭事件发生之前来不及处理完 前一个事件，这两个事件就会丢失，这就是死时间损失°PET出厂前 都要进行死时间损失测量：根据测量结果画出计数率一一药物强度曲 线，如图3所示图2临床PET数据校正的顺序图图3计数率一一药物强度曲线在强度低的时候，计数率随药物强度正比增加，呈直线上升，当药 物强度增加到某一限度后，曲线逐渐弯曲，它与直线的距离就是丢失 的计数率，可以据此计算与记录校正参数以便进行死时间校正。

死时 间校正是有范围的，例如当上述曲线随药物强度呈下降趋势时，无法 再进行校正事实上有效评估PET计数特性的是噪声等效计数(NEC) NEC定义为在无散射和偶然符合计数条件下达到同样的信噪比所需的 真符合计数［3］，由于散射和偶然符合的存在，使NEC先于计数率而 饱和，因此要注意死时间校正的有效范围1.4偶然符合校正是指两个或两个以上没有关联的光子被同时探测到而造成的符合计 数，也叫随机符合(random coincidence)，见图4，它与活度的平 方成反比，它增加图象的噪声，影响图象的对比度偶然符合校正硬 件方法是使用延迟符合电路只要延迟时间大于两倍的符合电路时间 窗宽度，就能保证该符合电路输出中没有真的湮灭符合事件而只有偶 然符合计数，然后再从总计数中减去该方法简明有效、实时、 速度快，易于实现，商用PET多采用这种方法图4偶然符合与散射符合形成的示意图偶然符合校正软件方法：用公式NAC=2 t N1N2估计偶然符合计数率， 然后在总符合计数率扣除，其中N1和N2分别是参与符合的两个探测 单元的单事例计数率，t是符合时间窗宽1.5散射符合校正散射符合：主要是指组织中正电子湮灭产生的两个光子在到达探测 器之前其中之一或全部发生了康普顿散射而偏移了原来的运行轨迹， 且无法用能量窗方法有效去除，造成错误的符合信息(如图4所示)。

散射符合影响图象的对比度散射校正有多种硬件与软件的校正方 法，如双能量窗法，三能量窗法、卷积扣除法、人工神经网络法、MONTE CARLO,模拟法等卷积扣除法(Convolution Subtraction)假设投影空间的散射符合 分布可以通过真实符合分布的积分变换来近似表述这种积分变换的 核(kernel)函数一般以指数分布函数或者高斯分布函数的形式出现 如果以T表示真实符合，S表示散射符合，R表示实际测量的符合分 布，则在投影空间S=T信h,因此真实符合就可以通过下式近似求取： T=R-S^r-r信h,或采用反卷积方法更精确地求解核函数h的求取 一般采用实验测量加函数拟合的方式具体作法是用线源或点源放置 一个模拟人体的水模中，在离中心轴线不同的距离上测量其符合投影 值而后对不同距离上的散射分布采取非线性最小二乘拟合的方法或 尾端拟合方法求取核函数通过对核函数积分求出散射分数(散射所 占比例)，可对散射做进一步校正双能窗法(Dual Energy Window):所有的符合事例都可以在两个相 邻的能量窗内获取无论散射的、还是不散射的符合事例都收集在高 窗(380〜850 keV)中，低窗(200〜380 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图5衰减校正前后桶源重建图象在PET中，某一条符合线上的符合计数A表示如下：Ac=A0e- #u (x)dx=A0・T，该公式表明，某条符合线上的衰减因子T与源点的位 置无关，即只要沿同一路径传播，不论湮灭点在哪里，测得的符合强 度都相等这样就可以用一条置于人体之外并与人体轴线平行的68Ge 线(棒)源，来测量过源点的各条符合路径的衰减情况。

线源绕人体 一周完成透扫(transmission scan)，就能取得沿所有符合线的衰减 结果Ac，与没有病人时的空扫结果A0 (blank scan)相比，就得到 衰减校正因子t =A0/Ac=e- J ba u dl,衰减量的测量精度取决于透扫 的计数密度，它的统计噪声会进入校正后的发射图象，因此透射图象 的计数一般要比被校正的发射图象的计数高十倍对透扫数据进行平 滑能够减少统计噪声，也能改善衰减校正的质量图5显示了均匀桶 源经衰减校正的效果此外常用的衰减校正方法还有区域分割(segmentation)法和使用137Cs源的单光子校正法等［4］1.7几何校正PET中的原始正弦(SINO)图是由探测器环上的探测器对通过事件 符合，探测器编码、角度换算而得到的由于探测器圆环结构，在某 一角度下相邻符合线间的实际距离从中心到两边逐渐减小，空间采样 间距是不等间距的，也就是说直接所得的正弦图是错位的因此应给 予几何弧度校正，才能用以图象重建，否则重建的图象是畸形的校 正方法是通过线性插值计算或其他插值运算等分坐标位置上的计数 值［5］,得到等物理间距的新的正弦图迭代法图象重建可通过修正 其系数矩阵而直接对原始正弦图进行重建，避免了线性插值计算，可 提高重建精度。

2 PET中的刻度与定标原始PET数据经过上述校正与重建处理之后形成断层图象，其图象 象素的绝对大小不反映相应组织中放射性核素浓度的大小在定量分 析时，要根据医师所选择的象素单位(如Bq/ml.uCi/ml等)，将重 建后的图象乘以刻度因子(scale factor )输出为计数率/象素 (count拟rate/pixel)值，再乘以定标因子(calibration factor), 则输出比活度值利用活度值和病人体重及用药量等可进一步计算 SUV值，这可用于肿瘤良恶的判断PET中的刻度因子与SINO图总测量计数、图象重建算法有关，是一 种SINO图计数与图象象素值总和的除运算，每种PET具体计算方法 不尽一样而定标因子的计算如下：CF二刻度圆柱体源体积活度•刻 度源分支比圆柱体重建图像感兴趣区的计数。