一种核磁共振多维谱的伪峰抑制方法

摘要

本发明公开了一种核磁共振多维谱的伪峰抑制方法，该方法首先进行直接维均匀采样和间接维随机采样，得到时域数据；接着时域数据依次进行直接维快速傅里叶变换和网格傅里叶变换，得到伪峰校正前的核磁共振谱图；然后对进行直接维快速傅里叶变换后的时域数据进行子集抽取，并进行网格傅里叶变换，得到子集核磁共振谱图；最后将随机子集核磁共振谱图与伪峰校正前的核磁共振谱图进行比较，得到伪峰校正前的核磁共振谱图中每点谱峰强度的相对标准偏差，并利用权重函数对伪峰校正前的核磁共振谱图进行校正，得到伪峰抑制后的核磁共振谱图。本发明的有益效果是：1、数据处理速度快；2、有效地区别信号峰与伪峰，并对伪峰有显著的抑制作用。

主权项

一种核磁共振多维谱的伪峰抑制方法，其特征在于，包括以下步骤：步骤1、对磁化矢量进动时引起的电平振荡进行直接维均匀采样和间接维随机采样，得到时域数据；随机采样具体的步骤为：步骤1.1、确定随机采样所需要的采样时间点坐标，这些采样时间点随机分布在由第一间接维和第二间接维构成的二维平面上，其中在第一间接维和第二间接维上的最大值即为第一间接维和第二间接维上的最大间接维演化时间；步骤1.2、根据步骤1.1确定的采样时间点坐标对HNCO脉冲序列的间接维演化时间点坐标修改；步骤1.3、根据步骤1.2修改的HNCO脉冲序列间接维演化时间点坐标进行时域数据采集；步骤2、对步骤1中的时域数据进行直接维快速傅里叶变换处理得到频域‑时域混合数据；步骤3、对步骤2中得到的频域‑时域混合数据进行网格傅里叶变换GFFT处理，得到伪峰抑制处理所需的核磁共振谱图；步骤4、对步骤2中得到的频域‑时域混合数据，在间接维上进行采样点的随机抽取，随机抽取的采样点组成N个随机子集，其中N为大于等于3的正整数；步骤5、对步骤4中得到的N个随机子集进行网格傅里叶变换GFFT处理，得到N个随机子集核磁共振谱图；步骤6、计算步骤5得到的N个随机子集核磁共振谱图中每点的谱峰强度与步骤3中的核磁共振谱图中相对应点的谱峰强度的相对标准偏差；步骤7、根据步骤6得到N个随机子集核磁共振谱图中每点谱峰强度的相对标准偏差，利用权重函数对步骤3得到的核磁共振谱图进行权重校正，即将N个随机子集核磁共振谱图中每点谱峰强度的相对标准偏差作为变量输入权重函数中，所得结果再与步骤3得到的核磁共振谱图中相应的该点谱峰强度相乘，以此得到伪峰抑制后的核磁共振谱图。