一种图像重建方法、装置及设备与流程

本申请涉及医疗设备技术领域，特别涉及一种图像重建方法、装置及设备。

背景技术：

正电子发射计算机断层显像(PET，positron emission tomography)设备是当前公认最先进的大型医用科研和临床诊断设备之一，在诊断肿瘤，心血管病等疾病方面有着卓越的性能。在PET设备采集病人数据的过程中，由于PET设备的多边形几何结构会导致晶体对的灵敏度(效率)不一致，使得采集的数据数目与病人体内实际发出的数据数目存在差别，进而重建出来的图像无法正确的反映病人体内放射性元素活度的分布情况。

为了解决该问题，现有技术中提出了对PET设备采集的数据进行几何校正。来保证晶体对灵敏度一致性，具体为：在PET设备视野中心放置一个(外置的)均匀的放射源装置，采集放射源(比如旋转棒源、平面源、环形源和移动线源等)发出的光子对，然后借助于采集的数据，设计相应的校正算法，计算每对晶体对的几何校正因子。

现有技术中，虽然能通过校正晶体对的几何因子，保证晶体对灵敏度的一致性来提高图像重建的准确性，但是，现有的这种方式中，一方面，需要放置发射源装置，而外置放射源装置设计复杂，放射源不容易存储、屏蔽，且价格昂贵；另一方面，校正过程需要工作人员操纵放射源装置，而放射源发出的射线会对工作人员的身体造成辐射危险。

技术实现要素：

有鉴于此，本申请提供一种图像重建方法、装置及设备，以通过降低PET设备中几何因子校正成本，以及校正几何因子所带来的辐射危险，来提高图像重建的准确性。

具体地，本申请是通过如下技术方案实现的：

第一方面，提供一种图像重建方法，所述方法包括：

确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

获取通过所述晶体对扫描患者的数据；

利用所述晶体对的灵敏度对所述数据进行校正；

根据校正后的数据进行图像重建。

第二方面提供一种图像重建装置，所述装置包括：

第一确定单元，用于确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

第二确定单元，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

获取单元，用于获取通过所述晶体对扫描患者的数据；

校正单元，用于利用所述晶体对的灵敏度对所述数据进行校正；

重建单元，用于根据校正后的数据进行图像重建。

第三方面提供一种医疗设备，包括：处理器；以及用于存储所述处理器可执行指令的存储器；其中，

所述处理器被配置为：

确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

获取通过所述晶体对扫描患者的数据；

利用所述晶体对的灵敏度对所述数据进行校正；

根据校正后的数据进行图像重建。

本申请实施例中，先确定晶体对中单个晶体的晶体效率因子和几何因子，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，然后，利用晶体对的灵敏度对获取到患者的数据进行校正，并根据校正后的数据进行图像重建。也就是说，本申请实施例中，根据单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，并利用该灵敏度对扫描的数据进行有效校正，提高了图像重建的准确性。另外，通过对晶体的本底辐射进行校正，无需购买放射源、制作放射源装置，简洁高效的计算出了晶体对的几何因子，不但降低了成本，还避免了操作人员的辐射危险。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。

附图说明

图1是本申请实施例提供的一种图像重建方法的流程图；

图2是本申请实施例提供的一种确定晶体对的灵敏度的流程图；

图3是本申请实施例提供的一种晶体本底符合事件的示意图；

图4是本申请实施例提供的一种计算晶体效率因子的示意图；

图5是本申请实施例提供的晶体对图像重建装置所在图像重建设备的一种硬件结构图；

图6是本申请实施例提供的一种图像重建装置的结构示意图；

图7是本申请实施例提供的一种图像重建装置的另一结构示意图；

图8是本申请实施例提供的一种图像重建装置的另一结构示意图；

图9是本申请实施例提供的一种图像重建装置的另一结构示意图；

图10是本申请实施例提供的一种图像重建装置的另一结构示意图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。

在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。

请参阅图1，是本申请实施例提供的一种图像重建方法的流程图；所述方法包括：

步骤101：确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

该步骤中，一种确定晶体对中单个晶体的晶体效率因子和几何因子的过程为：

先采集晶体对上的多个本底符合事件；将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件；然后，确定第一类本底符合事件的第一晶体效率因子，以及第二类本底符合事件的第二晶体效率因子；并利用所述第一晶体效率因子对所述第一类本底符合事件进行校正，得到校正后第一类本底符合事件，以及利用所述第二晶体效率因子对所述第二类本底符合事件进行校正，得到校正后第二类本底符合事件；最后，利用每对晶体对上校正后的第一类本底符合事件与第二类本底符合事件的和，确定所述晶体对上单个晶体对应的第一几何因子和第二几何因子。

其确定晶体对中单个晶体的晶体效率因子和几何因子的详细过程详见下述图2至图4的描述，在此不再赘述。

步骤102：根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

该步骤中，将步骤101中得到的所述第一晶体效率因子、第二晶体效率因子、第一几何因子和第一几何因子分别相乘，得到晶体对的灵敏度。也可以称为晶体对的灵敏度模型。

也就是说，晶体对i-j的灵敏度模型可以通过公式(1)建立：

η_ij＝ε_iε_jg_ig_j (1)

其中，η_ij表示晶体对i-j的灵敏度，i-j分别表示晶体对中单个晶体，ε_i(或ε_j)表示的单个晶体i(或j)的晶体效率因子，g_i(或g_j)表示单个晶体i(或j)的几何因子，反映的是单个晶体的倾斜角对其接受单个光子的影响。

步骤103：获取通过所述晶体对扫描患者的数据；

该步骤中，在使用PET设备扫描患者时，PET设备获取通过晶体对扫描患者的数据。

步骤104：利用所述晶体对的灵敏度对所述数据进行校正；

该步骤中，利用晶体对的灵敏度对所述数据进行校正的过程对于本领域技术人员来说，已是熟知技术，在此不再赘述。

步骤105：根据校正后的数据进行图像重建。

该步骤中，根据校正后的数据进行图像重建过程对于本领域技术人员来说，已是熟知技术，在此不再赘述。

本申请实施例中，先确定晶体对中单个晶体的晶体效率因子和几何因子，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，然后，利用晶体对的灵敏度对获取到患者的数据进行校正，并根据校正后的数据进行图像重建。也就是说，本申请实施例中，根据单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，并利用该灵敏度对扫描的数据进行有效校正，提高了图像重建的准确性。另外，本实施例通过对晶体的本底辐射进行校正，无需购买放射源、制作放射源装置，简洁高效的计算出了晶体对的几何因子，不但降低了成本，还避免了操作人员的辐射危险。

请参阅图2，是本申请实施例提供的一种晶体对几何因子的校正方法的流程图；所述方法包括：

步骤201：采集晶体对上的多个本底符合事件；

该步骤中，PET设备收集晶体对上的多个本底符合事件，每个本底符合事件包括：第一类本底符合事件和第二类本底符合事件。

该实施例中，PET设备中包括晶体环，其中，晶体环上的每个晶体可以是具有本底辐射的晶体，如含镥(Lu)晶体(比如LYSO晶体等)，其中的Lu元素在自然的状态下会发生β-衰变，同时伴随着三种不同能量级别的γ光子的发出，分别为307KeV，202KeV，88KeV，衰变发出的β粒子在晶体中发生电离效应,释放出590Kev左右的能量，被其所在的晶体探测到，即β粒子会被晶体“吸收”。同时，衰变发出的γ光子逃离该晶体，在视野中飞行一段时间，被PET设备中的环形探测器中另一个晶体检测到，如图3所示，图3为本申请实施例提供的一种晶体本底符合事件的示意图。如图所示，202表示晶体i中的Lu元素衰变时，β元素被晶体i“吸收”，γ光子被晶体j接收到；201表示与202相反的过程，m，n表示晶体环上的另外两个晶体。

为了简洁、方便的处理Lu元素衰变发出的三中不同能量级别的γ光子，通过设置能量阈值的方法去掉能量为202kev和88kev的γ光子，只收集307kev的γ光子。

该实施例中，预先给PET设备设置一个能量阈值ΔE，一个符合时间窗ΔT，满足下述条件的事件称为本地辐射符合事件。

如果晶体对测得的两个粒子的能量满足能量阈值ΔE，且这两个粒子被晶体探测到的时间间隔在符合事件窗ΔT内，则这两个粒子称为一个本底符合事件。每当晶体对i-j(其中，i和j分别表示晶体对的单个晶体)上有一个本底符合事件发生时，晶体对i-j在电路中对应的某个位置计数会加1，此时，PET设备收集上述本底符合事件，并将其保存。

步骤202：将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件；

该步骤中，一种划分方式为：按照能量阈值将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件，其中，能量阈值是预先设定的，按照能量阈值划分本底符合事件对于本领域技术人员来说，已是熟知技术，在此不再赘述。

步骤203：确定第一类本底符合事件的第一晶体效率因子，以及第二类本底符合事件的第二晶体效率因子；

该步骤中，由于采集到本底符合事件中一个是γ粒子，一个是β粒子，β粒子以概率1被单个晶体i(或j)吸收，具体如图2所示，而γ粒子从该单个晶体i(或j)中飞出打到对面单个晶体j(或i)上被单个晶体j(或i)接收，因此，本底符合事件反映的是单个晶体j(或i)接收γ粒子的效率(β粒子的作用在于借助于它确定出它对应的γ光子)。但是，在PET设备实际使用过程中，病人体内发出的是2个γ光子，其分别被两个晶体所吸收的符合事件，即，一个符合事件反映两个晶体接受光子的效率。

也就是说，任何一对晶体对在接收符合事件的时候，其效率既受到晶体固有效率的影响，也受到设备环形几何结构的影响，在计算晶体对的几何校正因子之前，先要去除晶体固有效率的影响。

因此，晶体效率因子的一种计算方式(本实施例以计算第一晶体效率因子为例)为：

由于图3展示了本底符合事件的特征可知，晶体对i-j采集到的本底符合事件(即数据，也就是说，数据是由多个本底符合事件组成的)有两种类型，即201和202，本文中，201称为第一类本底符合事件，202称为第二类本底符合事件。其中，201表示晶体i接收γ光子的本底符合事件，202表示晶体j接收γ光子的本底符合事件。假设晶体对i-j上采集到的本底符合事件总计数为n，则n中包含一定数目的201本底符合事件和一定数目的202本底符合事件。

由于β粒子与γ粒子能量不同，可以通过对晶体i探测到的粒子的能量阈值判断出晶体i探测到的粒子是β还是γ，即该本底符合事件是201还是202。假如，设n_ij为统计出的202本底符合事件的总数目，则201本底符合事件总数目为n_ij＝n-n_ji。

计算晶体i的固有效率，对晶体对i-j上收集的本底符合事件，本实施例以201这种类型的本底符合事件为例。根据采集到的201本底符合事件n_ij，利用扇束算法的来计算晶体对上的单个晶体i的效率校正因子(即第一晶体效率因子)，其扇束算法的公式为：

其中，N表示环形探测器上总的晶体数目，j_min表示扇束的起始晶体编号，j_max表示扇束的终止晶体编号，n_ij表示第一类本底符合事件的第一计数总数，ε_i表示第一晶体效率因子。

同理，统计第二本底符合事件的第二计数总数n_ji,,根据所述第二计数总数n_ji,利用扇束算法的来计算晶体对上的单个晶体j的效率校正因子(即第二晶体效率因子)，其扇束算法的公式为：

其中,N表示环形探测器上总的晶体数目，j_min表示扇束的起始晶体编号，j_max表示扇束的终止晶体编号，n_ji表示第二本底符合事件的第二计数总数，ε_j表示第二晶体效率因子。

其中，j_min和j_max具体如图4所示，图4为本申请实施例提供的一种计算晶体效率因子的示意图。

另外，晶体j的效率因子(即第二晶体效率因子)计算需要根据晶体对i-j上的202这种类型的本底符合事件，其它晶体效率因子的计算用的数据类型或者是与晶体i采用的本底符合事件类型相同或者是与晶体j的相同，其计算过程类似，在此不再赘述。

步骤204：利用所述第一晶体效率因子对所述第一类本底符合事件进行校正，得到校正后第一类本底符合事件，以及利用所述第二晶体效率因子对所述第二类本底符合事件进行校正，得到校正后第二类本底符合事件；

该步骤中，利用晶体效率因子分别对对应的本底符合事件进行校正，

由于得到晶体环上每个晶体的效率校正因子ε_i(i＝1～N)后，还需要利用效率校正因子ε_i对晶体对i-j上的本底符合事件进行晶体效率校正，其校正的目的是除去晶体固有效率带来的影响，然后再计算几何校正因子。

由于考虑到本底符合事件只受一个晶体的影响，对晶体对i-j上收集到本底符合事件的总计数n进行晶体效率校正时，需要分别对其包含的201和202事件用ε_i和ε_j单独校正，即是n'_ji＝n_ji*ε_j，其中n'_ji表示晶体效率校正后的202事件，n'_ji＝n_ji*ε_j，n’_ij＝n_ij*ε_i，此时晶体对i-j上的符合事件总数目为n'＝n’_ij+n'_ji。

步骤205：利用每对晶体对上校正后的第一类本底符合事件与第二类本底符合事件的和，确定所述晶体对上单个晶体对应的第一几何因子和第二几何因子；

该步骤中，与现有直接计算每对晶体对的几何因子的方法不同，本实施例考虑到本底符合事件的特征，提出了单独计算每个晶体的几何因子，然后再组合该几何因子，得到晶体对的几何因子。

由图3可知，假设沿着晶体对i-j的方向上共有M对晶体对(例如晶体对m-n是与晶体对i-j有相同方向的晶体对)，分别为：Xtal_i1Xtal_j1，Xtal_i2Xtal_j2，…Xtal_iMXtal_jM,每对晶体上晶体效率因子校正后的本底符合事件的总计数为n’₁，n'₂，…n'_M,则单个晶体Xtal_ik(或者Xtal_jk)(k＝1～N)的几何校正因子可以通过归一化算法得到，其中，一种归一化算法的公式为：

其中，由于晶体Xtal_jk(或者Xtal_ik)接收数据的时候受到的几何影响与Xtal_ik(或者Xtal_jk)相同，所以，g_jk的几何校正因子也可以通过归一化算法得到，其一种归一化算法的计算公式如公式(4)所示，其计算过程与g_ik的几何校正因子的计算过程类似，具体详见上述，在此不再赘述。

进一步，在上述实施例的基础上，还可以：组合第一几何因子和第二几何因子，得到晶体对的几何因子，也可以组合第一晶体效率因子和第二晶体效率因子，得到该晶体对的晶体效率因子，下面以组合第一几何因子和第二几何因子为例，其组合晶体效率因子与其类似。

该步骤中，一种组合方式是将所述第一几何因子乘以所述第二几何因子，得到晶体对的几何因子。

即，晶体对Xtal_ikXtal_jk(k＝1～N)的几何校正因子为：

g_ikjk＝g_ik*g_jk (5)

因此，根据公式(2A)、(2B)与(5)，晶体对i-j的灵敏度可以借助于上述建立的模型的公式(1)得到。

本申请实施例中，将采集到的多个本底符合事件划分为不同类型的本底符合事件，根据相应的本底符合事件分别计算相应的晶体效率因子，然后利用计算出的晶体效率因子校正之前对应的本底符合事件，并利用每对晶体对上校正后的本底符合事件计算单个晶体的几何因子，然后，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，利用晶体对的灵敏度对获取到患者的数据进行校正，并根据校正后的数据进行图像重建。也就是说，本申请实施例中，根据单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，并利用该灵敏度对扫描的数据进行有效校正，提高了图像重建的准确性。另外，本实施例通过对晶体的本底辐射进行校正，无需购买放射源、制作放射源装置，简洁高效的计算出了晶体对的几何因子，不但降低了成本，还避免了操作人员的辐射危险。

与前述图像重建方法的实施例相对应，本申请还提供了图像重建装置的实施例。

本申请提供的图像重建装置的实施例可以应用在图像重建设备上，该图像重建设备可以是上位机。如图5所示，为本申请实施例提供的图像重建装置所在图像重建设备的一种硬件结构图，该图像重建设备可以包括处理器501以及机器可读存储介质502，其中，处理器501和机器可读存储介质502通常借由内部总线503相互连接。在其他可能的实现方式中，所述设备还可能包括外部接口504，以能够与其他设备或者部件进行通信。进一步地，机器可读存储介质502上存储有图像重建的控制逻辑505，该控制逻辑505从功能上划分的逻辑模块，可以是图6所示的图像重建装置的结构。

请参考图6，为本申请实施例提供的一种图像重建装置的结构示意图，所述装置包括：第一确定单元61，第二确定单元62，获取单元63、校正单元64和重建单元65，其中，

第一确定单元61，用于确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

第二确定单元62，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

获取单元63，用于获取通过所述晶体对扫描患者的数据；

校正单元64，用于利用所述晶体对的灵敏度对所述数据进行校正；

重建单元65，用于根据所述校正单元校正后的数据进行图像重建。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述第一确定单元61包括：采集单元71，划分单元72，效率因子确定单元73，事件校正单元74和几何因子确定单元75，其结构示意图如图7所示，其中，

采集单元71，用于采集晶体对上的多个本底符合事件；

划分单元72，用于将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件；

效率因子确定单元73，用于确定第一类本底符合事件的第一晶体效率因子，以及第二类本底符合事件的第二晶体效率因子；

事件校正单元74，用于利用所述第一晶体效率因子对所述第一类本底符合事件进行校正，得到校正后第一类本底符合事件，以及利用所述第二晶体效率因子对所述第二类本底符合事件进行校正，得到校正后第二类本底符合事件；

几何因子确定单元75，用于利用每对晶体对上校正后的第一类本底符合事件与第二类本底符合事件的和，确定所述晶体对上单个晶体对应的第一几何因子和第二几何因子。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述划分单元72，具体用于按照能量阈值将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述效率因子确定单元73包括：统计单元731和计算单元732，其结构示意图如图8所示，其中，

统计单元731，用于统计第一类本底符合事件的第一计数总数和第二本底符合事件的第二计数总数；

计算单元732，用于根据所述第一计数总数和第二计数总数利用扇束算法分别计算第一晶体效率因子和第二晶体效率因子。

其中，可以根据第一计数总数利用公式计算第一晶体效率因子；以及根据所述第二计数总数利用公式计算第二晶体效率因子；

其中，N表示环形探测器上总的晶体数目，j_min表示扇束的起始晶体编号，j_max表示扇束的终止晶体编号,n_ij表示第一类本底符合事件的第一计数总数，n_ji表示第二本底符合事件的第二计数总数n_ji，ε_i表示第一晶体效率因子，ε_j表示第二晶体效率因子。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述事件校正单元74包括：第一校正子单元741和第二校正子单元742，其结构示意图如图9所示，其中，

第一校正子单元741，用于将所述第一晶体效率因子乘以所述第一类本底符合事件，得到校正后的第一类本底符合事件；

第二校正子单元742，用于将所述第二晶体效率因子乘以所述第二类本底符合事件，得到校正后的第二类本底符合事件。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述几何因子确定单元75包括：第一确定子单元751和第二确定子单元752，其结构示意图如图10所示，其中，

第一确定子单元751，用于根据校正后每对晶体对上的第一类本底符合事件与第二类本底符合事件的和，利用归一化算法计算所述晶体对上一个晶体的第一几何因子；

第二确定子单元752，用于根据校正后每对晶体对上的第一类本底符合事件与第二类本底符合事件的和，利用归一化算法计算所述晶体对上另一个晶体的第二几何因子。

可选的，在另一实施例中，该实施例在上述实施例的基础上，所述效率因子确定单元76还可以包括：第一组合单元和/或第二组合单元(图中未示)，其中，

第一组合单元，用于组合所述第一几何因子和第二几何因子，得到晶体对的几何因子；

第二组合单元，用于组合所述第一晶体效率因子和第二晶体效率因子，得到晶体对的晶体效率因子。

上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。

在不同的例子中，所述机器可读存储介质502可以是：RAM(Radom Access Memory，随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等)，或者类似的存储介质，或者它们的组合。

本申请实施例中，先确定晶体对中单个晶体的晶体效率因子和几何因子，然后，根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度，并利用晶体对的灵敏度对获取到患者的数据进行校正，最后利用校正后的数据进行图像重建，提高了图像重建的准确性。也就是说，本申请实施例中，根据单个晶体的晶体效率因子和几何因子确定晶体对的灵敏度，并利用该灵敏度对扫描的数据进行有效校正，提高了图像重建的准确性。另外，本实施例通过对晶体的本底辐射进行校正，无需购买放射源、制作放射源装置，简洁高效的计算出了晶体对的几何因子，不但降低了成本，还避免了操作人员的辐射危险。

另外，本申请实施例还提供一种医疗设备，包括：处理器；以及用于存储所述处理器可执行指令的存储器；该存储器可以是机器可读存储介质，其中，

所述处理器被配置为：

确定晶体对中单个晶体的晶体效率因子和几何因子，其中，所述晶体对为PET设备中晶体环上的晶体对；

根据所述晶体对中单个晶体的晶体效率因子和几何因子确定所述晶体对的灵敏度；

获取通过所述晶体对扫描患者的数据；

利用所述晶体对的灵敏度对所述数据进行校正；

根据校正后的数据进行图像重建。

其中，确定晶体对中单个晶体的晶体效率因子和几何因子，包括：

采集晶体对上的多个本底符合事件；

将所述多个本底符合事件划分成第一类本底符合事件和第二类本底符合事件；

确定第一类本底符合事件的第一晶体效率因子，以及第二类本底符合事件的第二晶体效率因子；

利用所述第一晶体效率因子对所述第一类本底符合事件进行校正，得到校正后第一类本底符合事件，以及利用所述第二晶体效率因子对所述第二类本底符合事件进行校正，得到校正后第二类本底符合事件；

利用每对晶体对上校正后的第一类本底符合事件与第二类本底符合事件的和，确定所述晶体对上每个晶体对应的第一几何因子和第二几何因子。

本申请实施例中，将采集到的多个本底符合事件划分为不同类型的本底符合事件，根据相应的本底符合事件分别计算相应的晶体效率因子，然后利用计算出的晶体效率因子校正之前对应的本底符合事件，最后，利用每对晶体对上校正后的本底符合事件计算单个晶体的几何因子，然后再组合单个晶体的几何因子，得到晶体对的集合因子。也就是说，本申请实施例中，根据设备晶体的特点，建立相应的晶体对灵敏度模型，无需购买放射源、制作放射源装置，简洁高效的计算出了晶体对的几何因子，不但降低了成本，还避免了操作人员的辐射危险。

对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。

以上所述仅为本申请的较佳实施例而已，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。