基于k-空间钱采样的动态MRI重建方法研究

发布时间：2020-06-01 23:15

【摘要】：磁共振成像(Magnetic Resonance Imaging,MRI)已成为医学临床诊断和科学研究的重要技术手段。受各种成像条件的限制,磁共振成像的检测扫描时间过长,扫描过程中由于病人需要长时间保持某种姿态,一方面会造成病人的不适感,另一方面也容易导致运动伪影的产生。动态磁共振成像(Dynamic MRI,dMRI)是一种在传统MRI设备上通过连续快速扫描获取人体组织器官随时间变化的运动影像技术。相比静态MRI,dMRI图像具有更高的空间分辨率和时间分辨率,但其成像条件要求更高、扫描时间更长。一段时间以来,由于无法克服奈奎斯特采样定律的限制,如何缩短扫描时间、提高成像速度成为困扰MRI研究人员的难题。压缩感知(Compressive Sensing,CS)理论的提出,使得从高度欠采样k空间数据中恢复重建MRI图像成为可能。利用MRI图像数据固有的稀疏性,将CS理论用于MRI成像应用(CS-MRI)可有效解决磁共振成像扫描时间长、成像速度慢的问题。但直接将CS-MRI技术用于dMRI图像重建,会面临因传统稀疏表示无法解决dMR丨图像在空间-时间域的稀疏问题,导致在自由呼吸状态下的动态心脏MRI或心肌灌注MRI成像的应用中表现出较差的重建效果。为解决CS-MRI技术在dMRI图像重建应用中的缺陷,本文利用dMRI图像在空域和时间域的稀疏性,提出了一种盲字典学习方法。该方法从k-t空间欠采样数据中学习生成空-时字典对dMRI图像进行稀疏表示,用于dMRI图像的加速重建。另外,在CS理论框架下,为降低字典学习训练中计算复杂性和字典中原子的过训练等问题,本文利用dMRI图像数据在空-时域的低秩性,基于低秩稀疏矩阵分解模型分别提出了一种凸优化和非凸优化方法,用于从k-t空间欠采样数据中加速dMRI图像的重建并改善重建dMRI图像的性能。本文主要工作与创新点如下:1)提出了一种基于盲字典学习的dMRI图像重建方法,直接从高度欠采样的k-t空间数据中自适应学习、生成字典并完成对图像的稀疏表示,并将字典学习与图像信号的稀疏表示组成联合优化问题,利用最小优化算法求解稀疏表示系数和时间基字典。该方法主要解决了经典CS方法和基于K-SVD的盲压缩感知方法对采样观测矩阵的限制。采用优化最小算法和共轭梯度算法求解dMRI图像的稀疏表示系数,克服了传统凸优化方法直接求解高维度矩阵稀疏解时,需要占用大量计算资源而导致计算收敛慢、图像重建时间过长等问题。实验测试结果表明,在相同采样加速因子下,本文使用的盲字典学习方法的图像重建性能优于对比方法。2)提出了一种利用截断核范数求解鲁棒主成分分析问题的dMRI图像重建方法(Truncated Nuclear Norm for Robust Principal Component Analysis,TNN-RPCA),用于解决在低秩稀疏矩阵分解模型中,传统凸优化方法直接使用核范数对非凸秩函数进行优化时,因忽略非零奇异值对秩函数贡献的差异而导致的重建误差。TNN-RPCA方法使用截断核范数代替传统核范数对秩函数凸优化,通过引入具有快速运算性能和高精度的非精确增广拉格朗日乘子法(Inexact Augmented Lagrange Multiplier,IALM)用于RPCA问题求解。通过仿真实验,并与两种用于求解RPCA问题重建dMRI图像的方法进行了对比,结果表明TNN-RPCA方法无论在重建图像的速度还是重建图像的质量方面,均具有明显的优势。3)提出了一种利用非凸优化奇异值收缩算法求解RPCA问题的dMRI图像重建方法(Optimal Singular value shrinkage for Robust Principal Component Analysis,OptShrink-RPCA),用于解决直接使用核范数对非凸秩函数进行凸优化无法获得最优解的应用缺陷。该方法首先对测量矩阵的奇异值进行截断,然后将截断剩下的奇异值通过门限阈值进行收缩处理,以逼近非凸秩函数的最优解。同时引入增广拉格朗日乘子加速RPCA问题的求解。通过实验与几种凸优化方法对dMRI图像重建的效果进行对比,发现使用OptShrink-RPCA方法重建dMRI图像的性能得到改善。
【图文】：

示意图,信号测量,示意图,变换矩阵

逦＝邋ＱＡ逦（２－５）逡逑２－５式中０也称信息算子，信号测量示意图如图２－１所示。ＣＳ理论指出，信息算子０或者测量逡逑矩阵０需要满足ＲＩＰ条件。Ｅ．邋Ｃａｎｄ红等人证明当Ｏ为随机矩阵时能以大概率满足ＲＩＰ条件，由逡逑于判断测量矩阵是否满足ＲＩＰ条件比较困难，Ｅ．邋Ｃａｒｎｅｓ在文献［９１］中提出相对简单的非相干逡逑性判别准则，测量矩阵０和稀疏变换矩阵Ｖ之间的相干性定义为：逡逑／／（＜！），ｖ｛／）邋＝邋Ｖ７７－逦（２－６）逡逑２－６式求解的相干性Ｗ表示稀疏信号的维度大小。／／（屯中）越小，表明测量逡逑矩阵￥与变换矩阵Ｖ的非相千性越大。文献［９１］实验结果给出：当采样矩阵为冲击函数，稀疏逡逑变换为傅里叶变换时，采样矩阵与变换矩阵之间的相干系数Ｍ（电Ｖ）邋＝邋ｌ，具有最大的非相干逡逑性。另外

脉冲序列,梯度线圈,扫描机,成像技术

ＤＯＣＴＯＲＡＬ邋ＤＩＳＳＥＲＴＡＴＩＯＮ逡逑２．２．邋１邋ＭＲ成像扫描方法逡逑为获得如图２－３所示的置身于外磁场中人体组织或器官影像的空间信息，需要由梯度线圈逡逑分别在Ｘ、７和２三个方向上产生梯度磁场，这三个梯度可分别作为层面选择梯度、相位编码梯逡逑度和频率编码（或读出）梯度。磁共振扫描过程中，可以根据成像需要任意选择并组合梯度线逡逑圈，分别在不同平面（横断、矢状和冠状）成像。为提高ＭＲＩ扫描成像效率，单扫描成像技术逡逑的应用比较流行［１３５］。扫描过程中通过调整ＭＲ成像参数，可获得组织器官的不同成像效果。由逡逑于影响ＭＲ成像效果的参数种类繁多，使得ＭＲ成像脉冲序列的设计变得非常复杂。目前普遍逡逑使用的扫描成像技术主要有快速自旋回波（ＦａｓｔＳｐｉｎＥｃｈｏ，ＦＳＥ）成像和回波平面成像（Ｅｃｈｏ－逡逑Ｐｌａｎａｒ邋Ｉｍａｇｉｎｇ，，ＥＰＩ）等方法。逡逑ｆＴＴＪ］逡逑图２－３指定扫描机内的梯度线圈方向逡逑其中ＦＳＥ是在自旋回波成像技术基础上通过单次激发同时填充多行ｋ空间数据
【学位授予单位】：华中师范大学
【学位级别】：博士
【学位授予年份】：2018
【分类号】：R445.2;TP391.41

【相似文献】