CT迭代重建算法的加速方法的研究.pdf

摘要研究生姓名:导师姓名:学校名称:�迭代重建算法的加速方法研究鲍旭东教授论文题目:黎晖东南大学�成像作为放射性断层成像的主要手段,在医学诊断、辅助治疗和工业无损检测等方面都有广泛应用。经过几十年的研究,传统的滤波反投影方法已渐趋成熟,目前的�技术对完整均匀的无噪声投影数据能以较快的速度完成重建。抵抗噪声的�成像重建方法主要是迭代重建法,它的缺陷是重建速度很慢,远远低于应用需求,对该类算法的加速非常值得探讨。本论文研究两种对�迭代重建算法进行加速的方法,并通过模拟实验来分析它们的主观效果和客观误差。第一部分工作是基于统计迭代投影估计的快速�重建算法。考虑到传统的图像空间迭代法中对庞大的系统矩阵的计算和处理都极为耗时,尝试在投影空间对噪声投影进行统计估计和迭代求解,获得所有方向上的理想投影数据,并结合具有速度优势的解析法来完最后的重建。首要的问题是对噪声建立统计模型,采用了针对校准后的数据建立的不同于传统泊松描述的高斯模型,在该模型下,用贝叶斯方法推导了最大后验概率���兰频哪芰亢����被雇频�了惩罚加权最小均方���估计的能量函数,在采用相同的二次惩罚项的情况下,得到了类似的结果。考虑到本文所处理的估计问题中,初始值可信度高,又是唯一可能的初始点,采用了迭代条件模式���呕�惴ǘ訫�目标函数进行迭代求解。完成投影估计步骤之后,再进行滤波反投影步骤的实验。实验中产生两组不同信噪比的随机噪声数据,分别进行�次迭代估计和�次图像重建,并对每次的投影估计和图像重建结果都进行主观比较和客观误差分析。同时分别选取一种常用的线性滤波器和非线性滤波器进行对照。实验的图像结果以及总的误差趋势分析和比较表明,在选择合理的系统模型与优化算法的情况下,该方法能获得理想的估计投影,信噪比得到有效提升,误差很快地减小,进而有效地改善了图像重建质量。由于不需要处理系统矩阵,并且选取了收敛相当快的优化算法,该统计方法的时间代价极小。从误差图中发现,选择较少的迭代次数可避免少数点的误差累积。正则因子��的选取对投影估计的迭代后期影响较大,正则因子较小时,估计结果较好,而对重建结果来说,由于选择了较少的迭代次数,其影响并不显著。第二部分工作是基于傅立叶方法的快速重投影。重投影是一般迭代法中必不可少、反复进行的步骤,也是迭代法的速度瓶颈所在,对这一过程实现加速很有意义。考虑到�成像的基本原理与直接傅里叶重建法有关,讨论了近年来非均匀傅里叶变换����际踉谥苯痈道镆�重建中的研究进展,引入该类技术对迭代算法的重投影步骤进行加速,并与标准的系统矩阵重投影法进行比较。同时用传统的叔线性插值法实现重投影作为傅立叶方法的对比。在����投影算子的实现中,利用�������砘�擞呕�矫娴男陆�菇�星跋蛲队八阕雍头聪蛲队�算子的加速,与最近的一种从最小化最坏近似误差��.��判据直接导出的卷积核的重投影实现进行比较,给出了实验结果和误差对比。东南大学硕士学位论文
由实验结果可以看到,第一,与传统的系统矩阵计算方式相比,基于���母道镆斗椒�关键词:抗噪性�成像,迭代重建,加速算法,投影估计,高斯统计模型,傅里叶重投影,在前向投影和反投影算子的实现上有明显的速度优势。第二,在采用优化���.���插值器以及��甅�插值器的情况下,其计算代价仅接近于传统的双线性插值傅里叶方法,而投影误差则明显降低,其归一化均方误差和归一化绝对误差均等于或小于标准的重投影方式�聪低�矩阵方法��又鞴弁枷窠峁�弦惨芽梢钥吹搅钊寺�獾谋硐帧R虼吮疚娜衔2捎没�贜��的傅里叶方法进行重投影来实现对迭代中的重投影步骤加速是可行的,能够以小于标准的系统矩阵重投影方法的时间代价来保障重投影精度,适用于不需要进行分块运算的任何�迭代重建算法。����逯涤呕�砘�恕�东南大学硕士学位论文�
甌���������������������,�������—����������������������������������������.�����东南大学硕士学位论文���������:������:��.��������:����瓹����������瑃����琭��������.������.���琻������琣����.���,��
篊���������,�����—�����东南大学硕士学位论文�����������������������������,�����,���������,�����.
东南大学学位论文独创性声明东南大学学位论文使用授权声明括刊登�畚牡娜ú炕虿糠帜谌荨B畚牡墓ú�包括刊登�谌ǘ�洗笱а芯可�喊炖怼�本人声明所呈交的学位论文是我个人在导师指导下进行的研究工作及取得的研究成果。尽我所知,除了文中特别加以标注和致谢的地方外,论文中不包含其他入已经发表或撰写过的研究成果,也不包含为