用于PET图像重建的线性衰减系数获取方法及重建方法与流程

用于pet图像重建的线性衰减系数获取方法及重建方法
技术领域
1.本发明涉及医学成像领域，尤其涉及一种正电子发射计算机断层成像系统中图像重建使用的线性衰减系数的获取方法、基于线性衰减系数的pet图像重建方法和pet探测系统。


背景技术：

2.正电子发射断层显像pet(positron emission tomography)是一种高端核医学影像诊断设备，在实际操作中利用放射性核素(如
18
f、
11
c等)对代谢物质进行标记并将核素注入人体，然后通过pet系统对患者进行功能代谢显像，反映生命代谢活动的情况，从而达到诊断的目的。在pet系统采集时，光子到达pet系统之前，会在人体中进行衰减，较之物体内部，物体表面的符合事件有着更大的探测效率。如果这种衰减因素没有被校正，重建图像中会产生物体边缘影像过亮，物体内部组织影像过暗的衰减伪影。
3.为了消除衰减伪影，准确定量分析放射性药物在患者体内的分布，通常会配套其他模态(比如ct、mri等)以得到患者的解剖结构成像。一方面可以对核素分布情况进行精确定位，提高了病灶定位的准确性；另一方面可以对应得到患者的组织密度分布，计算出组织对射线的衰减性质(线性衰减系数)，然后运用至图像重建中，对pet数据进行衰减校正，最终获得组织实际放射性分布的功能图像。两种图像最终将进行同机融合，兼容功能成像和解剖成像的优点，达到早期发现病灶和诊断疾病的目的，对肿瘤和心脏、脑部疾病的诊治指导更具优势。
4.然而在多模态采集应用中，有时无法准确获得与pet数据相匹配的衰减信息，以至于在pet图像上产生伪影：
5.首先，在pet多模态成像时，不同模态的图像位置可能存在相对偏差。以pet/ct系统为例，ct扫描通常可以在很短时间内完成，获得的图像几乎是某时刻的快照。但是pet扫描速度较慢，每个体位通常会花费几分钟时间，因此不可能在病人屏气状态下完成数据采集。受心脏搏动和呼吸运动的影响，pet图像和ct图像对相同病灶在位置及相位上存在一定程度的不匹配。一方面，pet采集将运动的病灶在不同位置时采集的数据叠加，对应的成像反映的是病灶位置的一个平均效应，这不可避免的带来了分辨率的下降，导致病灶成像模糊和定量分析suv值(标准化摄取值)准确性降低，与ct成像在形态上存在差异。另一方面，ct图像与pet图像配准融合产生偏差(特别在运动幅度最大的膈膜附近)，用瞬间的ct图像去对平均的pet图像做衰减校正必然出现位置偏差，会在pet图像上产生局部伪影，可能影响胸腹部肿瘤的准确诊断与治疗计划的制定。另外，在长时间的pet扫描中，病人身体有可能发生移动(如扫描时间较长时胳膊，头等部位会发生移动)，这同样会导致pet和ct图像不匹配，产生衰减伪影。
6.其次，pet临床采集上常用的
18
f-fdg是一种反映葡萄糖代谢的非特异性显像剂。相比正常细胞，恶性肿瘤细胞生长活跃，细胞异常增殖，对能量需求量大，大大增加了dna合成和氨基酸利用。因此通过探测葡萄糖聚集，pet能从代谢水平上评价肿瘤的活性。但是炎症
急性期和肉芽肿等良性病变也可因代谢增高而产生明显摄取，导致假阳性的发生。为了更好的鉴别良恶性病变，通常临床上采用延迟显像，通过多时间点采集，多时间点pet成像，肿瘤随间隔时间延长摄取fdg程度高于炎症，这样可以更好的鉴别良恶性病变。在实际应用中，将整个扫描时间分隔为几段，分别进行静态pet图像重建，最终获取多个时间点的静态pet图像并进行对比分析。通常延迟显像多时间点pet采集需要配套进行多时间点ct采集，以避免在扫描过程中由于病人中途离开病床或身体位置发生移动而引入衰减校正信息误差和衰减伪影，这不可避免的增加了病人的x射线辐射剂量。
7.再次，扫描过程中衰减图像存在明显的伪影，会导致pet图像衰减校正产生明显误差。例如体内含有金属物质患者的ct图像(比如心脏起搏器或金属牙套等)存在明显的高亮金属伪影，它使周围的组织难以准确分辨，从而使衰减图像产生明显伪影，严重影响pet图像的衰减校正。
8.再者，pet的扫描范围通常会大于其他模态(比如ct或mri)的扫描范围。在扫描体重比较大的患者时，其他模态成像很可能无法提供足够大的成像范围，这会导致衰减图像发生截断。这种不完整的衰减信息应用在pet重建中也会产生衰减伪影。
9.最后，当pet与其他模体联合成像时，有时并不能获得令人满意的衰减校正图像，比如pet/mr成像。与ct显像相比，mr主要利用磁自旋成像，而非组织密度分布成像，因此不能直接提供关于组织衰减性质的准确信息。目前依靠mr成像进行衰减校正的算法应用复杂，精度也不高，容易导致衰减伪影的产生。另外，mr无法对扫描床和mr线圈成像，这也会对后续的衰减校正产生影响。
10.同时，其他模态的应用条件也制约着pet成像的应用，比如带假牙或者心脏起搏器的患者无法做mr检查，也就影响了pet/mr的应用。另外ct成像需要极高的放射性防护要求，而mr成像需要进行严格的核磁共振屏蔽，这都导致多模态成像的扫描防护要求高，不易推广。


技术实现要素：

11.(一)要解决的技术问题
12.鉴于现有技术的上述缺点、不足，本发明提供一种用于pet图像重建的线性衰减系数获取方法及重建方法。
13.(二)技术方案
14.为了达到上述目的，本发明采用的主要技术方案包括：
15.第一方面，本发明实施例提供一种用于pet图像重建的线性衰减系数获取方法，其包括：
16.s10、基于采集的pet探测数据，构建估计无衰减校正的pet图像的对数似然函数l；
17.s20、根据预先定义的滤波信息和先验惩罚函数，对所述对数似然函数l进行调节并迭代优化，获得无衰减校正的pet图像；
18.s30、按照预先构建的映射关系将无衰减校正的pet图像映射为探测目标对应的线性衰减系数图像μ0；
19.s40、基于采集所述pet探测数据时对应的扫描床的位置参数变化量和所述探测目标对应的线性衰减系数图像μ0，获取用于pet图像重建的线性衰减系数μ1。
20.可选地，所述s10包括：
21.根据下述公式一重建，估计无衰减校正的pet图像的对数似然函数l；
22.公式一：
23.其中，r＝[r
1t
,r
2t
,
…
,r
nt
]
t
表示随机噪声的平均值，y＝[y
1t
,y
2t
,
…
,y
nt
]
t
表示采集的pet探测数据，m表示为pet图像空间的大小，a＝[a
ijt
]为系统矩阵，x＝[x1,x2,
…
,xj…
,xm]
t
表示未知pet图像，t为飞行时间。
[0024]
可选地，所述s20包括：
[0025]
根据下述公式二进行调节并迭代优化，获得无衰减校正的pet图像；
[0026]
公式二：
[0027]
其中，μ＝[μ1,μ2,
…
,μk]
t
表示未知的线性衰减系数分布，β为权重因子，f表示后滤波函数，r(x)为先验的标量惩罚函数。
[0028]
可选地，所述s30包括：
[0029]
预先构建的映射关系为预先训练的深度学习网络g，训练的深度学习网络g用于实现无衰减校正的pet图像与线性衰减系数图像μ0的映射；
[0030]
具体地，获取用于训练深度学习网络g的第一训练数据集，所述第一训练数据集包括：模拟仿真的训练数据和/或实际采集的训练数据；
[0031]
对所述第一训练数据集进行预处理并输入深度学习网络g，获取输出结果，并优化深度学习网络g的网络参数θ，使得损失函数l'最小，获得训练后的深度学习网络g，训练后的深度学习网络g中网络参数为的深度学习网络g中网络参数为
[0032]
基于训练的深度学习网络g，获取探测目标对应的线性衰减系数图像μ0，
[0033]
可选地，所述s30包括：
[0034]
预先构建的映射关系为预先训练的深度学习网络g’，训练的深度学习网络g’用于实现无衰减校正的pet图像与ct图像的映射；
[0035]
具体地，获取用于训练深度学习网络g’的第二训练数据集，所述第二训练数据集包括：模拟仿真的训练数据和/或实际采集的训练数据；
[0036]
对所述第二训练数据集进行预处理并输入深度学习网络g’，获取输出结果，并优化深度学习网络g’的网络参数θ，使得损失函数l'最小，获得训练后的深度学习网络g’，训练后的深度学习网络g’中网络参数为中网络参数为
[0037]
基于训练的深度学习网络g’，获取探测目标对应的ct图像；
[0038]
将所述ct图像转换为探测目标对应的线性衰减系数图像μ0。
[0039]
可选地，深度学习网络g、深度学习网络g’均为cnn网络、unet网络、gan网络或其他网络，本实施例不限定深度学习网络的具体结构，根据实际需要进行选择。
[0040]
可选地，所述s40包括：
[0041]
根据下述公式三对线性衰减系数图像μ0进行修正，获得获取用于pet图像重建的线性衰减系数图像μ1；
[0042]
公式三：μ1＝μ0+μ
bed
(δx,δy,δz)；
[0043]
μ
bed
(δx，δy，δz)为探测目标所在的扫描床的线性衰减系数图像，δx，δy，δz表示为三个方向的移动值。
[0044]
第二方面，本发明实施例还提供一种pet图像的重建方法，其包括：
[0045]
采用上述第一方面任一所述的线性衰减系数获取方法获得用于pet图像重建的线性衰减系数图像；
[0046]
将所述线性衰减系数图像值作为线性衰减系数的已知初始值；
[0047]
基于预先建立的pet图像的目标函数、已知初始值的线性衰减系数分布、pet放射性活度分布x，采用交替求解策略，获得最终重建的pet图像；
[0048]
其中，线性衰减系数分布μ和pet放射性活度分布x分别为pet图像的目标函数中的两个变量；
[0049]
所述交替求解策略为：在第一个变量为已知值时，获取第二个变量的估计值，将得到的第二个变量的估计值作为已知值，获取第一个变量的估计值，交替求解n次，n为大于1的自然数。
[0050]
可选地，基于预先建立的pet图像的目标函数、已知初始值的线性衰减系数分布、pet放射性活度分布x，采用交替求解策略，获得最终重建的pet图像，包括：
[0051]
线性衰减系数μ为已知值时，根据下述公式五，获取pet放射性活度分布x的估计值；
[0052][0053]
pet放射性活度分布x为已知值时，根据下述公式六，获取线性衰减系数分布μ的估计值；
[0054][0055]
其中，x＝[x1,x2,
…
,xm]
t
表示未知pet图像即pet放射性活度分布，μ＝[μ1,μ2,
…
,μk]
t
表示线性衰减系数分布，a＝[a
ijt
]为系统矩阵，r＝[r
1t
,r
2t
,
…
,r
nt
]
t
表示随机噪声和散射噪声的平均值，t表示飞行时间tof的维度，l＝[l
ik
]为线性衰减系数矩阵，j表示pet探测系统中空间位置点源，i表示响应线lor。
[0056]
第三方面，一种pet探测系统，包括：存储器和处理器；所述存储器中存储有计算机
程序指令，所述处理器执行所述存储器中存储的计算机程序指令，具体执行上述第二方面所述的pet图像的重建方法。
[0057]
(三)有益效果
[0058]
在本发明中，本技术中采用pet探测数据获取pet图像重建过程中的衰减校正信息，由此，在呼吸或心跳以及患者移动而导致pet多模态图像不匹配时，仍然能够对pet图像进行衰减校正，提高图像质量，为医生的分析和应用提供更准确的图像。
[0059]
另外，由于采用pet探测数据获取的线性衰减系数的初始值，进而迭代重建pet图像，针对衰减图像存在伪影的患者，如体内存在心脏起搏器或金属牙套的pet/ct扫描患者，ct图像存在明显金属伪影，可以进行准确的衰减校正，消除金属伪影的影响，在处理应用中有效避免了衰减图像截断的问题，便于医生对大体重患者进行扫描。
[0060]
可理解的是，pet图像重建过程中衰减校正迭代算法的初始值是通过深度学习网络从未衰减校正的pet图像中得到的，无论是定量性还是组织分布都比以前更加准确，因此大大提升了衰减校正算法的稳定性和迭代速度。
[0061]
在获取线性衰减系数图像时，通过图像阈中进行深度学习网络映射，处理速度快，额外增加的时间相对于重建过程可以忽略不计，保证了算法的可行性；深度学习网络生成的线性衰减系数图像经过衰减系数迭代算法依靠采集数据进行微调，解决了pet采集数据泛化的问题，也简化了深度学习网络的难度，提高了深度学习网络的稳定性；pet采集不依赖于其他模态，可以应用在单独pet扫描中，降低了扫描环境要求，拓展了应用场合。
附图说明
[0062]
图1和图2a均为本发明一实施例提供的用于pet图像重建的线性衰减系数获取方法的流程示意图；
[0063]
图2b为深度学习网络的训练过程示意图；
[0064]
图3为pet图像探测系统的框架示意图；
[0065]
图4为采用本发明的pet重建图像的方法和现有重建算法得到的pet图像的比对示意图。
具体实施方式
[0066]
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
[0067]
另外，现有技术中为了能够有效校正衰减伪影，拓宽pet成像应用场合，提供有两种衰减校正方法，其从pet采集到的飞行时间tof(time of flight)信息中迭代提取线性衰减系数分布图像(linear attenuation coefficient image)以用做pet重建的衰减校正，可以保证pet图像和衰减图像的严格匹配，有效地消除了运动伪影，提高图像质量，为医生的分析和应用提供更准确的图像。但是在实际应用时，由于很难获得准确的线性衰减系数初始分布，估计算法通常需要多次迭代：如果使用全成像视野均匀分布的线性衰减图像作为迭代初始，不准确的猜想会导致收敛过程变长，计算量增加；如果使用其他模态的图像作为先验迭代初始，受限于多模态图像间可能存在的不匹配，仍然需要多次迭代来消除伪影。
[0068]
因此，在实际应用中，为了解决衰减校正算法迭代收敛运算时间长的问题，通常需
要配套更高水平的计算资源(比如高性能gpu)，增加了成本。另外，如果迭代初始值与真实线性衰减系数存在较大差异，有可能会导致迭代算法收敛到局部最优，无法得到全局最优结果。特别地，现有技术中为了避免局部最优情况的发生，需要在迭代算法中加入很多限制和保护，也需要设置很多调节参数，这降低了算法的稳定性和鲁棒性。
[0069]
为了提高线性衰减系数迭代算法的运算速度和结果稳定性，本发明实施例中利用深度学习网络将没有进行衰减校正的pet图像映射为线性衰减系数图像作为后续算法的已知初始值，用来优化算法的收敛路径，达到尽快收敛到全局优化解的目的。相比于原先的初始值(全成像视野均匀线性衰减系数图像或其他模态图像转换得到的线性衰减系数图像)，此算法得到的线性衰减系数图像带有更准确的扫描物体衰减信息，作为真实衰减图像的良好近似，可以保证迭代过程快速收敛，增加算法的稳定性和定量性准确性。另外，深度学习模型提取的衰减信息来源于pet图像，不存在多模态图像之间的不匹配，避免了运动和截断伪影。
[0070]
需要说明的是，线性衰减系数图像和线性衰减系数在下述表示的是同一个含义，在不同的描述中使用了，其含义相同。
[0071]
实施例一
[0072]
如图1所示，本实施例提供一种用于pet图像重建的线性衰减系数获取方法，本实施例的方法可在任一电子设备上实现，优先在pet探测器关联的计算设备中实现，本实施例的方法可包括下述的步骤：
[0073]
s10、基于采集的pet探测数据，构建估计无衰减校正的pet图像的对数似然函数l。
[0074]
举例来说，无衰减校正的pet图像的对数似然函数l如下；
[0075]
公式一：
[0076]
其中，r＝[r
1t
,r
2t
,
…
,r
nt
]
t
表示随机噪声的平均值，y＝[y
1t
,y
2t
,
…
,y
nt
]
t
表示采集的pet探测数据，m表示为pet图像空间的大小，且空间大小不变，a＝[a
ijt
]为系统矩阵，x＝[x1,x2,
…
,xj…
,xm]
t
表示未知pet图像，t为飞行时间。
[0077]
s20、根据预先定义的滤波信息和先验惩罚函数，对所述对数似然函数l进行调节并迭代优化，获得无衰减校正的pet图像。
[0078]
在本实施例中，可根据下述公式二进行迭代优化，获得无衰减校正的pet图像；
[0079]
公式二：
[0080]
其中，μ＝[μ1,μ2,
…
,μk]
t
表示未知的线性衰减系数分布，β为权重因子，f表示后滤波函数，r(x)为先验的标量惩罚函数。
[0081]
s30、按照预先构建的映射关系将无衰减校正的pet图像映射为探测目标对应的线性衰减系数图像μ0。
[0082]
可理解的是，预先构建的映射关系为预先训练的深度学习网络g，训练的深度学习网络g用于实现无衰减校正的pet图像与线性衰减系数图像μ0的映射；
[0083]
具体地，获取用于训练深度学习网络g的第一训练数据集，所述第一训练数据集包括：模拟仿真的训练数据和/或实际采集的训练数据；
[0084]
对所述第一训练数据集进行预处理并输入深度学习网络g，获取输出结果，并优化
深度学习网络g的网络参数θ，使得损失函数l'最小，获得训练后的深度学习网络g，训练后的深度学习网络g中网络参数为的深度学习网络g中网络参数为如图2b所示。
[0085]
基于训练的深度学习网络g，获取探测目标对应的线性衰减系数图像μ0，
[0086]
在其他实施例中，预先构建的映射关系为预先训练的深度学习网络g’，训练的深度学习网络g’用于实现无衰减校正的pet图像与ct图像的映射。相应地，在训练之后，基于训练的深度学习网络g’，获取探测目标对应的ct图像；将所述ct图像转换为探测目标对应的线性衰减系数图像μ0。
[0087]
本实施例的深度学习网络g、深度学习网络g’均可为cnn网络、unet网络、gan网络或其他网络，本实施例中不限定深度学习网络g的具体结构，根据需要设置。即深度学习网络g实现将pet图像映射为线性衰减系数图像，深度学习网络g’实现将pet图像映射为ct图像，此后根据ct图像和线性衰减系数图像之间的一一对应关系，将ct图像转为线性衰减系数图像。
[0088]
s40、基于采集所述pet探测数据时对应的扫描床的位置参数变化量和所述探测目标对应的线性衰减系数图像μ0，获取用于pet图像重建的线性衰减系数μ1。
[0089]
举例来说，该步骤s40可根据下述公式三对线性衰减系数图像μ0进行修正，获得获取用于pet图像重建的线性衰减系数μ1；
[0090]
公式三：μ1＝μ0+μ
bed
(δx,δy,δz)；
[0091]
μ
bed
(δx，δy，δz)为探测目标所在的扫描床的线性衰减系数图像，δx，δy，δz表示为三个方向的移动值。
[0092]
在本实施例中，在获取线性衰减系数图像时，通过图像阈中进行深度学习网络映射，处理速度快，额外增加的时间相对于重建过程可以忽略不计，保证了算法的可行性；深度学习网络生成的线性衰减系数图像经过衰减系数迭代算法依靠采集数据进行微调，解决了pet采集数据泛化的问题，也简化了深度学习网络的难度，增加了深度学习网络的稳定性；pet探测数据的采集不依赖于其他模态，可以应用在单独pet扫描中，降低了扫描环境要求，拓展了应用场合。
[0093]
特别地，采用pet探测数据获取pet图像重建过程中的衰减校正信息，由此，在呼吸或心跳以及患者移动而导致pet多模态图像不匹配时，仍然能够对pet图像进行衰减校正，提高图像质量，为医生的分析和应用提供更准确的图像。
[0094]
实施例二
[0095]
本实施例提供一种利用深度学习网络，从未经衰减校正的pet图像中映射得到线性衰减系数图像，再将得到的线性衰减系数图像值作为线性衰减系数迭代算法的初始值，获取重建的pet图像，保证线性衰减系数估计算法快速稳定收敛。该方法可以在pet探测系统的计算设备中完成，其有效提高计算速度和计算时间，增加可靠性和稳定性，结合图2a至图4所示，具体步骤如下：
[0096]
以下的步骤201至步骤203均为现有的建模过程，由于后续步骤需要公式及说明，故在此列出。
expectation maximization，osem)。
[0112]
前述未做衰减校正的pet图像存在明显的衰减伪影，导致放射性活度分布定量性不准确：例如患者边缘影像过亮，患者内部组织影像过暗，肺部摄取过高等等。
[0113]
虽然pet图像中患者不同组织的放射性摄取值不准确，但是依然保留了不同组织的结构信息，比如虽然患者边缘高亮，其边缘范围依然可以确定；肺部摄取虽然对比度不正确，但是依然可以依靠图像勾划肺部轮廓。
[0114]
为此，在本实施例中可以通过未衰减校正的pet图像还原线性衰减系数图像，能够准确反映患者内部的组织分布。经过试验验证，决定于人体组织的密度，人体组织的线性衰减系数值在线性衰减系数图像中不会存在特别大的偏差，比如软组织区域大致为0.0975cm-1、脂肪区域大致为0.0864cm-1、肺部区域大致为0.0224cm-1。所以从未衰减校正的pet图像还原得到的线性衰减图像定量性准确性也比较容易保证。
[0115]
205、为了能够充分提取pet图像的特征，本实施例中选择利用深度学习网络g实现未衰减校正的pet图像与线性衰减系数图像的映射。
[0116]
深度学习网络g为预先训练的网络。在训练过程中，以pet/ct为例，将未衰减pet图像作为输入，线性衰减系数图像μ0作为输出，与ct扫描得到的真实线性衰减图像μ
ct
进行比较，通过优化训练网络参数θ使得损失函数l’最小，最终可以实现将未衰减校正的pet图像映射为准确的线性衰减系数图像，即：
[0117][0118]
其中，是训练处中优化得到的，训练数据集可以来自模拟仿真，也可以来自实际采集。训练数据集需要进行预处理，通过筛选以保证线性衰减系数图像和未衰减pet图像完全匹配，不存在截断或运动伪影。该处的预处理即为通过筛选以保证线性衰减系数图像和未衰减pet原始图像完全匹配，不存在截断或运动伪影。
[0119]
不失一般性，深度学习网络g可以选择cnn网络、unet网络、gan网络或其他网络，本实施例不对深度学习网络g的结构进行限定。
[0120]
206、由于未做衰减校正的pet图像仅属于进行功能成像，pet图像上无法显示扫描床，因此需要在深度学习网络g输出的线性衰减图像上再额外加上床的衰减信息。加上扫描床的衰减信息的线性衰减系数图像为μ1＝μ0+μ
bed
(δx,δy,δz)。
[0121]
在本实施例中，由于扫描床的形状已知且只能进行刚性运动，并且在扫描前pet探测系统进行过良好的机械调整，因此扫描床只需考虑三个方向的平移，扫描床的线性衰减系数分布可以表示为μ
bed
(δx，δy，δz)，δx，δy，δz表示为三个方向的移动值，μ
bed
(0，0，0)为床的初始位置。由于扫描时床不能水平移动，因此水平方向移动δx可以通过机械测量或校正过程来得到，在扫描过程中设定为常量；竖直方向移动δy可以通过读取床的升降高度值来确定；轴向方向移动δz可以通过读取扫描定位灯的轴向位置来决定。在其他实施例中，可以不依靠外部信号，也可以利用神经网络通过学习pet图像中患者和床的相对位置来得到。
[0122]
由于pet图像是功能性成像，无法对扫描床进行成像，但是因为患者是平躺在扫描床上的，因此身体下侧的活度分布是平的，而扫描床上的头托决定了头部的起始扫描位置，
因此pet图像本身虽然没有显示扫描床，但是其放射性活度分布含有扫描床的位置信息，可以通过图像形状识别来学习得到，具体为网络输入是衰减或无衰减的pet图像，输出为扫描床的位置，通过学习构建网络。
[0123]
需要说明的是，步骤205中借助于深度学习网络g将未衰减校正的pet图像映射为线性衰减系数图像μ1，而步骤206中扫描床的补偿也是线性衰减系数图。
[0124]
在另一可能的实现方式中，由于线性衰减系数与ct值符合双线性转换关系，满足一一对应的要求，且ct值的动态范围更大，因此步骤205和步骤206也可以采用深度学习网络g’将未衰减校正的pet图像映射为ct图像，然后对映射的ct图像加扫描床的补偿信息，再进一步将增加有扫描床的补偿信息的ct图像转换为线性衰减系数图像μ1。
[0125]
当然，若深度学习网络g’将未衰减校正的pet图像映射为ct图像，则在训练过程中，也是针对ct图像进行训练和补偿，最后床位补偿后的ct图像转换成线性衰减系数图像。ct图像转换成线性衰减系数中无需使用深度学习网络g’。
[0126]
在具体操作中，由于ct值的动态范围大，因此在优化的深度学习网络进行梯度计算时，会有更高的精度，使得优化结果会更加准确。
[0127]
另外，为更好的理解，以下对深度学习网络g的选择和训练说明如下。受不同场地、不同设备、不同扫描参数的影响，pet图像质量存在很大差异，因此在实际应用中很难保证当前pet图像(无衰减校正的pet图像)与训练pet图像质量相同，这会极大的影响学习网络结果的适用性，即由于pet探测数据泛化导致的问题。为了解决pet探测数据泛化的问题，通常需要为深度学习网络训练提供各种不同情况的图像数据，这种数据要求通常不现实，对于深度学习网络的构建也带来很大难度，训练时间和内存要求也很大。因此本实施例不直接在重建中应用μ1，而是将μ1作为线性衰减系数分布迭代算法的初始变量，这样一方面可以利用采集数据对训练得到的μ1进行微调，保证线性衰减系数结果与实际采集数据相符，变相降低了训练网络数据泛化的要求；另一方面，μ1作为真实线性衰减系数分布的一个良好估计，相比于原来算法设定的初始值，可以大大提升原算法的运算速度和定量性准确性，可以很快得到和实际采集更加匹配的线性衰减图像。即为了解决不匹配的数据泛化问题，采用实际采集数据对结果进行微调。
[0128]
207、由于公式(3)中的log-likelihood函数对于未知数x和μ是一个很复杂的函数，很难得到解析解，因此需要利用迭代算法逐渐逼近最优解。对未知的pet放射性活度分布x，保持线性衰减系数分布μ为常数，最大化log-likelihood函数，即pet图像重建通用的mlem算法：
[0129][0130]
公式(7)中，n表示当前迭代次数。
[0131]
208、保持pet放射性活度分布x为常数，针对未知衰减系数分布μ最大化log-likelihood函数，直接利用pet探测数据计算得到新的线性衰减系数分布μ，对应算法为：
[0132][0133]
在具体实施过程中，先保持线性衰减系数分布μ为常数，针对pet放射性活度分布x来最大化目标函数，采用上述步骤207，即传统mlem迭代重建算法；再选择保持pet放射性活度分布x为常数，针对未知线性衰减系数分布μ来最大化目标函数，采用步骤208。第一次迭代时pet图像重建所应用的线性衰减系数分布的初始值即为步骤206得到的线性衰减系数图μ1的值。这样交替进行操作，不断修正衰减校正使之逼近真实衰减情况，最终得到满足最大化目标函数要求的x和μ的估计值。
[0134]
以pet/ct系统为例，图3定义了多模态探测系统坐标系，图4(a)为未做衰减校正的pet图像，图4(b)为上述步骤208和步骤207迭代后得到的衰减校正的pet图像，图4(c)为利用深度学习网络得到的线性衰减系数分布，用作线性衰减系数迭代算法的初始值。
[0135]
相比传统的通过其他的模态图像进行衰减校正的方法，本实施例重建的pet图像质量更好，解决了pet采集数据泛化的问题，降低了扫描环境的要求，拓展了应用场合。
[0136]
实施例三
[0137]
本发明实施例提供一种pet图像的重建方法，本实施例的方法可在任一电子设备上实现，其包括：
[0138]
301、采用上述实施例一所述的线性衰减系数获取方法获得用于pet图像重建的线性衰减系数图像μ1；
[0139]
302、将所述线性衰减系数图像μ1值作为线性衰减系数分布μ的已知初始值；
[0140]
303、基于预先建立的pet图像的目标函数、已知初始值的线性衰减系数、pet放射性活度分布x，采用交替求解策略，获得最终重建的pet图像；
[0141]
其中，线性衰减系数分布μ和pet放射性活度分布x分别为pet图像的目标函数中的两个变量；
[0142]
所述交替求解策略为：在第一个变量为已知值时，获取第二个变量的估计值，将得到的第二个变量的估计值作为已知值，获取第一个变量的估计值，交替求解n次，n为大于1的自然数。
[0143]
举例来说，线性衰减系数μ为已知值时，根据下述公式五，获取pet放射性活度分布x的估计值；
[0144][0145]
pet放射性活度分布x为已知值时，根据下述公式六，获取线性衰减系数分布μ的估
计值；
[0146][0147]
其中，x＝[x1,x2,
…
,xm]
t
表示未知pet图像即pet放射性活度分布，μ＝[μ1,μ2,
…
,μk]
t
表示线性衰减系数分布，a＝[a
ijt
]为系统矩阵，r＝[r
1t
,r
2t
,
…
,r
nt
]
t
表示随机噪声和散射噪声的平均值，t表示飞行时间tof的维度，l＝[l
ik
]为线性衰减系数矩阵，j表示pet探测系统中空间位置点源，i表示响应线lor。
[0148]
另外，本发明实施例还提供一种pet探测系统，其包括：存储器和处理器；所述存储器中存储有计算机程序指令，所述处理器执行所述存储器中存储的计算机程序指令，具体执行上述任意实施例所述的pet图像的重建方法。
[0149]
应当注意的是，在权利要求中，不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中，这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用，仅是为了表述方便，而不表示任何顺序。可将这些词语理解为部件名称的一部分。
[0150]
此外，需要说明的是，在本说明书的描述中，术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述，是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0151]
尽管已描述了本发明的优选实施例，但本领域的技术人员在得知了基本创造性概念后，则可对这些实施例作出另外的变更和修改。所以，权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
[0152]
显然，本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也应该包含这些修改和变型在内。