定量磁化率成像多回波相位拟合算法研究.PDFVIP

波 谱 学 杂 志 第 33 卷第4 期 Vol. 33 No. 4 2016 年 12 月 Chinese Journal of Magnetic Resonance Dec. 2016 文章编号：1000-4556(2016)04-609-09 doi:10.11938/cjm 定量磁化率成像多回波相位拟合 算法研究 1 1 2 1,2,3 1 赵欣欣 ，薄斌仕 ，刘 田 ，王 乙 ，李建奇 （1. 华东师范大学 物理系，上海市磁共振重点实验室，上海 200062 ； 2. Department of Radiology, Weill Medical College, Cornell University, Ithaca, New York 10021, USA; 3. Department of Biomedical Engineering, Cornell University, Ithaca, New York 14853, USA ） 摘 要：定量磁化率成像（quantitative susceptibility mapping, QSM ）技术大多采用多回波梯 度回波序列采集相位数据，经加权最小二乘法（weighted linear least-square ，WLS ）拟合得 到局部磁场分布. 对于组织磁化率分布不均匀的区域，尤其是颅底部位，常规 WLS 算法拟 合得到的局部磁场误差较大，导致相应部位磁化率分布图信噪比较低. 针对常规 WLS 算 法的这一不足，该文提出了一种截断 WLS 算法. 对两种算法拟合得到的磁化率分布图对 比研究发现，截断WLS 算法可有效提高颅底部位定量磁化率分布图的图像质量，使其噪声 明显下降. 关键词：定量磁化率成像（QSM ）；多回波相位拟合；加权最小二乘法（WLS ）；阈值截断 中图分类号：O482.53 文献标识码：A 引言 定量磁化率成像（quantitative susceptibility mapping, QSM ）利用一般磁共振成像 （magnetic resonance imaging, MRI ）技术舍弃的相位信息得到局部磁场变化特性，再 通过复杂的场到源反演计算，可直接得到定量磁化率分布图[1–9]. 单回波或多回波梯度 回波（gradient echo, GRE ）序列可用来获取 QSM 反演计算需要的磁场分布图. 与单 收稿日期: 2016-04-04; 收修改稿日期: 2016-10-24 基金项目: 国家自然科学基金资助项目. 作者简介: 赵欣欣(1990-)，女，山东人，硕士研究生，无线电物理专业. *通讯联系人：李建奇，电话：021 E-mail ：jqli@. 610 波 谱 学 杂 志 第 33 卷 回波 GRE 序列相比，多回波序列一次激发可采集多个回波，这有利于相位解缠绕； 且对磁化率更加敏感，得到的磁化率分布图信噪比也更高[10–12]. 有研究表明采用 5~10 [13] 个回波可获得较好的颅内核团磁化率对比 . QSM 重建时，多回波相位数据首先需进行时间域相位解缠绕，然后通过线性拟 合得到磁场分布图，再进行后续处理.