基于异构平台的CT图像三维重建加速方法及其装置制造方法及图纸
【技术实现步骤摘要】
基于异构平台的CT图像三维重建加速方法及其装置
本专利技术属于X射线CT
，特别涉及一种基于异构平台的CT图像三维重建加速方法及其装置。
技术介绍
X射线计算机断层成像(ComputedTomography，CT)是一种由物体的X射线投影来反求其衰减分布的技术，涵盖核物理、数学、计算机、精密仪器等多个学科。由于CT可以在非接触、不破坏的条件下得到物体内部的高精度三维结构信息，因此自Hounsfield成功研制第一台CT以来，CT已在无损检测、医学诊断、材料分析等领域得到广泛应用。在实际应用中，锥束CT高分辨率三维重建所需的计算资源和存储资源都非常大，随着重建规模的增加，重建的存储需求和计算量迅速增加，在很多情况下，很难满足实际应用的需求，以重建算法中通用的反投影计算为例：设待重建三维图像的各维尺寸都为N，则相应反投影的计算复杂度将高达O(N4)，重建一个分辨率为10243的三维图像，计算将循环10995亿次，在普通PC机上完成如此大的计算量是非常耗时的，难以满足实际应用的要求。因此，锥束CT重建过程加速是当前工程人员亟需解决的问题，针对CT重建算法设计相应的...
【技术保护点】
一种基于异构平台的CT图像三维重建加速方法，其特征在于，异构平台包含主机及异构OpenCL计算设备，该加速方法包含如下内容：对FDK重建算法进行算粒分解，分析各个算粒的并行计算流程；通过异构平台中的主机及异构OpenCL计算设备对各个算粒进行加速优化处理。
【技术特征摘要】
1.一种基于异构平台的CT图像三维重建加速方法，其特征在于，异构平台包含主机及异构OpenCL计算设备，该加速方法包含如下内容：对FDK重建算法进行算粒分解，分析各个算粒的并行计算流程；通过异构平台中的主机及异构OpenCL计算设备对各个算粒进行加速优化处理。2.根据权利要求1所述的基于异构平台的CT图像三维重建加速方法，其特征在于，所述的主机为运行主程序的CPU，OpenCL计算设备包含运行内核程序的异构容器GPU和FPGA，CPU、GPU及FPGA之间通过PCI-E总线相通信，主程序通过定义上下文来管理内核程序的运行。3.根据权利要求2所述的基于异构平台的CT图像三维重建加速方法，其特征在于，对FDK重建算法进行算粒分解，包含：根据FDK算法内容，分解为：用于对投影数据进行加权的投影加权算粒，用于对加权后的投影数据进行滤波的滤波算粒，用于将滤波后的投影数据反投影到重建物体上的反投影算粒，及用于对反投影结果进行归约处理的归约算粒。4.根据权利要求3所述的基于异构平台的CT图像三维重建加速方法，其特征在于，根据FDK重建算法公式：通过对公式中积分的拆分及离散化，分为：投影加权算粒，表示为：其中，p'(θ,u,v)表示旋转分度为θ时对投影数据进行加权后的数据，为加权系数；滤波算粒，表示为：其中，df(θ,u,v)为滤波后的数据，h(u)为滤波算子的单位冲激响应，[-um,um]表示探测器每一行采集的2m个数据；反投影算粒，表示为：其中，f(x,y,z,θ)表示重建物体在旋转分度为θ时投影点对f(x,y,z)的贡献值；归约算粒，表示为：其中，φmax为重建物体旋转一周时离散采得的投影分度数。5.根据权利要求3所述的基于异构平台的CT图像三维重建加速方法，其特征在于，对各个算粒进行加速优化处理，包含：采用FPGA对投影加权算粒进行并行处理，通过异步传输至GPU，在传输过程中同时对滤波算粒进行处理；结合反投影时各体素点的数据并行操作，在GPU内按体素点对反投影算粒进行多线程并行反投影计算。6.根据权利要求5所述的基于异构平台的CT图像三维重建加速方法，其特征在于，根据FDK重建算法中重建区域在旋转轴方向上行的各层投影数据和探测器投影数据纵轴方向上的各行数据的投影对应关系，采用分块重建策略，沿旋转轴方向将待重建区域分成若干块，重建其中一块时从外部存储器取出对应的投影数据进行重建操作。7.根据权利要求5所述的基于异构平台的CT图像三维重建加速方法，其特征在于，采用FPGA对投影加权算粒进行并行处理，包含：将全局存储器分割为2块bank，通过负载分配实现随机存取存储器的访问均衡；通过常数存储器存储需要多次重复计算的中间变量。8.根据权利要求5所述的基于异构平台的CT图像三维重建加速方法，其特征在于，在GPU内按体素点对反投影算粒进行多线程并行反投影计算，包含：采用基于体素驱动方式，按...
【专利技术属性】
技术研发人员：闫镔，李磊，王林元，孙艳敏，路万里，蔡爱龙，张瀚铭，张文昆，
申请(专利权)人：中国人民解放军信息工程大学，
类型：发明
国别省市：河南,41
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。