[发明专利]时域核磁共振谱反演的方法

摘要

本发明涉及一种时域核磁共振谱反演的方法，首先，生成反演核并形成L1范数正则化问题先根据采样序列和采样参数，生成反演核，添加L1范数项进行惩罚，形成L1范数正则化问题；然后，数据拟合对L1范数正则化问题进行迭代求解，并得到反演谱。大幅提高了反演算法的执行效率与谱的分辨率，同时还具有很好的鲁棒性。

主权项

一种时域核磁共振谱反演的方法，其特征在于，具体包括如下步骤：1)生成反演核K并形成L1范数正则化问题：先根据采样序列和采样参数，生成反演核K，添加L1范数项进行惩罚，形成L1范数正则化问题；2)数据拟合：对L1范数正则化问题进行迭代求解，并得到反演谱；所述步骤1)中生成反演核K为单个二维矩阵，如采样序列为一维时域核磁共振的实验所得，那就是一维反演问题，生成的反演核K为单一的二维矩阵，一维反演包括横向弛豫时间T2谱反演、纵向弛豫时间T1谱反演、扩散系数D谱反演、内部梯度场G谱反演中任意一项；如采样序列为二维时域核磁共振的实验所得，存在两个任意一维反演核，对这两个反演核重排到一个二维矩阵中，作为最终的反演核K；所述步骤1)中形成L1范数正则化问题，为了得到反演谱时域谱s,确定采样参数之后，根据反演核K、采样数据m和正则化因子λ确定一个L1范数正则化问题：args>0min{F(s)≡||Ks-m||2+λ||s||1}]]>||s||1表示s的L1范数，其中m和s均为向量，s＞0表示向量s中的所有元素均不小于0。