CN105928965B - 一种核磁共振谱采样截断伪峰的抑制方法 - Google Patents

一种核磁共振谱采样截断伪峰的抑制方法 Download PDF

Info

Publication number
CN105928965B
CN105928965B CN201610310469.4A CN201610310469A CN105928965B CN 105928965 B CN105928965 B CN 105928965B CN 201610310469 A CN201610310469 A CN 201610310469A CN 105928965 B CN105928965 B CN 105928965B
Authority
CN
China
Prior art keywords
nmr spectra
time
peak
domain signal
threshold value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610310469.4A
Other languages
English (en)
Other versions
CN105928965A (zh
Inventor
蒋滨
孙鹏
罗凡
张许
杨运煌
李从刚
周欣
刘买利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Precision Measurement Science and Technology Innovation of CAS
Original Assignee
Wuhan Institute of Physics and Mathematics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Institute of Physics and Mathematics of CAS filed Critical Wuhan Institute of Physics and Mathematics of CAS
Priority to CN201610310469.4A priority Critical patent/CN105928965B/zh
Publication of CN105928965A publication Critical patent/CN105928965A/zh
Application granted granted Critical
Publication of CN105928965B publication Critical patent/CN105928965B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/087—Structure determination of a chemical compound, e.g. of a biomolecule such as a protein

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明公开了一种核磁共振谱采样截断伪峰的抑制方法，设定第一百分比阈值和第二百分比阈值；对时域信号a进行填零处理并进行傅立叶变换后获得NMR谱图A；根据第一比较阈值将NMR谱图A分为NMR谱图B和NMR谱图C；根据NMR谱图B和第二比较阈值判断是否进入循环迭代得到各个循环迭代中的NMR谱图B和NMR谱图C；各个循环迭代中获得的NMR谱图B以及最后一次循环迭代中获得的NMR谱图C叠加，得到NMR谱图D；对NMR谱图D作窗口平滑处理修正，得到最终NMR谱图E。本发明可以减少间接维采样点数，提高实验速度，避免采样点数不足带来的采样截断伪峰，避免常规窗函数处理方法带来的谱峰展宽和重叠，提高谱图分辨率。

Description

一种核磁共振谱采样截断伪峰的抑制方法
技术领域
本发明涉及核磁共振分析领域，具体涉及一种核磁共振谱采样截断伪峰的抑制方法，用于对化学结构的核磁共振分析。
背景技术
一、核磁共振技术简介
核自旋量子数不为零的原子核在磁场中会产生磁矩，当这些磁矩的方向与静磁场夹角不为零时，就会绕静磁场的方向产生进动。核磁共振(Nuclear Magnetic Resonance,NMR)就是利用射频脉冲使静磁场中的样品的磁化矢量与静磁场的方向偏离，并产生进动；磁化矢量进动时，其磁力线会切割样品旁的检测线圈，通过电磁感应在线圈中产生电平振荡，对电平振荡进行采样得到时域信号。对此时域信号作傅里叶变换，即可得到一维NMR谱图。图1为单脉冲一维实验的脉冲序列图。
然而，当样品组分复杂、或者分子量较大时，一维NMR谱图中有严重的谱峰重叠，这将严重地影响对谱峰的归属和定量，如图2为蛋白质ubiquitin的一维氢谱。
为了解决一维谱中严重的谱峰重叠问题，出现了多维谱技术，如图3所示。以COSY实验的脉冲序列为例(图4)，此脉冲序列中有两个演化时间，t1和t2，分别称为间接维和直接维演化时间。对于直接维演化时间t2，其实就是样品的磁化矢量在检测线圈中进动的时间，因此，只有在直接维演化时间t2内，探头才可以对样品进行直接采样，这也是直接维的由来。若第一次检测时，t1取值τ0，并采得信号s(τ0，t2)；在第二次采样中，将t1设定为τ0+Δt1，并采得信号s(τ0+Δt1，t2)；在第n次采样中，将t1设定为τ0+nΔt1，并采得信号s(τ0+nΔt1，t2)，以此类推。因此，如果需要在间接维上采集n个数据点，需要将实验重复n次(在实际中，由于其它因素的存在，通常需要重复次数是n的倍数)。
图3为蛋白质Ubiquitin的COSY谱图中的部分示意图。由于谱峰分布由一维扩展到二维，谱峰重叠现象得以解决，对谱峰的归属才得以实现。
图4为COSY实验的脉冲序列。其中，t1和t2分别为间接维和直接维的演化时间。
一维NMR实验采得的时域数据是一个数列，数列中的每个元素是不同时间的探头检测到的信号强度，因此，对一维实验作数据处理，仅需要对此数列作一次傅里叶变换即可得到频域谱图。而二维NMR实验采得的数据是一个二维数组(矩阵)，对其的处理步骤通常是先对矩阵中的每行(直接维)分别作傅里叶变换，然后再对矩阵中的每列(间接维)作傅里叶变换，才能得到二维NMR谱图。对于三维实验可以依次类推：其实验数据是个三维数组，需要沿三个维度各作一次傅里叶变换，共计三次傅里叶变换。
二、核磁共振谱图中的采样截断伪峰
在核磁共振谱图中，常常会出现采样截断伪峰。其产生的原因是：当采样结束时，信号还没衰减完毕，因此采样信号就相当于实际全部信号与一个矩形函数的乘积。当将采样信号傅立叶变换后，时域的矩形函数就会给信号峰带来一个振荡的卷积，使得谱图中出现振荡状的伪峰。
图5为采样截断伪峰的由来示意图。当时域信号衰减完毕时结束采样(A)，则谱图正常(B)；当时域信号还未衰减完毕就结束采样(即采样截断，C)，则会在谱图中引入振荡状的伪峰(D)。
一般而言，采样截断伪峰较多出现在多维谱的间接维，主要原因是为了节约实验时间，间接维的采样点通常被限制在较少数目，使得间接维采样时间不够，从而导致采样截断。
三、目前对核磁共振谱图中采样截断伪峰的常规处理方法
目前，核磁共振波谱学中一般利用窗函数来抑制采样截断伪峰。既然采样截断伪峰源于时域信号未衰减完毕，那么将一个衰减函数与采样信号相乘，强行使采样信号衰减完毕，那么自然就消除了采样截断伪峰。但是，窗函数方法也存在缺点，就是会导致谱峰增宽，降低谱图分辨率。如附图6所示，窗函数消除了伪峰，但展宽了信号峰，使得靠近的信号峰重叠在一起(图6C)，不利于识别。用于消除截断伪峰的常见窗函数，包括指数函数和余弦函数等。
图6为截断采样数据的三种方法处理结果对比示意图。(A)完全采样所得的谱图；(B)对采样截断信号，经FFT所得的谱图，其中采样截断造成的伪峰很明显，左侧强信号峰旁的弱信号被伪峰掩盖；(C)，利用常规的窗函数方法，可消除截断伪峰，但是的谱峰变宽，相互靠近的信号会重叠在一起。(D)利用本问所提出的新方法(TASI)在消除截断伪峰的同时，避免了谱峰展宽带来的重叠。
发明内容
本发明的目的是针对现有技术存在的上述问题，提供一种核磁共振谱采样截断伪峰的抑制方法，该方法对核磁共振谱图中的采样截断伪峰具有较好的抑制效果，且避免了其它方法常导致的谱峰展宽。
本发明的上述目的通过以下技术方案实现：
以一维谱为例，对于多维谱，其处理步骤与此相同。
步骤1：设定第一百分比阈值p和第二百分比阈值q，第一百分比阈值p的取值范围为95％～98％，第二百分比阈值q的取值范围为0.005％～0.015％。
步骤2：获得的时域信号a，时域信号a的采样数据点点数为N1，对时域信号a进行填零处理，填零后的时域信号a的采样数据点点数为N2，然后对进填零后的时域信号a进行傅立叶变换后获得NMR谱图A。
步骤3：求得NMR谱图A中最高信号峰的强度值I，以最高信号峰的强度值I与第一百分比阈值p的乘积pI为第一比较阈值，将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值分别减去第一比较阈值，并将NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值设定为0，得到NMR谱图B；将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值设定为第一比较阈值，且NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值不变，得到NMR谱图C。
步骤4、将步骤2获得的NMR谱图A中最高信号峰的强度值I与第二百分比阈值q的乘积qI为第二比较阈值，
若NMR谱图B中最高信号峰的强度值小于第二比较阈值，则进入步骤7；
若NMR谱图B中最高信号峰的强度值大于等于第二比较阈值，则进入步骤5。
步骤5：将NMR谱图B进行逆傅立叶变换得到时域信号b，并对时域信号b作截断处理，即将N2个采样数据点的时域信号b只保留前面的N1个采样数据点。
步骤6：将时域信号a减去截断处理后的时域信号b获得新的时域信号a，将新的时域信号a作为步骤2中的时域信号a，循环迭代进入步骤2。
步骤7：将各个循环迭代中获得的NMR谱图B，以及最后一次循环迭代中获得的NMR谱图C相叠加，即得到NMR谱图D。
步骤8：对NMR谱图D作窗口平滑处理修正，得到最终NMR谱图E。即为NMR谱图D中每个数据点设定一个点数窗口，点数窗口一般为5～10个点，这个点数窗口中所有数据点的平均值，作为该数据点修正后的强度值。经此步骤，得到最终NMR谱图E。
本发明是基于以下原理：
采样截断伪峰，来源于不完全采样。令采样时样品磁化强度演化的时域信号为f(t)，则将其经傅立叶变换得到核磁共振谱图然而实际采样获得的时域信号是原始时域信号与采样函数的成绩s(t)f(t)，得到的核磁共振谱图则是与的卷积，这里是s(t)的傅立叶变换。当信号未衰减完毕即停止采样，即采样截断(或称为不完全采样)时，s(t)实际上是一个离散矩形函数，其傅立叶变换是sinc函数。正是与sinc函数的卷积，给谱图中引入了震荡状的伪峰。
当对谱峰按某阈值，如0.95，作削峰处理后，由于峰尖远离震荡状截断伪峰所在的底部区域，因此被削掉的小峰不含截断伪峰，利用逆傅立叶变换(iFFT)将其变换到时域后，再经特意的截断处理，则相当于这部分时域信号与采样函数s(t)相乘，对应在频域里是小峰也在底部出现截断伪峰。
在时域中，将截断后的小峰时域信号，从原始时域信号中减去，相当于在谱图中减去小峰的同时，也减去了小峰在截断伪峰中的贡献。
以上步骤不断迭代重复，每次迭代中被从原谱图中移去的小峰，叠加在一起组成新的谱图，直至原谱图中谱峰高度被消减至噪声水平。经以上处理可有效消除截断伪峰。
本发明相对于现有技术具有以下有益效果：
1、对于多维实验，可以减少间接维采样点数，提高实验速度，同时能避免采样点数不足带来的采样截断伪峰，以及避免常规窗函数处理方法带来的谱峰展宽和重叠，提高谱图分辨率。
2、对于Pure Shift实验，可以减少伪间接维采样点数，提高实验速度，同时能避免采样点数不足带来的采样截断伪峰，以及避免常规窗函数处理方法带来的谱峰展宽重叠，提高谱图分辨率。
附图说明
图1为单脉冲一维NMR实验的脉冲序列。
图2为蛋白质Ubiquitin的一维氢谱。
图3为蛋白质Ubiquitin的COSY谱图中的部分示意图。
图4为COSY实验的脉冲序列。
图5为采样截断伪峰的由来示意图。当时域信号衰减完毕时结束采样(A)，则谱图正常(B)；当时域信号还未衰减完毕就结束采样，即采样截断(C)，则会在谱图中引入振荡状的伪峰(D)。
图6为截断采样数据的三种方法处理结果对比示意图。(A)完全采样所得的谱图；(B)对采样截断信号，经FFT所得的谱图，其中采样截断造成的伪峰很明显，左侧强信号峰旁的弱信号被伪峰掩盖；(C)，利用常规的窗函数方法，可消除截断伪峰，但是的谱峰变宽，相互靠近的信号会重叠在一起。(D)利用本问所提出的新方法(TASI)在消除截断伪峰的同时，避免了谱峰展宽带来的重叠。
图7为本发明的流程示意图，A为NMR谱图A，B为NMR谱图B，C为NMR谱图C，D为NMR谱图D，E为NMR谱图E，B1～Bn为各个迭代循环过程中得到的NMR谱图B，C1～Cn为各个迭代循环过程中得到的NMR谱图C。
图8为实施例2中对Pure Shift实验所获得的原始拼接时域信号示意图。
图9为实施例2中的NMR谱图A的示意图，谱峰的底部，有明显的截断伪峰存在。
图10为实施例2中的NMR谱图B的示意图。
图11为实施例2中经过iFFT以及截断处理后的时域信号b的示意图。
图12为实施例2中的步骤6获得的新的时域信号a的示意图。
图13为实施例2中最终NMR谱图E的示意图。与图9所示的NMR谱图A相比，谱峰底部的截断伪峰已经被有效抑制。
具体实施方式
以下结合附图对本发明的技术方案进一步详细说明：
实施例1：
步骤1、设定第一百分比阈值p和第二百分比阈值q，第一百分比阈值p的取值范围为95％～98％，第二百分比阈值q的取值范围为0.005％～0.015％。
步骤2、对获得的时域信号a作填零和傅立叶变换处理。时域信号a是一个长度为N1的复数数组复数数组中的各个元素的序号为M，复数数组中的各个元素是在采样时间为M·△t时的信号强度。对时域信号a的填零处理是依次将数字零添加在复数数组的尾部，将复数数组的长度由N1增加至N2。
傅立叶变换是通过如下公式，将填零处理后的时域信号a由时域变换到频域得到NMR谱图A，NMR谱图A为
其中xk为长度为N2的填零处理后的时序信号a中的第k个元素，j,k均为0到N2-1的整数，i为虚数单位，即e为自然对数；yj为NMR谱图A中的第j个元素，π为圆周率。
步骤3、令NMR谱图A中绝对值最大的数值为I。
将步骤2中所设定的第一百分比阈值p与I相乘，所得结果即为第一比较阈值，记第一比较阈值为T＝pI。
以T为阈值，将NMR谱图A截断为两部分得到NMR谱图B与NMR谱图C。
NMR谱图B为
其中，yj1为NMR谱图B的第j1个元素，j1的取值范围为0～N2-1；
NMR谱图C为
其中，yj2为NMR谱图C的第j2个元素，j2的取值范围为0～N2-1；
步骤4、将步骤2获得的NMR谱图A中最高信号峰的强度值I与第二百分比阈值q的乘积qI为第二比较阈值，
若NMR谱图B中最高信号峰的强度值小于第二比较阈值，则进入步骤7；
若NMR谱图B中最高信号峰的强度值大于等于第二比较阈值，则进入步骤5；
步骤5、将NMR谱图B进行逆傅立叶变换(iFFT)得到时域信号b，并对时域信号b作截断处理，即将N2个采样数据点的时域信号b只保留前面的N1个采样数据点。
将NMR谱图B通过逆傅立叶变换(iFFT)变换为时域信号b，时域信号b为逆傅立叶变换公式为
其中y'j1为NMR谱图B中的第j1个元素，x'k为长度为N2时序信号a中的第k个元素，j1,k均为为0到N2-1的整数，i为虚数单位，即e为自然对数；π为圆周率。
截断处理：将时域信号b只保留前面N1个数据点，后面(N2-N1)个数据点直接丢弃。
步骤6、将时域信号a与截断处理后的时域信号b相减，得到新的时域信号a并返回步骤2。
anew＝aold-b
步骤7：将各个循环迭代中获得的NMR谱图B，以及最后一次循环迭代中获得的NMR谱图C相叠加，即得到NMR谱图D。
上式中，CZ1为最后一次循环迭代中获得的NMR谱图C，BZ2为各个循环迭代中获得的NMR谱图B，Z1为循环迭代的次数。
步骤8：对NMR谱图D作窗口平滑处理修正，得到最终NMR谱图E。即为NMR谱图D中每个数据点设定一个点数窗口，点数窗口一般为5～10个点，这个点数窗口中所有数据点的平均值，作为该数据点修正后的强度值。经此步骤，得到最终NMR谱图E。
令NMR谱图E为NMR谱图D为则
其中w为平滑窗口参数，一般设定为3～5即可，s∈[0～N2-1]。
实施例2：
Pure Shift实验是一类直接维去偶实验，利用Pure Shift实验获得的谱图中，偶合裂分被抑制，仅保留化学位移(Chemical Shift)信息，显著减少了谱峰重叠，简化了谱图分析，因此被称为Pure Shift实验。
一维Pure Shift实验实际上一个伪二维实验，是将一系列一维实验的采样数据，提取各自的前若干个采样点，拼接成一个新的时域信号，再对该信号进行处理，得到最终谱图。
因此一维Pure Shift实验比常规一维实验要费时的多，但是如果减少间接维采样次数，固然能减少实验时间，但所拼接的时域信号就面临着不完全采样引入的截断效应，在谱图中导致截断伪峰的出现。
利用本发明所提出的处理方法，则可以消除截断效应。测试样品为蔗糖，完成PureShift实验及FID拼接之后，利用本发明方法处理过程如下：
步骤1、设定第一百分比阈值p＝95％，第二百分比阈值q＝1.0ⅹ10-5。
步骤2、对Pure Shift实验所获得的原始拼接时域信号作为时域信号a。时域信号a的采样数据点数为N1，对时域信号a进行填零处理，填零后的时域信号a的采样数据点数为N2，然后对进填零后的时域信号a行傅立叶变换后获得NMR谱图A；
步骤3、求得NMR谱图A中最高信号峰的强度值I，以最高信号峰的强度值I与第一百分比阈值p的乘积pI为第一比较阈值，将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值分别减去第一比较阈值，并将NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值设定为0，得到NMR谱图B；将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值设定为第一比较阈值，且NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值不变，得到NMR谱图C。
步骤4、将步骤2获得的NMR谱图A中最高信号峰的强度值I与第二百分比阈值q的乘积qI为第二比较阈值，
若NMR谱图B中最高信号峰的强度值小于第二比较阈值，则进入步骤7；
若NMR谱图B中最高信号峰的强度值大于等于第二比较阈值，则进入步骤5。
步骤5、将NMR谱图B进行逆傅立叶变换得到时域信号b，并对时域信号b作截断处理，即将N2个采样数据点的时域信号b只保留前面的N1个采样数据点。
步骤6、将时域信号a减去截断处理后的时域信号b获得新的时域信号a，将新的时域信号a作为步骤2中的时域信号a，循环迭代进入步骤2。
步骤7、将各个循环迭代中获得的NMR谱图B，以及最后一次循环迭代中获得的NMR谱图C相叠加，即得到NMR谱图D。
步骤8、对NMR谱图D作窗口平滑处理修正，得到最终NMR谱图E。即为NMR谱图D中每个数据点设定一个点数窗口，点数窗口一般为5～10个点，这个点数窗口中所有数据点的平均值，作为该数据点修正后的强度值。经此步骤，得到最终NMR谱图E。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (1)

1.一种核磁共振谱采样截断伪峰的抑制方法，其特征在于，包括以下步骤：
步骤1：设定第一百分比阈值p和第二百分比阈值q，第一百分比阈值p的取值范围为95％～98％，第二百分比阈值q的取值范围为0.005％～0.015％，
步骤2：获得的时域信号a，时域信号a的采样数据点点数为N1，对时域信号a进行填零处理，填零后的时域信号a的采样数据点点数为N2，然后对进填零后的时域信号a进行傅立叶变换后获得NMR谱图A，
步骤3：求得NMR谱图A中最高信号峰的强度值I，以最高信号峰的强度值I与第一百分比阈值p的乘积pI为第一比较阈值，将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值分别减去第一比较阈值，并将NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值设定为0，得到NMR谱图B；将NMR谱图A中所有强度值大于第一比较阈值的采样数据点对应的强度值设定为第一比较阈值，且NMR谱图A中所有强度值小于等于第一比较阈值的采样数据点对应的强度值不变，得到NMR谱图C，
步骤4、将步骤2获得的NMR谱图A中最高信号峰的强度值I与第二百分比阈值q的乘积qI为第二比较阈值，
若NMR谱图B中最高信号峰的强度值小于第二比较阈值，则进入步骤7；
若NMR谱图B中最高信号峰的强度值大于等于第二比较阈值，则进入步骤5，
步骤5：将NMR谱图B进行逆傅立叶变换得到时域信号b，并对时域信号b作截断处理，即将N2个采样数据点的时域信号b只保留前面的N1个采样数据点，
步骤6：将时域信号a减去截断处理后的时域信号b获得新的时域信号a，将新的时域信号a作为步骤2中的时域信号a，循环迭代进入步骤2，
步骤7：将各个循环迭代中获得的NMR谱图B，以及最后一次循环迭代中获得的NMR谱图C相叠加，即得到NMR谱图D，
步骤8：对NMR谱图D作窗口平滑处理修正，得到最终NMR谱图E。
CN201610310469.4A 2016-05-11 2016-05-11 一种核磁共振谱采样截断伪峰的抑制方法 Active CN105928965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610310469.4A CN105928965B (zh) 2016-05-11 2016-05-11 一种核磁共振谱采样截断伪峰的抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610310469.4A CN105928965B (zh) 2016-05-11 2016-05-11 一种核磁共振谱采样截断伪峰的抑制方法

Publications (2)

Publication Number Publication Date
CN105928965A CN105928965A (zh) 2016-09-07
CN105928965B true CN105928965B (zh) 2017-09-15

Family

ID=56834757

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610310469.4A Active CN105928965B (zh) 2016-05-11 2016-05-11 一种核磁共振谱采样截断伪峰的抑制方法

Country Status (1)

Country Link
CN (1) CN105928965B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111175335B (zh) * 2020-01-17 2023-05-02 武汉中科牛津波谱技术有限公司 一种核磁共振t1噪声的压制方法
CN111239657B (zh) * 2020-01-20 2022-05-06 上海东软医疗科技有限公司 谱图的相位校正方法、装置及设备
CN112395983B (zh) * 2020-11-18 2022-03-18 深圳市步锐生物科技有限公司 一种质谱数据峰位置对齐方法及装置
CN114487954B (zh) * 2022-04-14 2022-07-01 中国科学院精密测量科学与技术创新研究院 一种精确测量电磁体场强及分布的多通道收发nmr方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1116079A (zh) * 1994-05-31 1996-02-07 岛津制作所株式会社 磁共振成像装置
CN102305918A (zh) * 2011-05-24 2012-01-04 中国科学院武汉物理与数学研究所 一种核磁共振多维谱的伪峰抑制方法
CN103267768A (zh) * 2013-01-15 2013-08-28 长沙理工大学 可用于变压器故障诊断的混合相油核磁共振分析方法
WO2014024722A1 (ja) * 2012-08-04 2014-02-13 株式会社 東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6040972A (ja) * 1983-08-16 1985-03-04 Hitachi Ltd 核磁気共鳴装置
US6960913B2 (en) * 2003-01-14 2005-11-01 Schlumberger Technology Corporation Multi-measurement NMR analysis based on maximum entropy

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1116079A (zh) * 1994-05-31 1996-02-07 岛津制作所株式会社 磁共振成像装置
CN102305918A (zh) * 2011-05-24 2012-01-04 中国科学院武汉物理与数学研究所 一种核磁共振多维谱的伪峰抑制方法
WO2014024722A1 (ja) * 2012-08-04 2014-02-13 株式会社 東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
CN103267768A (zh) * 2013-01-15 2013-08-28 长沙理工大学 可用于变压器故障诊断的混合相油核磁共振分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NASR: An Effective Approach for Simultaneous Noise and Artifact Suppression in NMR Spectroscopy;Bin Jiang et al.;《American Chemical Society》;20130122;第85卷(第4期);2523-2528 *
基于迭代TSVD的NMR二维谱反演算法;周小龙 等;《波谱学杂志》;20131231;第30卷(第4期);541-551 *

Also Published As

Publication number Publication date
CN105928965A (zh) 2016-09-07

Similar Documents

Publication Publication Date Title
CN105928965B (zh) 一种核磁共振谱采样截断伪峰的抑制方法
CN101627910B (zh) 磁共振成像装置及磁共振成像方法
CN104237820B (zh) 一种单扫描获取磁共振二维j‑分解谱的方法
CN103885013B (zh) 一种不均匀磁场下获得核磁共振二维j分解谱的方法
CN106093099B (zh) 一种获得高分辨二维j分解谱的方法
CN102305918B (zh) 一种核磁共振多维谱的伪峰抑制方法
CN103860176B (zh) 减少成像时间的mr并行成像系统
CN106841270B (zh) 一种获得核磁共振二维相敏j谱的方法
CN103744042B (zh) 在不均匀磁场下获得核磁共振二维自旋回波相关谱的方法
Grissom et al. Additive angle method for fast large‐tip‐angle RF pulse design in parallel excitation
CN105232045A (zh) 基于双回波的单扫描定量磁共振扩散成像方法
CN103529411B (zh) 一种基于梯度编码的自动匀场方法
CN101598812A (zh) 去除数字检波器单点接收地震记录中的异常噪声方法
CN107144591A (zh) 一种测量独立同核自旋对间接偶合模式的方法
CN107423543A (zh) 一种超复数磁共振波谱的快速重建方法
WO2012127494A4 (en) A method and system for rapid mri acquisition using tailored signal excitation modules (rate)
Marion Processing of ND NMR spectra sampled in polar coordinates: a simple Fourier transform instead of a reconstruction
Shim et al. Two-dimensional NMR spectroscopy of 13C methanol at less than 5 μT
Karlsons et al. Identification of sampling patterns for high‐resolution compressed sensing MRI of porous materials:‘learning’from X‐ray microcomputed tomography data
CN106725480A (zh) 磁共振图像采集与重建方法及装置
WO2002079807A1 (en) Nmr spectroscopy data recovery method and apparatus
CN108680595A (zh) 一种低场磁共振食用油种类检测方法
CN111538944B (zh) 一种基于子空间的磁共振波谱快速重建方法
DE69029091T2 (de) Verfahren zur bilderzeugung mittels magnetischer kernresonanz
CN107728088A (zh) 一种新型快速核磁共振t1成像方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210223

Address after: 430071 Xiao Hong, Wuchang District, Wuhan District, Hubei, Shanxi, 30

Patentee after: Institute of precision measurement science and technology innovation, Chinese Academy of Sciences

Address before: 430071 Xiao Hong, Wuchang District, Wuhan District, Hubei, Shanxi, 30

Patentee before: WUHAN INSTITUTE OF PHYSICS AND MATHEMATICS, CHINESE ACADEMY OF SCIENCES

TR01 Transfer of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20160907

Assignee: Wuhan Lianying Life Science Instrument Co.,Ltd.

Assignor: Institute of precision measurement science and technology innovation, Chinese Academy of Sciences

Contract record no.: X2021980001566

Denomination of invention: A method to suppress false peaks in NMR spectra

Granted publication date: 20170915

License type: Common License

Record date: 20210310

EE01 Entry into force of recordation of patent licensing contract