[发明专利]一种基于LSQR的低场核磁共振二维谱反演算法

摘要

本发明涉及一种基于LSQR(LeastSquaresQRDecomposition，最小二乘QR分解)的低场核磁共振二维谱反演算法，属于核磁共振信号处理领域。本发明的优点表现为：本发明与目前国内外文献报道的其他低场核磁共振反演算法相比，其借助“L”曲线自适应选择合适的迭代次数，得到的解更能体现真实的谱分布；通过迭代修正来实现非负限制，计算精度高，鲁棒性好，在不同信噪比数据情况下，能够得到稳定的反演结果。

主权项

一种基于LSQR的低场核磁共振二维谱反演算法，其特征在于，包括如下步骤：a.读取低场核磁共振设备采集得到的CPMG原始数据文件；b.对原始数据进行预处理操作得到反演核心矩阵K和信号幅值m；c.根据核心矩阵K的线性关系预设一个LSQR迭代的最大迭代次数；d.进行LSQR迭代，得到每次迭代后的反演结果和残差的范数；e.根据不同迭代次数对应的反演结果的范数和残差的范数，绘出“L形”曲线；f.计算“L形”曲线曲率，定义拐角位置为LSQR迭代的最优解位置；g.使用迭代修正的方式进行约束，抑制负峰的产生，输出二维谱。