[發明專利]一種對擴散張量磁共振圖像進行恢復處理的方法

無效

• 申請號: 200810033965.5
• 申請日: 2008-02-28
• 公開（公告）號: CN101236241A
• 公開（公告）日: 2008-08-06
• 發明（設計）人: 張相芬;田蔚風;葉宏;李大志
• 申請（專利權）人: 上海師范大學
• 主分類號: G01R33/565
• 分類號: G01R33/565;A61B5/055
• 代理公司: 上海伯瑞杰知識產權代理有限公司
• 代理人: 吳澤群

1.一種對擴散張量磁共振圖像進行恢復處理的方法，其特征是，使用計算機，包括以下步驟：

第一步，獲取研究對象的圖像信息；

所述的獲取，是在磁共振成像系統上進行；

所述的圖像信息，包括未施加梯度脈沖獲得的DTI圖像信息、施加加權梯度脈沖獲得的DTI圖像信息即每層擴散加權的DWI圖像信息；

所述的每層擴散加權DWI圖像信息，是在每層中以6個以上的不共線方向對研究對象進行加權成像從而獲得該層的6組以上擴散加權DWI圖像信息；

第二步，把每層DWI圖像信息轉化為向量圖像信息，其表示式為：I＝(I₁(x，y)，...，I_m(x，y))^T；

其中I_n，n＝1，2，...，m表示第n加權脈沖加權圖像，(x，y)為二維圖像中的橫、縱坐標，T表示矩陣轉置；

第三步，按自然層數為序或者按另行編排的先后層數為序進行排列，并令最前層為第1層、最后層為第P層；

令層數計數器為1；

設置步長為dt，設置在每一層的循環次數為Q；

第四步，令循環次數計數器為1；

第五步，對當前層數、當前循環次數的H和J進行計算；

所述的H和J是兩個相互獨立的仿射不變算子；

H=fxxTfyy-fxyTfxy]]>

=Σi=1m(fixxfiyy-fixyfixy)]]>

=Σi=1mHi]]>

上式中，f_jk＝[f_1jk，...，f_mjk]^T?j，k＝x，y

T表示矩陣轉置；當m＝1時，f_jk是標量；

對任意尺度的向量，J由下式構造：

J=fyT,-fxTdiag(fxx),diag(fxy)diag(fxy),diag(fyy)fy-fx]]>

=Σi=1mJi]]>

由上式可得，

Ji=fiyT,-fixTdiag(fixx),diag(fixy)diag(fixy),diag(fiyy)fiy-fix]]>

=fiy,-fixfixx,fixyfixy,fiyyfiy-fix]]>

第六步，計算仿射不變梯度G_aff和歐幾里得曲率κ；

Gaff=HJ2+1]]>

κ＝J?|▽f|³

上式中，f_i(x，y，0)＝I_i(x，y)(-1，...，m)

|▿f|2=Σi=1m|▿fi|2]]>

“▽”表示梯度；

第七步，確定向量圖像仿射不變各向異性平滑的數值

f_it＝G_affκ^1/3|▽f_i|；

第八步，根據以下公式計算得到恢復后的圖像信息并存儲該信息；

f(x，y，t)＝f(x，y，t-1)+dt×f_t-1(t＝1，...，N)；

上式中，f(x，y，t)＝(f₁(x，y，t)，...，f_m(x，y，t))^T，f_t＝(f_1t，...，f_mt)^T；

第九步，檢查循環次數計數器的計數值是否為Q，如果為Q則跳過第十步、直接進行第十一步，如果不為Q則進行下一步；

第十步，當第九步檢查循環次數計數器的計數值不為Q時，循環次數計數器加1后返回第五步；

第十一步，檢查層數計數器的計數值是否為P，如果為P則跳過第十二步、直接進行第十三步，如果不為P則進行下一步；

第十二步，當第十一步檢查層數計數器的計數值不為P時，層數計數器加1后返回第四步；

第十三步，存儲圖像信息；進入結束。

2.根據權利要求1所述的一種對擴散張量磁共振圖像進行恢復處理的方法，其特征是，所述的設置在每一層的循環次數為20。