基于GPU并行计算的X射线动态显微CT的快速重构

来源 :中国科学院大学(中国科学院上海应用物理研究所) | 被引量 : 0次 | 上传用户:jimmycjriyue
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
近年来,随着第三代同步辐射装置的广泛使用,X射线成像线站的光子通量密度得到了极大的改善,为实现X射线动态显微CT成像提供了条件,X射线动态显微CT成像技术趋向于更高的空间分辨率和更快的时间分辨率。上海光源成像组基于长工作距离显微镜透镜系统和互补金属氧化物半导体(Complementary Metal-Oxide-Semiconductor,简称CMOS)数字高速相机开发了快速X射线成像探测器,利用这个探测器可实现时间分辨率25Hz的X射线动态显微CT成像。一天用户实验产生的CT数据量大约500套,上海光源X射线成像线站的X射线动态显微CT用户实验在短时间内将会产生海量的投影数据。当前成像线站采用基于CPU重构程序进行CT重构,由于CPU重构程序是串行计算,重构耗时较长。因此,图像重构速度远比原始投影数据的产生速度要慢,实验原始数据将大量积压,如何实现快速显微CT重构是当前X射线动态显微CT实现广泛应用所面临的挑战。图形处理器(Graphics Processing Unit,简称GPU)因其多核架构、可大规模并行化特点在高性能计算领域得到广泛流行,将其用于CT重构可实现X射线动态显微CT的快速重构。本文采用CUDA编程标准,将滤波反投影(Filtered Back Projection,FBP)CT重构算法的反投影运算实现了并行化,并在NVDIA RTX2080的GPU上成功实现了CT重构的并行计算。论文取得的结果如下。1)基于GPU并行计算实现了FBP重构算法的加速,相比于CPU串行运算的CT重构算法,GPU并行重构算法在重构时间上提高了近200倍,将一套CT数据(800帧)的重构时间由原来的分钟级减少到了秒级,显著提升了海量CT实验数据的图像重构效率。2)发现在GPU并行运算和CPU串行运算的CT重构中采用图像反相算法对微弱信号有一定的放大作用。CT重构过程中,对于样品中细致微小的结构,重构结果信号比较微弱,重构图像的灰度值较小。借助图像反相算法,对于原本灰度值较小的像素点在一定程度上可以得到放大,进而达到增强微弱信号的目的。经过实验验证,采用图像反相算法可以提高重构切片的衬度,三维重构图像明暗变化更大,细节信息更明显。3)研究了反投影重构过程中的星状伪影现象。基于GPU并行计算重构分别选择不同的滤波器处理了同一套投影数据,对比了不同滤波器的性能差异,实验结果发现,采用Hamming滤波器滤波器处理得到的三维CT重构结果伪影较少,图像重构质量较好。4)设计了动态显微CT实验并搭建了10 Hz时间分辨率的动态显微CT实验平台对基于GPU并行计算和Hamming滤波器的CT重构算法进行了验证。实验结果证明,50套黄粉虫体内气泡CT数据完成重构总共耗时5.93 min,而基于CPU串行重构耗时20.4 h。通过选取不同时间点的气泡CT重构数据,在时间维度上观察到气泡的动态变化过程,并对运动气泡做了定量化分析,经过计算,通过计算,黄粉虫体内的气泡半径为0.624 mm,气泡平均运动速度为1.08 mm/s。5)在对CT重构加速的同时,结合相位恢复算法,获得了样品的相位信息,从而准确获取了样品的三维结构信息。实验结果证明,通过相位恢复算法提高了CT重构结果的图像衬度,可以明显区分开来黄粉虫的组织结构,图像局部区域的像素值非常接近,通过普通的图像分割算法就可以将内部组织结构进行提取和划分,这对后期图像信息的分析有着重要价值。综上所述,本文基于GPU并行计算对FBP算法中的反投影过程做了并行程序的开发,并结合了图像反相算法、滤波函数以及相位恢复算法对图像的重构质量做了进一步的优化。最后通过设计动态显微CT实验和搭建快速成像实验平台对研究内容做了实验验证。实验结果证明,本文实现了动态显微CT的快速重构,图像重构质量得到了一定程度的提高,有效缓解了上海光源X射线动态显微CT实验用户面临的原始实验数据大量积压的问题。
其他文献
"10以内数的认识和加减法"是数概念和计算教学中最基础的内容,是小学生学习数学的开始,也是后续知识学习的基础。作为教材的重要组成部分——习题,它设计的优劣会对学生的学习产生很大的影响,也需要我们不断加以改进和完善。为此,笔者作了一些研究和思考。
中国石油天然气股份有限公司辽河油田分公司紧紧抓住中国石油人事劳动分配三项制度改革和"油公司"模式改革的有利时机,持续创新政策机制,统筹优化生产经营要素,调整业务结构,优化组织机构,用好用活用新人力资源,大力实施人力资源创效工程,在"战严冬"中充分释放企业高质量发展动能。
目的了解北京市昌平区外来务工人员免疫预防相关传染病的健康素养状况和影响因素,为辖区制定相应的干预策略提供依据。方法采用多阶段整群随机抽样的方法抽取调查对象,对调查对象基本情况,免疫预防相关传染病健康素养水平等进行问卷调查。用Epidata 3.1软件建立数据库进行数据录入,采用SPSS 18.0软件进行数据分析,率的比较采用χ2检验,采用非条件Logistic回归模型对外来务工人员免疫预防相关传染
超导转变边缘探测器是硬X射线自由电子激光装置(Shanghai HIgh repetitio N rate XFEL and Extreme light facility,简称“SHINE”)中部署的探测器之一,是灵敏度极高的一类低温超导探测器,广泛应用于引力波探测、毫米波观测、X射线能谱测量,暗物质测量等领域。超导转变边缘探测器的核心是超导转变边缘传感器(Transition Edge Sens
脱氧核糖核酸(DNA)不仅能作为自然界中生物的遗传信息,近几十年来还作为一种新型的材料用于构建各种二维、三维的纳米结构。具有可编程性、可预测性、可寻址性等优点。DNA纳米材料近年的发展不仅仅局限于早年间结构的设计与合成,而是积极地向各个领域应用,与其他领域结合以另一种视角解决科学问题。DNA纳米技术与分子筛领域结合的相关工作较少,有少量将DNA纳米结构作为DNA分子筛的研究。DNA由于其良好的生物
全凭静脉麻醉是围手术期常见的麻醉方式,与吸入麻醉相比具有一定优势,在一些儿科手术中也更加适用。临床小儿麻醉的主要原则是提高麻醉安全性,合理控制麻醉深度。由于儿童的身体器官组织处于发育阶段,药动学与药效学在不同年龄阶段有较大差异,需要更精准地掌握麻醉用药,其中镇静、镇痛药物的选择和使用至关重要。除合理用药方案外,输注方式对提高小儿静脉麻醉安全性也不容忽视。对此,根据不同年龄段儿童的特点,不断优化静脉
相对论重离子碰撞可以产生高温、高密、解禁闭且手征对称性恢复的夸克胶子等离子体(QGP),同时也会产生迄今为止最强的电磁场,其中电磁场的量级可以达到1019 Gauss。非平庸的高温QCD真空会频繁地发生拓扑涨落和跃迁,因而会引起拓扑荷Qw的改变,并造成强相互作用QCD理论发生局域宇称(P)对称性和电荷共轭-宇称(CP)对称性破缺。解禁闭的QGP介质在QCD真空拓扑跃迁的条件下会产生手征不平衡和净电
熔盐堆作为第四代先进核反应堆的六种候选堆型之一,采用高温液态熔盐作为冷却剂、慢化剂和燃料,也是候选堆型中唯一使用液态燃料的反应堆系统。熔盐堆在高温常压下运行,其系统在固有安全性、中子经济性、防核扩散和燃料循环等方面具有显著的优势。当今研究现状和趋势表明,第四代反应堆已经逐渐成为解决世界新的能源和环境问题的关键,而熔盐堆具有较好的技术积累和显著的优点,逐渐在第四代堆中脱颖而出,成为当下研究焦点。放射
FLiBe熔盐作为钍基熔盐堆项目的燃料盐载体,高温下熔融会蒸发大量蒸气。熔盐蒸气冷却后凝聚形成熔盐颗粒物,随氩气吹扫进入尾气处理装置。这些含铍颗粒因其毒性需要准确监测、高效去除和达标排放,而粒径分布是解决这些问题的关键。本文主要探究了高温FLiBe和FLiBeZr熔盐颗粒的粒径分布特征。基于熔盐减压蒸馏实验装置和低压冲击采样器,搭建了高温熔盐颗粒粒径分布实验平台。在温度设置、样品采集、样品处理及测