螺旋锥束工业CT的Katsevich精确重建算法研究

来源 :重庆大学 | 被引量 : 0次 | 上传用户:gsdx2009
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
CT(computed Tomography),即计算机断层成像,它是数学、物理、光电、机械和计算机等多学科相结合的产物。CT按其应用领域主要划分为医学CT和工业CT两大类。在近四十年的时间里,工业CT的扫描方式已由最初的二维平行束扫描发展到今天的三维螺旋锥束扫描。螺旋锥束工业CT不仅扫描时间短、射线利用率高、轴向分布均匀,而且它的扫描轨迹满足三维精确重建的完全条件,目前已成为无损检测领域的研究热点。螺旋锥束CT重建算法大致分为迭代法和解析法两大类。迭代重建需要多次迭代运算,计算比较耗时。解析法又分为精确重建和近似重建,两者比较而言,精确重建不受锥角的影响,而且重建的图像理论上没有伪影。精确重建算法主要有Grangeat算法和Katsevich算法两类。由于螺旋锥束Katsevich算法具有算法实现效率较高、重建图像质量较好等优点,目前它受到越来越多的关注。精确重建要求投影数据齐全、不含噪声。可在实际的螺旋锥束工业CT中,由于部分探测器失调、散射、串扰等原因,获得的实际投影数据含有大量的噪声。因此需要改进精确重建算法或对投影数据降噪,最大可能地满足精确重建的条件。本文对含噪的螺旋锥束CT投影数据进行Katsevich算法重建时,采取了两种策略减小噪声对重建图像的影响。第一种策略主要是对Katsevich算法中的求导方法进行改进;第二种策略主要是利用改进的自适应混合滤波对投影数据降噪。计算机仿真实验证实,第一种策略对投影数据含有少量的异常噪声和随机噪声具有较好的抑制作用;第二种策略更适用于含有大量噪声的投影数据。另外,在实际的螺旋锥束工业CT中,由于射线源焦点与探测器中心的水平高度不一致以及转心偏移等原因,可使获得的实际投影数据发生不同程度的偏移。如果直接利用偏移的投影数据进行图像重建,得到的重建图像与理想的重建图像相差较大。本文通过修正螺旋锥束Katsevich算法,以适应这种偏心的CT成像系统。首先将投影数据的偏移转换为探测器中心的偏移;在重建过程中,根据探测器中心的偏移量,对滤波线方程、滤波前后的重排和校正以及反投影过程分别进行了修正。其中,在反投影过程中,将PI线的传统迭代求根转变为一元函数的迭代求根,并利用高次迭代求根法对其进行了优化。计算机仿真结果证实,这种偏心修正方法可以有效地减小投影数据的偏移对重建图像的影响。同时,改进的PI线求根方法缩短了重建时间。图像的边缘往往包含了图像的重要特征,在图像的识别、分割、测量过程中具有重要的作用。通过工业CT获得物体三维边缘的传统方法一般分为两步:首先利用实际扫描的投影数据进行图像重建,然后对重建的图像进行三维边缘提取。本文对螺旋锥束CT的Katsevich重建公式进行了修改,研究了一种利用螺旋锥束CT获得的投影数据直接重建物体三维增强边缘的方法。通过该方法得到的三维增强边缘,可以初步判断物体内部结构的形状及类型,而且三维增强边缘重建时间与螺旋锥束Katsevich算法重建时间相当,节省了三维边缘提取的时间。
其他文献
明确课程的总体思路,增加实训内容,提出了具体的实训项目。 Clear the general idea of ​​the course, increase training content, put forward a specific training pro
反问题的研究是数学物理中一个较新的研究领域.声波反散射问题是一类典型的反问题,它在雷达和声纳等领域有广泛的应用前景.声波反散射问题被证明是不适定的,其求解起源于20世纪
本文主要讨论了在左截断假设下,样本为α-混合的情况下,密度函数的非线性小波估计的大样本性质.首先,在适当的条件下,给出了非线性小波密度估计的MISE渐近展开式,并且证明了非线性
随着时间序列相关技术日益广泛的应用,协整分析作为处理非平稳时间序列的一种技术逐渐受到人们的重视,引发了深刻的理论研究,展现出了广阔的应用前景。 协整向量的估计以及协
辅导员管理工作是高校日常管理的重要组成,精细化的管理模式是高校教育发展重要的新模式,能够促进辅导员管理工作的有效开展.当前我国高校的辅导员管理工作理念落后、工作团
本文研究一类Jacobi矩阵特征值反问题的敏感性分析。 第一章和第二章主要介绍了代数特征值反问题以及 Jacobi矩阵代数特征值反问题的基本问题和研究成果; 第三章和第四
n-李代数是李代数的推广,它是乘法运算为n元运算的一种多元李代数.我们知道n-李代数在物理及几何上都有它的背景,因此研究n-李代数的结构及应用是非常有必要的.本文主要研究Z2