[發(fā)明專利]磁共振圖像重建方法及裝置

有效

• 申請(qǐng)?zhí)枺?/td>: 201310280457.8
• 申請(qǐng)日: 2013-07-05
• 公開（公告）號(hào): CN104181486A
• 公開（公告）日: 2014-12-03
• 發(fā)明（設(shè)計(jì)）人: 翟人寬
• 申請(qǐng)（專利權(quán)）人: 上海聯(lián)影醫(yī)療科技有限公司
• 主分類號(hào): G01R33/56
• 分類號(hào): G01R33/56;G01R33/48
• 代理公司: 無
• 代理人: 無

【技術(shù)領(lǐng)域】

本發(fā)明涉及磁共振成像領(lǐng)域，尤其是涉及磁共振圖像重建方法及裝置。

【背景技術(shù)】

目前，在磁共振多通道采集圖像重建過程中，一般采用SOS等方法對(duì)各個(gè)通道采集到的數(shù)據(jù)進(jìn)行簡(jiǎn)單的合并。此方法的缺陷在于：信噪比較低的通道在合并過程中的權(quán)重沒有很好的被抑制，會(huì)影響整體的成像質(zhì)量。因此，有必要進(jìn)一步研究提高圖像信噪比的方法。

【發(fā)明內(nèi)容】

為了解決上述中圖像噪聲抑制不明顯，獲取的圖像質(zhì)量不夠理想的問題，本發(fā)明提供了一種磁共振圖像重建方法及裝置。

一種磁共振圖像重建方法，所述方法包括：

a1)采集獲得若干通道全采樣的原始K空間；或采集獲得若干通道欠采樣的K空間，并利用相關(guān)方法獲得全采樣的原始K空間；

a2)對(duì)各通道原始K空間按照相同的規(guī)則進(jìn)行預(yù)處理分組，獲得若干分組K空間，所述分組K空間為與所述原始K空間相同大小的矩陣，所述每一分組K空間都包括原始K空間中的部分已采樣數(shù)據(jù)點(diǎn)，以及待填補(bǔ)數(shù)據(jù)；

a3)利用原始K空間數(shù)據(jù)作為校準(zhǔn)數(shù)據(jù)，選取卷積核K0求得合并系數(shù)；根據(jù)所述合并系數(shù)使用并行采集重建方法對(duì)所述每一分組K空間中的每一待填補(bǔ)數(shù)據(jù)進(jìn)行填補(bǔ)，將填補(bǔ)完成后的各分組K空間數(shù)據(jù)進(jìn)行合并，并更新所述各通道原始K空間數(shù)據(jù)；

a4)利用所述更新之后的K空間數(shù)據(jù)作為校準(zhǔn)參數(shù)，選取卷積核Ki(i為迭代處理的次數(shù)，i≥1)求得新的合并系數(shù)；根據(jù)所述新的合并系數(shù)使用并行采集重建方法對(duì)所述每一分組K空間中的每一待填補(bǔ)數(shù)據(jù)進(jìn)行填補(bǔ)，將填補(bǔ)完成后的各分組K空間數(shù)據(jù)進(jìn)行合并，并再次更新所述各通道K空間數(shù)據(jù)；

a5)判斷迭代處理是否完成，若否，返回執(zhí)行步驟a4)；

a6)若是，將最后一次更新的各通道K空間數(shù)據(jù)轉(zhuǎn)換至圖像域，獲得各通道圖像，合并所述各通道圖像獲得最終圖像。

優(yōu)選地，所述迭代次數(shù)i大于等于3時(shí)，判定迭代計(jì)算已經(jīng)完成。

優(yōu)選地，利用兩次相鄰迭代處理之后的K空間數(shù)據(jù)之差的絕對(duì)值小于等于閾值時(shí)，判定迭代處理已經(jīng)完成。

優(yōu)選地，當(dāng)兩次相鄰迭代處理之后的K空間數(shù)據(jù)之差絕對(duì)值中的最大值小于等于閾值A(chǔ)時(shí)，判定迭代處理已經(jīng)完成。

優(yōu)選地，當(dāng)兩次相鄰迭代處理之后的K空間數(shù)據(jù)之差絕對(duì)值中所有的值之和小于等于閾值B時(shí)，判定迭代處理已經(jīng)完成。

優(yōu)選地，進(jìn)行不同次迭代處理的所選取的卷積核相同或者不同。

優(yōu)選地，所述卷積核K0大小為4×3(選取PE方向4個(gè)數(shù)據(jù)點(diǎn)，RO方向3個(gè)數(shù)據(jù)點(diǎn))或者2×5(選取PE方向2個(gè)數(shù)據(jù)點(diǎn)，RO方向5個(gè)數(shù)據(jù)點(diǎn))。

優(yōu)選地，使用GRAPPA方法對(duì)所述每一分組K空間中的待填補(bǔ)數(shù)據(jù)點(diǎn)進(jìn)行填補(bǔ)，具體包括以下步驟：

根據(jù)原始K空間數(shù)據(jù)計(jì)算GRAPPA方法的合并系數(shù)；

利用所述合并系數(shù)，對(duì)每一分組K空間中的每一待填補(bǔ)數(shù)據(jù)進(jìn)行填補(bǔ)，所述待填補(bǔ)數(shù)據(jù)具體可通過如下公式獲得：

DcTm=(G-I)DTa,]]>

其中，G為GRAPPA方法的合并系數(shù)，所述D^Ta為每一通道相應(yīng)分組K空間的已采樣數(shù)據(jù)，為待填補(bǔ)數(shù)據(jù)。

優(yōu)選地，步驟a1)中所述相關(guān)方法為并行采集重建方法或部分傅里葉重建方法。

優(yōu)選地，將填補(bǔ)完成后的各分組K空間數(shù)據(jù)進(jìn)行合并，具體為：將填補(bǔ)完成后的各分組K空間數(shù)據(jù)相加或者將填補(bǔ)完成后的各分組數(shù)據(jù)分別取絕對(duì)值后相加。

本發(fā)明還提供了另一種磁共振圖像的重建方法，所述方法包括如下步驟：

b1)采集獲得若干通道欠采樣的原始K空間，所述原始K空間中至少包括一個(gè)全采樣區(qū)域，所述全采樣區(qū)域中每一數(shù)據(jù)點(diǎn)都已采樣；

b2)對(duì)各通道原始K空間的已采樣數(shù)據(jù)按照相同的規(guī)則進(jìn)行預(yù)處理分組，獲得若干分組K空間，所述每一分組K空間都包括部分已采樣數(shù)據(jù)點(diǎn)，以及待填補(bǔ)數(shù)據(jù)；