基于GPU的三维医学图像处理算法研究

论文摘要

医学图像处理与分析技术不仅极大地提高了现代医学临床诊断水平,而且为医学培训、教学与研究、手术导航等提供了数字实现手段,成为推动现代医学发展的重要力量。因此,关于医学图像处理与分析技术的研究一直是国内外广大学者所关注的焦点,这对保证人们身体健康和提高人们身体素质具有重大意义。根据当前实际应用的需求,采用直观的包含丰富信息的3D数据代替简单的2D图像已成为医学图像处理领域的发展趋势。但是,在目前的个人电脑中进行三维数据处理仍是一个非常耗时的计算过程,不能满足实际应用中的实时性要求,这主要归因于三维图像包含数据量较大,而且各种医学图像处理算法相对复杂、计算量较大。因此,如何进一步提高三维医学图像处理算法的计算效率,使其满足实时应用的需求,是当前该领域内一个亟待解决的问题。近年来,计算机图形处理器（Graphics Processing Unit,GPU）得到了极大地发展,从最初局限于图形渲染的图形卡,发展为如今可编程的并行计算平台。针对GPU的高级编程语言如Cg、HLSL等使得人们更加容易对其进行编程控制,并应用在渲染之外的通用计算方法,如偏微分方程求解、矩阵运算等。与CPU的串行计算模式不同,GPU是一种高度并行的流处理器,具有更强的浮点计算能力,在物理模拟、信号分析等许多领域人们将高强度的计算任务以适当的方式转化为流数据计算模式,通过编程在GPU中进行加速计算,通常能获得一个数量级的速度提高,这也是当今的热点研究问题之一。因此,通过对GPU编程发挥其高强度计算能力来解决三维医学图像处理的速度瓶颈,提高算法的执行效率,将是一个十分有应用价值的研究内容。本文首先研究了可编程图形处理器的体系结构及其流式编程模型。GPU的体系结构以图形流水线为基础,提供了顶点处理器和片段处理器的可编程特性,通过OoenGL等编程接口与主程序进行通信。GPU将用户编写的程序作为运算核对数据流中的多个元素同时进行处理,实现高度的并行计算。但根据流计算的特点,GPU程序对存储器的访问将受到一定的限制。在此基础上,本文研究了高效的三维医学图像处理算法,主要包括基于3D纹理的精确体绘制方法、CPU-GPU联合工作模式加速基于互信息的3D图像刚性配准、GPU加速的圆弧轨迹Katsevich锥束CT重建算法,具体内容为:（1）基于3D纹理的体绘制方法是一种高效的三维医学体数据可视化手段,但由于受低采样率的影响,在最终绘制图像中往往会出现切片痕迹,造成图像质量较差。分段积分是提高绘制图像质量的有效方法,它首先对每个采样段应用体绘制积分方程,然后沿视线方向将所有采用段进行累加,得到较好质量的绘制图像。采样段内的颜色值通常采用线性表达式代入到体绘制积分方程中,但由于人体各组织之间关联复杂,线性插值不能很好地近似采样段内的颜色变化,本文通过计算每个采样段中点采样位置得到更加精确的抛物线近似表达,并进一步求解体绘制积分方程,将其化简为易于在GPU中实现的形式。通过顶点处理器计算采样段三个控制点位置的纹理坐标,利用片段处理器计算采样段上的积分方程。实验结果表明该方法在保证实时绘制速度的同时取得了更好的绘制图像质量。在此基础上,本文进一步改进了基于剪切几何体的剖切绘制方法,降低内存使用,提高剖切绘制的速度。（2）针对基于互信息的3D医学图像刚性配准,提出一种基于CPU-GPU联合工作模式的加速方法。最外层的搜索算法需要复杂的逻辑控制,但计算强度小,适合CPU计算。在搜索过程中调用互信息计算函数,首先将浮动图像进行空间变换和三线性插值,然后统计与参考图像之间的联合直方图,并根据联合直方图计算熵和互信息值。空间变换和三线性插值约消耗互信息计算函数的85%时间,而且适合并行计算,因此将其移植至GPU中执行。采用平面三维纹理、渲染到纹理等技术编写顶点程序和片段程序,在GPU中执行三线性插值,并将结果读入系统内存,继续由CPU完成互信息值计算。实验结果表明,与传统方法相比,该算法在获得相似配准精度的情况下,将配准速度提高了一个数量级。（3）改进滤波反投影方法在GPU中的实现技术,并应用于圆弧轨迹Katsevich精确锥束重建算法中。Katsevich算法首先对投影数据进行滤波处理,包括偏导数计算、傅立叶变换、Hilbert变换等,这个过程由CPU完成,并将滤波后的投影数据传递到GPU内存,进行反投影计算。本文利用帧缓冲对象（FBO）的多个绑定点,将上一次的输出作为下一次渲染的输入,在GPU中完成所有射源转角下的反投影数据累加;采用浮点纹理映射保证重建图像的质量;设置RGBA四通道纹理实现数据级的并行计算;采用四叉树编码的矩形网格加速绘制过程。本文GPU加速的重建方法所得图像质量已非常接近于传统CPU重建方法,而重建速度却提高了10倍以上。最后,本文总结了三维医学图像处理算法在GPU平台上运行时的编程特点,主要包括GPU的存储器模型、三维数据结构、并行策略和程序优化等。

论文目录

• 摘要
• Abstract
• 1 绪论
• 1.1 研究背景及意义
• 1.2 GPU在通用计算领域中的应用
• 1.3 本文研究内容及主要贡献
• 1.4 论文结构安排
• 2 GPU通用计算性能及其开发环境
• 2.1 引言
• 2.2 可编程图形处理器的发展简介
• 2.2.1 计算机图形硬件的发展
• 2.2.2 图形硬件流水线
• 2.2.3 可编程图形管线
• 2.3 GPU通用计算模型
• 2.4 GPU通用计算特点及发展趋势
• 2.4.1 GPU通用计算的特点
• 2.4.2 GPU通用计算的发展趋势
• 2.5 可编程图形处理器的开发环境
• 3 基于三维纹理的精确体绘制积分
• 3.1 引言
• 3.2 体绘制方法的主要步骤
• 3.2.1 数据模型
• 3.2.2 滤波与分类
• 3.2.3 光照计算
• 3.2.4 颜色合成
• 3.3 精确分段体绘制积分方程
• 3.3.1 传统体绘制积分方程
• 3.3.2 抛物线插值的分段积分
• 3.4 三维纹理加速的体绘制
• 3.4.1 设置数据纹理
• 3.4.2 生成代理几何体
• 3.4.3 渲染代理几何体
• 3.4.4 硬件颜色合成
• 3.5 剖切绘制
• 3.5.1 基于剪切几何体的体绘制剖切方法
• 3.5.2 结合传递函数的剖切方法
• 3.6 实验结果
• 3.6.1 实验数据
• 3.6.2 三维数据体绘制结果
• 3.7 小结
• 4 GPU加速基于互信息的三维医学图像刚性配准
• 4.1 引言
• 4.2 医学图像配准算法概述
• 4.2.1 图像配准的基本概念
• 4.2.2 图像配准的分类
• 4.2.3 配准算法的基本要素
• 4.2.4 图像的离散化特性
• 4.3 基于最大互信息值的图像配准
• 4.3.1 互信息值作为相似性测度
• 4.3.2 优化策略
• 4.3.3 加速方法
• 4.4 基于GPU的快速配准算法
• 4.4.1 GPU在图像配准中的应用
• 4.4.2 基于CPU-GPU联合平台的配准算法结构
• 4.4.3 设置平面三维纹理
• 4.4.4 执行GPU程序
• 4.4.5 渲染到纹理与数据传递
• 4.4.6 基于CPU的运算过程
• 4.5 实验结果
• 4.5.1 GPU实现体数据3D空间变换
• 4.5.2 GPU加速配准过程
• 4.6 小结
• 5 GPU加速的圆弧轨迹Katsevich锥束CT重建
• 5.1 引言
• 5.2 圆弧轨迹Katsevich锥束重建
• 5.2.1 CT成像的基本原理
• 5.2.2 滤波反投影基础
• 5.2.3 圆弧轨迹Katsevich锥束重建算法
• 5.2.4 算法实现结构与CPU滤波
• 5.3 基于透视纹理映射的反投影计算
• 5.3.1 GPU在锥束CT重建算法中的应用
• 5.3.2 基于纹理映射的反投影方法
• 5.4 GPU加速锥束重建算法
• 5.4.1 设置RGBA四通道纹理
• 5.4.2 细化矩形网格与纹理映射
• 5.4.3 FBO渲染到纹理与多次渲染传递
• 5.4.4 顶点程序与片段程序
• 5.5 实验结果与分析
• 5.6 小结
• 6 三维医学图像处理中GPU程序优化方法
• 6.1 GPU存储器模型
• 6.2 GPU三维数据结构
• 6.2.1 3D纹理
• 6.2.2 2D纹理堆栈
• 6.2.3 平面3D纹理
• 6.3 数据并行计算
• 6.3.1 指令级并行
• 6.3.2 数据级并行
• 6.4 计算性能
• 6.5 小结
• 7 总结与展望
• 7.1 全文工作总结
• 7.2 工作展望
• 参考文献
• 攻读博士学位期间发表学术论文情况
• 致谢
• 创新点摘要