锥束CT迭代重建方法的改进与加速

来源 :东南大学 | 被引量 : 0次 | 上传用户:tangguoxun3726
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
锥束CT是三维断层成像技术的代表,它具有扫描速率快,辐射利用率高,可获得高精度的空间分辨率等诸多优势,在医学和工业检测等诸多领域具有广泛的应用。  X射线对人体健康有一定的危害,X射线的照射时间越长,照射剂量越大,对人体的危害越大,且在CT扫描期间病人自主或非自主的运动会造成伪影影响重建结果,因此减少X射线照射的时间和剂量十分有意义。对待测成像物体进行稀疏角度投影数据测量,可以有效减少数据采集的时间,减少X射线的辐射剂量。CT重建的算法包括解析重建算法和迭代重建算法两类。解析算法对噪声比较敏感且要求投影数据完备,无法解决数据缺失的情况。迭代重建算法对投影数据的完备性要求较低,可以根据具体成像条件引入先验约束信息,可以抑制噪声,使得重建图像更为准确,效果更好。但是迭代重建算法的时间和空间复杂度高,速度较慢,因此实际应用并不广泛。  在迭代重建中,投影和反投影的计算时间最为耗时,远远超过其他步骤,因此本文使用Nvidia公司推出的统一计算设备架构(CUDA)环境,在GPU上实现投影和反投影的快速计算。对于投影的快速计算,对传统射线驱动的投影算法进行改进,采用固定的采样点个数,使用一个线程块计算一条射线的投影值,线程块中的每个线程计算一个采样点的值,并使用共享存储器实现采样点的规约相加。将体数据绑定为三维纹理,利用纹理拾取加速体数据的滤波采样。对于反投影的快速计算,改进传统的基于体素驱动的反投影算法,并使用二维纹理绑定投影数据。经实验证明,在保证全浮点运算精度的情况下,本文提出的算法与传统GPU加速投影和反投影算法相比,速度有大幅提升。  MAP重建算法基于贝叶斯估计,既考虑投影数据的统计特性,更加符合实际的数据采集过程,又充分利用了先验知识,利用先验对图像中的噪声施加一定的约束。压缩感知理论证明了,当采样信号不满足Nyquist采样定理时,稀疏信号仍可以被恢复出来,它为稀疏角度投影数据的图像重建问题提供了一条解决途径。根据压缩感知理论,本文采用了图像稀疏变换后的l0范数作为先验项,通过一系列势函数对如范数进行逼近,并设计了一个渐进的优化策略,采用自适应的参数,根据图像局部方差自适应计算先验项,可以在平滑图像的同时加强图像的边界。通过计算机仿真实验证明,该算法较TV算法更能保持图像的细节,更接近原始图像。
其他文献
网格是下一代Internet上的计算平台,其核心任务是管理分布在Internet广域环境中的各种类型的数据与服务资源,并为基于Internet的分布式应用提供一个统一的、虚拟的共享资源的计
随着计算机网络技术的飞速发展,各种网络设备、网络协议和网络服务推陈出新。网络的维护和管理面临着新的挑战,传统的集中式网络管理已经不能满足这种大型复杂网络的需要。移动
随着计算机和互联网络技术的普及应用和飞速发展,业界的分析家估计,对存储容量的需求以后会以每年50%—100%的速度递增;对众多企业而言,由于其工作性质的要求,致使数据必须在
数控加工是CAM的主要组成部分,数控技术及计算机辅助设计的发展使传统的机械设计与制造发生了根本性的变化。在开发基于PMAC卡的数控加工上位机监控系统,完成对加工过程的实时
MMC卡(多媒体存储卡)是一种广泛应用的存储设备,它可以作为单片机系统的外存储设备。本文所讨论的用于单片机的MMC卡的文件系统,可以提供一组文件操作函数,通过调用这些函数,单片
语义复合研究是一项面向计算机的以解决短语语义理解问题为目标的工作,直接地涉及到语言学和计算机科学两个学科,是对信息内容进行语义理解的重要环节。本文以概念网络模型为基
数字水印技术作为知识产权与信息安全保护的有利武器,成为当前研究热点。本论文主要对数字图像水印技术及其应用进行了研究,目的是解决其中两个问题:盲水印鲁棒性和不可见性
本文对互联网信息搜索技术进行研究,并在此基础之上建立一个基于分布式的互联网信息深度搜索系统。 首先分析当前普遍采用的互联网信息搜索引擎的优缺点,鉴于互联网信息搜索
一般说来,离群点是远离其他数据点的数据,但很可能包含着极其重要的信息。对离群点进行识别是相当重要的,它不但能够提高分析的质量,而且能够发现和研究离群点。本文对算法如
外包镍丝石墨条是应用于航空航天及核电核能领域的特种密封产品,该产品生产工艺独特,受专利保护及商业性限制,该产品的研发生产在国内尚不多见,其破损检测更属空白,目前对该产品的