[發明專利]一種磁共振成像多回波磁場分布圖擬合方法

有效

• 申請號: 201910831438.7
• 申請日: 2019-09-04
• 公開（公告）號: CN110687489B
• 公開（公告）日: 2021-08-20
• 發明（設計）人: 陳佳林;李建奇
• 申請（專利權）人: 華東師范大學
• 主分類號: G01R33/565
• 分類號: G01R33/565;G01R33/48;A61B5/055
• 代理公司: 上海藍迪專利商標事務所(普通合伙) 31215
• 代理人: 徐筱梅;張翔

1.一種磁共振成像多回波磁場分布圖擬合方法，其特征在于，該方法包括以下步驟：

步驟1：利用多回波梯度回波序列采集得到磁共振原始數據；

步驟2：對所有回波的相位進行線性度檢查，舍棄線性度差的數據點的部分或者全部回波信號；

步驟3：將經過線性度檢查后的相位圖利用非線性最小二乘法進行擬合得到初步的磁場分布圖；

步驟4：利用周圍可靠數據點即相位線性度好的點的局部磁場的平均值來替代回波信號全部舍去的點的局部磁場，最終得到完整的磁場分布圖；其中：

步驟2中所述的相位進行線性度檢查，具體步驟為：

第一步，對前三個回波進行線性擬合得到初始相位近似值和角頻率近似值ω′，具體表述為如下線性方程：

其中，為初始相位理論值，ω為角頻率理論值，TE_j為第j個回波的回波時間，為第j個回波的相位；

第二步，奇異點判斷

采用擬合得到的初始相位近似值和角頻率近似值計算圖像上每一個像素點前三個回波的擬合相位與測量相位之間的平均差值如下式：

設置一個誤差閾值th，如果則認為該像素點相位誤差嚴重，在場圖上為奇異點，對奇異點之后再做進一步處理；

第三步，非奇異點回波可靠性判斷

利用下式計算第j(N≥j＞3)個回波的擬合相位與測量相位的差值δ_j：

如果δ_j＞th，則認為該回波相位可靠性低，舍棄該回波及其后面的所有回波信號。

2.根據權利要求1所述的磁共振成像多回波磁場分布圖擬合方法，其特征在于，步驟3中所述利用非線性最小二乘法進行擬合得到初步的磁場分布圖是利用最小二乘法進行多回波場圖擬合，根據擬合值與測量值之間的誤差不斷進行迭代，求解出局部磁場即磁場分布圖的最優解，如下式：

其中，s(r，TE_j)為信號測量值，TE_j為回波時間，N為回波個數，經過相位檢查后線性度不同的數據點所使用的回波個數不同，a(r，TE_j)為每個回波信號的幅值，γ為旋磁比，γΔB(r)為所求場圖，單位為頻率(Hz)。