图像重建块迭代算法中松弛参数选取的研究.PDFVIP

中国科技论文在线 第16 卷 第2 期 CT 理论与应用研究 Vol. 16 No. 2 2007 年6 月 14-19 CT Theory and Applications Jun., 2007 文章编号 100402-0014-06 图像重建块迭代算法中松弛参数选取的研究 孙杰 渠刚荣 北京交通大学理学院 北京100044 摘要 在图像重建中 Landweber迭代算法是图像重建算法中的重要方法 本文将针对Landweber 分块迭代算法中松弛参数的选取进行研究 在重建过程中采取对投影矩阵按投影角度分块的方 法 选取特定的松弛参数 通过数值实验得出结论 对于按角度分块的块迭代算法 松弛参数 选取为 乘以块矩阵与其共轭转置矩阵乘积的最大特征值分之一 当采集完全投影数据 且 接近(1/6) (1/7)时效果最好 另外 本文按角度分块的做法和松弛参数的选取方法对于有限 角度图像重建问题也是可行的 关键字 代数重建算法 松弛参数 块矩阵 投影矩阵 有限角 中图分类号 TP391.41文献标识码 A 目前 图像重建算法中主要的两种是解析重建算法和基于迭代的重建算法 解析重 建算法有滤波反投影重建算法和 Fourier重建算法等 滤波反投影算法对数据的采集有一 定的要求 如果投影数据不是均匀分布的 或者无法得到足够多的数据等情况下 滤波反 [1] 投影重建算法就不如代数重建算法具有优势 文献[2]研究了迭代松弛参数的选取在一般 Landweber迭代算法 包括同时迭代算法 简称SART算法 和块迭代算法中的收敛性问题 文献[3]用矩阵的奇异值定理研究了 SART迭代算法对更广泛松弛参数选取的收敛性 怎样 选取松弛参数加快收敛速度 并改进图像的质量是人们所关注的 但是 现在我们没有好 的方法去选择松弛参数 可靠的方法是通过实验验证 这个问题在理论方面需要进一步的 [4] 研究 本文在所选松弛参数满足文献[2]和[3]中块迭代算法和SART算法收敛性的条件下 通过数值实验 研究了投影矩阵按角度分块的块迭代算法中松弛参数的不同选取对迭代算 法的收敛速度和成像质量的影响 实验表明当按本文采集模式采集完全投影数据时 对于 标准shepp-logan图或其他图像 松弛参数选取为 (0 1 R)乘以块矩阵与其共 轭转置矩阵乘积的最大特征值分之一的重建效果比直接选取为块矩阵与其共轭转置矩阵乘 积的最大特征值分之一要好 且 接近(1/6) (1/7)时效果最好 重建图像与原始图像之 间的平均误差最小 按角度分块的方法和松弛参数的选取方法 对有限角度图像重建问题 有启发作用 我们也做了相关数字实验 并且得到可以接受的结果 收稿日期 2007-04-20 基金项目 国家自然科学基金项目助 转载 中国科技论文在线 2 期 孙 杰等 图像重建块迭代算法中松弛参数选取的研究 15 1代数重建算法以及按角度分块 开始 将图像离散化 并固定J 个基图像 b b b 1 2 J 1 bj x , y