一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法

doi:10.11938/cjmr20160404

波谱学杂志 ›› 2016, Vol. 33 ›› Issue (4): 549-558.doi: 10.11938/cjmr20160404

一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法

庄孝星, 马崚嶒, 蔡聪波, 陈忠

厦门大学电子科学系, 福建省等离子体与磁共振研究重点实验室, 福建厦门 361005

收稿日期:2016-04-08 修回日期:2016-10-24 出版日期:2016-12-05 发布日期:2016-12-05
通讯作者: 陈忠,电话:0592-2181712,E-mail:chenz@xmu.edu.cn E-mail:chenz@xmu.edu.cn
作者简介:庄孝星(1990-),男,福建连江人,硕士研究生,主要从事核磁共振成像算法研究
基金资助:
国家自然科学基金资助项目（11375147，81171331）.

An Image Reconstruction Algorithm for Spiral MRI Based on Spatio-Temporal Transform and Compressed Sensing

ZHUANG Xiao-xing, MA Ling-ceng, CAI Cong-bo, CHEN Zhong

Department of Electronic Science, Fujian Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen 361005, China

Received:2016-04-08 Revised:2016-10-24 Online:2016-12-05 Published:2016-12-05

摘要/Abstract

摘要： 螺旋采样磁共振快速成像在功能性成像、并行成像和动态成像等领域发挥着越来越重要的作用．螺旋采样图像重建的传统算法是用核函数将螺旋状分布的k空间数据插值到均匀网格中，再利用傅里叶变换和最小二乘法进行重建．但是基于网格化的算法对核函数过于依赖，在网格化过程中产生难以避免的误差．该文提出了基于时空变换和压缩感知的l₁范数的最优化模型和重建算法．时空变换矩阵描述了空间上的磁共振图像与采集到的时域信号间的关系，使得算法直接使用采集到的数据作为保真约束项，避免了网格化过程产生的误差．此外，基于图像处理单元的并行计算被用来提高时空变换矩阵的运算速度，使得算法具有较强的应用价值．

关键词: 时空变换, 螺旋采样磁共振成像, 压缩感知, 并行计算, 非均匀傅里叶变换

Abstract: Spiral magnetic resonance imaging (MRI) plays an important role in functional imaging, parallel imaging and dynamic imaging. Most of the traditional image reconstruction algorithms for spiral imaging are based on kernel functions which interpolate the k-space data acquired with spiral trajectories onto a uniform grid in Cartesian space. Non-uniform fast Fourier transform (NUFFT) and least square method could then be applied after gridding to reconstruct the images. With these algorithms, the results are dependent on the choice of kernel function, and reconstruction errors are unavoidable during the gridding process. In this study, a spatio-temporal transform (STT) matrix, representing the relation between the image and sampled k-space data, is introduced in l₁ norm to optimize the problem based on spatio-temporal transform and compressed sensing (CS). The k-space data, rather than the after-gridding data, are used as the fidelity term, such that the gridding errors can be avoided. In addition, parallel computing on GPU can be applied to reduce the computation time for the STT matrix, making the algorithm more efficient.

Key words: spatio-temporal transformation(STT), compressed sensing(CS), parallel computing, spiral MRI, Non-uniform fast Fourier transform(NUFFT)

中图分类号:

O482.53

庄孝星, 马崚嶒, 蔡聪波, 陈忠. 一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法[J]. 波谱学杂志, 2016, 33(4): 549-558.

ZHUANG Xiao-xing, MA Ling-ceng, CAI Cong-bo, CHEN Zhong. An Image Reconstruction Algorithm for Spiral MRI Based on Spatio-Temporal Transform and Compressed Sensing[J]. Chinese Journal of Magnetic Resonance, 2016, 33(4): 549-558.

参考文献

[1] Yang Bao-lian(杨保联). Future of ultra high field MRI in basic research and clinical applications(超高场磁共振人体成像应用研究和医学前景)[J]. Chinese J Magn Reson(波谱学杂志), 2015, 32(4):707-714.
[2] Haacke E M, Brown R W, Thompson M R, et al. Magnetic Resonance Imaging:Physical Principles and Sequence Design[M]. New York:Wiley-Liss, 1999.
[3] Sachs T S, Meyer C H, Hu B S, et al. Real-time motion detection in spiral MRI using navigators[J]. Magn Reson Med, 1994, 32(5):639-645.
[4] Meyer C H, Hu B S, Nishimura D G, et al. Fast spiral coronary artery imaging[J]. Magn Reson Med, 1992, 28(2):202-213.
[5] Tolouee A, Alirezaie J, Babyn P. Compressed sensing reconstruction of cardiac cine MRI using golden angle spiral trajectories[J]. J Magn Reson Imaging, 2015, 260:10-19.
[6] O'Sullivan J D. A fast sinc function gridding algorithm for Fourier inversion in computer tomography[J]. IEEE T Med Imaging, 1985, 4(4):200-207.
[7] Jackson J, Meyer C H, Nishimura D G, et al. Selection of a convolution function for Fourier inversion using gridding[J]. IEEE T Med Imaging, 1991, 10(3):473-478.
[8] Sha L, Guo H, Song A W. An improved gridding method for spiral MRI using nonuniform fast Fourier transform[J]. J Magn Reson Imaging, 2003, 162(2):250-258.
[9] Fessler J A. On NUFFT-based gridding for non-Cartesian MRI[J]. J Magn Reson Imaging, 2007, 188(2):191-195.
[10] Candès E J, Wakin M B. An introduction to compressive sampling[J]. IEEE Signal Proc Mag, 2008, 25(2):21-30.
[11] Donoho D L. Compressed sensing[J]. IEEE T Inform Theory, 2006, 52(4):1289-1306.
[12] Li Zhi-min(李智敏), Xie Hai-bin(谢海滨), Zhou Min-xiong(周敏雄), et al. Spike noise removal for magnetic resonance imaging based on sparse reconstruction(基于稀疏重建的磁共振图像尖峰噪声消除方法)[J]. Chinese J Magn Reson(波谱学杂志), 2015, 32(01):41-50.
[13] Lustig M, Donoho D, Pauly J M. Sparse MRI:The application of compressed sensing for rapid MR imaging[J]. Magn Reson in Med, 2007, 58(6):1182-1195.
[14] Owens J D, Houston M, Luebke D, et al. GPU computing[J]. P IEEE, 2008, 96(5):879-899.
[15] Sanders J, Kandrot E. CUDA by Example:An Introduction to General-purpose GPU Programming[M]. Boston:Addison-Wesley, 2010.
[16] Yang J, Zhang Y, Yin W. A fast alternating direction method for TVL1-L2 signal reconstruction from partial fourier data[J]. IEEE J STSP, 2010, 4(2):288-297.
[17] Zhan Z, Cai J, Guo D, et al. Fast multi-class dictionaries learning with geometrical directions in MRI reconstruction[J]. IEEE T Bio-Med Eng, DOI:10.1109/TBME.2015.2503756.

一种基于时空变换和压缩感知的磁共振螺旋采样的图像重建算法

An Image Reconstruction Algorithm for Spiral MRI Based on Spatio-Temporal Transform and Compressed Sensing

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 13

编辑推荐

Metrics

本文评价

[1]	柴青焕, 苏冠群, 聂生东. 双树小波变换与小波树稀疏联合的低场CS-MRI算法[J]. 波谱学杂志, 2018, 35(4): 486-497.
[2]	徐嘉文, 徐健, 周晓东, 张聪, 陈群. 基于Gadgetron平台的多GPU分布式磁共振图像重建[J]. 波谱学杂志, 2018, 35(3): 303-317.
[3]	黄丽洁, 宋阳, 赵献策, 谢海滨, 吴东梅, 杨光. 一种结合并行成像和压缩感知的快速磁共振成像新方法[J]. 波谱学杂志, 2018, 35(1): 31-39.
[4]	罗庆, 张成秀, 李文静, 郑慧, 谢海滨, 杨光. 一种新的磁共振图像处理流水线的设计与实现[J]. 波谱学杂志, 2018, 35(1): 40-51.
[5]	宋阳, 谢海滨, 杨光. 用于压缩感知磁共振成像的分割字典学习算法[J]. 波谱学杂志, 2016, 33(4): 559-569.
[6]	聂莉莎, 蒋滨, 张许, 刘买利. 基于压缩感知/重采样的NMR噪声抑制新方法[J]. 波谱学杂志, 2016, 33(2): 244-256.
[7]	高芒, 谢海滨, 李智敏, 张成秀, 奚伟, 姜小平, 杨光. 压缩感知同步扫描重建及其采样方案的研究[J]. 波谱学杂志, 2016, 33(2): 257-268.
[8]	韩明月，郑慧，胡炳文，杨光. 迭代软阈值法压缩感知重建谱峰较宽的二维固体谱 [J]. 波谱学杂志, 2015, 32(4): 551-562.
[9]	肖洒1,2，吕植成1，孙献平1，叶朝辉1，周欣1. CS-MRI* 中评价随机欠采样矩阵的新方法 [J]. 波谱学杂志, 2015, 32(4): 584-595.
[10]	董芳，王前锋，郑慧，李智敏，杨光，李建奇*. 结合导航回波与压缩感知进行运动伪影矫正[J]. 波谱学杂志, 2015, 32(3): 419-429.
[11]	李智敏1，谢海滨1，周敏雄1，张成秀2，奚伟2，姜小平2，杨光1. 基于稀疏重建的磁共振图像尖峰噪声消除方法[J]. 波谱学杂志, 2015, 32(1): 40-50.
[12]	郑慧，韩明月，胡炳文，杨光. 不同采样模式的固体DQ-SQ 实验的压缩感知重建比较[J]. 波谱学杂志, 2014, 31(4): 535-547.
[13]	高明生，谢海滨，严序，周敏雄，李智敏，杨光*. 选择性双向顺序压缩感知重建动态磁共振成像[J]. 波谱学杂志, 2013, 30(2): 194-203.