一种磁共振定量磁化率成像多回波相位拟合方法
【专利摘要】本发明公开了一种磁共振定量磁化率成像多回波相位拟合方法，包括以下步骤：利用多回波梯度回波序列采集得到磁共振原始数据；对其中被试第一回波模图噪声做统计，得到噪声强度水平；根据模图噪声水平设置截断信号的阈值；从第五回波起，对磁共振多回波信号进行选择性截断，将信号幅值低于既定阈值的回波截断，只允许高于或等于既定阈值的回波参与相位拟合；对符合阈值要求的回波采用加权最小二乘法进行相位拟合，得到磁场分布图；对磁场分布图进行时间域解缠绕、去背景场等后续处理，最终由场到源反演出磁化率分布图。本发明磁共振定量磁化率成像多回波相位拟合方法具有快速、精确的优点。
【专利说明】
一种磁共振定量磁化率成像多回波相位拟合方法
技术领域
[0001 ]本发明涉及磁共振成像的技术领域，具体涉及一种定量磁化率成像多回波相位拟 合方法。
【背景技术】
[0002] 定量磁化率成像利用一般磁共振成像技术舍弃的相位信息得到局部磁场变化特 性，再通过复杂的场到源的反演计算，可直接得到定量磁化率分布图。单回波或多回波梯度 回波序列可用来获取磁化率分布图反演计算需要的磁场分布图。与单回波梯度回波序列相 比，多回波序列一次激发可采集多个回波，这有利于相位解缠绕，且对磁化率更加敏感，得 到的磁化率分布图信噪比也更高。有研究表明采用5~10个回波可获得较好的颅内核团磁 化率对比。
[0003] 磁化率分布图重建时，多回波相位数据首先需进行时间域相位解缠绕，然后通过 线性拟合得到磁场分布图，再进行后续处理。常用的线性拟合方法是加权最小二乘法，其原 理是按照一定规则设置加权系数，以调整各回波在拟合中的加权，使拟合出的相位逐步逼 近真实值。加权系数的设置可根据磁共振信号随回波时间的衰减特性，从模图提取信息，将 回波时间较短、信号较强的回波赋予较大加权;反之，回波时间较长、信号较弱的回波赋予 较小加权。这种加权系数结合了模图中的组织结构信息，一定程度上突出了组织间的对比。
[0004] 常规加权最小二乘方法对于所有组织都采用相同的回波数，无法兼顾不同组织获 取较高信噪比。大脑的某些组织，尤其颅底部位，其磁化率分布很不均匀，局部磁场变化较 大，造成相位迅速离散，若回波时间较长，信号甚至可能衰减至噪声水平。因此对于这些核 团，采用所有回波反而将噪声信号也一并拟合到局部磁场中，降低了磁化率分布图的信噪 比。常规加权最小二乘法得到的磁场分布图在颅底等部位噪声较大，造成最终的磁化率分 布图信噪比较低，无法准确测量颅底组织的磁化率值。
【发明内容】
[0005] 本发明的目的是为了克服上述常规加权最小二乘法的不足，而提供的一种磁共振 定量磁化率成像多回波相位拟合方法，以期减小磁场分布图的拟合误差，提高颅底的磁化 率分布图的图像质量，准确测量大脑各组织的磁化率值。
[0006] 本发明提出一种磁共振定量磁化率成像多回波相位拟合方法。该方法首先对多回 波梯度回波序列采集的第一回波的模图噪声做统计，从而得到磁共振图像的噪声强度水 平。然后在采用常规加权最小二乘法拟合之前引入阈值截断条件，根据不同组织的信号衰 减程度，合理选择回波数目。将幅值小于既定阈值的回波截断，只允许高于或等于阈值的回 波参与相位拟合，从而减小磁场分布图的拟合误差，提高颅底部位磁化率分布图的信噪比。
[0007] 实现本发明上述目的所采用的技术方案是：
[0008] -种磁共振定量磁化率成像多回波相位拟合方法，特点是所述方法包括以下步 骤：
[0009] A)利用多回波梯度回波序列采集得到磁共振原始数据；
[0010] B)对其中第一回波的模图噪声进行统计，得到噪声强度水平；
[0011] C)根据模图噪声水平设置信号幅值的截断阈值；
[0012] D)设置加权系数，将信号幅值低于既定阈值的回波截断，只允许高于或等于既定 阈值的回波参与相位拟合；
[0013] E)对符合阈值要求的回波采用加权最小二乘法进行相位拟合，得到磁场分布图；
[0014] F)对磁场分布图进行时间域解缠绕、去除背景场等后续处理，最终由场到源反演 出磁化率分布图。
[0015] 所述步骤B)中的噪声强度水平是第一回波模图四个角内感兴趣区中体素的幅值 标准差作为噪声，表示为公式[1 ]
[0016]
[X]
[0017] 其中，SiS第一回波模图背景图像中的某体素信号强度，？为感兴趣区内体素信号 强度的平均值，I是感兴趣区内的体素个数。
[0018] 所述步骤C)中截断阈值设置为模图噪声的m倍，表示为公式[2]
[0019] a = 〇*m [2]
[0020] 其中，σ为步骤B)获取的噪声强度水平，m为经验值。
[0021]所述步骤D)中设置的加权系数表示为公式[3]
[0022]
[3]
[0023]其中，η取值为1，2,…，N，N为最大回波数目；为保证组织对比度，从第五回波开始 进行判断，Mn为某体素第η个回波的信号幅值，a为步骤C)中获取的阈值。
[0024] 所述步骤E)中加权最小二乘法进行多回波相位拟合，根据拟合值与测量值之间的 误差不断进行迭代加权，可求解出初始相位<Po(r)和局部磁场变化ΔΒ( Γ)的最优解，ΔΒ(Γ) 即为所求磁场分布图，表示为公式[4]。
[0025]
[0026] 其中，表示为某体素第η个回波的相位测量值，γ为氢质子的旋磁比，Wn 为该体素第η个回波的加权系数，TEn为第η个回波的回波时间。
[0027] 本发明中，"加权最小二乘法"表示拟合方法中的一种数学术语。
[0028] 与现有技术相比，采用本发明所述技术方案，可达到以下技术效果：
[0029] 1、白质等磁化率均匀的区域保持较高信噪比，组织对比明显。
[0030] 2、黑质、红核等颅底部位信噪比明显提高，所测得的磁化率值更精确。
[0031 ] 3、本发明方法适用于放射科临床上对脑铁等磁敏感源的定量测量。
【附图说明】
[0032]图1为本发明流程图；
[0033] 图2为本发明磁共振多回波梯度回波采集得到的数据图例；
[0034] 图3为某一像素的信号幅值随回波时间衰减示意图；
[0035] 图4为利用本发明得到的磁场分布图；
[0036] 图5为利用本发明得到的磁化率分布图。
【具体实施方式】
[0037] 结合以下具体实施例和附图，对本发明作进一步详细说明。实施本发明的过程、条 件、实验方案方法等，除以下专门提及的内容之外，均为本领域的普遍知识和公知常识，本 发明没有特别限制内容。
[0038] 本发明定量磁化率多回波相位拟合方法，在采用常规加权最小二乘法拟合之前引 入阈值截断条件。根据不同组织的信号衰减程度，合理选择回波数目，将幅值小于既定阈值 的回波截断，只允许幅值高于或等于既定阈值的回波参与相位拟合，从而达到降低多回波 拟合误差的目的。该方法一定程度上能够提高颅底部位磁化率分布图的信噪比。本发明中， 噪声通过对磁共振图像中被试第一回波模图的四个角感兴趣区中体素幅值的标准差统计 得到。
[0039] 以下分步介绍本发明对多回波梯度回波序列采集得到的原始数据进行截断-加权 最小二乘法拟合，生成磁场分布图的具体实施过程。其中，采集的磁共振图像数据由多回波 梯度回波序列获取，数据来源于3T磁共振成像设备系统（西门子MAGNETOM Trio a Tim 3T)，本实施例中所采用的回波个数为8。
[0040] Α)在3Τ磁共振成像设备系统上对被试人脑进行多回波梯度回波序列扫描，获得多 回波的原始数据，包括原始相位图与原始模图。
[0041] 其中，利用多回波梯度回波序列扫描被试人脑的具体实施为磁共振扫描的一般过 程，所获得的原始图像为不同回波时间ΤΕ采集得到的多个三维图像。本实施例中，8个回波 时间 ΤΕ 分别为 6.8、13.6、20.4、27.2、34.0、40.8、47.6、54.4(单位:1^)。8个回波模图强度随 回波时间ΤΕ的增大而衰减，如附图2所示。
[0042] Β)对其中的第一回波模图噪声做统计，得到噪声强度水平。
[0043]选取被试第一回波模图四个角内感兴趣区中体素的幅值标准差作为噪声，每个感 兴趣区内包含1000个体素，具体可表示为公式[1 ]。
[0044]
[1]
[0045] 本实施例中，第一模图背景噪声水平为2.28。
[0046] C)根据模图噪声水平设置信号的截断阈值。本实施例中，为得到磁共振信号的截 断阈值，提取步骤Β)获得的模图噪声，再与经验值相乘，得到截断阈值为68.4,具体由公式
[2]计算得出。
[0047] a = 〇*m [2]
[0048] 其中〇是模图背景噪声，此处取经验值m为30。
[0049] D)设置加权系数，将低于既定阈值的回波截断，只允许高于或等于既定阈值的回 波参与相位拟合，具体加权矩阵可表示为公式[3]。
[0050]
[3]
[0051] 其中，为保证组织对比度，从第五回波开始进行阈值判断。本实施例中，选取被试 靠近颅底区域某一像素点，其8个回波的信号幅值分别为106.93、98.53、93.16、88.0、 84.29、82.95、65.31、49.38，信号幅值随回波时间TE的衰减趋势如附图3所示。经判断，第7、 8回波的模图幅值小于步骤C)获得的阈值68.4,因此，从第七回波开始截断，只允许前6个回 波参与相位拟合。
[0052] E)对保留的所有回波采用加权最小二乘法进行相位拟合，得到磁场分布图。对于 多回波梯度回波序列，空间某个像素第η回波的相位可表示为：
[0053]
[5]
[0054] 其中，<pQ(r)表示初相位，γ表示氢质子的旋磁比，r表示像素的空间坐标，Ν表示最 大回波数目，AB(r)表示局部磁场变化。对于等间隔回波时间，TEnS：表示为TEfTE^n-l) Δ TE，△ TE表示回波间隔。利用加权最小二乘法可拟合出初相位φσ0-)与局部磁场变化Δ B ㈦，即：
[0055]
[0056] 其中，Wn为各回波的加权系数。根据拟合值与测量值之间的误差不断进行迭代加 权，可求解出物的和Δ B (r)的最优解。
[0057]本实施例中，将步骤D)中获取的加权系数％带入公式[4]中，最终拟合出Δ B(r)。 如附图4所示为本实施例利用本发明得到的磁场分布图。
[0058] F)对磁场分布图进行时间域解缠绕、去除背景场等后续处理，最终由场到源反演 出磁化率分布图。如附图5所示为本实施例利用本发明方法得到的磁场分布图采用常规反 演方法得到的磁化率分布图。
【主权项】
1. 一种磁共振定量磁化率成像多回波相位拟合方法，其特征在于，所述方法包括以下 步骤： A) 利用多回波梯度回波序列采集得到磁共振原始数据； B) 对其中第一回波的模图噪声进行统计，得到噪声强度水平； C) 根据模图噪声水平设置信号幅值的截断阈值； D) 设置加权系数，将信号幅值低于既定阈值的回波截断，只允许高于或等于既定阈值 的回波参与相位拟合； E) 对符合阈值要求的回波采用加权最小二乘法进行相位拟合，得到磁场分布图； F) 对磁场分布图进行时间域解缠绕、去除背景场等后续处理，最终由场到源反演出磁 化率分布图。2. 根据权利要求1所述磁共振定量磁化率成像多回波相位拟合方法，其特征在于，所述 步骤B)中的噪声强度水平是第一回波模图四个角内感兴趣区中体素的幅值标准差作为噪 声，表示为公式[1][1] 其中，S1S第一回波模图背景图像中的某体素信号强度，？为感兴趣区内体素信号强度 的平均值，I是感兴趣区内的体素个数。3. 根据权利要求1所述磁共振定量磁化率成像多回波相位拟合方法，其特征在于，所述 步骤C)中截断阈值设置为模图噪声的m倍，表示为公式[2] α = σ*ηι [2] 其中，σ为步骤B)获取的噪声强度水平，m为经验值。4. 根据权利要求1所述磁共振定量磁化率成像多回波相位拟合方法，其特征在于，所述 步骤D)中设置的加权系数表示为公式[3][3] 其中，η取值为1，2,…，N，N为最大回波数目；为保证组织对比度，从第五回波开始进行 判断，Mn为某体素第η个回波的信号幅值，α为步骤C)中获取的阈值。
【文档编号】A61B5/00GK105997074SQ201610273836
【公开日】2016年10月12日
【申请日】2016年4月28日
【发明人】赵欣欣, 李建奇, 薄斌仕, 王乙
【申请人】华东师范大学