[發明專利]一種磁共振非笛卡爾采樣的快速重建方法

無效

• 申請號: 201310311478.1
• 申請日: 2013-07-23
• 公開（公告）號: CN103399285A
• 公開（公告）日: 2013-11-20
• 發明（設計）人: 郭紅宇
• 申請（專利權）人: 沈陽工業大學
• 主分類號: G01R33/56
• 分類號: G01R33/56;A61B5/055
• 代理公司: 沈陽智龍專利事務所(普通合伙) 21115
• 代理人: 宋鐵軍;周楠

1.一種磁共振非笛卡爾采樣的快速重建方法，其特征在于：該方法按下列步驟進行：

a、密度補償函數的快速計算：

用Jackson密度補償函數的后補償網格重建的過程表示為式（1）：

Mswcs(kx,ky)=[Ms(kx,ky)*C(kx,ky)]·III(kx,ky)[S(kx,ky)*C(kx,ky)]·III(kx,ky)---(1)]]>

由式（1）得到，后補償方法中的密度補償函數為式（2），

W(kx,ky)=1[S(kx,ky)*C(kx,ky)]·III(kx,ky)---(2)]]>

從式（2）的分子部分看出，在后補償網格重建中，密度補償實際上是對采樣值全部等于1的采樣數據進行網格化重建計算；在網格重建中對每個采樣點計算它對于處在卷積窗內的網格點的貢獻，也就是說此時的密度補償函數只與采集點的分布即采樣軌跡有關；

b、根據對稱性快速計算采樣密度函數：

在非笛卡爾采樣中，很多采樣軌跡的分布是對稱分布的；第一象限內的網格點(k_x,k_y)的密度補償函數W(k_x,k_y)與其他象限內與它對稱的網格點的采樣密度函數相同，即如式（3）：

W(k_x,k_y)＝W(-k_x,k_y)＝W(-k_x,-k_y)＝W(k_x,-k_y)???（3），

在計算密度補償函數時，只需要計算位于直角坐標系第一象限內的采樣點的密度補償函數，考慮卷積窗寬度的影響，實際實現過程中需要計算密度補償函數的采樣點為第一象限的點加上半個卷積窗寬；

c、進行網格化計算：

利用已經計算完成的密度補償函數，對采樣數據按照式（4）進行網格化重建，

MSWCS(kx,ky)=[MS(kx,ky)*C(kx,ky)]·III(kx,ky)W(kx,ky)---(4).]]>