Pet系统符合探测响应的生成方法
【专利摘要】本发明提供了一种PET系统符合探测响应的生成方法。S1.利用一曲线方程拟合得到不同入射角度下单光子响应的拟合曲线；S2.根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射角度的方程，从而获取不同入射角度下单光子入射的方程；S3.获取多个沿多个指定响应线入射的两个γ光子的符合响应；S4.利用模拟得到的符合响应以及不同入射角度下单光子入射的方程二维曲面拟合得到所述多个指定响应线的曲面方程；S5.对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响应线入射的两个γ光子入射角度的二维曲面方程，从而获取单光子响应与符合探测响应的关系方程。本发明可以提升PET系统成像质量。
【专利说明】PET系统符合探测响应的生成方法
【技术领域】
[0001] 本公开涉及核医学成像【技术领域】，具体涉及一种由单光子入射响应符合生成PET 系统符合探测响应的方法。
【背景技术】
[0002] 核医学成像技术是核科学与核技术在医学诊断中最为典型和直接的应用。核医学 成像技术的成像真实，可以完整直接的显示细胞或分子水平的生理和病理过程，能够对肿 瘤、心血管、神经系统等许多疾病进行分子水平的早期诊断，不仅代表了医学影像学的发展 方向，而且具有广阔的市场前景。
[0003] PET (Positron Emission Computed Tomography，正电子发射断层扫描）系统是核 医学成像【技术领域】的一种重要设备，其可有效的探测到注入生物体内的标记了放射性示踪 剂的时空分布。由于生物体的病灶组织（如肿瘤）的代谢相对于正常细胞较旺盛，示踪剂 在这些组织的分布也较正常细胞多；所以PET系统探测到的这些示踪剂的分布可以提供携 带生物体功能信息的图像（如肿瘤的相关信息），为疾病的早期诊断和研究提供有力的参 考。
[0004] PET系统可以具有各种不同的几何形状，例如平板PET系统、环形PET系统等等。而 其形状主要由PET系统中每个探测器模块的不同排列而形成，但不同形状PET系统的探测 原理基本相同。例如，均是在进行扫描前，为生物体注射含有正电子放射性核素的示踪剂， 示踪剂中每一个正电子可与生物体内的一个负电子发生湮没生成两个背对背的Y光子即 一个Y光子对，这两个Y光子穿过生物体组织，几乎同时打在正电子扫描仪的一对探测器 模块上；探测器模块主要由闪烁晶体、光电倍增管以及前端逻辑电路部分组成。当Y光子 首先进入闪烁晶体后，闪烁晶体将高能光子信号转化为低能光子，低能光子通过光电倍增 管转化为电信号并进行放大，由前端逻辑电路输出至计算机，从而被记录下来；两个探测器 模块同时输出的电信号记录下有关上述Y光子对就是一次淹没事件的信息。这样所有的 Y光子对都被记录下来，然后经过图像重建，即可获得反映放射性示踪剂分布的图像。上述 过程即被称作探测符合。
[0005] 现有技术中，通常用一些简单的模型来描述上述探测符合响应情况，可以方便地 分析PET系统的响应。例如，由于湮没事件产生的两个γ光子是背对背的，可以近似为一条 直线，于是可以将上述所记录的这两个探测器模块输出的信号决定的一次湮没事件假想为 一条连接这两个探测器模块的闪烁晶体的直线，这条线被称为响应线（Line of Response, L0R)。
[0006] 在PET系统采集过程中，γ光子入射到探测器模块的闪烁晶体10上会产生晶体 间的穿透作用和散射作用，具体可以如图1以及图2中所示。这样所造成的深度效应严重 影响了 PET系统的图像质量；具体表现为图像分辨率、图像对比度以及信噪比等都会明显 的下降。
[0007] 目前，针对上述问题的主要解决方案是尽可能精确的得到PET系统探测符合响 应，然后引入后续的重建算法来较精确的恢复图像。获取PET系统探测符合响应的方法主 要有解析计算法、蒙特卡罗模拟法以及实验测量法。其中蒙特卡罗模拟和实验测量方法可 以较为准确的得到PET系统的响应函数。现有技术中，有些方法利用蒙特卡罗模拟法或者 实验测量法直接得到PET系统符合探测后的特殊点的响应，然后根据一些数学方法和PET 系统的对称性来得到整个视野内的PET系统响应情况；本公开所提供的方法中通过蒙特卡 罗模拟实验或实际实验测得的数据可重复利用率低，而且测量过程费时费力，对于不同形 状的PET系统需要重复测量。
[0008] 另有一些方法，利用不同角度入射闪烁晶体的单个Y光子入射闪烁晶体的响应 情况，将符合探测信号涉及的两个Y光子的这种单光子响应信息符合成PET系统所要的符 合探测响应信息，再根据探测器的对称性，得到整个视野中的PET系统响应情况。本公开所 提供的方法的关键在于将两个单光子响应信息符合生成符合探测的响应信息的过程是否 准确合理。
[0009] 上述现有方案的方法，对于连续的闪烁晶体而言，将两个Y光子在每个微分点的 响应进行乘积的积分。而对于常见的像素闪烁晶体而言，则是将上述积分离散化运算，即两 个Y光子入射处产生响应的闪烁晶体的响应进行两两相乘；但是本公开所提供的方法的 解析性太强，而PET系统合探测的过程又非常复杂，受到物理因素的影响较多，所以这种方 法的准确性欠佳。
【发明内容】
[0010] 针对现有技术中的部分或者全部问题，本公开提供一种PET系统符合探测响应的 生成方法，从而使得符合过程更加细致、精确，符合结果更加准确、可信，进而提高最终PET 系统成像的图像质量。
[0011] 本公开的其他特性和优点将通过下面的详细描述变得显然，或部分地通过本公开 的实践而习得。
[0012] 根据本公开的一个方面，一种PET系统符合探测响应的生成方法，包括：
[0013] 一种PET系统符合探测响应的生成方法，其特征在于，包括：
[0014] S1.利用一曲线方程拟合得到不同入射角度下单光子响应的拟合曲线；
[0015] S2.根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射 角度的方程，从而获取不同入射角度下单光子入射的方程；
[0016] S3.获取多个沿多个指定响应线入射的两个γ光子的符合响应；
[0017] S4.利用模拟得到的符合响应以及不同入射角度下单光子入射的方程二维曲面拟 合得到所述多个指定响应线的曲面方程；
[0018] S5.对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响 应线入射的两个Y光子入射角度的二维曲面方程，从而获取单光子响应与符合探测响应 的关系方程。
[0019] 在本公开的一种示例实施方式中，所述步骤S1包括：
[0020] 对所述PET系统中探测器模块的闪烁晶体进行编号；
[0021] 设所述曲线方程为Ch, a2, a3, a4......an ;x);其中a2, a3, a4…an为所述曲线方 程的系数，X为单光子入射的闪烁晶体编号；
[0022] 根据所述曲线方程，对不同入射角度下单光子响应拟合得到拟合曲线C(an，a 21， a3i，a41......anl ;x)......C(alm, a2m, a3m, a4m......;x);其中，m 为拟合曲线的个数。
[0023] 在本公开的一种示例实施方式中，所述拟合曲线的个数多于所述曲线方程系数的 个数。
[0024] 在本公开的一种示例实施方式中，所述步骤S2包括：
[0025] 根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射角 度 Θ 的方程 4(0),4( Θ)，Α3(Θ)，Α4(Θ)，……Αη(θ);
[0026] 得到不同入射角度下单光子入射的方程…… Αη( Θ ) ;x]，记作 C[ θ ;χ]。
[0027] 在本公开的一种示例实施方式中，所述步骤S3包括：
[0028] 所述步骤S3中，利用蒙特卡罗仿真模拟获取所述多个沿多个指定响应线入射的 两个Υ光子的符合响应。
[0029] 在本公开的一种示例实施方式中，所述多个指定响应线为处于不同位置的响应 线。
[0030] 在本公开的一种示例实施方式中，所述步骤S4包括：
[0031] 对于任一所述指定响应线，设所述曲面方程为：
[0032] SIX, b2，b3, b4......bn ;& ( Θ i ;xa)，C2 ( Θ 2 ;xb)];其中 b2，b3, b4-bn 为所述曲面方 程的系数；所述xa，xb为所述指定响应线对应的两个闪烁晶体的编号；
[0033] 将所述步骤S2中得到的方程C[ θ ;x]代入所述曲面方程得到所述指定响应线的 曲面方程为S(xa，xb);
[0034] 对于每一所述指定响应线，用所述曲面方程S(xa，xb)以^，x b变量进行二维曲面拟 合，得到各指定响应线的曲面方程SlXi, b2i，b3i, bybni ;& ( Θ i ;xai), c2( Θ 2 ;xbi)];其中i 为各所述指定响应线的编号。
[0035] 在本公开的一种示例实施方式中，所述步骤S5包括：
[0036] 对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响应 线入射的两个Y光子入射角度0 a，的曲面方程也（93, 0b)，B2(0a，0b)，B3(0 a，0b)， B4(ea，0b)，……B n(0a，0b);
[0037] 得到单光子响应与符合探测响应的关系方程为：
[0038] Θ^,Β^Θ^ 0^,63(0^ θ2)……θ2) 4(0。 X)，C2( θ 2 ;χ2)]
[0039] 在本公开的一种示例实施方式中，所述曲线拟合以及曲面拟合通过Matlab软件 中的拟合数据包实现。
[0040] 在本公开的实施例所提供PET系统符合探测响应的生成方法中，通过精准的蒙卡 模拟实验，引入PET系统符合探测的多个的先验数据，使得符合过程更加细致、精确，符合 结果更加准确、可信，从而提高了最终PET系统成像的图像质量。
【专利附图】
【附图说明】
[0041] 通过参照附图详细描述其示例实施方式，本公开的上述和其它特征及优点将变得 更加明显。
[0042] 图1是入射Y光子在闪烁晶体间的穿透作用示意图；
[0043] 图2是入射Υ光子在闪烁晶体间的散射作用示意图；
[0044] 图3是本公开一种实施方式中PET系统符合探测响应的生成方法的流程示意图；
[0045] 图4A是不同角度下Y光子的入射响应示意图；
[0046] 图4B是不同角度下的单光子响应示意图；
[0047] 图5是Y光子沿不同响应线入射的示意图；
[0048] 图6是Y光子沿图5中某一响应线入射的放大示意图；
[0049] 图7是单光子入射响应在各闪烁晶体中的不对称分布示意图；
[0050] 图8是高斯分段分布函数曲线示意图；
[0051] 图9是响应线组织方式示意图；
[0052] 图10是Y光子沿所述指定响应线两端的闪烁晶体入射角度示意图。
[0053] 附图标记说明：
[0054] 10:闪烁晶体
[0055] 20:响应线
[0056] S1-S5 :步骤
【具体实施方式】
[0057] 现在将参考附图更全面地描述示例实施方式。然而，示例实施方式能够以多种形 式实施，且不应被理解为限于在此阐述的实施方式；相反，提供这些实施方式使得本公开将 全面和完整，并将示例实施方式的构思全面地传达给本领域的技术人员。在图中相同的附 图标记表示相同或类似的结构，因而将省略它们的详细描述。
[0058] 此外，所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施 例中。在下面的描述中，提供许多具体细节从而给出对本公开的实施例的充分理解。然而， 本领域技术人员将意识到，可以实践本公开的技术方案而没有所述特定细节中的一个或更 多，或者可以采用其它的方法、组元、材料等。在其它情况下，不详细示出或描述公知结构、 材料或者操作以避免模糊本公开的各方面。
[0059] 本示例实施方式中提供了一种由单光子入射响应符合生成PET系统符合探测响 应的方法。如图3中所示，该方法可以包括以下步骤：
[0060] S1.利用一曲线方程拟合得到不同入射角度下单光子响应的拟合曲线；
[0061] S2.根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射 角度的方程，从而获取不同入射角度下单光子入射的方程；
[0062] S3.获取多个沿多个指定响应线入射的两个Y光子的符合响应；
[0063] S4.利用模拟得到的符合响应以及不同入射角度下单光子入射的方程二维曲面拟 合得到所述多个指定响应线的曲面方程；
[0064] S5.对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响 应线入射的两个Y光子入射角度的二维曲面方程，从而获取单光子响应与符合探测响应 的关系方程。
[0065] 在该PET系统符合探测响应的生成方法中，通过引入PET系统符合探测的多个的 先验数据，使得符合过程更加细致、精确，符合结果更加准确、可信，从而提高了最终PET系 统成像的图像质量。
[0066] 下面对上述步骤S1-S5分别加以详细的说明。
[0067] S1.利用一维曲线方程拟合得到不同入射角度下单光子响应的拟合曲线；本示例 实施方式中，该步骤可以如下：
[0068] 该步骤中主要输入为已经得到的单光子在不同角度下入射到闪烁晶体阵列的响 应。在PET系统探测器模块的闪烁晶体大小和材料一定的情况下，γ光子在闪烁晶体上穿 透的数目和响应情况主要由入射角度决定，例如图4A中所示，入射角度不同，得到的响应 也不同。一般而言，入射角度越大，其穿透的闪烁晶体10数目越多，响应情况越复杂。图4B 为三个入射角度下的响应曲线20 ;为方便计算，本示例实施方式中对所述PET系统中探测 器模块的闪烁晶体10进行了编号；在图4B中，编号为7的闪烁晶体即为入射闪烁晶体。 [0069] 用相同的一维曲线方程对不同入射角度的这种单光子响应进行拟合，可以得出不 同入射角度下的多个拟合曲线。设曲线方程为：
[0070] C(a1； a2, a3, a4......an ；x)；
[0071] 其中ai，a2, a3，a4?an为所述曲线方程所涉及至I撤系数，x为单光子入射的闪烁晶 体10编号，其为该曲线方程的变量。
[0072] 根据所述曲线方程，对于不同入射角度下单光子响应可以拟合得到多个拟合曲 线·
[0073] C (Sin，&21，&31，&41 ······ &nl ;X)，
[0074] C (a12, a22，已32，已42......an2 ;x)，
[0075] C (a13, a23，a33, a43......an3 ;x)，
[0076] C (&14，&24，&34，&44^ · " · · &n4 ;X)，
[0077] ......
[0078] C (alm, a2m，a3m, a4m......a· ;x);
[0079] 其中，m为拟合曲线的个数；为了方便计算，一般要求所述拟合曲线的个数多于所 述曲线方程系数的个数。
[0080] S2.根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射 角度的方程，从而获取不同入射角度下单光子入射的方程；本示例实施方式中，该步骤可以 如下：
[0081] 对于步骤1中得到的一系列曲线方程C(an, a21，a31, a41......anl ;x), C(a12, a22, a32> a42......an2 ；x)?C (a13, a23? a33, a43......an3 ；x)?C (a14, a24j a34, a44......an4 ；x),......C (alm, a2m? a3m, a4m......;x)，其中的系数ap a2，a3, a4......an可以看作是入射角Θ的离散分布，即ai? Θ , a2 ?Θ , a3 ?Θ , a4 ?Θ ,......, an ?Θ。
[0082] 对所述曲线方程的各系数进行拟合得到各系数关于入射角度Θ的方程：
[0083] Ai( θ ),Α2( θ ),Α3( θ ),Α4( θ ),……Αη(θ);
[0084] 这样，则可以得到不同入射角度下单光子入射的方程：
[0085] 0^(0),4(0)^3(0)^4( Θ ),……Αη(θ) ;x];
[0086] 所以C还可以记作：
[0087] C[ θ ；χ] 0
[0088] S3.获取多个沿多个指定响应线20入射的两个Υ光子的符合响应；本示例实施 方式中，该步骤可以如下：
[0089] 本方法的关键技术之一在于能够引入符合探测先验性信息，使得符合过程更加可 靠。因此本示例实施方式中可以利用蒙特卡罗仿真模拟得到沿不同响应线20入射的两个 Y光子在PET系统中的符合响应的情况，如图5中所示。本示例实施方式中，仿真得到的符 合响应是以两端闪烁晶体编号进行记录。图6显示了沿其中一条响应线20入射的放大图。 所以Y光子沿每一条响应线20方向入射的响应均可看作为关于两端闪烁晶体编号x a，xb 的散点分布。在选择所述指定响应线时，引入的先验数据越多越好，但是考虑到实验条件， 可以尽量能选取差别较大的具有代表性的响应线20,即可以取处于不同位置的响应线20 来进行仿真模拟。
[0090] S4.利用模拟得到的符合响应以及不同入射角度下单光子入射的方程二维曲面拟 合得到所述多个指定响应线的曲面方程；本示例实施方式中，该步骤可以如下：
[0091] 对于任一所述指定响应线，将拟合的二维曲面方程设为以两个单光子入射角度响 应为变量的方程：
[0092] S [b" b2, b3, b4......bn ;& ( Θ i ;xa)，C2 ( Θ 2 ;xb)];
[0093] 其中汕1，132，133，13 4?比是相应的参数列，即为所述曲面方程的系数；(：1(01; xa)，C2( Θ 2 ;Xb)即步骤S3中得到的方程，是两个对应的变量。对于每个所述指定响应线，Θ i 和Θ 2都是确定值，Xa，Xb分别为Q C2的变量。在步骤S3中提到，蒙特卡罗仿真模拟得出的 符合响应是以两端闪烁晶体即闪烁晶体对的编号记录的，所以Y光子沿每一条响应线20 方向入射的响应均可看作关于两端闪烁晶体编号x a，xb的散点分布。
[0094] 将步骤S3中得到的方程代入所述曲面方程，可以得到xa，xb为变量的方程，即
[0095] S (xa, xb)
[0096] 此时由于对于每个指定响应线，Θ i和Θ 2都是给定的值，所以 〇[^(0)，4(0)^3(0),^(0),……Α η(θ) ;χ]中的 4(0),4( Θ)，Α3(Θ)，Α4(Θ),…… Αη( Θ )都是确定的值，因此，此时拟合的方程中只有t^，b2, b3, b4*"bn是未知量。
[0097] 对于每一所述指定响应线，通过上述的曲面方程SIX, b2, b3, b4......bn ；〇! ( Θ j ； xa)，C2(02;Xb)]以Xa，Xb为变量进行二维曲面拟合，可以得到一列各指定响应线的曲 面方矛王：S [bn, b21，b31, b41 …bnl S[b12, b22，b32, b42 …;C^ ( Θ 丄； xi)，C2 ( Θ 2 ;x2) ]，S [b13, b23，b33, b43…bn3 ;C^ ( Θ i ，C2 ( Θ 2 ;x2)]......。
[0098] S5.对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响 应线入射的两个Y光子入射角度的二维曲面方程，从而获取单光子响应与符合探测响应 的关系方程；本示例实施方式中，该步骤可以如下：
[0099] 在步骤 S4 得到的一列曲面方程 S[bn, b21，b31, b41 …bnl ;& ( Θ i ;Xl), C2( θ 2 ;χ2)]， S [b12, b22，b32, b42…bn2 ;C^ ( Θ i :叉），C2 ( θ 2 ;χ2) ]，S [b13, b23，b33, b43…bn3 ;C^ ( Θ i :叉），C2 ( θ 2 ; x2) 1......中系数bp b2，b3, b4......bn都各不相同，但是对于每个具体的指定响应线，两个γ光 子的入射角度Θ^Ρ θ2都是给定的，于是bi，b2，b3，b4……b n可以看作为Θ^Ρ 02的散点分 布；即 h ?（θ Θ 2)，b2 ?（θ Θ 2)，b3 ?（θ Θ 2)，b4 ?（θ Θ 2)，......bn?（θ Θ 2)。 对各系数进行二维曲面拟合得到各系数关于所述沿指定响应线入射的两个Y光子入射 角度 ea，二维曲面方程：Bl(ea，0b)，B 2(0a，0b)，B3(0a，0 b)，B4(0a，0b)，…… Bn(9a，0b)。这样就找到了单光子响应与符合探测响应的关系方程为：
[0100] 0,),62(0^ Θ2), 63(0^ 0,),64(0^ Θ2)……θ2) 4(0。 X)，C2( θ 2 ;χ2)]
[0101] 这个关系方程可以直接作为每次符合过程的计算关系式，继而代入后续的图像重 建过程。
[0102] 值得注意的是，本公开所提供的方法中多次用到了拟合的过程，拟合方程的选取 是关键点之一。拟合方程选取的越精确，符合过程也就越精确。本示例实施方式中，曲线拟 合以及曲面拟合可以通过Matlab软件中的拟合数据包进行实现。
[0103] 上述方法过程在具有相同规格闪烁晶体的各种PET系统中只需进行一次即可，得 到的关系方程与PET系统状无关。另外，本公开所提供的方法在引入符合先验性条件的时 候，只需用特定的一个PET系统进行即可，因为此方法只在于得出符合过程相对于两端入 射角的关系，而不针对PET系统的形状。
[0104] 进一步的，本示例实施方式中还提供了一种上述方法的具体实施例；其具体可以 如下所述：
[0105] 结合单光子入射在各闪烁晶体的分布特点，如图7中所示，当入射角度越大时，其 两边呈不对称分布。所以本示例实施方式中可以采用分段拟合的方法，当然其也可以采用 其他的曲线公式进行拟合。以采用分段一维高斯曲线进行拟合为例，其中以峰值为分界线 的分段拟合曲线示意图如图8所示，其公式如下，：
[0106]
【权利要求】
1. 一种PET系统符合探测响应的生成方法，其特征在于，包括：
51. 利用一曲线方程拟合得到不同入射角度下单光子响应的拟合曲线；
52. 根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射角度 的方程，从而获取不同入射角度下单光子入射的方程；
53. 获取多个沿多个指定响应线入射的两个γ光子的符合响应；
54. 利用模拟得到的符合响应以及不同入射角度下单光子入射的方程二维曲面拟合得 到所述多个指定响应线的曲面方程；
55. 对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响应线 入射的两个Y光子入射角度的二维曲面方程，从而获取单光子响应与符合探测响应的关 系方程。
2. 根据权利要求1所述的PET系统符合探测响应的生成方法，其特征在于，所述步骤 S1包括： 对所述PET系统中探测器模块的闪烁晶体进行编号； 设所述曲线方程为C(ai，a2, a3，a4……an ;x);其中&1，a2, a3, a^·· an为所述曲线方程的 系数，x为单光子入射的闪烁晶体编号； 根据所述曲线方程，对不同入射角度下单光子响应拟合得到拟合曲线C(an，a21， a3i，a41......anl ;x)......C(alm, a2m, a3m, a4m......;x);其中，m 为拟合曲线的个数。
3. 根据权利要求2所述的PET系统符合探测响应的生成方法，其特征在于，所述拟合曲 线的个数多于所述曲线方程系数的个数。
4. 根据权利要求2所述的PET系统符合探测响应的生成方法，其特征在于，所述步骤 S2包括： 根据各所述拟合曲线对所述曲线方程的各系数进行拟合得到各系数关于入射角度Θ 的方程 AJ θ )，Α2( θ )，Α3( θ )，Α4( Θ )，......Αη( θ ); 得到不同入射角度下单光子入射的方程： 〇[^(0)，4(0)^3(0),^(0),……Αη(θ) ;x],记作 C[0 ;x]。
5. 根据权利要求4所述的PET系统符合探测响应的生成方法，其特征在于，所述步骤 S3包括： 所述步骤S3中，利用蒙特卡罗仿真模拟获取所述多个沿多个指定响应线入射的两个 Y光子的符合响应。
6. 根据权利要求1所述的PET系统符合探测响应的生成方法，其特征在于，所述多个指 定响应线为处于不同位置的响应线。
7. 根据权利要求4所述的PET系统符合探测响应的生成方法，其特征在于，所述步骤 S4包括： 对于任一所述指定响应线，设所述曲面方程为： SIX, b2，b3, b4......bn A ( Θ i ;xa), c2( Θ 2 ;xb)];其中 bi, b2，b3, V"bn为所述曲面方程的 系数；所述xa，xb为所述指定响应线对应的两个闪烁晶体的编号； 将所述步骤S2中得到的方程C[ θ ;x]代入所述曲面方程得到所述指定响应线的曲面 方程为S (xa，xb); 对于每一所述指定响应线，利用所述曲面方程S(xa，xb)以所述xa，x b变量进行二维曲面 拟合，得到各指定响应线的曲面方程S[bn, b2i，b3i, bybni ;&( Θ i ;Xai)，C2( θ 2 ;Xbi)];其中 i为各所述指定响应线的编号。
8. 根据权利要求7所述的PET系统符合探测响应的生成方法，其特征在于，所述步骤 S5包括： 对各所述曲面方程的各系数进行二维曲面拟合得到各系数关于所述沿指定响应线入 射的两个Y光子入射角度9a，0b的曲面方程：0b)，B2(0 a，0b)，B3(0a，0b)， B 4(ea，0b)，……Bn(0a，0 b); 得到单光子响应与符合探测响应的关系方程为： Θ^,Β^Θ^ 0^,3(9^ θ2)……Βη( Θ Θ 2) ;Ci( Θ i ; Χ1)，。2 ( 9 2 ;X2)]。
9. 根据权利要求1-8任意一项所述的PET系统符合探测响应的生成方法，其特征在于， 所述曲线拟合以及曲面拟合通过Matlab软件中的拟合数据包实现。
【文档编号】A61B6/03GK104146726SQ201410409253
【公开日】2014年11月19日 申请日期:2014年8月19日 优先权日:2014年8月19日
【发明者】樊馨, 曹学香, 刘双全, 孙校丽, 贠明凯, 孙翠丽, 周小林, 王海鹏, 魏存峰, 刘宝东, 李道武, 魏龙 申请人:中国科学院高能物理研究所