基于EM算法的动态PET参数图像分部重建算法制造技术

本发明专利技术公开了基于EM算法的动态PET参数图像分部重建算法，涉及医学成像技术领域，将耦合血流生理学模型与成像线性方程得到的参数图像成像线性方程进行拆解，将重建目标图像分解为两个目标图像，并通过交替进行的EM迭代图像重建算法重建斜率图像κ与截距图像b，本发明专利技术基于EM迭代图像重建算法和Patlak数学模型的分部直接动态PET参数图像重建方法，将Patlak曲线的斜率图像与截距图像在线性方程中分部分交替迭代重建，极大程度减小斜率图像与截距图像在迭代计算中的互相干扰，使重建图像具有更高的图像质量，极大程度减少参数图像重建过程中的误差引入，使所得参数图像更加有利于后续的诊断与分析。利于后续的诊断与分析。利于后续的诊断与分析。

【技术实现步骤摘要】
基于EM算法的动态PET参数图像分部重建算法


[0001]本专利技术涉及医学成像
，具体涉及基于EM算法的动态PET参数图像分部重建算法。

技术介绍

[0002]动态PET成像技术应用的关键是对病人生理代谢情况的定量分析。基于Patlak生理学模型对图像序列进行后处理是目前主流的动态PET数据分析方法。间接法重建参数图像指，使用图像采集得到的投影数据首先重建出所有动态PET图像序列，从动态图像中得到病人各组织的时间活度曲线(TAC)，结合时间活度曲线与成像过程中的血液输入函数，基于合适的方式模型和Patlak生理学模型，计算出病人各个组织的显像剂流入率Ki，即Patlak曲线有效区段的斜率。各像素点κ即构成动态PET参数图像。在间接法之外，动态PET学科前沿科学家提出从图像采集所得投影数据直接计算参数图像的方法，即直接法动态PET参数图像重建。这种方法由于不经过两次图像估算，相比较间接法参数图像重建技术，极大减少了参数图像计算结果中的误差和干扰信息。
[0003]间接法重建参数图像由于经过两次图像估算的过程，即动态PET图像重建与参数图像重建，参数图像重建结果中存在较大的误差与干扰信息。相关技术提出的使用最大后验直接重建参数图像的方法很大程度减少了重建参数图像中的误差与干扰信号。但是，由于这种方法将Patlak曲线的斜率图像κ与截距图像b整合为一个目标图像一次性求解，迭代中所要求解的未知数数量较大，导致重建参数图像质量仍然存在限制。

技术实现思路

[0004]为了解决现有技术中的问题，本专利技术提出基于EM算法的动态PET参数图像分部重建算法，将Patlak曲线的斜率图像与截距图像在线性方程中分部分交替迭代重建的方法，极大程度减小斜率图像与截距图像在迭代计算中的互相干扰，使重建图像具有更高的图像质量，更加有利于后续的诊断与分析。
[0005]为了实现以上目的，本专利技术提供了基于EM算法的动态PET参数图像分部重建算法，包括：将耦合血流生理学模型与成像线性方程得到的参数图像成像线性方程进行拆解，将重建目标图像分解为两个目标图像，并通过交替进行的EM迭代图像重建算法重建斜率图像κ与截距图像b。
[0006]进一步地，所述参数图像成像线性方程拆解为下式：
[0007][0008]其中，与分别为当前动态PET成像任务对应的血液输入函数C
p
(t)的时间积分函数和C
p
(t)的衰减积分函数；P为成像系统矩阵；r为投影数据中的散射与随机符合时间；为克罗内克积运算；k和b分别为待重建的Patlak曲线斜率图像与截距图像。
[0009]进一步地，所述时间积分函数表示为：
[0010][0011]其中，n表征某个时间帧，t
s,n
为时间帧n的起始时间；t
e,n
为时间帧n的结束时间；τ为积分时间变量；ξ为内层函数积分时间变量；λ为标记同位素的衰减常数。
[0012]进一步地，所述衰减积分函数表示为：
[0013][0014]进一步地，所述EM迭代图像重建算法包括如下斜率图像κ迭代公式：
[0015][0016]其中，N为目标参数图像像素数，1
N
为长度为N的元素值全为1的列向量；T表示矩阵转置。
[0017]进一步地，所述EM迭代图像重建算法包括如下截距图像b迭代公式：
[0018][0019]进一步地，所述EM迭代图像重建算法采用迭代公式交替迭代方式，迭代计算达到稳定收敛之后，即获得直接重建的Patlak曲线重建斜率图像κ和截距图像b。
[0020]进一步地，所述算法包括人体脑部、胸部、肝部、肾脏或全身动态PET参数图像。
[0021]进一步地，所述算法重建的斜率图像κ的PSNR值达到28.0636，截距图像b的PSNR值达到22.2039。
[0022]与现有技术相比，本专利技术将耦合血流生理学模型与成像线性方程得到的参数图像成像线性方程进行拆解，将重建目标图像{κ,b}分解为两个目标图像{κ}和{b}，并通过交替进行的EM迭代图像重建算法重建斜率图像κ与截距图像b，本专利技术基于EM迭代图像重建算法和Patlak数学模型的分部直接动态PET参数图像重建方法，将Patlak曲线的斜率图像与截距图像在线性方程中分部分交替迭代重建，极大程度减小斜率图像与截距图像在迭代计算中的互相干扰，使重建图像具有更高的图像质量，极大程度减少参数图像重建过程中的误差引入，使所得参数图像更加有利于后续的诊断与分析。
附图说明
[0023]图1为本专利技术的实验验证结果图。
具体实施方式
[0024]为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本专利技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本专利技术实施例的组件可以以各种不同的配置来布置和设计。因此，以下对本专利技术的实施例的详细描述并非旨在限制要求保护的本专利技术的范围，而是仅仅表示本专利技术的选定实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0025]本专利技术的实施例提供了基于EM算法的动态PET参数图像分部重建算法，包括：将耦合血流生理学模型与成像线性方程得到的参数图像成像线性方程进行拆解，将重建目标图像{κ,b}分解为两个目标图像{κ}和{b}，并通过交替进行的EM迭代图像重建算法重建斜率图像κ与截距图像b。
[0026]参数图像成像线性方程拆解为下式：
[0027][0028]其中，与分别为当前动态PET成像任务对应的血液输入函数C
p
(t)的时间积分函数和C
p
(t)的衰减积分函数；P为成像系统矩阵；r为投影数据中的散射与随机符合时间；为克罗内克积运算；κ和b分别为待重建的Patlak曲线斜率图像与截距图像。
[0029]时间积分函数表示为：
[0030][0031]其中，n表征某个时间帧，t
s,n
为时间帧n的起始时间；t
e,n
为时间帧n的结束时间；τ为积分时间变量；ξ为内层函数积分时间变量；λ为标记同位素的衰减常数。
[0032]衰减积分函数表示为：
[0033][0034]EM迭代图像重建算法包括如下斜率图像κ迭代公式：
[0035][0036]其中，N为目标参数图像像素数，1
N
为长度为N的元素值全为1的列向量；T表示矩阵转置。
[0037]EM迭代图像重建算法包括如下截距图像b迭代公式：
[0038][0039]EM迭代图像重建算法采用迭代公式交替迭代方式，迭代计算达到稳定收敛之后，即获得直接重建的Patlak曲线重建斜率图像κ和截距图像b。
[0040]本专利技术包括人体脑部、胸部、肝部、肾脏或全身动态PET参数图像，重建的斜率图像κ本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于EM算法的动态PET参数图像分部重建算法，其特征在于，包括：将耦合血流生理学模型与成像线性方程得到的参数图像成像线性方程进行拆解，将重建目标图像分解为两个目标图像，并通过交替进行的EM迭代图像重建算法重建斜率图像κ与截距图像b。2.根据权利要求1所述的基于EM算法的动态PET参数图像分部重建算法，其特征在于，所述参数图像成像线性方程拆解为下式：其中，与分别为当前动态PET成像任务对应的血液输入函数C
p
(t)的时间积分函数和C
p
(t)的衰减积分函数；P为成像系统矩阵；r为投影数据中的散射与随机符合时间；为克罗内克积运算；κ和b分别为待重建的Patlak曲线斜率图像与截距图像。3.根据权利要求2所述的基于EM算法的动态PET参数图像分部重建算法，其特征在于，所述时间积分函数表示为：其中，n表征某个时间帧，t
s，n
为时间帧n的起始时间；t
e，n
为时间帧n的结束时间；τ为积分时间变量；ξ为内层函数积分时间变量；λ为标记同位素的衰减常数。4.根据权利要求2所述的基于EM算法的动态...

【专利技术属性】
技术研发人员：陈子翔，胡战利，郑海荣，梁栋，刘新，杨永峰，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：