[發(fā)明專利]基于復(fù)奇異譜分析的磁共振部分K數(shù)據(jù)圖像重建方法

無效

• 申請?zhí)枺?/td>: 200710039881.8
• 申請日: 2007-04-24
• 公開（公告）號: CN101051075A
• 公開（公告）日: 2007-10-10
• 發(fā)明（設(shè)計）人: 駱建華;朱躍敏
• 申請（專利權(quán)）人: 駱建華
• 主分類號: G01R33/56
• 分類號: G01R33/56
• 代理公司: 上海智信專利代理有限公司
• 代理人: 王潔

1.一種基于復(fù)奇異譜分析的磁共振部分K數(shù)據(jù)圖像重建方法，其特征在于，所述的方法包括以下步驟：

(1)從磁共振成像掃描儀中預(yù)設(shè)的相位編碼范圍中采集部分K數(shù)據(jù)G(k)；

(2)根據(jù)該部分K數(shù)據(jù)信息進行模型參數(shù)估計；包括以下步驟：

(a)對該部分K數(shù)據(jù)G(k)的缺失部分補零，并根據(jù)以下公式得到缺失數(shù)據(jù)補零后的付里葉譜數(shù)據(jù)G_z(k)：

G_z(k)＝G(k)R_s-e(k)；

其中，G(k)為信號g(x)，x＝0，1，...，N-1的付里葉譜數(shù)據(jù)，k＝-N/2-1，...，N/2-1，為矩形函數(shù)，其中s為截斷上限頻率，e為截斷下限頻率，N為完整的K空間數(shù)據(jù)的相位編碼數(shù)，g(x)，x＝0，1，...，N-1為復(fù)數(shù)字信號，{b₁，b₂，...，b_Q}為g(x)上的Q個奇異點，為分別以{b₁，b₂，...，b_Q}為奇異點的Q個奇異函數(shù)，a₁，a₂，...，a_Q是該Q個奇異點{b₁，b₂，...，b_Q}上的復(fù)奇異值，G(k)＝DFT(g(x))，b_q＝0，1，...，Q，DFT(·)為離散付里葉變換算子；

(b)根據(jù)以下公式計算d_z(x)：

d_z(x)＝g_z(x)-g_z(x-1)；

其中，g_z(x)＝DFT^-1(G_z(k))，k＝-N/2-1，...，N/2-1，DFT^-1(·)表示離散付里葉反變換算子；

(c)將所得到的d_z(x)的模按照從大到小排序，并取前L個點作為預(yù)選奇異點集{b₁，b₂，...，b_L}，其中L為矩形函數(shù)R_s-e(k)的寬度，即L＝e-s，且已知頻譜為{G(k₁)，G(k₂)，...，G(k_L)}；

(d)根據(jù)以下公式構(gòu)造奇異譜方程：

(e)用偽逆矩陣法解出所述的奇異譜方程，得到一個最小誤差解，獲得L個復(fù)數(shù)奇異值{a₁，a₂，...，a_L}；

(f)將{a₁，a₂，...，a_L}作為模型參數(shù)估計的結(jié)果返回；

(3)根據(jù)模型參數(shù)估計的結(jié)果，利用復(fù)系數(shù)加權(quán)奇異函數(shù)的磁共振圖像的數(shù)學(xué)模型或者復(fù)奇異譜分析模型進行磁共振復(fù)圖像的重構(gòu)。

2.根據(jù)權(quán)利要求1所述的基于復(fù)奇異譜分析的磁共振部分K數(shù)據(jù)圖像重建方法，其特征在于，所述的系統(tǒng)預(yù)設(shè)的相位編碼范圍為-N/16～N/2-1，其中N為完整的K空間數(shù)據(jù)的相位編碼數(shù)。