基于低秩稀疏分解的心肌灌注动态PET图像恢复-生物医学工程专业论文.docxVIP

优秀毕业论文 精品参考文献资料 为简单的图像恢复技术为高斯滤波，其滤波核在空间上是移不变的，无法保持 为简单的图像恢复技术为高斯滤波，其滤波核在空间上是移不变的，无法保持 组织结构的边缘细节信息，因此容易导致图像的模糊效应。对于边缘保持滤波 或非局部均值滤波等方法，其恢复效果更多的依赖于滤波过程中超参数的选取， 而如何选取参数目前仍是一个公开的问题。引入CT／MRI解剖信息的PET图像恢 复方法受PET图像与解剖图像之间匹配程度的影响较大，要求患者在两次扫描 过程中不能存在明显的运动幅度，实施难度较大。另外一种较为主流的动态PET 图像恢复策略HYPR考虑了合成图像的特性，利用合成图像的高信噪比改善时间 序列图像的信噪比。HYPR策略要求序列图像之间的空间位置必须一一对应。上 述几种图像恢复方法主要针对二维图像滤波，即依次对单帧PET图像进行恢复。 由于未考虑到动态序列中图像与图像间存在的关联，因此用于动态PET图像恢 复并不十分合适，恢复效果有限。鉴于此，有学者提出基于稀疏特性的图像恢 复技术，如引入总变分正则化项的图像恢复算法或是基于小波分解的图像恢复 方法。相比于单个时间帧的滤波方法，基于稀疏特性的图像恢复方法明显表现 更优。 目前来说，针对心肌灌注动态PET图像恢复的研究较少，临床上通常采用 传统的图像恢复技术进行恢复，尚无专用于心肌灌注动态PET图像的恢复方法。 然而，由于心肌灌注扫描区域的特殊性，传统的图像恢复方法并不完全适用。 考虑到心肌灌注动态PET图像的特点，本文重点开展基于心肌灌注动态PET特 性的图像恢复模型研究工作，并开展心肌灌注动态PET图像分析与评价指标的 研究，主要包括：1、基于不同评价指标分析心肌灌注动态PET图像的恢复效果： 2、应用动力学模型对动力学参数进行拟合估计并评价准确性。本文在经典评价 指标的基础上利用靶心图进行更为全面的评估。靶心图，又称为牛眼图，作为 一种心肌断层图像的定量分析方法，可定量评估心肌活力和心肌缺血程度，作 为心肌灌注动态PET图像的评价指标更具针对性。因此，利用靶心图对心肌灌 注动态PET图像及动力学参数图像(如K1参数图，k2参数图)进行评价更为全 面直观。归纳起来，本文的主要工作有： (1)针对心肌灌注动态PET图像特点，提出基于低秩稀疏分解的心肌灌注 动态PET图像恢复算法。由于动态PET序列图像不同时间帧的图像活度变化将 导致图像噪声的变化，使得动态PET图像的噪声分布变得更为复杂，因此针对 万方数据 单帧图像的图像恢复方法难以应用于整个动态PET 单帧图像的图像恢复方法难以应用于整个动态PET序列图像。本文从矩阵的观 点出发，将心肌灌注动态PET图像恢复问题作为一个矩阵恢复问题。矩阵的行 对应单个像素点的时间活度曲线(Time Activity Curve，TAC)，矩阵的列对应 单个时间帧的PET图像。结合心肌灌注PET显像独有的特点及低秩稀疏分解理 论，该矩阵可以分解为两个矩阵的叠加：第一个矩阵表示示踪剂缓慢变化的区 域，也称为背景区域(如肌肉和肺区)；第二个矩阵表示示踪剂快速变化的区域， 也称为目标区域(如心肌组织和血池)。相比于背景区域，目标区域的帧间活度 变化更为明显。连续的序列图像之间具有高度相关的冗余信息，因此图像矩阵 具有低秩性。而在序列图像的某些区域，由于含有血流灌注信息因此具有稀疏 性。也就是说，上述背景区域和目标区域分别满足矩阵的低秩特性和稀疏特性。 鉴于此，本文将低秩稀疏分解引入到PET图像恢复模型，具体研究包括：1、基 于低秩稀疏分解的心肌灌注动态PET图像恢复模型的构建及优化求解；2、评估 本文所提恢复模型的可行性和有效性。实验结果表明，基于本文恢复模型的动 态PET恢复图像噪声水平大幅度降低，组织结构清晰可见，展示了较为优质的 恢复效果。 (2)对心肌灌注动态PET图像进行动力学参数的拟合估计，并以靶心图作 为评价指标进行恢复效果评估。鉴于动力学参数在临床诊断上具有重大意义， 本文设计实现Rb一82 MP PET心肌灌注成像仿真实验，基于单房室模型、己知的 血输入函数和动力学参数Kl、k2对心肌、肺区等组织的时间活度曲线进行仿真。 在仿真图像数据基础上，根据恢复的动态PET图像，提取感兴趣区(如心肌组 织)的时间活度曲线。本文采用经典的动力学参数估计方法一一加权最小二乘 估计对感兴趣区动力学参数K1，k2进行拟合估计，并以靶心图作为评价指标评 估动力学参数估计的准确性，从而验证本文恢复模型的可靠性。最终结果表明， 基于本文恢复模型的动态PET恢复图像能对动力学参数进行较好的拟合估计， 可为临床医生的诊断提供更为精确的定量参数依据。 关键词：PET成像心肌灌注显像低秩稀疏分解 图像恢复动力学参数 万方数据 一 L