一种高维核磁共振时域信号补全方法
【技术领域】
[0001] 本发明涉及核磁共振高维谱信号处理，特别是涉及一种高维核磁共振时域信号补 全方法。
【背景技术】
[0002] 核磁共振波谱（NMR)，简称为核磁共振谱，已经成为物理学、化学以及生命科学等 诸多学科研宄物质组成的成分、组织形态及其变化、脑功能、分子结构和动力学强有力的手 段。特别在化学结构分析领域，核磁共振谱有着重要的应用。
[0003] 随着检测对象和设备的不断发展，高维NMR(多2维）在应用领域发挥着越来越 重要的作用。例如，2维、3维和4维NMR广泛用于解析复杂蛋白结构（D. Sakakibara，A. Sasaki, T. Ikeyaj J. Hamatsuj T. Hanashimaj M. Mishimaj et al.，''Protein structure determination in living cells by in-cell NMR spectroscopy，"Nature，vol. 458， pp. 102-105, 2009)。但数据采集时间随着维数的增加而指数级增长。2维谱到4维谱的采 集时间从几分钟不断上升到几十天（M. Mobli and J.C. Hoch, 〃Nonuniform sampling and non-Fourier signal processing methods in multidimensional NMR, ''Progress in Nuclear Magnetic Resonance Spectroscopy, vol. 83, pp. 21-41，2014)。冗长的高维蛋白 质NMR实验时间，不但使得实验必须耗费大量的谱仪机时，而且提高了不稳定蛋白质样品 的实验难度，从而限制了高维NMR技术的在蛋白质研宄中的应用。为了缩短高维NMR实验 时间，一个常见的做法是采用非均匀采样的方式来减少采样的点数，但同时也造成了信号 的丢失，需要信号补全。
[0004] 在NMR实验中，采样得到的时域信号不可避免会受到噪声的影响。由于核磁共振 信号末尾处的信号强度小造成这些数据点的信噪比低，截掉这些数据点有利于提高核磁共 振谱的信噪比，但这也会造成信号丢失，需要信号补全。
[0005] 频谱宽度和分辨率是评价核磁共振谱质量的两个重要指标。当采样所得的时间信 号未达到预期的频谱宽度和分辨率时，需要对时域信号进行信号插值。信号插值也是一种 信号补全。
[0006] 因此信号补全技术在NMR实验中应用很广泛，设计有效的信号补全技术意义重 大。
【发明内容】
[0007] 本发明的目的在于针对高维核磁共振波谱实验存在的上述问题，提供一种高维核 磁共振时域信号补全方法。
[0008] 本发明包括以下步骤：
[0009] 1)设计模板确定信号丢失位置：根据实验中的采样方式和给定的核磁共振谱频 谱宽度和分辨率，确定需要补全的信号位置并设计成模板，模板大小与预期完整的信号大 小相同；模板由数字0和1构成，数字0表示该位置对应的信号需要补全，数字1表示该位 置对应的信号未丢失；
[0010] 2)构建信号重建模型：在m维核磁共振实验中，预期完整的信号X是一个m维矩 阵，其中第一维用J1表示，这一维的信号长度是J1，第二维用j2表示，这一维的信号长度是 J2，以此类推；第一维是直接维，其余是间接维，当给定除jk以外的其它维度时，可以得到一 个一维数据，称为一维阵列，用表示；根据核磁共振信号具有指数函数形式的 特点，信号中沿任意维度的一个一维阵列排列成一个汉克尔矩阵，该矩阵的秩比行数小很 多，这种矩阵称为低秩矩阵；构造以下的模型来求解完整的核磁共振信号X:
【主权项】
1. 一种高维核磁共振时域信号补全方法，其特征在于包括w下步骤： 1) 设计模板确定信号丢失位置：根据实验中的采样方式和给定的核磁共振谱频谱宽 度和分辨率，确定需要补全的信号位置并设计成模板，模板大小与预期完整的信号大小相 同；模板由数字0和1构成，数字0表示该位置对应的信号需要补全，数字1表示该位置对 应的倍号未丢失； 2) 构建信号重建模型；在m维核磁共振实验中，预期完整的信号X是一个m维矩阵，其 中第一维用jl表示，该一维的信号长度是J1，第二维用J'2表示，该一维的信号长度是J2,W 此类推；第一维是直接维，其余是间接维，当给定除jkW外的其它维度时，可W得到一个一 维数据，称为一维阵列，用X,,,:...心jw..,，。表示；根据核磁共振信号具有指数函数形式的特点， 信号中沿任意维度的一个一维阵列排列成一个汉克尔矩阵，该矩阵的秩比行数小很多，该 种矩阵称为低秩矩阵；构造W下的模型来求解完整的核磁共振信号X:
在公式（1)中，算子A根据步骤1)中设计的模板确定X中信号丢失的位置，并把未 丢失的信号排列成一维向量巧表示采集到并且没有丢失的数据；1_F-AX|g表示（y-AX) 中的元素的平方和；KX,,,，...A_,;,w..,,。,表示将一维阵列Xw;...心,:ky...^,排列成一个汉克尔矩阵 ?化..山,加...,。,；13化..山也,..,1表示矩阵11,沁.,,,_,山,.,,。的核范数，即对矩阵的奇异值求和；参 数e与实际存在的噪声大小有关，e> 0 ; 求解公式（1)中的模型时，可W把约束优化问题转化为W下的非约束优化问题，
(4) 其中A是正则化参数（A〉〇)，用于权衡
两项的重要性； 3) 求解最优化问题：求解最优化问题即公式（4)，即可得到完整的m维核磁共振信号。
2. 如权利要求1所述一种高维核磁共振时域信号补全方法，其特征在于在步骤3)中， 所述求解最优化问题即公式（4)是一个核范数最优化问题，采用数值计算方法求解，得到 完整的m维核磁共振时间信号X;然后对X进行傅立叶变换得到核磁共振谱。
3. 如权利要求2所述一种高维核磁共振时域信号补全方法，其特征在于所述数值计算 方法包括但不限于奇异值阔值化法或交替方向乘子法。
【专利摘要】一种高维核磁共振时域信号补全方法，涉及核磁共振高维谱信号处理。先根据采集得到的数据和给定的核磁共振谱频谱宽度和分辨率，确定需要补全的时域信号的位置并设计模板；再利用提出的高维核磁共振时域信号补全方法来构建重建模型，然后通过最优化算法求解出完整的高维核磁共振时域信号；最后对补全后的时域信号做傅立叶变换得到核磁共振谱。实现了在高维核磁共振实验中对有丢失的核磁共振时域信号进行信号补全，进而得到完整的核磁共振时域信号。由于利用了高维核磁共振信号自身的特征，可以补全任意高维核磁共振时域信号。可以达到降低采样时间、提高信噪比、达到给定核磁共振谱频谱宽度和分辨率的目的。
【IPC分类】G01R33-46
【公开号】CN104793159
【申请号】CN201510235929
【发明人】屈小波, 应佳熙, 郭迪, 陈忠
【申请人】厦门大学
【公开日】2015年7月22日
【申请日】2015年5月11日