1.一种核磁共振成像方法,其特征在于,包括:
在同一扫描部位采集多组不同对比度图像的数据;
在k空间,分别对每组对比度图像的数据在相位方向外围进行欠采样且在自校准信号线处进行全采样,获得采样数据;
将线圈灵敏度信息作为多组不同对比度图像中每组图像的共享信息,将所述线圈灵敏度信息和多组不同对比度图像的待重建信息作为核磁共振成像的未知数,将所述线圈灵敏度信息和多组不同对比度图像的待重建信息形成重建数据矩阵;
将所述采样数据和线圈灵敏度平滑信息作为输入数据,将所述重建数据矩阵作为输出,形成用于求解所述线圈灵敏度信息和多组不同对比度图像的待重建信息的目标函数;
采用非线性迭代算法求解所述目标函数,同时重建多组不同对比度图像,计算出所述线圈灵敏度信息;
通过以下公式将所述线圈灵敏度信息和多组不同对比度图像的待重建信息形成重建数据矩阵:其中,F(X)是重建数据矩阵,X表示线圈灵敏度信息和多组不同对比度图像的待重建信息组成的未知数,P是k空间数据欠采样模型,F是傅里叶变换,.*是点乘运算,C1...CN是一组经线圈灵敏度平滑信息处理过的线圈灵敏度信息,ρ1是第一组图像的待重建信息,ρ2是第二组图像的待重建信息,ρn是第n组图像的待重建信息,n是不同对比度图像的组数;
利用所述线圈灵敏度平滑信息通过以下公式处理所述线圈灵敏度信息:
Ci=F-1(W·*Si)
其中,Ci是k空间的线圈灵敏度信息经所述线圈灵敏度平滑信息转换成的图像域的线圈灵敏度信息,Si是线圈灵敏度信息在k空间的原始数据,i表示线圈灵敏度信息中的第i个数值, W为索伯列夫式平滑矩阵,K是k空间各数据点到k空间中心的距离,I是索伯列夫指数,F-1是傅里叶逆变换;
通过以下公式表示所述目标函数:
其中,Y表示所述目标函数,d表示所述采样数据, 表示第一组图像的采样数据, 表示第二组图像的采样数据, 表示第n组图像的采样数据。
2.一种核磁共振成像装置,其特征在于,包括:
数据采集模块,用于在同一扫描部位采集多组不同对比度图像的数据;
数据采样模块,用于在k空间,分别对每组对比度图像的数据在相位方向外围进行欠采样且在自校准信号线处进行全采样,获得采样数据;
重建数据生成模块,用于将线圈灵敏度信息作为多组不同对比度图像中每组图像的共享信息,将所述线圈灵敏度信息和多组不同对比度图像的待重建信息作为核磁共振成像的未知数,将所述线圈灵敏度信息和多组不同对比度图像的待重建信息形成重建数据矩阵;
目标函数生成模块,用于将所述采样数据和线圈灵敏度平滑信息作为输入数据,将所述重建数据矩阵作为输出,形成用于求解所述线圈灵敏度信息和多组不同对比度图像的待重建信息的目标函数;
求解模块,用于采用非线性迭代算法求解所述目标函数,同时重建多组不同对比度图像,计算出所述线圈灵敏度信息;
所述重建数据生成模块通过以下公式将所述线圈灵敏度信息和多组不同对比度图像的待重建信息形成重建数据矩阵:其中,F(X)是重建数据矩阵,X表示线圈灵敏度信息和多组不同对比度图像的待重建信息组成的未知数,P是k空间数据欠采样模型,F是傅里叶变换,.*是点乘运算,C1...CN是一组经线圈灵敏度平滑信息处理过的线圈灵敏度信息,ρ1是第一组图像的待重建信息,ρ2是第二组图像的待重建信息,ρn是第n组图像的待重建信息,n是不同对比度图像的组数;
所述目标函数生成模块利用所述线圈灵敏度平滑信息通过以下公式处理所述线圈灵敏度信息:Ci=F-1(W·*Si)
其中,Ci是k空间的线圈灵敏度信息经所述线圈灵敏度平滑信息转换成的图像域的线圈灵敏度信息,Si是线圈灵敏度信息在k空间的原始数据,i表示线圈灵敏度信息中的第i个数值, W为索伯列夫式平滑矩阵,K是k空间各数据点到k空间中心的距离,I是索伯列夫指数,F-1是傅里叶逆变换;
所述目标函数生成模块通过以下公式表示所述目标函数:
其中,Y表示所述目标函数,d表示所述采样数据, 表示第一组图像的采样数据, 表示第二组图像的采样数据, 表示第n组图像的采样数据。