学位论文 > 优秀研究生学位论文题录展示
磁共振稀疏重建算法与直接离散傅里叶变换重建算法并行处理加速方法的研究
作 者: 王聪
导 师: 陈武凡;冯衍秋
学 校: 南方医科大学
专 业: 计算机应用技术
关键词: 磁共振成像 非笛卡尔采样 稀疏采样 并行运算 多核处理器 通用图形处理器
分类号: TP391.41
类 型: 硕士论文
年 份: 2011年
下 载: 66次
引 用: 0次
阅 读: 论文下载
内容摘要
磁共振成像具有无电离辐射、多参数成像、功能成像、可任意方向断层成像等优点。在临床上磁共振已经得到广泛的应用,成为最重要的成像方法之一。具有诸多优点的同时,扫描速度慢是磁共振成像的一大瓶颈,不但给病人造成长时间等待的痛苦,而且使磁共振难以应用到动态成像等需要快速成像的方面。自磁共振成像出现以来,提高成像速度一直是一个研究课题,非笛卡尔采样(如螺旋轨迹、放射状轨迹、propeller等)已经得到了较好的应用,新近出现的稀疏磁共振成像方法也正在作为一个热点被研究。对于非笛卡尔采样的到的数据一般采用各种插值方法重采样到均匀的笛卡尔坐标点上,然后用经典的fft方法重建图像。插值方法已经比较成熟,但不插值而直接进行离散傅里叶变换的共轭相位重建算法仍被认为是最准确的方法。共轭相位算法最大的问题是时间复杂度太高,因此一般不用于临床而只在科研中用来生成对照图像以评价各种插值算法的准确性。新近出现的稀疏磁共振成像的重建算法要解决的问题是从病态程度很高的欠定方程组的无穷解中找到一个具有稀疏性或变换稀疏性的解,这是一个是迭代算法,若采用非线性共轭梯度下降算法,迭代过程中要进行多次空间变换、梯度等复杂的矩阵操作,因此也是一个非常耗时的算法。解决密集运算耗时多的方法是并行运算,即把一项大的任务分解为多个可同时进行的子任务,然后将子任务分发到位于同一计算机或不同计算机的多个处理器上同时执行,子任务执行完毕再将各子任务的结果综合为宿主任务的最终结果。经典的并行计算使用的是多CPU计算机或者计算机集群,这些超级计算机的占地与功耗都非常大,应用到临床成本很大;具有多核CPU的桌面计算机虽然有一定的提速效果但不够明显,新近出现的GPGPU则将网格计算所需的硬件集成到一块板卡上,很大程度上降低了体积和功耗。从另一个角度来看,GPGPU是对传统GPU的改进,传统GPU只能按照固定的流水线进行图形处理,虽然流水线的每一级都具有高性能的并行运算能力,但使用很不方便,GPGPU则将GPU改进为可用C/C++或Fortran语言编程的处理器,可以比较方便的应用到各种密集计算任务中。本文首先将原理上比较简单的共轭相位算法进行两种粒度的任务分解,然后分别利用较容易理解的共享式存储多核桌面计算机和稍复杂的GPGPU对其进行加速,与单核串行运算相比分别得到了5倍和72倍的加速比。稀疏磁共振重建算法过程比较复杂,因此本文将其子过程分别做并行分解,子过程由GPU来处理,对整体迭代过程的控制计算量很小,作为串行程序由CPU来执行。加速效果比较明显与Michael Lustig提供的SparseMRI0.2相比有76倍的加速。
|
全文目录
摘要 3-5 ABSTRACT 5-10 第一章 引言 10-14 1.1 磁共振成像技术简介 10-11 1.2 研究背景 11-12 1.3 本文的主要工作 12-13 1.4 本文的结构 13-14 第二章 磁共振成像的基本原理 14-27 2.1 核磁共振现象 14-17 2.1.1 自旋核及其磁矩 14-15 2.1.2 静磁场对自旋核的作用 15-17 2.1.3 核磁共振现象及其产生条件 17 2.2 磁共振信号的产生和检测 17-21 2.2.1 弛豫 18-19 2.2.2 自由进动 19-20 2.2.3 磁共振信号的检测 20-21 2.3 层面选择和空间信息编码技术 21-25 2.3.1 层面选择技术 21-23 2.3.2 空间信息编码技术 23-25 2.4 K空间 25-27 2.4.1 K空间定义 25-26 2.4.2 K空间数据的性质 26-27 第三章 并行处理 27-42 3.1 什么是并行处理 27-28 3.2 概念和术语 28-30 3.2.1 冯诺依曼结构 28 3.2.2 弗林分类 28-30 3.3 并行计算机的存储结构 30-33 3.3.1 共享式存储 30-31 3.3.2 分布式存储 31-32 3.3.3 分布-共享混合式存储 32-33 3.4 并行编程模型 33-36 3.5 OpenMP 36-37 3.6 NVIDIA GPU与CUDA 37-42 3.6.1 GPU的历史 37-39 3.6.2 支持CUDA的GPU 39-40 3.6.3 CUDA架构下的编程模型 40-42 第四章 共轭相位重建算法的并行化 42-49 4.1 什么是共轭相位重建算法 42-43 4.2 算法的时间复杂度和并行化分解 43-44 4.3 DFT在多核CPU上的实现 44-45 4.4 DFT在GPGPU上的实现 45-46 4.5 实验结果 46-49 第五章 稀疏磁共振重建及其GPGPU实现 49-63 5.1 什么是稀疏磁共振成像 49-51 5.2 稀疏磁共振成像的重建方法 51-52 5.3 重建方法的GPGPU实现 52-63 5.3.1 需要的数据结构 52-56 5.3.2 二维小波变换的实现 56-58 5.3.3 其它子过程的实现 58-60 5.3.4 实验结果与讨论 60-63 第六章 总结与展望 63-65 参考文献 65-70 致谢 70-72 统计合格证明 72
|
相似论文
- MR导向冷冻消融兔坐骨神经的影像与病理对照实验研究,R445.2
- 背景抑制磁共振全身弥散加权成像(DWIBS)与hPET/CT临床价值对比研究,R445.2
- 3.0T高场磁共振扩散张量成像对特发性癫痫的研究,R445.2
- MRS联合DWI对前列腺癌的诊断价值,R737.25
- 颅脑外伤及上呼吸道感染后嗅球MR体积测量与嗅觉功能相关性研究,R651.1
- MRA在颅内动脉瘤栓塞术后随访应用的可行性分析,R651.1
- 健康国人嗅觉相关中枢结构MRI体积与年龄及嗅觉功能相关性研究,R445.2
- 3.0T BOLD MR技术对大鼠急性单侧输尿管梗阻肾损害的研究,R693.2
- 进行性核上性麻痹的临床及磁共振影像学特征的研究,R741
- 3.0T高场磁共振对中央型肺癌的初步研究,R734.2
- MR弥散加权成像在股骨头坏死的应用价值研究,R681.8
- 磁共振弥散加权及动态增强成像量化分析在乳腺病变诊断中的应用研究,R445.2
- MRI在多系统萎缩、进行性核上性麻痹和原性帕金森病诊断及鉴别诊断中的应用研究,R445.2
- STIR-FLAIR与T2W-SPIR序列在眼眶病变诊断中的对比研究,R777
- 功能磁共振成像技术在前列腺癌诊断中的应用研究,R737.25
- 颞下颌关节弹响治疗前后的磁共振对比研究,R445.2
- TD-SCDMA无线链路控制协议实现研究,TN929.533
- 两种新型磁共振造影剂的合成和结构表征,R914
- 3D-MRI诊断胎儿体表畸形的研究及和2D-MRI,3D-US的对比,R714.5
- 3.0T磁共振灌注加权成像与扩散张量成像在脑胶质瘤分级诊断中的应用价值,R445.2
- MRI在肱骨外上髁炎中的探索和应用,R681.8
中图分类: > 工业技术 > 自动化技术、计算机技术 > 计算技术、计算机技术 > 计算机的应用 > 信息处理(信息加工) > 模式识别与装置 > 图像识别及其装置
© 2012 www.xueweilunwen.com
|