基于不同先验获取的PET图像优质重建新方法研究

摘要

作为一种有效的功能医学图像成像方法，正电子发射成像(Positron EmissionTomography，PET)能够在分子水平上利用影像技术反映人体心脑代谢和功能，已经在肿瘤学，心血管疾病学和神经系统疾病学研究中，以及新医药学开发研究等邻域中显示出它卓越的性能。可是在实际中由于受到低计数率和一些物理噪声的影响，PET、图像的重建问题在理论上是一个病态的问题。传统的滤波反投影(Filtered Back-Projection，FBP)重建方法虽然具有成像速度快的优点，其重建图像却含有大量噪声，图像质量较差。 最大似然一期望最大法(Maximum-Likelihood Expectation-Maximization，ML-EM)能够针对系统模型的物理效应和探测数据和噪声的统计泊松特性建立数学模型，其重建的图像质量要优于传统的FBP方法。然而，单纯的传统ML-EM方法依然有两个缺点：1，在迭代过程中，会产生质量退化的图像而导致的棋盘效应，从而导致非收敛的迭代过程；2，收敛速度较慢，需要大量迭代次数才能重建出可以接受的图像。 近二十年来，国际上出现了很多解决以上两个问题的方法。一方面，基于马尔可夫随机场(Markov Random Fields，MRF)的贝叶斯(Bayesian)重建方法或者最大化后验估计(Maximum A Posteriori，MAP)的方法通过引入正则化项来引进目标同位素密度数据在空间上的概率分布的先验信息，能够明显改善重建图像质量以及迭代过程的收敛性，该方法已被证明了其在理论上的正确性和实际上的有效性。另一方面，很多科学家也提出一些有效的方法和算法来解决ML-EM方法收敛速度慢的问题，1994年，H.Malcolm Houdson和Richard s.Larkin提出使用有序子集(Ordered Subsets，OS)的方法将完整探测数据sinogram分割成有序部分数据的方法来减小每次迭代的运算量，从而达到了加快收敛速度的目的；密歇根大学的J.A.Fessler小组提出了能够提高收敛速度的SAGE(Space-Alternating Generalized EM)收敛重建算法和PSCD(ParaboloidalSurrogate CoordinateAscent)，收敛重建算法；另外，Erkan U Mumcuoglu，David S．Lalush及Fessler也分别提出将共轭梯度法(Conjugate Gradient，CG)应用于PET重建算法中以加快收敛速度。对于重建速度的问题，虽然能够反映探测数据统计特性的迭代重建能够重建出较传统的FBP重建更好的图像，然而其应用一直受到重建速度慢的缺点的限制；在重建图像质量方面，作者发现，sinogram探测值中大量的噪声一直对重建有着很大的负面影响，而且此种负面影响持续贯穿于整个在迭代过程中。虽然引入了图像先验信息的Bayesian能够在很大程度上改善迭代重建，但是依赖于传统的局部邻域信息的，Bayesian方法只能为重建提供有限的局部先验信息，一方面，传统的平滑二次QM(Quadratic Membrane)先验易在重建的结果图像中产生过平滑效果，另一方面，具有边缘保持作用的非二次先验则会给重建图像带来阶梯状的边缘伪影。 本文对于PET重建算法的研究工作同样也是基于如何进一步加快重建算法的收敛速度和进一步提高重建图像的质量。作者做了以下四项PET重建算法方面的工作： 1，提出将有序子集的思想和共轭梯度法相结合的新的快速子集共轭梯度重建算法(OSCG Ordered Subsets Conjugate Gradient)来加快PET重建的收敛速度。 2，提出基于修改sinogram探测值的新的耦合反馈(Coupled Feedback，CF)迭代模型来降低在重建迭代中探测数据sinogram中的噪声对PET重建的负面影响。 3，提出一种新的有效的综合了QM先验和QP(Quadratic Plate)先验的MRF二次Membrane-Plate混合先验模型，实现了在Bayesian重建中根据目标图像中各点的性质来自适应的选择QM先验和OP先验。 4，提出一个新的且具有二次先验能量形式的MRF非局部(nonlocal)先验模型，实现了利用图像中的全局的信息来为PET重建提供更加有效的正则化作用。 试验中将以上四种方法应用于相应的PET断层发射成像(EmissionTomography)和断层透射成像(Transmission Tomography)，相关试验分析表明：本文所提出的四种基于MRF和优化理论的四种新算法均能够分别在不同程度上改善PET重建。

目录

文摘

英文文摘

声明

第一章引言

1.1 PET 简介

1.2 PET成像中的问题

1.3作者的主要工作

1.4主要创新

参考文献

第二章PET成像原理及常用重建算法

2.1 PET结构与成像

2.1.1 PET的历史

2.1.2 PET成像的原理

2.1.3 PET的结构及数据的获取

2.1.4衰减、散射及随机计数

2.2 Emission断层成像和transmission断层成像

2.2.1 Emission断层成像

2.2.2衰减信息和transmission断层成像

2.3解析法PET重建-FBP和CBP

2.4基于统计模型的迭代法PET重建

2.4.1符号和基本概念

2.4.2 ML-EM迭代重建算法

参考文献

第三章基于快速子集共轭梯度法的PET图像Bayesian重建

3.1应用于Bayesian图像重建的共轭梯度算法

3.2有序子集共轭梯度(OSCG)重建算法

3.3模拟实验

3.4结论

参考文献

第四章基于纠正正弦图探测值的PET图像Bayesian重建

4.1 CF重建模型的理论

4.2 CF重建模型在Emission断层图像重建中的应用

4.3模拟实验

4.4结论

参考文献

第五章混合MRF Membrane-Plate先验PET重建

5.1 MRF混合多阶先验模型

5.2基于QHM先验的图像Bayesian重建

5.2.1联合Bayesian估计算法

5.2.2图像f的收敛Bayesian重建算法及参数ρQP的选择

5.2.3基于ML-EM算法的Gamma分布参数估计算法

5.2.4阈值KG的确定

5.3在PET图像重建中的应用及相关重建试验

5.3.1同FBP重建、QM先验和QP先验的实验比较

5.3.2同TV正则化和MRP先验的实验比较

5.4结论

参考文献

第六章基于非局部二次MRF先验的PET图像重建

6.1非局部MRF先验模型理论

6.1.1传统的局部先验模型

6.1.2非局部二次先验模型

6.2基于非局部MRF先验模型的PET图像重建算法

6.3试验与分析

6.3.1模拟体摸数据重建

6.3.2真实数据重建

6.3.3平均标准差(average standard deviations,astd)、偏差(bias)和对比度(contrast)实验

6.4结论

参考文献

第七章结论与进一步工作的展望

攻读学位期间成果

致谢