基于NLTV和NESTA的MRI/MRA图像感兴趣区域的压缩感知重构

doi:10.11896/j.issn.1002-137X.2017.09.058

摘要/Abstract

摘要： 压缩感知理论借助信号内在的稀疏性或可压缩性,利用随机投影实现在远低于奈奎斯特频率的采样频率下对压缩数据进行采集。将该技术应用于医学成像领域可以加快MRI/MRA的扫描速度,提高扫描效率,减少患者的不适感。以NLTV(Nonlocal Total Variation)正则化来改善传统TV导致的边缘模糊、阶梯效应等缺点,提出改进的NESTA算法(简称NLTV-ROI-NESTA算法)实现MRI/MRA图像感兴趣区域(Region of Interests,ROIs)的精确重构,增强低对比度血管的细节信息,以峰值信噪比、结构化相似度、相对误差3个指标来定性、定量地评价算法的性能。实验结果表明,与传统的压缩感知重构算法相比,NLTV-ROI-NESTA算法在重构精度和细节保留方面均具有明显优势,能较好地保持低对比度血管或其他感兴趣区域的细节特征,在快速医学成像领域具有广阔的应用前景。

关键词: 压缩感知,图像重构,非局部TV,NESTA算法,感兴趣区域

Abstract: Compressed sensing (CS) is a novel developed theoretical framework for information acquisition and proces-sing.Taking advantages of the inherent sparsity or compressibility in real world signals,CS can collect compressed data at the sampling rate much lower than that needed in Shannon’s theorem.CS is used to medical imaging techniques to accelerate the scanning speed of MRI/MRA,improve the scanning efficiency and alleviate the patients’ suffering.The flow chart of our NLTV-ROI-NESTA algorithm is as follows.The nonlocal total variation (NLTV) is applied to overcome the disadvantages of the traditional TV for its edge blurring and stepstairs effects.An improved NESTA algorithm is used to reconstruct the region of interests (ROIs) for MRI/MRA images accurately and fast to maintain or enhance the details of the low-contrast vessels.Three indices:peak signal to noise rate (PSNR),structural similarity index (SSIM) and relative l2 norm error (RLNE) are adopted to compare the reconstruction performances and clinical value of the ordinary CS-MRI algorithms and our ROI-CS-MRI algorithm qualitatively and quantitatively.Experimental results demonstrate that the proposed NLTV-ROI-NESTA algorithm is superior in reconstruction accuracy and detail features when the undersampled ratio changed from 10% to 50%,which can be extended and widely used in rapid medical imaging technology.

Key words: Compressed sensing,Image recovery,Nonlocal total variation,NESTA algorithm,Region of interests

赵扬,王伟,董蓉,王敬时,汤敏. 基于NLTV和NESTA的MRI/MRA图像感兴趣区域的压缩感知重构[J]. 计算机科学, 2017, 44(9): 308-314. https://doi.org/10.11896/j.issn.1002-137X.2017.09.058

ZHAO Yang, WANG Wei, DONG Rong, WANG Jing-shi and TANG Min. Compressed Sensing Recovery Algorithm for Region of Interests of MRI/MRA Images Based on NLTV and NESTA[J]. Computer Science, 2017, 44(9): 308-314. https://doi.org/10.11896/j.issn.1002-137X.2017.09.058

参考文献

[1] DONOHO D.Compressed Sensing[J].IEEE Transactions onInformation Theory,2006,52(4):1289-1306.
[2] CANDES E J,ROMBERG J,TAO T.Stable Signal Recoveryfrom Incomplete and Inaccurate Measurements [J].Communications on Pure and Applied Mathematics,2006,59(8):1207-1223.
[3] CANDES E J,ROMBERG J,TAO T.Robust Uncertainty Principles:Exact Signal Reconstruction from Highly Incomplete Fourier Information[J].IEEE Transactions on Information Theo-ry,2006,52(2):489-509.
[4] ELDAR Y,KUTYNIOK G.Compressed Sensing:Theory andApplications [M].UK:Cambridge University Press,2012:1-68.
[5] 杨正汉,冯逢,王霄英.磁共振成像技术指南——检查规范、临床策略及新技术应用[M].北京:人民军医出版社,2013:162-250.
[6] 陈武凡,康立丽.MRI原理与技术[M].北京:科学出版社,2014:50-150.
[7] RUNGE V M.Current Technological Advances in Magnetic Re-sonance with Critical Impact for Clinical Diagnosis and The-rapy[J].Investigative Radiology,2013,48(12):869-877.
[8] ZHANG G S,XIAO G,DAI Z Z,et al.Compressed SensingTechnology and its Application in MRI [J].Chinese Journal of Magnetic Resonance Imaging,2013,4(4):314-320.(in Chinese) 张桂珊,肖刚,戴卓智,等.压缩感知技术及其在MRI上的应用[J].磁共振成像,2013,4(4):314-320.
[9] LUSTIG M,DONOHO D,PAULY J.Sparse MRI:The Application of Compressed Sensing for Rapid MR Imaging [J].MagneticResonance in Medicine,2007,58(6):1182-1195.
[10] JIAO P F,LI L,ZHAO J I.New Advance of Compressed Sen-sing in Medical Image Reconstruction [J].CT Theory and Applications,2012,21(1):133-147.(in Chinese) 焦鹏飞,李亮,赵骥.压缩感知在医学图像重建中的最新进展[J].CT理论与应用研究,2012,21(1):133-147.
[11] LV Y,WU W T,LI P.Simulation of Application of Compressive Sensing to Medical Ultrasound Imaging[J].Technical Acoustics,2013,32(2):106-110.(in Chinese) 吕燚,吴文焘,李平.压缩感知在医学超声成像中的仿真应用研究[J].声学技术,2013,32(2):106-110.
[12] WANG L Y,LIU H K,LI L,et al.Review of Sparse Optimization-based Computed Tomography Image Reconstruction from Few-view Projections[J].Acta Physics Sinica,2014,63(20):208702.(in Chinese) 王林元,刘宏奎,李磊,等.基于稀疏优化的计算机断层成像图像不完全角度重建综述[J].物理学报,2014,63(20):208702.
[13] LIU Y S,ZHAN Z F,CAI J F,et al.Projected Iterative Soft-thresholding Algorithm for Tight Frames in Compressed Sen-sing Magnetic Resonance Imaging [J].IEEE Transactions on Medical Imaging,2016,5(9):2130-2140.
[14] ZHANG Y D,DONG Z C,PHILLIPS P,et al.ExponentialWavelet Iterative Shrinkage Thresholding Algorithm for Compressed Sensing Magnetic Resonance Imaging [J].Information Sciences,2015,322(1):115-132.
[15] HOT E,SEKULIC P.Compressed Sensing MRI using Masked DCT and DFT Measurements .http://mailer.oalib.net/paper/4072265.
[16] PEJOSKI S,KAFEDZISKI V,GLEICH D.Compressed Sensing MRI Using Discrete Nonseparable Shearlet Transform and FISTA [J].IEEE Signal Processing Letters,2015,22(10):1566-1570.
[17] JIANG M F,LIU Y,XU W L,et al.Efficient Compressed Sen-sing Reconstruction Using Group Sparse Total Variation Regularization [J].Journal of Medical Imaging and Health Informa-tics,2015,5(5):907-917.
[18] ZHANG H M,CHEN J L,YAN B,et al.Reconstruction of CT Sparse Angle Based on Nonlocal TV Regularization [C]∥National Symposium on Digital Radiography and CT in 2014.Xiamen,China,2014.(in Chinese) 张瀚铭,陈建林,闫镔,等.基于非局部TV正则化的CT稀疏角度重建[C]∥2014年全国射线数字成像与CT新技术研讨会论文集.厦门,中国,2014.
[19] AMARESHA S K,SHIVARAJ A,S HASHIKALA H C,et al.Application of Region of Interest Compressed Sensing to Acce-lerate Magnetic Resonance Angiography[C]∥2014 36th Annual InternationalConference of the IEEE Engineering in Medicine and Biology Society (EMBC).IEEE,2014:2428-2431.
[20] ZHANG C G,VNAN DE G G, EISEMANN E,et al.User-guided Compressed Sensing for Magnetic Resonance Angiography [C]∥2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).IEEE,2014:2416-2419.
[21] OH H,LEE S.Visually Weighted Reconstruction of Compressive Sensing MRI[J].Magnetic Resonance Imaging,2014,32(3):270-280.
[22] GILBOA G,OSHER S.Nonlocal Operators with Applications to Image Processing[J].Multiscale Modeling & Simulation,2008,7(3):1005-1028.
[23] ELMOATAZ A,LEZORAY O, BOUGLEUXOUGLEUX S.Nonlocal Discrete Regularization on Weighted Graphs:A Framework for Image and Manifold Processing [J].IEEE Transaction on Image Processing,2008,17(7):1047-1060.
[24] ZHANG X Q,BURGER M,BRESSON X,et al.BregmanizedNonlocal Regularization for Deconvolution and Sparse Reconstruction[J].SIAM Journal on Imaging Science,2010,3(3):253-276.
[25] NESTEROV Y.A Method for Unconstrained Convex Minimization Problem with the Rate of Convergence O(1k²)[J].Doklady AN USSR (translated as Soviet Math.Docl.),1983,269(3):543-547.
[26] NESTEROV Y.Introductory Lectures on Convex Optimization:A Basic Course[M].Boston:Kluwer Academic Publishers,2004:50-150.
[27] NESTEROV Y.Smooth Minimization of Non-Smooth Functions[J].Mathematical Programming,2005,103(1):127-152.
[28] NESTEROV Y.Gradient Methods for Minimizing Composite Objec-tive Function [R].Universite Catholique de Louvain,2007.
[29] BECKER S,BOBIN S J,CANDES E J.NESTA:A Fast and Accurate First-Order Method for Sparse Recovery [J].SIAM Journal on Imaging Sciences,2011,4(1):1-39.
[30] XIAO C,SALERNO M,YANG Y,et al.Motion-Compensated Compressed Sensing for Dynamic Contrast-Enhanced MRI Using Regional Spatiotemporal Sparsity and Region Tracking:Block LOw-rank Sparsity with Motion-guidance (BLOSM) [J].Magnetic Resonance in Medicine,2014,72(4):1028-1038.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed