[发明专利]一种动态脑功能磁共振成像方法及系统

权利要求书

1.一种动态脑功能磁共振成像方法，其特征在于，该方法包括：

在(k,t)空间采集高时间低空间分辨率的导航数据S^NAV与高空间低时间分辨率的动态图像数据S^IMG；

采用奇异值对所述导航数据S^NAV进行分解获得模型阶数L和时间基函数再结合所述图像数据S^IMG进行Tikhonov正则化约束求解出空间基函数c_l(k)；

利用所述模型阶数L、时间基函数及空间基函数c_l(k)对所述(k,t)空间进行高时间高空间覆盖密度的插值恢复，再进行傅里叶逆变换，从而获得高时间分辨率高空间分辨率的动态脑功能磁共振图像。

2.根据权利要求1所述的方法，其特征在于，所述采用奇异值对所述导航数据S^NAV进行分解获得模型阶数L和时间基函数的公式包括：

SNAV=Σl=1min(M,N)λlulvl;]]>

其中，λ_l是导航数据S^NAV的奇异值按降序排列，u_l和v_l分别为左奇异向量和右奇异向量，M表示(k,t)空间中沿时间轴时间点总数，N表示(k,t)空间中每一个时刻的k空间点数目；

所述时间基函数由所述左奇异向量ul计算得到，所述模型阶数L由所述导航数据S^NAV的最佳估计得到。

3.根据权利要求1或2所述的方法，其特征在于，所述结合所述图像数据S^IMG进行Tikhonov正则化约束求解出空间基函数c_l(k)包括：

其中，α表示正则约束项参数，ε表示噪声因子，||c_l(k)||₂表示空间基函数c_l(k)的2范数，t_m为特定时刻；

利用非线性迭代方法求解上述公式，获得空间基函数c_l(k)。

4.根据权利要求1所述的方法，其特征在于，所述利用所述模型阶数L、时间基函数及空间基函数c_l(k)对所述(k,t)空间进行高时间高空间覆盖密度的插值恢复包括：

利用下述公式，并结合模型阶数L、时间基函数及空间基函数c_l(k)对所述(k,t)空间进行高时间高空间覆盖密度的插值恢复，其公式为：

5.根据权利要求1所述的方法，其特征在于，采集所述导航数据S^NAV与图像数据S^IMG包括：

基于预设的相位编码步数和回波链长度ETL，产生对应的激发次数；

其中，每次射频脉冲激发后填充K空间相位编码线条数，且在同一射频周期内将导航数据S^NAV与图像数据S^IMG进行同步采集。