低剂量PET图像重建方法、装置、设备及存储介质与流程

本发明属于pet成像技术领域，尤其涉及一种低剂量pet图像重建方法、装置、设备及存储介质。

背景技术：

正电子发射断层成像(positronemissiontomography，简称pet)是一种发射型成像技术(emissiontomography，简称et)，它通过把放射性药物注入体内的方法来显示不同组织的新陈代谢情况。pet技术是继计算机断层成像(computedtomography，简称ct)和磁共振成像(magneticresonanceimaging，简称mri)之后应用于临床的一种新型影像技术。

在pet成像中，放射性药物实际上是个分子载体，它依附于特定的生理组织或病理过程。放射性物质在药物的带领下在人体内有目的的分布。pet成像的目的实际上就是得到放射性物质在人体内部的分布图，它的工作原理是：将一些放射性核元素，如o-15、c-11、n-13和f-18等标记在人体代谢所需的化合物上，然后通过手臂静脉血管注射等方式输入受检者体内。标记化合物在参与体内代谢的过程中，放射性核元素发生衰变，释放出正电子(带一个正电荷的电子)，正电子与其周围的(带负电)电子发生湮灭，产生两个能量为511kev的伽马光子。这对光子在一条直线上朝相反的方向射出，利用体外的伽马照相机可以探测到特定区域放射的所有光子，然后设计一定的算法，就可以近似得到放射性物质在人体内部的分布情况。

由于在pet检查中使用的放射性药物会对近距离接触该药物的人员产生辐射，而受到辐射的人员患癌的几率会远高于正常人，同时放射性药物的消耗在pet检查的成本中占有一定比重。因此，根据国际放射防护委员会(internationalcommissiononradiologicalprotection，简称icrp)提出的合理使用低剂量(aslowasreasonablyachievable，简称alara)原则，在pet临床诊断时，以期用最小的剂量获得满足临床需求的图像，尽量降低对患者的辐射剂量。

然而，在对低剂量采样得到的测量数据进行pet图像重建时，现有传统的pet图像重建算法重建图像的速度慢，进而使得重建图像产生运动伪影，这些伪影将会直接影响医生的诊断行为。

技术实现要素：

本发明的目的在于提供一种低剂量pet图像重建方法、装置、设备及存储介质，旨在解决由于现有技术无法提供一种有效的低剂量pet图像重建方法，导致低剂量pet图像重建速度慢、且重建图像质量差的问题。

一方面，本发明提供了一种低剂量pet图像重建方法，所述方法包括下述步骤：

当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备采集到的投影数据，并获取所述pet设备的系统矩阵；

根据所述投影数据以及所述系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建，获得初始重建pet图像；

对所述初始重建pet图像进行梯度运算，获得所述初始重建pet图像对应的水平梯度图像和垂直梯度图像；

对所述水平梯度图像和所述垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示；

通过最小二乘法将所述水平梯度图像稀疏表示和所述垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像；

判断所述重建pet图像是否满足预设的目标函数，是则，输出所述重建pet图像，否则，将所述重建pet图像设置为所述待重建pet图像，并跳转到通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建的步骤。

另一方面，本发明提供了一种低剂量pet图像重建装置，所述装置包括：

参数获取单元，用于当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备采集到的投影数据，并获取所述pet设备的系统矩阵；

初始重建单元，用于根据所述投影数据以及所述系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建，获得初始重建pet图像；

梯度图像获取单元，用于对所述初始重建pet图像进行梯度运算，获得所述初始重建pet图像对应的水平梯度图像和垂直梯度图像；

字典学习单元，用于对所述水平梯度图像和所述垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示；

重建图像获取单元，用于通过最小二乘法将所述水平梯度图像稀疏表示和所述垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像；以及

重建图像判断单元，用于判断所述重建pet图像是否满足预设的目标函数，是则，输出所述重建pet图像，否则，将所述重建pet图像设置为所述待重建pet图像，并触发所述初始重建单元执行通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建的步骤。

另一方面，本发明还提供了一种计算设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述低剂量pet图像重建方法所述的步骤。

另一方面，本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上述低剂量pet图像重建方法所述的步骤。

本发明根据采集到的投影数据和pet设备的系统矩阵，通过pet图像重建算法对待重建pet图像进行重建，获得初始重建pet图像，通过梯度运算和字典学习，分别获得初始重建pet图像对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像，判断重建pet图像是否满足目标函数，是则输出重建pet图像，否则，继续对待重建pet图像执行图像重建、梯度运算以及字典学习的操作，从而提高了低剂量pet图像的重建速度，且降低重建图像的伪影程度，进而提高了低剂量pet图像重建的图像质量。

附图说明

图1是本发明实施例一提供的低剂量pet图像重建方法的实现流程图；

图2是本发明实施例二提供的低剂量pet图像重建装置的结构示意图；

图3是本发明实施例二提供的低剂量pet图像重建装置的优选结构示意图；以及

图4是本发明实施例三提供的计算设备的结构示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

以下结合具体实施例对本发明的具体实现进行详细描述：

实施例一：

图1示出了本发明实施例一提供的低剂量pet图像重建方法的实现流程，为了便于说明，仅示出了与本发明实施例相关的部分，详述如下：

在步骤s101中，当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备采集到的投影数据，并获取pet设备的系统矩阵。

本发明实施例适用于医学图像处理平台、系统或设备，例如个人计算机、服务器等。当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备在低剂量条件下采集到的投影数据，并获取pet设备的系统矩阵，该系统矩阵是根据pet设备的几何结构信息计算得到的。

在步骤s102中，根据投影数据以及系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建，获得初始重建pet图像。

在本发明实施例中，根据投影数据以及系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行预设次数的迭代操作，以对待重建pet图像进行图像重建，获得初始重建pet图像，其中，待重建pet图像是二维图像，预设的pet图像重建算法为最大似然期望最大算法(maximumlikelihoodexpectationmaximized，简称mlem)或者有序子集期望值最大算法(orderedsubsetexpectationmaximization，简称osem)或者最大后验概率算法(maximumaposterior，map)。

在初始化待重建pet图像时，作为示例地，将待重建pet图像的像素值都初始化为零。

在步骤s103中，对初始重建pet图像进行梯度运算，获得初始重建pet图像对应的水平梯度图像和垂直梯度图像。

在本发明实施例中，在对初始重建pet图像进行梯度运算时，优选地，根据预设的水平梯度图像公式对初始重建pet图像进行梯度运算，获得该初始重建pet图像对应的水平梯度图像，根据预设的垂直梯度图像公式对初始重建pet图像进行梯度运算，获得该初始重建pet图像对应的垂直梯度图像，其中，水平梯度图像公式为垂直梯度图像公式为为初始重建pet图像，j和k分别为中像素的坐标，从而通过梯度运算将初始重建pet图像从图像域转换到梯度域，提高了初始重建pet图像的稀疏性，进而提高后续pet图像重建的成像效果。

在对初始重建pet图像进行梯度运算之前，优选地，对重建得到的初始重建pet图像进行正约束，从而提高后续pet图像重建的成像效果。

在步骤s104中，对水平梯度图像和垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示。

在本发明实施例中，对水平梯度图像和垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，由于医学图像的近邻像素值具有相关性，即转换至梯度域中的图像比图像域中的图像更稀疏，因此，对梯度域中的图像进行字典学习比对图像域中的图像进行字典学习，可以获得更稀疏的图像稀疏表示，从而提高后续pet图像重建的成像效果。

在对水平梯度图像和垂直梯度图像分别进行字典学习时，优选地，根据水平梯度图像和垂直梯度图像，通过预设的字典训练算法(例如，k-svd算法)进行字典学习，再根据学习结束的字典，对水平梯度图像稀疏表示和垂直梯度图像稀疏表示进行更新，从而提高后续pet图像重建的成像效果。

在对水平梯度图像和垂直梯度图像分别进行字典学习时，进一步优选地，根据预设的图像块数量对水平梯度图像和垂直梯度图像进行图像块提取，得到对应的水平梯度图像块和垂直梯度图像块，对水平梯度图像块和垂直梯度图像块进行字典学习，获得水平梯度图像块对应的稀疏表示和垂直梯度图像块对应的稀疏表示，由水平梯度图像块对应的稀疏表示构成水平梯度图像稀疏表示，由垂直梯度图像块对应的稀疏表示构成垂直梯度图像稀疏表示，从而通过自适应字典将全局图像的稀疏表示转换成基于图像块的稀疏表述，有效地捕捉初始重建pet图像的局部特征，并去除该图像的伪影和块效应，提高后续pet图像重建的成像效果。

在步骤s105中，通过最小二乘法将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像。

在本发明实施例中，通过最小二乘法将水平梯度图像稀疏表示和垂直梯度图像稀疏表示进行融合，以将初始重建pet图像从梯度域恢复到图像域，获得重建pet图像。

在步骤s106中，判断重建pet图像是否满足预设的目标函数。

在步骤s107中，输出重建pet图像。

在步骤s108中，将重建pet图像设置为待重建pet图像。

在本发明实施例中，当重建pet图像满足预设的目标函数时，执行步骤s107，结束pet图像的重建过程，当重建pet图像不满足预设的目标函数时，执行步骤s108，并跳转到步骤s102执行通过预设的pet图像重建算法对待重建pet图像进行图像重建。

在本发明实施例中，优选地，目标函数表示为：

其中，g为投影数据，m为系统矩阵，f为待重建pet图像，v1为预设大于0的惩罚参数，和表示水平方向和垂直方向的差分算子，表示f的梯度图像，和分别表示f的水平梯度图像和垂直梯度图像，为从中提取出的第l个图像块的向量表示，d为预设的过完备字典，αl为第l个图像块所对应的稀疏系数，t0为每个图像块的稀疏度，通过自适应稀疏变换算法结合迭代重建算法，来实现重建图像对真实图像的不断逼近，从而降低在pet检查中使用的放射性药物的剂量，并提高了pet图像重建的成像效果。

在本发明实施例中，根据采集到的投影数据和pet设备的系统矩阵，通过pet图像重建算法对待重建pet图像进行重建，获得初始重建pet图像，通过梯度运算和字典学习，分别获得初始重建pet图像对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像，判断重建pet图像是否满足目标函数，是则输出重建pet图像，否则，继续对待重建pet图像执行图像重建、梯度运算以及字典学习的操作，从而提高了低剂量pet图像的重建速度，且降低重建图像的伪影程度，进而提高了低剂量pet图像重建的图像质量。

实施例二：

图2示出了本发明实施例二提供的低剂量pet图像重建装置的结构，为了便于说明，仅示出了与本发明实施例相关的部分，其中包括：

参数获取单元21，用于当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备采集到的投影数据，并获取pet设备的系统矩阵。

本发明实施例适用于医学图像处理平台、系统或设备，例如个人计算机、服务器等。当接收到对低剂量pet图像进行重建的请求时，获取通过pet设备在低剂量条件下采集到的投影数据，并获取pet设备的系统矩阵，该系统矩阵是根据pet设备的几何结构信息计算得到的。

初始重建单元22，用于根据投影数据以及系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行图像重建，获得初始重建pet图像。

在本发明实施例中，根据投影数据以及系统矩阵，通过预设的pet图像重建算法对预先初始化的待重建pet图像进行预设次数的迭代操作，以对待重建pet图像进行图像重建，获得初始重建pet图像，其中，待重建pet图像是二维图像，预设的pet图像重建算法为最大似然期望最大算法(maximumlikelihoodexpectationmaximized，简称mlem)或者有序子集期望值最大算法(orderedsubsetexpectationmaximization，简称osem)或者最大后验概率算法(maximumaposterior，map)。

在初始化待重建pet图像时，作为示例地，将待重建pet图像的像素值都初始化为零。

梯度图像获取单元23，用于对初始重建pet图像进行梯度运算，获得初始重建pet图像对应的水平梯度图像和垂直梯度图像。

在本发明实施例中，在对初始重建pet图像进行梯度运算时，优选地，根据预设的水平梯度图像公式对初始重建pet图像进行梯度运算，获得该初始重建pet图像对应的水平梯度图像，根据预设的垂直梯度图像公式对初始重建pet图像进行梯度运算，获得该初始重建pet图像对应的垂直梯度图像，其中，水平梯度图像公式为垂直梯度图像公式为为初始重建pet图像，j和k分别为中像素的坐标，从而通过梯度运算将初始重建pet图像从图像域转换到梯度域，提高了初始重建pet图像的稀疏性，进而提高后续pet图像重建的成像效果。

在对初始重建pet图像进行梯度运算之前，优选地，对重建得到的初始重建pet图像进行正约束，从而提高后续pet图像重建的成像效果。

字典学习单元24，用于对水平梯度图像和垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示。

在本发明实施例中，对水平梯度图像和垂直梯度图像分别进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，由于医学图像的近邻像素值具有相关性，即转换至梯度域中的图像比图像域中的图像更稀疏，因此，对梯度域中的图像进行字典学习比对图像域中的图像进行字典学习，可以获得更稀疏的图像稀疏表示，从而提高后续pet图像重建的成像效果。

在对水平梯度图像和垂直梯度图像分别进行字典学习时，优选地，根据水平梯度图像和垂直梯度图像，通过预设的字典训练算法(例如，k-svd算法)进行字典学习，再根据学习结束的字典，对水平梯度图像稀疏表示和垂直梯度图像稀疏表示进行更新，从而提高后续pet图像重建的成像效果。

在对水平梯度图像和垂直梯度图像分别进行字典学习时，进一步优选地，根据预设的图像块数量对水平梯度图像和垂直梯度图像进行图像块提取，得到对应的水平梯度图像块和垂直梯度图像块，对水平梯度图像块和垂直梯度图像块进行字典学习，获得水平梯度图像块对应的稀疏表示和垂直梯度图像块对应的稀疏表示，由水平梯度图像块对应的稀疏表示构成水平梯度图像稀疏表示，由垂直梯度图像块对应的稀疏表示构成垂直梯度图像稀疏表示，从而通过自适应字典将全局图像的稀疏表示转换成基于图像块的稀疏表述，有效地捕捉初始重建pet图像的局部特征，并去除该图像的伪影和块效应，提高后续pet图像重建的成像效果。

重建图像获取单元25，用于通过最小二乘法将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像。

在本发明实施例中，通过最小二乘法将水平梯度图像稀疏表示和垂直梯度图像稀疏表示进行融合，以将初始重建pet图像从梯度域恢复到图像域，获得重建pet图像。

重建图像判断单元26，用于判断重建pet图像是否满足预设的目标函数，是则，输出该重建pet图像，否则，将该重建pet图像设置为待重建pet图像，并触发初始重建单元22执行通过预设的pet图像重建算法对待重建pet图像进行图像重建步骤。

在本发明实施例中，优选地，目标函数表示为：

其中，g为投影数据，m为系统矩阵，f为待重建pet图像，v1为预设大于0的惩罚参数，和表示水平方向和垂直方向的差分算子，表示f的梯度图像，和分别表示f的水平梯度图像和垂直梯度图像，为从中提取出的第l个图像块的向量表示，d为预设的过完备字典，αl为第l个图像块所对应的稀疏系数，t0为每个图像块的稀疏度，通过自适应稀疏变换算法结合迭代重建算法，来实现重建图像对真实图像的不断逼近，从而降低在pet检查中使用的放射性药物的剂量，并提高了pet图像重建的成像效果。

因此，如图3所示，优选地，梯度图像获取单元23包括：

梯度图获取子单元231，用于根据预设的水平梯度图像公式和预设垂直梯度图像公式对初始重建pet图像进行梯度运算，获得初始重建pet图像对应的水平梯度图像和垂直梯度图像，其中，水平梯度图像公式为垂直梯度图像公式为为初始重建pet图像，j和k分别为中像素的坐标。

字典学习单元24包括：

图像块提取单元241，用于根据预设的图像块数量对水平梯度图像和垂直梯度图像进行图像块提取，得到对应的水平梯度图像块和垂直梯度图像块；以及

字典学习子单元242，用于对水平梯度图像块和垂直梯度图像块进行字典学习，获得对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示。

在本发明实施例中，低剂量pet图像重建装置的各单元可由相应的硬件或软件单元实现，各单元可以为独立的软、硬件单元，也可以集成为一个软、硬件单元，在此不用以限制本发明。

实施例三：

图4示出了本发明实施例三提供的计算设备的结构，为了便于说明，仅示出了与本发明实施例相关的部分。

本发明实施例的计算设备4包括处理器40、存储器41以及存储在存储器41中并可在处理器40上运行的计算机程序42。该处理器40执行计算机程序42时实现上述低剂量pet图像重建方法实施例中的步骤，例如图1所示的步骤s101至s108。或者，处理器40执行计算机程序42时实现上述各装置实施例中各单元的功能，例如图2所示单元21至26的功能。

在本发明实施例中，根据采集到的投影数据和pet设备的系统矩阵，通过pet图像重建算法对待重建pet图像进行重建，获得初始重建pet图像，通过梯度运算和字典学习，分别获得初始重建pet图像对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像，判断重建pet图像是否满足目标函数，是则输出重建pet图像，否则，继续对待重建pet图像执行图像重建、梯度运算以及字典学习的操作，从而提高了低剂量pet图像的重建速度，且降低重建图像的伪影程度，进而提高了低剂量pet图像重建的图像质量。

本发明实施例的计算设备可以为个人计算机、服务器。该计算设备4中处理器40执行计算机程序42时实现低剂量pet图像重建方法时实现的步骤可参考前述方法实施例的描述，在此不再赘述。

实施例四：

在本发明实施例中，提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述低剂量pet图像重建方法实施例中的步骤，例如，图1所示的步骤s101至s108。或者，该计算机程序被处理器执行时实现上述各装置实施例中各单元的功能，例如图2所示单元21至26的功能。

在本发明实施例中，根据采集到的投影数据和pet设备的系统矩阵，通过pet图像重建算法对待重建pet图像进行重建，获得初始重建pet图像，通过梯度运算和字典学习，分别获得初始重建pet图像对应的水平梯度图像稀疏表示和垂直梯度图像稀疏表示，将水平梯度图像稀疏表示和垂直梯度图像稀疏表示从梯度域恢复到图像域，获得重建pet图像，判断重建pet图像是否满足目标函数，是则输出重建pet图像，否则，继续对待重建pet图像执行图像重建、梯度运算以及字典学习的操作，从而提高了低剂量pet图像的重建速度，且降低重建图像的伪影程度，进而提高了低剂量pet图像重建的图像质量。

本发明实施例的计算机可读存储介质可以包括能够携带计算机程序代码的任何实体或装置、记录介质，例如，rom/ram、磁盘、光盘、闪存等存储器。

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