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

無效

• 申請號: 200710039881.8
• 申請日: 2007-04-24
• 公開（公告）號: CN101051075A
• 公開（公告）日: 2007-10-10
• 發明（設計）人: 駱建華;朱躍敏
• 申請（專利權）人: 駱建華
• 主分類號: G01R33/56
• 分類號: G01R33/56
• 代理公司: 上海智信專利代理有限公司
• 代理人: 王潔

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

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

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

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

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

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

(b)根據以下公式計算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個點作為預選奇異點集{b₁，b₂，...，b_L}，其中L為矩形函數R_s-e(k)的寬度，即L＝e-s，且已知頻譜為{G(k₁)，G(k₂)，...，G(k_L)}；

(d)根據以下公式構造奇異譜方程：

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

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

(3)根據模型參數估計的結果，利用復系數加權奇異函數的磁共振圖像的數學模型或者復奇異譜分析模型進行磁共振復圖像的重構。

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