[發(fā)明專利]基于迭代投影的MRI圖像重構(gòu)方法

無(wú)效

• 申請(qǐng)?zhí)枺?/td>: 201210007708.0
• 申請(qǐng)日: 2012-01-02
• 公開(kāi)（公告）號(hào): CN102609905A
• 公開(kāi)（公告）日: 2012-07-25
• 發(fā)明（設(shè)計(jì)）人: 張小華;焦李成;王尚禮;王爽;鐘樺;張揚(yáng);田小林
• 申請(qǐng)（專利權(quán)）人: 西安電子科技大學(xué)
• 主分類號(hào): G06T5/00
• 分類號(hào): G06T5/00
• 代理公司: 陜西電子工業(yè)專利中心 61205
• 代理人: 王品華;朱紅星

1.一種基于迭代投影的核磁共振MRI圖像重構(gòu)方法，包括如下步驟：

(1)對(duì)MRI圖像x進(jìn)行傅里葉變換，得到變換系數(shù)：x_f＝DFT(x)，用變密度欠采樣觀測(cè)矩陣Φ對(duì)圖像x的傅里葉變換系數(shù)x_f進(jìn)行觀測(cè)采樣，獲得傅里葉域的觀測(cè)值y＝Φx_f，對(duì)觀測(cè)值y進(jìn)行傅里葉反變換得到MRI圖像x的初始值：x⁰＝DFT^-1(y)；

(2)結(jié)合MRI圖像局部光滑特性和小波域稀疏的先驗(yàn)知識(shí)，將MRI圖像x的初始值x⁰用全變差方法進(jìn)行濾波，然后用小波域雙變量閾值方法對(duì)該濾波結(jié)果x^TV進(jìn)行優(yōu)化，得到MRI圖像x的優(yōu)化結(jié)果x′；

(3)對(duì)優(yōu)化結(jié)果x′進(jìn)行投影，得到投影結(jié)果x^*＝x′+Φ^T(y-Φx′)，并判斷x^*是否滿足迭代終止條件：|D⁽ⁱ⁺¹⁾-D⁽ⁱ⁾|＜10^-5，其中N是圖像的像素個(gè)數(shù)，i表示第i次迭代，若滿足條件，將投影結(jié)果x^*作為MRI圖像x的重構(gòu)結(jié)果并輸出，否則返回步驟(2)繼續(xù)執(zhí)行。

2.根據(jù)權(quán)利要求1所述的基于迭代投影的MRI圖像重構(gòu)方法，其中步驟(2)所述的用小波域雙變量閾值方法對(duì)全變差濾波結(jié)果x^TV進(jìn)行優(yōu)化，按如下步驟進(jìn)行：

2.1)選取雙變量閾值為：Threshold(f,λ)=(f2+fP2-λ3σ(i)σξ)f2+fP2·f,]]>λ是控制算法收斂性的常數(shù)，取λ＝20，f＝Ψx^TV是x^TV的小波變換系數(shù)，Ψ表示小波變換，i表示低i次迭代，f_p是f的父系數(shù)，如果則取閾值為0，是系數(shù)ξ3×3鄰域的邊緣方差估計(jì)，σ⁽ⁱ⁾是中值估計(jì)函數(shù)，

2.2)根據(jù)確定的雙變量閾值Threshold(f，λ)選擇比閾值大的小波變換系數(shù)f_s保留，其余的系數(shù)置0，再對(duì)保留系數(shù)f_s進(jìn)行小波反變換就得到對(duì)x^TV的優(yōu)化結(jié)果x′＝Ψ^-1f_s。