基于邻域共享压缩感知的动态对比增强磁共振快速成像方法
【技术领域】
[0001] 本发明是一种集成压缩感知（compressed sensing, CS)技术的邻域共享（view sharing)K空间采样轨迹优化采集方法，属于磁共振医学成像技术领域，能够提供高时间分 辨率和高信噪比的脏器动态对比增强（DCE-MRI)图像，可以用于腹部脏器等的疾病诊断、 功能参数定量测量和预后评价等。
【背景技术】
[0002] 慢性肾脏病和肝癌已经成为全球性公共健康问题，在发达国家，普通人群的慢性 肾病患病率高达6. 5 %~16 %，给全球卫生财政带来了沉重负担。根据2012年一项《中国慢 性肾脏病流行病学调查》结果显示，我国成年人群中慢性肾脏病的患病率已经高达10. 8%， 预计全国现有成年慢性肾脏病患者1. 2亿人。此外，发达国家普通人群的肝癌发病率高达 ASR 2. 7~8. 6,且发病率在逐年提高，患者死亡率高达ASR 2. 5~7.1。肝癌已经成为世界 范围内男性癌症致死的第二大疾病，仅2012年全球新增肝癌患者78万，死亡74万，其中， 中国肝癌患病人数和死亡人数约占到50%。腹部临床诊断中，基于磁共振成像技术的动态 对比增强检查是诊断肾脏疾病和肝脏疾病的有效手段。因此，改善腹部动态对比增强磁共 振成像的时间分辨率和图像质量，以达到更加准确的脏器功能定量测量，对于腹部脏器疾 病的临床诊断及预后评价具有非常重要的意义。
[0003] 在腹部动态对比增强DCE-MRI成像中，由于三维成像技术能够更好的显示整个肾 脏、肝脏等器官的形态，并提供更多的图像信息及灌注曲线，从而被临床广泛应用。目前，在 临床中主要采用三维快速扰相梯度回波3D FSPGR序列以及三维随机轨迹时间分辨TWIST 序列进行腹部动态对比增强磁共振成像。
[0004] 三维快速扰相梯度回波3D FSPGR序列采用小角度射频脉冲激发信号，然后进行相 位编码采集信号，最后在梯度编码方向上加入扰相位梯度将残留的横向磁化矢量散相，以 缩短序列采集时间。3D FSPGR序列通常需要结合并行成像技术，实现较快的腹部动态对比 增强扫描，临床上由于腹部成像覆盖范围大，相位编码步数多，成像时间分辨率受到很大限 制。
[0005] TWIST序列采用邻域共享（view sharing)的三维伪随机轨迹序列，数据采集时K 空间中心区域（A)进行密集采样而外周区域（B)进行亚采样，并在数据重建时将相邻K空 间数据进行合并共享（view-shared)重建，从而缩短每一帧图像的采集时间。在三维K空 间中，TWIST采样机制在相位编码梯度k y和层选相位编码梯度k z上执行，采用极坐标方式 表达，采样轨迹随着辐射半径I以及方位角Θ (〇〈 Θ〈2 π )增加逐渐向外周进行，设定一个 辐射半径的阈值k。，在阈值k。之内K空间密集采样，k。之外K空间伪随机亚采样。K空间外 周亚采样方式首先将外周K空间区域B随机分成相互不重叠的N个子集B 1Q = Ρ··Ν)，在 进行数据采集时，首先采集一个K空间点为全采样的数据集#0,然后依次采集A，B1, A，B2, A，B1…区域数据点获得#1，#2, #3, #4, #5, #6…数据集。在重建时，将奇数次采集的#1，#3 和#5前后两个偶数K空间数据集合并到当前数据集，组合成一个完整的全采样K空间进行 图像重建，如图1所示。所以，TWIST重建图像的图像帧率取决于区域A采集的帧数。
[0006] 但是，临床上腹部三维成像时为了实现大范围覆盖并获得足够高分辨率的增强图 像，通常设定较多层数（彡32层）、较大采样矩阵（相位编码k y> 128,频率编码k 128) 和多次激发（Nex多1)，序列时间分辨率在10~60秒之间，使得采集得到的时间-信号强 度曲线分辨率低，远不能满足功能参数的精确测量，或者以牺牲空间分辨率来进行补偿，导 致图像空间分辨率下降。临床3D FSPGR序列扫描是全采样方式，通常需要结合并行成像技 术（parallel imaging)进行图像信号采集，然而，并行成像技术导致图像信噪比SNR降低， 且加速倍数通常限制在三倍以内，加速倍速越高，信噪比越低，阻碍了时间分辨率的进一步 提高。另一方面，TWIST序列扫描方式是间隔采样K空间中心区域和外周区域，重建过程中 需要将相邻K空间数据集进行合并，保证了图像较高的空间分辨率，但是使得图像帧率降 低，限制了其时间分辨率的提高。所以，研究一种既提高时间分辨率又能保证足够高的空间 分辨率的磁共振快速成像方法是提高疾病诊断和功能参数精确测量的重要前提。
【发明内容】
[0007] 针对上述问题，本发明提出了一种将压缩感知技术集成到K空间邻域共享采样轨 迹中（view-shared compressed sensing, VCS)的K空间数据快速亚采样方法，提高了腹部 动态对比增强磁共振成像的时间分辨率，同时保证了较高的图像信噪比和图像质量。
[0008] 为实现上述目的，本发明采取以下技术方案：一种动态对比增强磁共振成像的快 速三维K空间亚采样填充方法，包括K空间的随机划分、CS采样矩阵设计、邻域共享压缩感 知采样矩阵的获取、基于邻域共享压缩感知方法的数据采集、基于邻域共享压缩感知方法 的数据重建五个部分。具体内容为：首先，采用相邻K空间外周数据共享采样轨迹，一个K 空间中同时采集中心区域和部分外周区域数据，邻域共享的K空间外周亚采样轨迹相互间 无重叠，所有相邻的亚采样K空间数据组合可以构成一个完整的K空间。其次，将压缩感知 加速采样技术嵌入到相邻K空间的共享采样轨迹中，每一个K空间填充轨迹由压缩感知方 法设计确定，其中K空间的中心区域保留不共享，外周区域经压缩感知采样矩阵设计后随 机分成2~5个共享子集，以实现连续多个K空间采样轨迹的邻域共享压缩感知采样。邻 域共享压缩感知采样在一次K空间采样时，同时采集压缩感知K空间的中心区域和1个外 周共享子集。K空间数据在重建时，连续多个不相重叠的共享子集合并到当前K空间，重组 成符合压缩感知采样矩阵轨迹的数据集用于进一步压缩感知图像重建。集成到邻域共享采 样轨迹的压缩感知加速倍数可以设置为2~8倍。压缩感知采样矩阵设计采用非均匀变密 度函数生成，且符合高斯分布，压缩感知的稀疏转换可以为傅里叶变换、小波变换、有限差 分变换等。最后，将邻域共享重组的K空间亚采样数据经I 1范式最小化的非线性算法重建 得到图像，重建公式如下：
[0011] 式中S为目标图像，λ Jp λ 2表示Φ稀疏表达和有限差分稀疏表达TV之间的权 重，F表示与亚采样对应的傅立叶变换，Cl1表示VCS 1采集到的k-space数据，ε通常设置在 背景噪声水平之下来保证重建图像的保真度。并按照DICOM格式重写成临床标准化医学图 像。
[0012] 本发明由于采取以上技术方案，具有如下优点：本发明是一种高时间分辨率