基于深度学习的低剂量PET数据三维迭代更新重建方法与流程

基于深度学习的低剂量pet数据三维迭代更新重建方法
技术领域
[0001]
本发明属于医学影像技术领域，具体地涉及一种基于深度学习的低剂量pet数据三维迭代更新重建方法。


背景技术：

[0002]
正电子发射断层扫描（positron emission tomography，pet）是一种能提供活体生化和定量信息的医学影像，由于可以在分子水平成像，并且能同时体现人体生物学功能代谢信息和形态解剖结构信息，使其在肿瘤学、心脏学、神经学以及精神疾病等方面具有重要应用。pet的成像过程包括扫描前对病人注入放射性示踪剂，示踪剂在人体内参与生理代谢时发生衰变产生正电子，正电子和邻近组织中的电子发生湮灭效应，产生逆向运动的511kev 光子对。用两个位置相对的探测器对光子对进行符合探测可记录一条符合响应线。采集一定数量的符合响应线保存成三维pet原始数据sinogram，对sinogram进行校正后再用基于sinogram统计信息的迭代重建方法对其进行三维迭代重建，即可获得表征人体内各组织代谢强度的三维pet图像。
[0003]
随着大型医疗设备pet/ct的应用推广，越来越多的人开始接触和接受用pet/ct进行肿瘤筛查，甚至健康体检。然而，由于大众接受到更多的辐射安全知识，pet/ct检查的辐射量问题开始获得广泛关注，尤其在pet脑检查对于低龄儿童的辐射危害问题上，来自患者家属、医疗机构以及国际社会的担忧使得pet脑检查在儿童肿瘤和精神疾病的诊断与治疗中一直没有很大进展。我国《辐射防护规定》（gb8703-88）明确指出一名在医院放射科工作的医务人员每年辐射的接受剂量不能超过20毫西弗（msv），而研究结果表明患者进行一次pet/ct扫描接收到的辐射剂量大约是10msv，可见pec/ct扫描中降低受检者辐射剂量，其实是减少注入患者体内的放射性核素剂量，用低于标准剂量的放射性示踪剂获取相同质量的pet图像是提升pet/ct检查安全性亟需解决的问题。然而，随着低剂量pet扫描中给病人注射的示踪剂剂量的降低，扫描仪接收到的符合响应线数量严重不足，从而使传统基于统计学的迭代重建算法获取的pet图像中存在大量噪声，严重影响了临床医生的读片与诊断。
[0004]
为了抑制低剂量pet重建图像中的噪声并保留病灶信息，基于神经网络的深度学习方法已经被引入。目前，此类方法主要集中在图像后处理范畴，即在图像域建立低剂量pet图像到标准剂量的pet图像映射。由于低剂量图像的重建过程受噪声影响严重，原始pet数据中的有效高频信息并不能被完整保存，只对重建图像进行后处理无法利用这些有效信息，从而导致由深度学习泛化能力不足所产生的图像伪影和定量误差成为了此类技术路线的性能瓶颈。


技术实现要素：

[0005]
本发明的目的在于针对现有技术的不足，提供一种基于深度学习的低剂量pet数据三维迭代更新重建方法。该方法通过获取pet原始数据在图像域的无损三维反投影laminogram，利用三维深度神经网络拟合laminogram与pet重建图像之间的正、逆映射，从
训练样本中学习标准剂量先验知识并将其应用于迭代更新标准剂量laminogram和标准剂量pet三维重建图像，以获取比传统重建方法以及图像域降噪处理噪声更低，分辨率更高，无伪影的低剂量pet重建图像。
[0006]
本发明的技术方案是通过如下技术方案实现的：一种基于深度学习的低剂量pet数据三维迭代更新重建方法，具体包括以下步骤：（1）将低剂量sinogram和标准剂量sinogram，分别经过衰减、随机、散射校正、反投影及归一化，获取低剂量pet数据在图像域的3d表征和标准剂量pet数据在图像域的3d表征；（2）将步骤（1）获得的低剂量pet数据在图像域的3d表征和标准剂量pet数据在图像域的3d表征作为输入，以标准剂量pet重建图像作为神经网络映射的标签，训练神经网络，使用adam优化算法最小化所述神经网络的损失函数，获得从pet数据到重建图像的第一映射网络估计；（3）随后以所述标准剂量pet重建图像作为输入，所述标准剂量pet数据在图像域的3d表征作为神经网络映射的标签，训练神经网络，使用adam优化算法最小化所述神经网络的损失函数，获得从pet重建图像到数据的第二映射网络估计；（4）将待重建低剂量pet数据在图像域的3d表征作为第一映射网络估计的输入，经网络权值计算，获得第一标准剂量pet重建图像的估计；（5）将步骤（4）获得的第一标准剂量pet重建图像的估计作为第二映射网络估计的输入，经网络权值计算，获得标准剂量pet原始数据图像域3d表征的估计；（6）将待重建剂量pet数据在图像域的3d表征和步骤（5）获得的标准剂量pet原始数据图像域3d表征的估计作为第一映射网络估计的输入，经网络权值计算，获得第二标准剂量pet重建图像的估计；（7）交替执行步骤（5）-（6）n次，直至第n次获得的第n标准剂量pet重建图像的估计与第n-1次获得的第n-1标准剂量pet重建图像的估计的差值的绝对值小于，即获取最终标准剂量pet重建图像的估计。
[0007]
进一步地，所述反投影是将sinogram反投影至图像域获得高度模糊的pet图像laminogram；所述高度模糊的pet图像laminogram与pet重建图像存在以下关系：其中，和分别表示三维pet重建图像和高度模糊的pet图像laminogram上某像素点处的活度值，是旋转对称的点扩散函数的三维傅里叶变换。
[0008]
进一步地，所述神经网络的损失函数为：其中，、和分别表示laminogram或pet图像在水平、竖直和轴向的像素点总个数，表示标准剂量laminogram，表示低剂量laminogram，表
示准剂量pet重建图像，表示3d深度神经网络拟拟合的从和到的映射，是从到的逆映射。
[0009]
进一步地，所述神经网络由3d卷积层、3d反卷积层，以及它们之间的跳跃连接构成。
[0010]
与现有技术相比，本发明的有益效果是：本发明在迭代重建过程中利用低剂量pet原始数据中的有用高频信息，抑制了由于神经网络泛化能力不足所导致的伪影和定量误差的产生。此外，本发明使用神经网络实现petlaminogram与重建图像之间的正、逆映射，较传统迭代方法中基于系统矩阵的正、反投影过程，具有计算复杂度低，耗时少的特点，使用本发明所提出的重建方法实现一次低剂量pet数据重建所需时间远低于传统迭代重建方法的耗时。
附图说明
[0011]
图1是基于低剂量、标准剂量pet原始数据的正、逆映射网络训练方法流程图。
[0012]
图2是网络训练完成后，低剂量pet原始数据迭代更新重建方法流程图。
[0013]
图3是传统迭代重建算法与本发明基于深度学习的迭代重建算法对低剂量pet数据的重建结果对比图，其中，（a）为传统迭代重建算法对低剂量pet数据的重建结果图,（b）为本发明低剂量pet数据的重建结果图。
具体实施方式
[0014]
本发明提供了一种基于深度学习的低剂量pet数据三维迭代更新重建方法，具体包括以下步骤：如图1为基于低剂量、标准剂量pet原始数据的正、逆映射网络训练方法流程图：（1）将低剂量sinogram和标准剂量sinogram，分别经过衰减、随机、散射校正、反投影及归一化，获取低剂量pet数据在图像域的3d表征和标准剂量pet数据在图像域的3d表征；所述反投影是将sinogram反投影至图像域获得高度模糊的pet图像laminogram；所述高度模糊的pet图像laminogram与pet重建图像存在以下关系：其中，和分别表示三维pet重建图像和高度模糊的pet图像
laminogram上某像素点处的活度值，是用球面坐标表示的旋转对称的点扩散函数的三维傅里叶变换，为该旋转对称的点扩散函数的图像域表征，与卷积使得低频分量被放大，高频分量被衰减，从而得到平滑模糊的laminogram，为ramp型图像域三维滤波器，laminogram与之卷积可恢复出高分辨率pet图像。
[0015]
本发明提出用3d深度神经网络分别拟合图像域滤波器和点扩散函数。该深度卷积网络由3d下采样卷积层，3d上采样反卷积层，以及它们之间的跳跃连接（shortcuts）构成。其中下采样卷积层用来编码输入图像，提取多尺度特征，上采样反卷积层再对所提取的特征进行解码，生成输出图像。多尺度特征的提取使得神经网络拟合的正、逆映射网络具有空间变异性，与传统的具有空间不变性的点扩散函数和ramp型滤波器相比，可以更好地对pet扫描系统的空间特性以及pet光子探测过程的物理特性进行矫正。另外，通过shortcuts连接将反卷积层的输出与相应卷积层的输出进行叠加，实现在不增加网络参数的基础上，改善网络训练效果并有效防止退化问题。
[0016]
（2）将步骤1获得的低剂量pet数据在图像域的3d表征和标准剂量pet数据在图像域的3d表征作为输入，以标准剂量pet重建图像作为神经网络映射的标签，训练神经网络，使用adam优化算法最小化所述神经网络的损失函数，获得从pet数据到重建图像的第一映射网络估计；（3）随后以所述标准剂量pet重建图像作为输入，所述标准剂量pet数据在图像域的3d表征作为神经网络映射的标签，训练神经网络，使用adam优化算法最小化所述神经网络的损失函数，获得从pet重建图像到数据的第二映射网络估计。
[0017]
所述神经网络训练损失函数为：其中，、和分别表示laminogram或pet图像在水平、竖直和轴向的像素点总个数，表示标准剂量laminogram，表示低剂量laminogram，表示准剂量pet重建图像，表示3d深度神经网络拟拟合的从和到
的映射，是从到的逆映射。对网络的训练实现了用标准剂量pet数据辅助重建低剂量pet数据，而训练好的网络所学得的先验知识可在测试过程中实现对低剂量pet laminogram的补偿。
[0018]
如图2是网络训练完成后，低剂量pet原始数据迭代更新重建方法流程图：（4）训练完成后，输入已知测试对象的低剂量pet数据和训练好的第一映射网络估计、第二映射网络估计，通过联合迭代更新重建方法，将待重建低剂量pet数据在图像域的3d表征作为第一映射网络估计的输入，经网络权值计算，获得第一标准剂量pet重建图像的估计；（5）将步骤4获得的第一标准剂量pet重建图像的估计作为第二映射网络估计的输入，经网络权值计算，获得标准剂量pet原始数据图像域3d表征的估计；（6）将待重建剂量pet数据在图像域的3d表征和步骤5获得的标准剂量pet原始数据图像域3d表征的估计作为第一映射网络估计的输入，经网络权值计算，获得第二标准剂量pet重建图像的估计；（7）交替执行步骤（5）-（6）n次，直至第n次获得的第n标准剂量pet重建图像的估计与第n-1次获得的第n-1标准剂量pet重建图像的估计的差值的绝对值小于，即获取最终标准剂量pet重建图像的估计。
[0019]
所述联合迭代方法的目标函数为：首先用低剂量laminogram初始化标准剂量laminogram，，将输入网络得到pet图像的初始估计，然后（a）用作为输入带入网络更新，（b）再将输入更新得，依次重复步骤（a）、（b）直到。
[0020]
下面证明上述联合迭代方法得到的解与是目标函数（4）的极小解。根据，可推出将式（5）代入目标函数（4）得到
又因为所以，得证与是目标函数q的极小解。
[0021]
传统迭代重建算法对低剂量pet数据的重建结果如图3中的（a）所示，迭代重建过程包括多次投影、反投影操作，耗时长，重建图像噪声大，对于病灶的诊断受到噪声的严重干扰；使用本发明所提出的基于深度学习的迭代重建算法可获取如图3中的（b）所示的低剂量pet重建图像，图3中的（b）的信噪比明显高于图3中的（a），另外基于神经网络的正、逆映射的实施大大降低了传统迭代所需时间，使得重建算法更加符合临床实时性要求。