基于最大后验概率的PET图像重建算法研究
The Research of PET Image Reconstruction Algorithm Based on Maximum A Posterior
正电子发射断层成像(Positron Emission Tomography,PET)是继计算机断层成像(Computed Tomography,CT)和磁共振成像(Magnetic Resonance Imaging,MRI)之后应用于临床的一种新型影像技术,现已广泛地使用于肿瘤细胞的探测、心脏病的诊断、神经和精神类疾病的诊断以及新药物的开发等领域.PET成像的目的是得到一个放射性物质在人体内部的分布图,因此,如何根据扫描数据来重建出高质量的图像,一直是PET领域的一个重要研究课题.总体来说,PET图像重建算法可以分为解析法和迭代法两类.解析法的代表是以中心切片定理和傅立叶变换为基础的滤波反投影算法,它具有计算简单,成像速度快等优点.但是当投影数据中含有大量噪声时,解析法很难重建出令人满意的图像,从而会影响临床诊断的效果.迭代法可以分为代数迭代法和统计迭代法两类.其中,代数迭代法的工作原理与解析法类似,由于其在重建图像的过程中较难引入各种物理成像条件及统计模型,因此很难重建出高质量的图像,故此方法在PET图像重建中使用较少.统计迭代法建立在观测数据的统计模型基础上,能够重建出高精度的重建图像,是目前PET图像重建中使用最广泛的一种方法.统计迭代法中经典的PET图像重建算法有很多,如最大似然期望最大算法(Maximum Likelihood Expectation Maximized,MLEM)、有序子集期望值最大算法(Ordered Subset Expectation Maximization,OSEM)和最大后验概率算法(Maximum A Posterior,MAP)等.本文的主要研究内容是MAP算法,又称为惩罚最大似然算法或Bayesian算法.本学位论文共五章,具体内容安排如下:在第一章中,我们首先介绍了PET成像技术的背景及意义,然后回顾了PET图像重建算法的历史与发展概况,最后简单地概述了本文的主要研究内容与结构安排.第二章是PET图像重建算法的基础理论知识部分,主要介绍了PET中的一些经典重建算法及它们的优缺点.剩下三章是本文的主要工作,所取得的研究成果如下所述.在第三章中,我们通过把各向异性中值扩散滤波器(Anisotropic Median-Diffusion,AMD)融合到中值根先验算法(Median Root Prior,MRP)中,提出了一种新的Bayesian图像重建算法.由于中值滤波器对高斯和Poisson两种噪声的抑制效果不明显,所以MRP算法很难取得令人满意的重建结果.新算法通过把AMD滤波器融合到MRP算法中,有效地抑制了重建图像中的所有噪声.模拟仿真实验结果表明,新算法在抑制噪声和保护边缘两方面取得了良好的折中,较大程度地提高了重建图像的质量.在第四章中,我们通过把AMD模型引入到PET图像重建算法中,提出了一种基于惩罚最大似然的PET图像重建算法.通过实验比较可知,新算法能取得较好的重建结果.此外,跟类似算法相比(如MLEM-PDE),新算法由于吸收了AMD模型的优点,参数设置简单,迭代过程中对梯度阈值和扩散次数的值不敏感,实用性较强.在第五章中,我们通过把一种修正的全变分模型(Total Variation,TV)和MLEM算法结合起来,提出了一种新的基于惩罚最大似然的PET图像重建算法.PET图像中的噪声主要是以Poisson噪声为主,一些传统的PET图像重建算法,如MLEM、MRP和MAP等,它们对一般的加性噪声有较好的抑制效果,但是对与信号相关的Poisson噪声的抑制效果却很不理想.新算法通过把PMTV模型(Poisson-modified Total Variation,PMTV)引入到MLEM算法中,有效地抑制了重建图像中的Poisson噪声,较大程度地提高了重建图像的质量.
- 作者：
- 何骞
- 学位授予单位：
- 湖南大学
- 专业名称：
- 应用数学
- 授予学位：
- 博士
- 学位年度：
- 2016年
- 导师姓名：
- 黄立宏
- 中图分类号：
- TP391.41
- 关键词：
- 图像重建；正电子发射断层成像；最大后验概率；惩罚最大似然；Bayesian算法；各向异性中值扩散；全变分模型
- Image Reconstruction; Positron Emission Tomography; Maximum A Posterior; Penalized Maximum Likelihood; Bayesian algorithm; Anisotropic Median-Diffusion; Total Variation Model;