磁共振成像仿真平台设计和实现.docVIP

磁共振成像仿真平台设计和实现 　　作者简介作者简介：水力（1987-），男，中南民族大学生物医学工程学院硕士研究生，研究方向为磁共振成像。0引言 　　磁共振成像（magnetic resonance imaging，MRI）是一种利用核磁共振原理进行人体断层成像的技术［1］。近年来，随着磁共振成像技术的不断进步，其成像设备的研发得到了飞速发展，对MRI系统成像方法以及临床应用等方面的研究工作也日益增多。但磁共振成像仪是一种昂贵的医疗仪器，研究人员常常需要用仿真方法来检验他们的研究结果，此时，磁共振成像仿真平台的建立尤为必要。MRI仿真平台可方便研究人员比较和评价数据采集和图像重建的新算法。在研究中，一般先借助仿真模型，通过人工合成数据验证算法的可行性和准确性。模型数据测试可行后，才用MRI扫描仪得到的试验数据进行验证。利用仿真平台进行图像重建算法研究的优点是: 模型是已知的，重建结果有可比较的参考对象，而且可以按照成像原理获得不同采集模式的数据，算法研究前期可节省用扫描仪采集数据的花费。 　　本文利用Matlab图形用户界面仿真MRI成像过程，模拟K空间数据采集和图像重建。该平台模拟了线性、螺旋及并行扫描3种成像方式，可以导入不同模型图像，设计K空间扫描轨迹，模拟K空间数据采集过程，并显示采集信号的波形。还对3种成像方式的图像重建方法进行了实现，对重建误差进行了分析。界面简明友好，为后续MRI成像算法的研究提供了一种有效工具。 　　1软件系统设计 　　本文介绍了Matlab图形用户界面开发环境GUIDE［2］，实现了磁共振成像仿真模拟平台的结构和功能。该软件可以实现MRI模型图像导入、模型切换、K空间轨迹设计［3］、K空间数据采集、图像重建以及误差分析等功能。是学习以及进一步研究磁共振成像方法的有力平台工具。磁共振成像仿真平台软件由图像导入系统、K空间轨迹设计和数据采集系统、图像重建系统三大功能模块组成，如图1所示。 　　图1软件系统模块 　　2软件系统功能 　　 采用Core2-3.0G CPU，1.99G内存和120G硬盘的PC机在Matlab2007上实现，软件界面如图2所示。 　　以下分别对软件各个模块的功能进行介绍。 　　2.1图像导入系统 　　通过选项框有5个图像模型可选，如图3所示，分别为不同分辨率的Shepp_logan模型和分辨率为128×128和256×256的头部模型，点击显示图像按键后导入图像到指定轴。 　　2.2K空间数据采集系统 　　K空间数据采集系统分为K空间轨迹设计和模拟数据采集过程的信号波形显示两个部分。 　　目前,本仿真平台研究了3种采样轨迹的设计和数据采集：①等间隔采样轨迹，即线性采集；②均匀螺旋轨迹设计和采集；③并行轨迹设计和采集。根据扫描轨迹设置相关参数设计出K空间填充轨迹后，再利用模拟数据采集的过程，得到图像重建需要的K空间数据，即采集到的信号。其中线性、并行扫描显示的是一行行采集波形，而螺旋扫描显示的是一支支采集波形。 　　图2软件界面 　　图3导入模型 　　在数据采集波形显示部分，设置了多种扫描方式：单行（支）数据采集、连续动态采集，以及轨迹扫描和波形显示的动态同显。本界面还设计了良好的排错系统，使得界面更加人性化。 　　2.2.1线性采集 　　 等间隔采样又称笛卡尔网格采样，它在K空间中是等间隔均匀采样，即基于线性填充轨迹的数据采样。在数据采集期间，所得数据沿Kx方向依次存放，Ky不变，采集的数据从左至右地顺序写入K平面的第一行，本行数据写完后，Ky加一再写第二行，直至写完。标准的2D K空间的标记轨迹为K空间内的一组平行线，即线性轨迹。如图4所示。 　　 其频域信号为：S(kx,ky)=??Mxy(x,y,t)× 　　exp(－j(kx(t)x+ky(t)y))dxdy（1）2.2.2螺旋采集 　　螺旋轨迹的形式很多，本文采用的是均匀螺旋K轨迹（dK/dθ为常数）阿基米德螺线（Archimedean spiral）。在笛卡尔坐标系中，阿基米德螺线K轨迹可表示为： 　　图4等间隔K-空间采样轨迹Kx=Kcosθ=Aθ(t)cosθ 　　Ky=Ksinθ=Aθ(t)sinθ （2）螺旋轨迹一般从K空间中央出发，然后呈螺旋状向外扩展，轨迹如图5所示。其中，（a）是采集单只螺旋轨迹，（b）是采集20支螺旋轨迹。 　　图5均匀螺旋采样轨迹 　　2.2.3并行采集 　　并行磁共振成像的数据用相控阵线圈，以加速因子（acceleration factor）R 进行采集，也就是每隔R个相位编码步做一次相位编码步数据采集，如图6所示。 　　(a)(b) 　　图6线性和并行采集的K空间 　　(a)为正常采集的K空间数据示意图，每一条横