第３１卷２０１１年第１Ｏ期 １０月 核电子学与探测技术 Ｎｕｃｌｅａｒ Ｅｌｅｃｔｒｏｎｉｃｓ＆Ｄｅｔｅｃｔｉｏｎ Ｔｅｃｈｎｏｌｏｇｙ ＶｏＩ＿３ｌ Ｎｏ．１０ ０ｃｔ． ２０ｌ１ ＰｅｎｅｌｏＰＥＴ：一种ＰＥＴ专用的蒙特卡罗仿真工具 刘豪佳，赵书俊，张 斌，周山虎 （郑州大学，河南郑州４５０００１） 摘要：ＰｅｎｅｌｏＰＥＴ是一种基于ＰＥＮＥＬＯＰＥ的ＰＥＴ专用蒙特卡罗仿真工具。文章介绍了ＰｅｎｅｌｏＰＥＴ 的主要特点及使用ＰｅｎｅｌｏＰＥＴ进行ＰＥＴ仿真的一般过程，并以ＧＥ Ｈｅａｈｈｃａｒｅ的双环型ｅＸｐｌｏｒｅ Ｖｉｓｔａ ｍｉ— ｅｒｏＰＥＴ为原型验证了ＰｅｎｅｌｏＰＥＴ仿真的有效性。仿真结果显示，ＰｅｎｅｌｏＰＥＴ仿真数据和实际实验结果及 ＧＡＴＥ仿真结果之间均具有良好的一致性，证明了ＰｅｎｅｌｏＰＥＴ是一种精确的ＰＥＴ仿真工具。 关键词：ＰｅｎｅｌｏＰＥＴ；ＰＥＴ；蒙特卡罗仿真 中图分类号：ＴＬ ９９ 文献标识码：Ａ 文章编号：０２５８－０９３４（２０１１）１０—１１４３－０４ 正电子发射断层成像（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ 何结构建模工具及高效的可视化程序。然而， Ｔｏｍｏｇｒａｐｈｙ，ＰＥＴ）是一种无创的分子影像技 若想使这些工具包用于ＰＥＴ和ＳＰＥＣＴ，则需要 术，它在疾病诊断和医药学研究中均具有独特 用户具备丰富的发射断层成像设备建模的专业 的价值。蒙特卡罗（Ｍｏｎｔｅ Ｃａｒｌｏ）仿真方法又称 知识。此外，ＰＥＴ和ＳＰＥＣＴ专用的蒙特卡罗仿 为统计实验方法，它是通过大量的计算机模拟 真工具有ＳｉｍＳＥＴ 、ＰＥＴＳＩＭ 、Ｅｉｄｏｌｏｎ 及 来检验系统的动态特性并归纳得到统计结果的 ＧＡＴＥ 等。ＳｉｍＳＥＴ是一个很强大的ＰＥＴ和 一种随机分析方法。近几十年来，随着高能计 ＳＰＥＣＴ仿真专用代码，可以对物理过程和设备 算机的发展，蒙特卡罗仿真方法已逐渐成为当 本身进行精确高效的建模，但它在探测器几何 前及未来发射断层技术发展的基本工具。蒙特 结构上有一些。比如，不能对探测器的晶 卡罗仿真方法在正电子发射断层成像和单光子 体间隙和晶体阵列进行建模。ＧＡＴＥ基于 发射计算机断层成像（Ｓｉｎｇｌｅ—ｐｈｏｔｏｎ Ｅｍｉｓｓｉｏｎ Ｇｅａｎｔ４，在继承Ｇｅａｎｔ４优点的基础上还增加了 Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ，ＳＰＥＣＴ）的分析计算及实 对时间相关过程的精确仿真，是目前广泛使用 验研究的各方面得到了广泛应用，如设计新的 的仿真平台，但它的缺点在于仿真速度比较慢。 成像设备，优化采集协议，发展新的图像重建算 本文介绍的ＰｅｎｅｌｏＰＥＴ ，是一种基于 法和校正技术等Ｈ Ｊ。 ＰＥＮＥＬＯＰＥ（ＰＥＮｅｔｒａｔｉｏｎ ａｎｄ Ｅｎｅｒｙｇ ＬＯＳＳ ｏｆ 目前，可用于ＰＥＴ和ＳＰＥＣＴ仿真的蒙特卡 Ｐｏｓｉｔｒｏｎｓ ａｎｄ Ｅｌｅｃｔｒｏｎｓ）的ＰＥＴ专用仿真工具０ 罗工具包有很多，每种都有各自的优缺点和不 ＰＥＮＥＬＯＰＥ采用“混合（ｍｉｘｅｄ）仿真法” ， 同的可靠度。通用的蒙特卡罗仿真工具包如 是一种能够模拟能量范围从几百电子伏到１ Ｇ ＥＧＳ４ 、ＭＣＮＰ 、Ｇｅａｎｔ３及Ｇｅａｎｔ４ ，这些工 电子伏的电子、正电子和光子在物质中的传输 具包具有精确有效的物理过程模型，丰富的几 过程的蒙特卡罗代码。它的通用性比ＧＥＡＮＴ４ 差，但它健壮、快速、精确，更适合ＰＥＴ的需要。 此外，它还广泛用于剂量学和放射治疗。 收稿日期：２０１０—１１—０５ 作者简介：刘豪佳（１９８５一），女，河南平顶山人，硕士 １ ＰｅｎｅｌｏＰＥＴ的特点 研究生，研究方向为核医学图像处理。Ｐｅｎｅｌ。ＰＥＴ是由马德里康普顿斯大学核物 １ ｌ４３ 理组研发的，目前的最新版本是２００７年６月发 布的１．１版本。ＰｅｎｅｌｏＰＥＴ具有以下特点： Ｍａｃ ＯＳ ｘ平台下运行ＰｅｎｅｌｏＰＥＴ，均不会出现 移植问题。而另一仿真软件ＧＡＴＥ，则不能在 Ｗｉｎｄｏｗｓ系统下运行 Ｊ，这对于广大Ｗｉｎｄｏｗｓ （１）ＰｅｎｅｌｏＰＥＴ仿真速度快，可以在计算 机集群上并行运行。ＰＥＮＥＬＯＰＥ是用Ｆｏｒｔｒａｎ 编写的，Ｆｏｒｔｒａｎ语言高度标准化，被认为是科 用户来说是很不方便的。 此外，ＰｅｎｅｌｏＰＥＴ还拥有高效的可视化工 具（ｇｖｉｅｗ２ｄ，ｇｖｉｅｗ３ｄ），用于在仿真中显示和测 试几何结构，这在ＰＥＴ的设计阶段是非常有用 学计算的专用语言。Ｆｏｔｒｒａｎ编译器的特点是 高速度和高精度，ＰＥＮＥＬＯＰＥ用Ｆｏｔｒｒａｎ语言编 写，因此必然快速便捷。同时，ＰｅｎｅｌｏＰＥＴ支持 Ｐｙｌｏｎ脚本语言，这使得用户可以在使用不同 操作系统的多台计算机上实现仿真程序的并行 运行，仿真时间会随所使用的ＣＰＵ数量线性递 减，因此大大加快了仿真过程。 ＧＡ１’Ｅ作为目前使用非常广泛的发射断层 仿真软件，在探测器配置上的灵活性是其主要 优点，但这一点是以增加计算资源为代价的。 同其他ＰＥＴ仿真工具相比，ＰｅｎｅｌｏＰＥＴ设置仿 真过程所需的时间适中，计算量适中，而程序运 行效率非常高，速度很快。 （２）ＰｅｎｅｌｏＰＥＴ的主要目的是优化小动物 ＰＥＴ的设计。目前，各种小动物ＰＥＴ已广泛应 用于生物医药学研究，成为实验科学和临床科 学的重要桥梁 。在优化ＰＥＴ设计的过程中， 我们通常要在放射源相同的情况下进行多次相 似的仿真实验，即改变ＰＥＴ的一个或多个设置 （如探测器环的直径、晶体尺寸、电子学电路的 参数等）分别进行仿真，以比对仿真结果。 ＰｅｎｅｌｏＰＥＴ快速稳定，并且能够同时在计算机 集群上执行多个仿真实例，这一点正是优化 ＰＥＴ设计所需要的。所以，ＰｅｎｅｌｏＰＥＴ是目前 优化ＰＥＴ设计的一种非常好的工具。 （３）ＰｅｎｅｌｏＰＥＴ使用方便，易学易用。 ＰｅｎｅｌｏＰＥＴ目标之一就是实现易用的仿真程 序。由于缺乏面向对象的结构，Ｆｏｒｔｒａｎ语言在 年轻的科研人员中并不是很流行，因此对很多 用户来说，ＰＥＮＥＬＯＰＥ的使用并不方便。正是 基于这个原因，ＰｅｎｅｌｏＰＥＴ研发组对ＰＥＮＥＬＯＰＥ 进行封装，使得用户能够很容易地使用ＰＥＮＥＬ－ ＯＰＥ定义一个完整的ＰＥＴ仿真系统，这就是 ＰｅｎｅｌｏＰＥＴ。使用ＰｅｎｅｌｏＰＥＴ，用户不需要任何 的Ｆｏｒｔｒａｎ或其他编程语言知识，只需要编制几 个简单的文本文件就可以对ＰＥＴ系统进行全 面精确的仿真。 ＰｅｎｅｌｏＰＥＴ可以在任何拥有Ｆｏｔｒｒａｎ编译器 的操作系统下运行，在Ｗｉｎｄｏｗｓ，Ｌｉｎｕｘ，Ｕｎｉｘ和 １１４４ 的。图１就是ＰｅｎｅｌｏＰＥＴ仿真中得到的ＧＥ Ｈｅａｌｔｈｃａｒｅ双环型（ＤＲ）ｅＸｐｌｏｒｅ Ｖｉｓｔａ ｍｉｃｒｏＰＥＴ 的结构图。 图１ 双环型ｅＸｐｌｏｒｅ Ｖｉｓｔａ ｍｉｃｒｏＰＥＴ结构示意图 （４）ＰｅｎｅｌｏＰＥＴ提供了多种数据输出格 式，便于数据的后期分析处理。ＰｅｎｅｌｏＰＥＴ能 够生成多种形式的输出，如投影正弦图，ＬＯＲ （１ｉｎｅ ｏｆ ｒｅｓｐｏｎｓｅ）直方图，单光子事件和符合事 件分布图及能谱图。ＰｅｎｅｌｏＰＥＴ的输出可设置 为三个不同的详细度等级。详细度最高时，记 录所有交互作用的相关信息。详细度二级，只 记录单光子（ｓｉｎｇｌｅ）事件及其相关信息。详细 度三级，也就是最低级，只以表模式（１ｉｓｔ ｍｏｄｅ） 记录符合事件和脉冲堆积事件，这里的符合事 件包括散射符合，随机符合和自符合（ｓｅｌｆ—ＣＯ— ｉｎｃｉｄｅｎｃｅ） Ｊ。此外，表模式下得到的输出文件 可以直接转换为ＲＯＯＴ文件，便于后期使用 ＲＯＯＴ １刮这个面向对象的数据分析工具对仿真 结果进行分析处理。 为了简化图像重建过程，ＰｅｎｅｌｏＰＥＴ输出 的正弦图符合ＳＴＩＲ（Ｓｏｆｔｗａｒｅ ｆｏｒ Ｔｏｍｏｇｒａｐｈｉｃ Ｉｍａｇｅ Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ） ， 的格式要求，可以直 接使用ＳＴＩＲ得到图像重建结果。 ２仿真过程 ＰｅｎｅｌｏＰＥＴ由两个主要模块组成。第一个 模块用于处理ＰＥＮＥＬＯＰＥ的规则，第二个模块 负责处理仿真环境下ＰＥＴ系统的具体特征，如 探测器和电子学电路。使用ＰｅｎｅｌｏＰＥＴ进行 ＰＥＴ仿真通常包括以下几个过程 引： （１）定义输人文件ｍａｉｎ．ｉｎｐ。ｍａｉｎ．ｉｎｐ中 ＬＹＳＯ晶体，长度为７ ｍｍ，外层ＧＳＯ晶体，长度 为８ ｍｍ。每个晶体体素的横截面积是１．４５ ｍｍ×１．４５ ｍｍ，它们之间有一个０．１ ｍｍ厚的 设置仿真的一般参数，如数据采集协议，采集时 间，能量窗和符合窗，死时间，输出格式等。如 果需要的话，也对次级粒子、正电子的射程和非 线性进行设置。 白色反射层。下面以ＤＲ ｅＸｐｌｏｒｅ Ｖｉｓｔａ为原型， 分别给出ＰｅｎｅｌｏＰＥＴ和ＧＡＴＥ的仿真结果，并 与实验数据进行比对。对灵敏度和噪声等效计 数率（ＮＥＣＲ）的测试方法符合ＮＥＭＡ＿】 标准。 在轴向灵敏度测试中，选择直径为１ ｍｍ的培Ｆ 点源，能量窗设为２５０ ｋｅＶ～７００ ｋｅＶ，仿真结果 （２）定义输入文件ｓｃａｎｎｅｒ．ｉｎｐ。在ｓｃａｎ— ｎｅｒ．ｉｎｐ中完成对探测器结构的定义，包括探测 器的直径、环数，相邻环间距，每环的模块数，每 个探测器模块包含的晶体阵列大小，闪烁晶体 的层数和各层的材料等。通常ＰＥＴ都由多环 的长方体形探测器模块组成，每个探测器模块 包含一层或多层的闪烁晶体。ＰｅｎｅｌｏＰＥＴ中已 定义了一些晶体的材料，用户也可以自行定义 新的材料文件，需要定义新材料的名称，能量分 辨率，时间分辨率，光脉冲的上升和降落时间 等，然后将新文件名添加到ｍａｔ—ｎａｍｅｓ．ｉｎｐ文 件中，并将新文件复制到ＰｅｎｅｌｏＰＥＴ目录下即 可。 （３）定义输入文件ｏｂｊｅｃｔ．ｉｎｐ。ｏｂｊｅｃｔ．ｉｎｐ 中定义其他的对象，如幻影（ｐｈａｎｔｏｍ）和屏蔽材 料（ｓｈｉｅｌｄｉｎｇ）。 （４）定义输入文件ｓｏｕｒｃｅ．ｉｎｐ。ｓｏｕｒｃｅ．ｉｎｐ 中定义放射源的形状，位置，活度，同位素类型， 如果需要的话，还可以对发射方向角进行定义。 （５）输入文件定义完成后，运行可执行文 件ｐｅｎｅｌｏｐｅｔ．ｅｘｅ即可开始仿真过程，当屏幕上 显示“ＣＲＥＡＴＩＮＧ ＰＥＮＥＬＯＰＥ ＣＲＯＳＳ ＳＥＣＴＩＯＮ ＤＡＴＡＢＡＳＥ…”时停止仿真，此时已在工作目 录下生成ｓｃａｎｎｅｒ．ｇｅｏ，ｈｉｔｓ．１ｉｓｔ，ｓｉｎｇｌｅｓ．１ｉｓｔ和 ｃｏｉｎｃｉｄｅｎｃｅｓ．１ｉｓｔ等文件（具体生成哪些文件跟 ｍａｉｎ．ｉｎｐ中的设置有关），用户可以使用程序 ｇｖｉｅｗ２ｄ．ｅｘｅ或ｇｖｉｅｗ３ｄ．ｅｘｅ打开ｓｃａｎｎｅｒ．ｇｅｏ 以显示探测器的几何结构，也可以使用其他软 件，如ｇｎｕｐｌｏｔ、ＲＯＯＴ来分析其他输出文件。 ３仿真实例 ｅＸｐｌｏｒｅ Ｖｉｓｔａ ｍｉｃｒｏＰＥＴ是ＧＥ公司为满足 临床前研究的需要而开发的小动物ＰＥＴ，包括 双环（ＤＲ）和单环（ＳＲ）两个型号。ＤＲ ｅＸｐｌｏｒｅ Ｖｉｓｔａ包括两个环，环直径为１１．８ ｃｍ，孔径为８ ｃｍ。每环由１８个ＰＭＴ探测模块组成，每个模 块连接着一个１３×１３双层闪烁晶体阵列，内层 如图２。ＮＥＣＲ测试中，选择低活度的 Ｆ线源， 以减少死时间、随机符合和堆积的影响，仿真结 果如图３。图２和图３中的实验数据来自Ｙｕ． ｃｈｕａｎ Ｗａｎｇ等的工作 。 实验数据—— 芝趟鲻 ．‘、 ４ ４ ３ ３ ２ ２ ｌ ｌ Ｏ Ｏ ＰｅｎｅｌｏＰＥＴ… ５ ５ ５ ５ ５ ＧＡＴＥ …一 ∥ ∥ ｜｛ 嚣 ．２５ —２Ｏ一１５－１０－５ ０ ５ １０ １５ ２Ｏ ２５ 轴向位置／ｃｍ 图２轴向灵敏度 图３噪声等效计数率 ４ 结论 ＰｅｎｅｌｏＰＥＴ作为ＰＥＴ专用的仿真软件，可 用于系统矩阵的计算，新扫描仪的开发，验证散 射校正、随机符合校正、堆积校正，衰减校正等 校正技术的效果，发展更优的图像重建算法等。 同时，ＰｅｎｅｌｏＰＥＴ具有运行速度快、界面友好、 仿真结果易于分析处理等优点。仿真结果表 明，ＰｅｎｅｌｏＰＥＴ仿真数据和实验数据之间的差 异很小且其仿真速度优于ＧＡＴＥ。因此，我们 可以说ＰｅｎｅｌｏＰＥＴ是一个精确高效的ＰＥＴ仿真 １１４５ 工具。未来的工作是将ＰｅｎｅｌｏＰＥＴ用于ＳＰＥＣＴ 的仿真： 。 参考文献 ［１］Ｓ Ｊａｎ，Ｇ Ｓａｎｔｉｎ，Ｄ．Ｓｔｍｌ，ｅｔ ａ１．ＧＡＴＥ：ａ ｓｉｍｕｌａｔｉｏｎ ｔｏｏｌｋｉｔ ｆｏｒ ＰＥＴ ａｎｄ ＳＰＥＣＴ［Ｊ：．Ｐｈｙｓｉｃｓ ｉｎ Ｍｅｄｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ，２００４，４９（１９）：４５４３—４５６１． ［２］Ｗ．Ｒ．Ｎｅｌｓｏｎ，Ｈ＿Ｈｉｒａｙａｍａ，Ｄ．Ｗ．Ｏ．Ｒｏｇｅｒｓ． Ｓｔａｎｆｏｒｄ Ｌｉｎｅａｒ Ａｃｃｅｌｅｒａｔｏｒ Ｃｅｎｔｅｒ Ｒｅｐｏｒｔ［Ｒ］．Ｓｔａｎ— ｆｏｒｄ，１９８５． ［３］Ｊｕｄｉｔｈ Ｆ．Ｂｒｉｅｓｍｅｉｓｔｅｒ．ＭＣＮＰ一一Ａ Ｇｅｎｅｒａｌ Ｍｏｎｔｅ Ｃａｒｌｏ Ｎ—Ｐａｒｔｉｃｌｅ Ｔｒａｎｓｐｏｒｔ Ｃｏｄｅ． Ｌｏｓ Ａ１ａｍｏｓ Ｎａ— ｔｉｏｎａｌ Ｌａｂｏｒａｔｏｒｙ ｒｅｐｏｒｔ［Ｒ］．Ｌｏｓ Ａ］ａｍｏｓ，１９９３． ［４］Ｇｅａｎｔ４ Ｃｏｌｌａｂｏｒａｔｉｏｎ，Ａｇｏｓｔｉｎｅｌｌｉ Ｓ，Ａｌｌｉｓｏｎ Ｊ，ｅｔ ａ１． ＧＥＡＮＴ４一ａ ｓｉｍｕｌａｔｉｏｎ ｔｏｏｌｋｉｔ［Ｊ］．Ｎｕｃｌｅａｒ Ｉｎｓｔｒｕ— ｍｅｎｔｓ ａｎｄ Ｍｅｔｈｏｄｓ ｉｎ Ｐｈｙｓｉｃｓ Ｒｅｓｅａｒｃｈ Ａ，２００３，５０６ （３）：２５０—３０３． ［５］Ｈａｒｒｉｓｏｎ Ｒ Ｌ．Ａｃｃｅｌｅｒａｔｉｏｎ ｏｆ ＳｉｍＳＥＴ ｐｈｏｔｏｎ ｈｉｓｔｏｒｙ ｇｅｎｅｒａｔｏｒ ｍｏｄｕｌｅ ｏｆ ａ ｐｕｂｌｉｃ ｄｏｍａｉｎ ｓｉｍｕｌａｔｉｏｎ ｓｙｓｔｅｍ ｏｆｒ ｅｍｉｓｓｉｏｎ ｔｏｍｏｇｒａｐｈｙ［Ｃ］．Ｃｏｎｆ．Ｒｅｃ．ＩＥＥＥ Ｎｕ— ｃ１．Ｓｃｉ．Ｓｙｍｐ ａｎｄ Ｍｅｄ．Ｉｍａｇｉｎｇ Ｃｏｎｆ．３：１８３５— １８３８． ［６］Ｔｈｏｍｓｏｎ Ｃ，Ｃａｎｔｕ Ｊ，Ｐｉｃａｒｄ Ｙ．ＰＥＴＳＩＭ：Ｍｏｎｔｅ Ｃａｒｌｏ ｐｒｏｇｒａｍ ｓｉｍｕｌａｔｉｏｎ ｏｆ ａｌｌ ｓｅｎｓｉｔｉｖｉｔｙ ａｎｄ ｒｅｓｏｌｕ— ｔｉｏｎ ｐａｒａｍｅｔｅｒｓ ｏｆ ｃｙｌｉｎｄｒｉｃａｌ ｐｏｓｉｔｒｏｎ ｉｍａｇｉｎｇ ｓｙｓｔｅｍｓ ［Ｊ］．Ｐｈｙｓｉｃｓ ｉｎ Ｍｅｄｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ，１９９２，３７（３）： ７３ｌ一７４９． ［７］Ｚａｉｄｉ Ｈ，Ｓｃｈｅｕｒｅｒ Ａ，Ｍｏｒｅｌ Ｃ．Ａｎ ｏｂｊｅｃｔ—ｏｒｉｅｎｔｅｄ Ｍｏｎｔｅ Ｃａｒｌｏ ｓｉｍｕｌａｔｏｒ ｆｏｒ ３ Ｄ ｃｙｌｉｎｄｒｉｃａｌ ｐｏｓｉｔｒｏｎ ｔｏ— ｍｏｇｒａｐｈｓ［Ｊ］．Ｃｏｍｐｕｔｅ Ｍｅｔｈｏｄｓ Ｐｒｏｇ Ｂｉｏｍｅｄ， １９９９，５８（２）：１３３—１４５． ［８］Ｄ．Ｓｔｒｌｕ，Ｇ．Ｓａｎｔｉｎ，Ｄ．Ｌａｚａｒｏ，ｅｔ ａ１．ＧＡＴＥ （Ｇｅａｎｔ４ Ａｐｐｌｉｃａｔｉｏｎ ｆｏｒ Ｔｏｍｏｇｒａｐｈｉｃ Ｅｍｉｓｓｉｏｎ）：ａ ＰＥＴ／ＳＰＥＣＴ ｇｅｎｅｒａｌ—ｐｕｒｐｏｓｅ ｓｉｍｕｌａｔｉｏｎ ｐｌａｔｆｏｒｍ ［Ｊ］．Ｎｕｃｌｅａｒ Ｐｈｙｓｉｃｓ Ｂ（Ｐｒｏｃ．Ｓｕｐｐ１．），２００３，１２５： ７５—７９． （９ Ｓ Ｅｓｐａｎａ，Ｊ Ｌ Ｈｅ￣ａｉｚ，Ｅ Ｖｉｃｅｎｔｅ，ｅｔ ａ１．ＰｅｎｅｌｏＰ— ＥＴ．ｆｌ Ｍｏｎｔｅ Ｃａｒｌｏ ＰＥＴ ｓｉｍｕｌａｔｉｏｎ ｔｏｏ１ ｂａｓｅｄ ｏｎ ＰＥ— ＮＥＬＯＰＥ：ｆｅａｔｕｒｅｓ ａｎｄ ｖａｌｉｄａｔｉｏｎ［Ｊ］．Ｐｈｙｓｉｃｓ ｉｎ Ｍｅｄｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ，２００９，５４（６）：１７２３—１７４２． １ １４６ ［１Ｏ Ｓ．Ｅｓｐａｎａ，Ｊ．Ｌ．Ｈｅｒｒａｉｚ，Ｊ．Ｊ．Ｖａｑｕｅｒｏ，ｅｔ ａ１． ＰｅｎｅｌｏＰＥＴ：Ａ Ｍｏｎｔｅ Ｃａｒｌｏ ＰＥＴ ｓｉｍｕｌａｔｉｏｎ ｔｏｏｌ ｂａｓｅｄ ｏｎ ＰＥＮＥＬＯＰＥ『Ｃ］．ＩＥＥＥ Ｎｕｃｌｅａｒ Ｓｃｉｅｎｃｅ Ｓｙｍｐｏｓｉｕｍ ａｎｄ Ｍｅｄｉｃａｌ Ｉｍａｇｉｎｇ Ｃｏｎｆｅｒｅｎｃｅ，２００６． ［１１］ｈｔｔｐ：／／ｎｕｃｌｅａｒ．ｆｉｓ．ｕｃｍ．ｅｓ ［１２ Ｊ．Ｂａｒｏ，Ｊ．Ｓｅｍｐａｕ，Ｊ．Ｍ．Ｆｅｒｎａｎｄｅｚ—Ｖａｒｅａ，ｅｔ １ａ．ＰＥＮＥＬＯＰＥ：Ａｎ ａｌｇｏｒｉｔｈｍ ｆｏｒ Ｍｏｎｔｅ Ｃａｒｌｏ ｓｉｍｕｌａ— ｔｉｏｎ ｏｆ ｔｈｅ ｐｅｎｅｔｒａｔｉｏｎ ａｎｄ ｅｎｅｒｇｙ ｌｏｓｓ ｏｆ ｅｌｅｃｔｒｏｎｓ ａｎｄ ｐｏｓｉｔｒｏｎｓ ｉｎ ｍａｔｔｅｒ［Ｊ］．Ｎｕｃｌｅａｒ Ｉｎｓｔｒｕｍｅｎｔｓ ａｎｄ Ｍｅｔｈｏｄｓ ｉｎ Ｐｈｙｓｉｃｓ Ｒｅｓｅａｒｃｈ Ｂ，１９９５，１００（１）：３１— ４６． ［１３］Ｊ．Ｓｅｍｐａｕ，Ｅ．Ａｃｏｓｔａ，Ｊ．Ｂａｒｏ，ｅｔ １ａ．Ａｎ ｌａｇｏｒｉｔｈｍ ｏｆｒ Ｍｏｎｔｅ Ｃａｒｌｏ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｃｏｕｐｌｅｄ ｅｌｅｃｔｒｏｎ—ｐｈｏ－ ｔｏｎ ｔｒａｎｓｐｏｒｔ［Ｊ］．Ｎｕｃｌｅａｒ Ｉｎｓｔｒｕｍｅｎｔｓ ａｎｄ Ｍｅｔｈｏｄｓ ｉｎ Ｐｈｙｓｉｃｓ Ｒｅｓｅａｒｃｈ Ｂ，１９９７，１３２（３）：３７７—３９０． ［１４］唐刚华．小动物ＰＥＴ及其在医药学研究中的应用 ［Ｊ］．中国药理学通报，２００４，２０（５）：４８６—４９０． ［１５］Ｓ Ｊａｎ，Ｇ Ｓａｍｉｎ，Ｄ Ｓｔｒｕｌ，ｅｔ ａ１．ＧＡＴＥ：ａ ｓｉｍｕｌａｔｉｏｎ ｔｏｏｌｋｉｔ ｏｆｒ ＰＥＴ ａｎｄ ＳＰＥＣＴ［Ｊ］．Ｐｈｙｓｉｃｓ ｉｎ Ｍｅｄｃｉｎｅ ａｎｄ Ｂｉｏｌｏｇｙ，２００４，４９（１９）：４５４３—４５６１． ［１６］Ｒｅｎｅ Ｂｒｕｎ，Ｆｏｎｓ Ｒａｄｅｍａｋｅｒｓ．ＲＯＯＴ—Ａｎ ｏｂｊｅｃｔ ｏｒｉｅｎｔｅｄ ｄａｔａ ａｎａｌｙｓｉｓ ｆｒｍａｅｗｏｒｋ［Ｊ］．Ｎｕｃｌｅａｒ Ｉｎｓｔｒｕ— ｍｅｒｉｔｓ ａｎｄ Ｍｅｔｈｏｄｓ ｉｎ Ｐｈｙｓｉｃｓ Ｒｅｓｅａｒｃｈ Ａ，１９９７，３８９ （１—２）：８１—８６． ［１７］ＳＴＩＲ Ｕｓｅｒ’ｓ Ｇｕｉｄｅ Ｖｅｒｓｉｏｎ ２．０． 『１８１ ＰｅｎｅｌｏＰＥＴ Ｖ １．１ Ｕｓｅｒ Ｍａｎｕａ１． ［１９］Ｎａｔｉｏｎａｌ Ｅｌｅｃｔｉｒｃａｌ Ｍａｎｕｆａｃｔｕｒｅｓｒ Ａｓｓｏｃｉａｔｉｏｎ ２００１ ＮＥＭＡ Ｓｔｎａｄａｒｄｓ Ｐｕｂｌｉｃａｔｉｏｎ ＮＵ２—２００１：Ｐｅｒｆｏｒｍ— ａｎｃｅ Ｍｅａｓｕｒｅｍｅｎｔｓ ｏｆ Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｓ （Ｗａｓｈｉｎｇｔｏｎ ＤＣ：ＮＥＭＡ） ［２０］Ｙｕｅｈｕａｎ Ｗａｎｇ，Ｓｅｉｄｅｌ Ｊ，Ｔｓｕｉ Ｂ．Ｍ．Ｗ，ｅｔ １ａ．Ｐｅｒ－ ｆｏｒｍａｎｃｅ Ｅｖａｌｕａｔｉｏｎ ｏｆ ｔｈｅ ＧＥ Ｈｅｌｔｈｃａｒｅ ｅＸｐｌｏｒｅ ＶＩＳ— ＴＡ Ｄｕａｌ—Ｒｉｎｇ Ｓｍａｌｌ—Ａｎｉｍａｌ ＰＥＴ Ｓｃａｎｎｅｒ［Ｊ］，Ｊ． Ｎｕｃ１．Ｍｅｄ．，２００６，４７（１１）：１８９１—１９００． ［２１］Ｓ．Ｅｓｐａｎａ，Ｊ．Ｌ．Ｈｅｒｒａｉｚ，ｅｔ ａｌ，Ｖａｌｉｄａｔｉｏｎ ｏｆ ＰｅｎｅｌｏＰＥＴ Ａｇａｉｎｓｔ Ｔｗｏ Ｓｍａｌｌ Ａｎｉｍａｌ ＰＥＴ Ｓｃａｎｎｅｒｓ ［Ｃ］．２００７ ＩＥＥＥ Ｎｕｃｌｅａｒ Ｓｃｉｅｎｃｅ Ｓｙｍｐｏｓｉｕｍ Ｃｏｆｎｅｒ— ｅｎｃｅ Ｒｅｃｏｒｄ，２００７． （下转第１１６２页，Ｃｏｎｔｉｎｕｅｄ ｏｎ ｐａｇｅ １１６２）
