[發明專利]磁共振成像方法、裝置及計算機存儲介質

有效

• 申請號: 202010782357.5
• 申請日: 2020-08-06
• 公開（公告）號: CN112014782B
• 公開（公告）日: 2023-07-28
• 發明（設計）人: 王海峰;梁棟;鄭海榮;劉新;蘇適;丘志浪
• 申請（專利權）人: 深圳先進技術研究院
• 主分類號: G01R33/565
• 分類號: G01R33/565;A61B5/055
• 代理公司: 深圳市威世博知識產權代理事務所(普通合伙) 44280
• 代理人: 黎堅怡

1.一種磁共振成像方法，其特征在于，所述方法包括：

獲取基于第一三維磁共振成像序列采集的目標對象的三維欠采樣數據；其中，所述第一三維磁共振成像序列在相位方向具有第一正弦梯度場，在選層方向具有第二正弦梯度場，所述第一正弦梯度場和所述第二正弦梯度場的持續時間在讀出方向施加的讀出梯度場的讀出平臺持續時間內，且所述第一正弦梯度場和所述第二正弦梯度場的0階矩均為0，所述第一正弦梯度場與所述第二正弦梯度場的相位相差π/2；

基于所述三維欠采樣數據或所述目標對象的二維映射數據計算得到三維點擴散函數；其中，所述二維映射數據與所述三維欠采樣數據的視野相同；

基于所述目標對象的三維低分辨率全采樣數據計算得到所述目標對象的靈敏度圖；

使用所述三維點擴散函數、所述靈敏度圖對所述三維欠采樣數據進行圖像重建，得到重建磁共振圖像；

其中，第一正弦梯度場公式表示如下：

第二正弦梯度場公式表示如下：

其中，t為時間，A為正弦梯度場幅度值，D_C為一個正弦周期的持續時間，D_R為讀出梯度場的讀出平臺持續時間，a為正奇數，b為正自然數，c為正自然數，d為正自然數，那么bD_C為第一持續時間，dD_C為第二持續時間，bD_C與cD_C均小于D_R，t＝0為讀出平臺的開始時間，此處表示原點。

2.根據權利要求1所述的方法，其特征在于，所述第一正弦梯度場與所述第二正弦梯度場的正弦頻率相同。

3.根據權利要求1所述的方法，其特征在于，所述二維映射數據包括基于第二三維磁共振成像序列采集到的所述目標對象的第一二維映射數據、基于第三三維磁共振成像序列采集到的所述目標對象的第二二維映射數據、基于第四三維成像序列采集到的所述目標對象的第三二維映射數據，以及基于第五三維成像序列采集到的所述目標對象的第四二維映射數據；其中，所述第二三維磁共振成像序列的相位方向包括所述第一正弦梯度場和第一梯度場，所述第三三維磁共振成像序列的相位方向包括所述第一梯度場，所述第四三維磁共振成像序列的選層方向包括所述第二正弦梯度場和所述第一梯度場，所述第五三維磁共振成像序列的選層方向包括所述第一梯度場；

所述基于所述目標對象的二維映射數據計算得到三維點擴散函數，包括：

將所述第一二維映射數據除以所述第二二維映射數據，得到相位方向的二維點擴散函數；

將所述第三二維映射數據除以所述第四二維映射數據，得到選層方向的二維點擴散函數；

基于所述相位方向的二維點擴散函數和所述選層方向的二維點擴散函數，得到所述三維點擴散函數。

4.根據權利要求1所述的方法，其特征在于，所述獲取所述目標對象的靈敏度圖，包括：

獲取基于第六三維磁共振成像序列采集的所述目標對象的低分辨率三維全采樣數據，所述低分辨率三維全采樣數據與所述三維欠采樣數據的所述視野相同；

基于所述低分辨率三維全采樣數據計算得到所述靈敏度圖。

5.根據權利要求4所述的方法，其特征在于，所述基于所述低分辨率三維全采樣數據計算得到所述靈敏度圖，包括：

使用特征值迭代自穩定并行成像重建技術對所述低分辨率三維全采樣數據進行計算，得到所述靈敏度圖。

6.根據權利要求4所述的方法，其特征在于，所述第六三維磁共振成像序列為三維梯度回波序列、快速自旋回波序列或三維平衡穩態自由旋進序列。

7.根據權利要求1所述的方法，其特征在于，所述三維欠采樣數據是使用數據欠采樣策略對所述目標對象基于所述第一三維磁共振成像序列產生的磁共振信號進行采集的，所述數據欠采樣策略為可控混疊的快速并行成像技術、并行成像中靈敏度編碼技術或整體自動校準部分并行采集技術。