理学ct统计重建论文Word下载.docx

1、 研究了基于罚似然（Penalized Likelihood，PL）的统计重建算法:OSL（One Step Late）- EM 算法和 PL-SPS 算法。 重点讨论了基于 Gibbs 分布的罚函数，从理论上分析了势函 数需要满足的条件以及势函数的选取对图像的影响，着重分析了二次势函数和 Huber 势函数的优缺点，通过仿真实验对它们的重建结果进行误差分析。 提出了罚似然重建中正则项参数的自适应选取的方法，该方法充分利用每一次迭 代的重建结果信息，不断对正则项参数进行更新。 并用于仿真实验和重建 CT 实际投影 数据，重建结果表明该方法降低了噪声，提高了图像质量。 将 OS-OSL-EM 算

2、法推广到三维锥束轨迹下的图像重建，取得了较好的仿真结果。 关键词：最大似然估计，CT 重建，罚似然估计，有序子集 The Study of CT Statistical Reconstruction Algorithm Based on Maximum Likelihood and Penalized Likelihood Estimates Abstract When the CT data has serious noise in the process of acquisition system and projection data is incomplete , analytic r

3、econstruction algorithm get images with artifact. The statistical reconstruction algorithm with a accurately physical model isnt sensitive to noise and easy to add constraints etc.Therefore, the reconstructed image quality is superior to conventional FBP methods. For the CT Statistics reconstruction

4、, This paper mainly studies the following content: The maximum likelihood estimate for the statistical reconstruction algorithms is studied and mainly includes expectation maximum algorithm （ML-EM）,separable paraboloidal surrogate（ML-SPS） algorithm and the both algorithms with ordered subset（OS-ML-E

5、M and OS-ML-SPS）. In simulation experiments, it shows that the initial convergence rate of ML- SPS is faster than the ML-EM ,however, the both algorithms is slower than OS-ML-EM and OS-ML-SPS algorithms. At last, OS-ML-SPS and OS-ML-EM algorithms are used to reconstruct actual CT projection data and

6、 get better image. Penalized maximum likelihood（PL） reconstruction for the statistical algorithms is based on the penalized likelihood estimates and adds a penalty term to suppress noise, which mainly includes OSL-EM algorithm and PL-SPS algorithm. It is focused on the penalty function which is base

7、d on Gibbs distribution and analyzes the need of potential functions to meet the conditions, the influence on image for the selection of potential functions and the advantages and disadvantages of the quadratic function and the Huber function. At last, simulation experiments give the error analysis

8、of reconstructed image. This method which is an adaptively regularied CT image reconstruction can adaptively choose regularization parameters with making full use of the results of every iteration to update regularization parameters. Simulation results show that this method reduces the noise and imp

9、rove image quality. The method is used to reconstruct the actual CT projection data to demonstrate the feasibility and practicality. OS-OSL-EM algorithm is applied to three-dimensional cone-beam image reconstruction and gets the better image. Key words: maximum likelihood estimate, CT Reconstruction

10、, Penalized maximum likelihood, ordered subsets I 目目 录录 第一章 绪论 1.1 课题研究背景及意义 1 1.2 国内外研究现状 1 1.3 论文的主要工作及结构安排 4 第二章 CT 统计重建的基础 2.1 CT 成像理论.5 2.2 统计重建的数学模型 7 2.3 优化变换原理 8 2.4 图像重建算法评价.9 第三章 基于最大似然准则的 CT 统计重建算法 3.1 最大似然准则的基本理论 11 3.2 ML-EM 算法.12 3.3 ML-SPS 算法14 3.4 有序子集的统计重建方法.17 3.4.1 有序子集 ML-EM17 3.

11、4.2 有序子集 ML-SPS19 3.5 实验结果及分析 20 3.5.1 ML-EM 算法与 ML-SPS 算法的重建结果21 3.5.2 OS-ML-EM 算法与 OS-ML-SPS 算法的重建结果.23 第四章 基于罚似然的 CT 统计重建算法 4.1 罚似然的基本理论.25 4.2 OSL-EM 算法 .31 II 4.3 OS-PL-SPS 算法.33 4.3.1 PL-SPS 算法.33 4.3.2 有序子集 PL-SPS 算法.34 4.4 自适应选择正则项参数 35 4.5 实验结果与分析.36 4.5.1 OSL-EM 算法实验结果分析.37 4.5.2 OS-PL-SPS

12、 算法实验结果分析 38 4.5.3 自适应选择正则项参数的实验结果分析40 第五章 锥束轨迹下有序子集罚似然统计（PL）重建 5.1 锥束扫描下的投影矩阵的生成 42 5.2 锥束下的 OS-OSL-EM 算法重建结果及分析.44 第六章 总结与展望 6.1 研究主要内容及成果 47 6.2 存在的问题及以后的工作展望 48 参考文献 攻读硕士学位期间发表的论文 致 谢 1 第一章 绪论 1.1 课题研究背景及意义 计算机层析成像技术是利用有某种能量的射线源对物体进行断层扫描，根据所获 得的某种物理量（指物质对射线的衰减系数），运用特定的重建算法，重建出物体某个断 层的无影像重叠的二维图像。

13、 自 1967 年 CT（Computed Tomography）在英国问世以来， 技术迅速发展，被公认为是 20 世纪影响人类发展的十大技术之一，其应用涉及到医学、 工业、石油、物探、材料、生物、安全等众多领域。 在国防工业中，CT 已成为航天器、 固体火箭发动机、高新技术等产品关键部件检测的必不可少的工具。 由于 CT 技术是发 展国防的基础技术之一，被发达国家列为敏感技术，作为国家机密进行严格的控制。 因此，CT 技术被视为国家科技实力的标志之一。 从数学上说，CT是从投影到重建图像的反问题1，具有普遍性，在数学界已经引 起了广泛的重视。 重建算法是工业CT技术中较为关键的一部分，CT重

14、建算法主要可以 分为滤波反投影重建算法和迭代重建算法。 当投影数据完备、噪声不很严重时，解析 法（滤波反投影（FBP）23等）可以得到很好的重建图像，但当CT系统在数据采集过程中 受到严重的噪声干扰，或者投影数据采集不完全时，用解析重建算法得到的图像有伪 迹。 而迭代重建算法（统计迭代算法和代数迭代算法等）可以用来重建不完全数据、动态 数据和噪声数据。 其中统计迭代重建方法的优势是其对噪声的鲁棒性，通过加入先验 模型和统计规律可以有效的抑制噪声，提高图像质量。 本课题主要研究统计算法在 CT 重建中的应用，通过研究有序子集的方法来提高收 敛速度，加入惩罚项，约束项来降低噪声的影响，提高图像质量

15、，为国内三维工业 CT 的研制提供理论基础和技术支持。 1.2 国内外研究现状 2 近十几年来， CT 广泛的应用在医学和工业上，相应的 CT 算法也迅速发展。 其 中解析算法重建速度较快，易于实现，是目前 CT 图像重建技术中最常用的算法。 但是 该算法通常要求采样数据是完全的，而且对噪声比较敏感，这在一定程度上影响了重 建图像的质量，限制了它在实际中的应用。 相对解析算法、代数迭代算法而言，统计 迭代算法能准确描述投影数据的物理模型，对误差不敏感，易于加入约束条件等优点 逐渐引起人们的重视，然而它与代数重建算法4一样重建速度慢，运算时间长，这些缺 点极大地限制了它的应用，但是随着加速算法和

16、计算机硬件技术的发展，统计重建将 在实际中得到了广泛地应用。 一个完整的统计重建框架包括图像的离散化系统的物理几何模型测量的统 计模型（泊松分布）优化准则（最大似然准则和罚似然准则）求解方法。 第步是考虑 如何对重建的图像进行建模，将待重建的图像离散成二维或者三维的像素矩阵，第 步考虑如何确定投影矩阵，与系统的几何结构和建模方法相关5-10，第步 ij Aa A 确定投影数据的模型，在统计重建中，一般认为投影数据服从泊松分布。 第步是引 入一定的优化准则。 第步在的基础之上进行问题的求解，该问题的求解实际上是 一个参数估计的数学优化问题。 统计重建的优化方法中涉及到优化变换（optimizat

17、ion transfer，OT），增量方法 （incremental method），有序子集加速（ordered subset，OS），这三种优化方法是统计重建 方法发展的主线11。 在统计重建中，最大似然估计（Maximum Likelihood，ML）是一种常用方法。 1977 年，A. P. Dempster 等人把期望值最大化（Expectation Maximization, EM）算法引入图像 重建之中。 1982年，Shepp12和Vardi13首次把期望最大迭代算法14应用到CT算法中， 使似然函数最大化，之后ML-EM算法得到了广泛的应用。 1997年，Fessler15等

18、人提出 了基于替代函数（surrogate function）的优化理论，其基本思想是在M步（maximization step）直 接对E步（expectation step）中所得到目标函数的替代函数作最大化处理，这一突破性的改 进在很大程度上应当归功于De Pierro的凸性（convexity）算法1617。 当然替代函数的选取 需要满足一定的条件，其它替代函数算法还包括Fessler18-21，Zheng22提出的算法以及 全局收敛的增量优化传递算法。 3 1994年，Hudson和Larkin提出了有序子集（Ordered Subsets） 23EM算法（OS-EM），大 大提高了

19、收敛速度。 OS-EM算法将投影数据分解成有限个有序子集，每次迭代时只处 理其中一个子集的数据，由于加速效果明显，OS-EM算法及其各种变形算法此后便被 广泛地应用。 其中RAMLA（row-action maximun likelihood algorithm）24重建算法的思想 来源于Herman的代数重建算法和Hudson等人的OS-EM算法，该算法也是将投影数据分 成一系列不正交的投影数据子集，并引入一个松弛因子，该算法可以从理论上证明其 收敛到ML的收敛点，并且其收敛速度和OS-EM算法差不多。 虽然子集类算法收敛速度 快，但不能全局收敛，迭代到一定次数会出现极限环，为了解决该问题，

20、S. Ahn 和J. A. Fessler提出了改进的块迭代EM算法25，E .Quan和D .S Lalush等人提出了基于快速子 集的凸算法在CT中的应用26。 最大后验概率重建算法（maximum a posteriori，MAP）2728是一种贝叶斯重建算法, 也可以称为罚似然重建（Penalized Likelihood，PL）29-32，该算法和 ML 重建算法可以说 是一对关系亲密的姊妹，因为 ML 算法是寻求合适的图像的估计值使得图像得到的已 知投影的概率最大；而贝叶斯重建算法是从已知的投影出发，要求在给定投影情况下 所求图像的概率最大，通过选择合适的先验分布模型来提高重建结果

21、的质量，其作用 等同于优化理论中的正则项（惩罚项） 。 正则化技术最早由 Tikhonov 在著作中提出， 他研究了求解不适定问题稳定解的基本理论，并提出了著名的 Tikhonov 正则化技术。 先验模型的引入虽然能改善重建效果，但同时增加了估计的难度。 就泊松观测模型而 言，它本身属于非线性模型，因此它的 ML 估计和 MAP 估计均没有闭型解析式。 针对 这一问题，Herbert 最先提出 GEM（generalized EM）算法33，Green 给出了经典的 OSL（one-step-late）算法34，Fessler35提出了罚似然的 SAGE（space-alternating g

22、eneralized EM）算法，另外 De Pierro 的3637凸性算法也涉及到了最大后验估计的问题。 在以后的 工作中，De Pierro 和 Yamagishi38将 RA（row-action）思想与最大后验估计的思想相结合， 提出了 BSREM（block sequential regularizedexpectation maximization）算法。 Chung Chan, Roger Fulton39等人提出了自适应结构先验在 CT 中的应用，根据图像结构块中灰度的 不同自适应的选择滤波。 Xiaochuan Pan40阐述了最小化 TV 正则项（total variat

23、ion）在图 像重建中的应用。 Jing Wang, Tianfang Li41等人把基于最小二乘的罚函数方法用于投影 数据的去噪，提高了图像重建的质量。 4 1.3 论文的主要工作及结构安排 本文 CT 统计迭代重建算法为研究课题，主要做了以下工作： 第一章主要介绍了课题的研究背景、意义，综述了 CT 及其发展历史，CT 技术的 发展意义及统计重建算法的研究现状。 最后提出了本文的研究内容，确定了本文的研 究方向。 第二章首先给出了 CT 统计重建的物理和数学基础，然后介绍了优化变换原理，最 后介绍了图像重建质量的评价标准。 第三章首先详细介绍了最大似然估计类的重建算法：ML-EM 算法和可

24、分离抛物面 型替代（Separable Paraboloidal Surrogate，ML-SPS）算法，它们的推导过程和实现步骤， 然后对有序子集的方法进行阐述，并详细介绍了 OS-ML-EM 算法和 OS-ML-SPS 算法 的实现步骤。 实验中对比了 ML-EM 算法和 ML-SPS 算法的实验结果，并用 OS-ML- EM 算法和 OS-ML-SPS 算法对 CT 采集的实际投影数据进行重建，得到了较好的重建 结果。 第四章先介绍了罚似然估计的理论知识，重点介绍了基于 Gibbs 分布的罚函数， 从理论上分析了势函数需要满足的条件以及势函数的选取对图像的影响。 然后介绍了 OSL （O

25、ne Step Late）-EM 算法和 PL-SPS 算法的推导和实现，接着提出了自适应选择正 则项参数的罚似然重建，通过对正则项参数的改进来提高图像的质量。 算法实现中选 用了二次势函数和 Huber 势函数来重建图像，对它们的重建结果进行误差分析，并用 提出的自适应选取正则项参数的方法来重建 CT 实际投影数据，取得了较好的结果，验 证了方法的实用性。 第五章介绍了 OS-OSL-EM 算法在三维锥束轨迹下的应用，首先介绍了锥束扫描 结构下的投影矩阵的计算，然后用 OS-OSL-EM 算法对仿真投影数据进行重建，并与 OS-OSL-EM 算进行比较分析。 第六章对本课题的研究内容进行了总

26、结，指出本文的创新之处和解决的问题，并 提出了尚未解决的问题，指明了今后研究的重点。 5 第二章 CT 统计重建的基础 在过去的几十年里，统计迭代法一直是断层重建研究的热点，其重建图像的过程就 是利用参数估计，优化变换原理等来逼近最优解一个过程。 在重建过程中，首先确定 能够正确反映成像系统的物理模型和统计模型以及输入输出关系，然后利用已知的条 件确定求解方法。 而且统计重建可以利用与物体相关的约束条件和其他的与特定系统 相关的边界条件，利用精确的物理模型和适当的统计模型，可以在数据缺失和不规则 采样的情况下重建出满意的图像。 2.1 CT 成像理论 CT 图像重建的原理是用穿透力强的射线扫描

27、被测物体，当一束强度大致均匀的射 线投照到物体上时，射线一部分被吸收和散射，另一部分透过物体沿原方向传播。 由 于物体的各种结构在密度、厚度等方面存在不同，它们对射线的吸收量也不相同，从 而使透过物体的射线强度分布发生变化，通过 CT 采集、转换得到投影数据，并运用相 应的数学模型，进行图像重建出二维或三维图像，最后通过算法将射线强度分布转换 成图像上的灰度分布，形成人眼可视的图像。 取一理想的 X 射线源，它发出极细的笔束 X 射线，在被测物体对面置一探测器， 如图 2.1（a）所示。 假设强度为的 X 射线穿过均匀分布衰减系数为的物体，行进了 0 Iu 的距离，强度变为，按 Beer 定理

28、有xI 00 ln（/ ） ux ii euxII 或 （2.1） 如图 2.1（b）所示，若物体是分段均匀的，衰减系数分别是,相应的长度为 1,23 ,.u u u 则下式成立： 1,23 ,.x x x 1230 .ln（/ ）uxuxuxII （2.2） 6 更一般，若物体在平面内为不均匀介质，则在某一方向 ，射线沿某一路径的总衰xylL 减值为 00 （ , ）ln（/ ）（ ）exp（ , ） LL x y dlIII LIx y dl （2.3） 其中表示物体在二维平面上的衰减系数。 是检测到的强度；是射线经（ , ）x y （ ）I LL 过的直线，是直线的弧长。 dl 若没有具

29、体路径，只说了沿着某一个方向，那么是投影，由等式（2.3）可知，dl 测得与，即可知道。 在 CT 中，入射的射线强度与出射的射线强度之比经过 0 IIdl 取对数运算后，则成了沿射线路径上衰减系数的线积分。 CT 重建问题实际上是给定一 个待重建物体的被测量的投影即线积分，然后运用算法去计算物体的衰减分布即dl 被积函数。 x o I 射线源 I 探测器 o I 射线源 I 探测器 1 x 2 x 3 x N x 1 2 3 N （a）均匀物体 （b）非均匀物体 图 2.1 单能 X 射线束在物体中的衰减图示 1917 年，丹麦数学家雷当（JRadon）研究已经为 CT 技术建立了数学基础，

30、从数学 理论上证明了某种物理参量的二维分布函数，由该函数在其定义域内的所有线积分确 定。 该结论指出了如果没有无穷多个且积分路径互不完全重叠的线积分，只能确定实 际分布的一个估计近似值。 有了上述数学、物理基础后，为了实现工程技术上的应用，还需要解决两个主要 问题。 首先是如何采集到检测断层衰减系数线积分的数据集，其次是如何充分利用该 数据集确定出衰减系数的二维分布。 解决第一个问题可采用不同的扫描方式，即用射 7 线束穿过被测体所检测断层并相应进行射线强度测量的规律性，可采用各种不同的扫 描检测模式围绕提高扫描检测效率。 解决第二个问题则是选择合适的图像重建算法， 一类是解析法，以滤波反投影算法为代表。 另一类是迭代重建法，有代数重建和统计 重建。 统计重建是本论文研究的主要方向，下面具体介绍该算法。 2.2 统计重建的数学模型 在 CT 图像重建中，图像区域是重建区域，首先应当明确，图像区域是一个正方形， 其中心在坐标原点；图像函数是一个二元函数，其值在图像区域之外为零。 一个具有 个元素的网格把这个图像区域分成 2 n个相等的正方形，每个小正方形（像素）内图像函n 数的值是均匀相等的。 一幅图