CN112731234B - 一种基于圆弧切线模型的核磁共振波谱基线校正方法 - Google Patents

一种基于圆弧切线模型的核磁共振波谱基线校正方法 Download PDF

Info

Publication number
CN112731234B
CN112731234B CN202011594412.4A CN202011594412A CN112731234B CN 112731234 B CN112731234 B CN 112731234B CN 202011594412 A CN202011594412 A CN 202011594412A CN 112731234 B CN112731234 B CN 112731234B
Authority
CN
China
Prior art keywords
spectrum
baseline
data
arc
magnetic resonance
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.)
Expired - Fee Related
Application number
CN202011594412.4A
Other languages
English (en)
Other versions
CN112731234A (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.)
Xiamen University
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN202011594412.4A priority Critical patent/CN112731234B/zh
Publication of CN112731234A publication Critical patent/CN112731234A/zh
Application granted granted Critical
Publication of CN112731234B publication Critical patent/CN112731234B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/46NMR spectroscopy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

一种基于圆弧切线模型的核磁共振波谱基线校正方法,涉及核磁共振波谱基线校正。1)将所采集到的fid数据末尾补零至2的整次幂个数据点,进行快速傅里叶变换、相位校正,得到样品的核磁共振波谱数据,并进行数据预处理;2)对样品核磁共振波谱进行谱峰识别,计算每个谱峰的半高宽,设置圆弧半径R为谱峰半高宽的中值,构建圆弧方程;3)滑动圆弧的圆心位置,计算样品核磁共振波谱与圆弧的切线点,得到圆弧滑动轨迹,根据圆弧滑动轨迹计算得到波谱基线点;4)对样品核磁共振波谱基线点进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正。可以校正基线失真较大的波谱。避免校正后出现负峰,可推广至二维波谱的基线校正。

Description

一种基于圆弧切线模型的核磁共振波谱基线校正方法
技术领域
本发明涉及核磁共振波谱基线校正,尤其是涉及一种基于圆弧切线模型的核磁共振波谱基线校正方法。
背景技术
核磁共振波谱是研究分子和物质结构的重要工具,基线失真会对波谱解析以及波谱定量分析造成很大影响,消除基线失真是波谱数据处理的一个必要步骤。目前,常用的基线校正方法以频域基线校正方法为主,包括基线点识别和构建基线模型两个具体步骤,其主要原理是通过有效的数据后处理,构建一条平滑的基线模型,以实现基线校正。
根据基线模型的不同,可以将基线校正方法分为无参数和含参数两类方法:多项式拟合是无参数基线校正的典型方法之一,它假设波谱所有位置上的基线都服从同样的非线性关系,并利用特定的函数来拟合这种非线性关系;但实际样品的波谱基线往往不满足这一假设,因此基于多项式拟合的方法很难给出理想的结果;含参数的基线校正算法方法,利用连续小波变换、滑动窗口法等基线识别方法进行基线点识别,并利用平滑算法如Whittaker滤波等构建基线模型,能取得更好的效果,但在处理复杂数据的时候,尤其是包含宽峰和多个重叠峰的数据时,可能会使谱峰信号发生畸变。
目前,许多标准谱仪软件允许用户手动设置基线点,并使用分析函数进行插值得到基线模型,这个方法虽然能取得较理想的效果,但需要手动干预,无法用于批处理,并且不适用于二维波普数据。目前,核磁共振波谱数据,特别是二维波谱数据的基线校正,仍然是一个没有完全解决的问题。
发明内容
本发明的目的在于克服现有技术的上述不足之处,提供计算速度快、可自适应得出波谱不同位置上的基线,且避免基线校正引入的谱峰畸变的一种基于圆弧切线模型的核磁共振波谱基线校正方法。
本发明包括以下步骤:
1)数据收集和预处理:将所采集到的fid数据末尾补零至2的整次幂个数据点,进行快速傅里叶变换、相位校正,得到样品的核磁共振波谱数据,并进行数据预处理;
2)圆弧方程:对样品核磁共振波谱进行谱峰识别,计算每个谱峰的半高宽,设置圆弧半径R为谱峰半高宽的中值,构建圆弧方程;
3)圆弧切线模型:滑动圆弧的圆心位置,计算样品核磁共振波谱与圆弧的切线点,得到圆弧滑动轨迹,根据圆弧滑动轨迹计算得到波谱基线点;
4)基线校正:对样品核磁共振波谱基线点进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正。
在步骤1)中,数据预处理方法包括:
(1)高斯平滑:选择所述波谱数据的一段无谱峰区间,计算该区间数据点的标准差σ作为噪声水平,并对所述波谱数据进行窗宽为2σ的高斯平滑,得到高斯平滑后的波谱数据;
(2)尺度归一化:采用公式
Figure BDA0002869917870000021
对高斯平滑后的波谱数据进行尺度变换,其中,y表示尺度变换后的波谱数据,ysmooth表示高斯平滑后的波谱数据,n表示高斯平滑后的波谱数据点数,ymax表示高斯平滑后的波谱数据的峰强度最大值。
在步骤2)中,所述构建圆弧方程的具体方法为:
(1)以预处理后的波谱数据维度为横坐标,谱峰强度为纵坐标构建坐标系并绘制波谱;
(2)以参数xo、yo分别表示圆心横坐标、纵坐标,以R为半径构建含参数的圆弧方程:
Figure BDA0002869917870000022
在步骤3)中,得到波谱基线点的具体方法为:
(1)滑动圆弧中心横坐标,滑动范围从1~n,滑动步长为1,滑动后纵向平移圆弧,使得所述圆弧与所述波谱相切;
(2)计算样品核磁共振波谱与圆弧的切线点,确定yo的值,得到确定的圆弧方程,并记录圆弧滑动轨迹;
(3)对于波谱的每一数据点,选择该数据点所对应的圆弧滑动轨迹的最大值,作为该数据点的基线值。
在步骤4)中,所述对基线进行尺度逆变换,采用公式
Figure BDA0002869917870000023
对预处理后的波谱数据的基线进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正,其中,b为波谱数据的基线,btemp为预处理后的波谱数据的基线。
与现有技术相比,本发明具有如下突出的优点:
1.本发明的基线校正方法,可以校正基线失真较大的波谱。
2.本发明的基线校正方法,可以避免校正后出现负峰,同时对正峰的峰强度的影响较小,可以避免在复杂波谱中由于基线校正引入的谱峰失真。
3.本发明的基线校正方法,不仅可以用于核磁共振波谱的基线校正,也可以推广至近红外光谱、质谱等其他波谱数据的基线校正。
4.本发明的基线校正方法,可以将圆弧切线模型替换为圆球切面模型,可以很好的推广至二维波谱的基线校正。
附图说明
图1为本发明实施例基线校正效果示意图。其中,(a)为样品的核磁共振波谱数据;
(b)为应用本发明方法构建的基线模型;虚线曲线(·)表示应用本方法构建的基线模型;
(c)为应用本发明方法进行基线校正后的样品核磁共振波谱数据。
具体实施方式
以下实施例将结合附图对本发明作进一步的说明。
本发明所述的一种基于圆弧切线模型的核磁共振波谱基线校正方法,包括以下步骤:采集核磁共振波谱数据;对该数据进行预处理;对预处理后的波谱数据进行谱峰识别,并计算谱峰半高宽;设置圆弧半径R为谱峰半高宽的中值,构建圆弧方程;构建圆弧切线模型,通过滑动圆弧中心位置得到波谱基线点;对得到的基线点进行平滑和尺度逆变换,得到波谱数据的基线模型;用原始波谱数据减去最终的基线,实现核磁共振波谱的基线校正。
本发明实施例具体包括以下步骤:
1、数据采集
本发明使用在某医学院采集的44个类风湿性关节炎大鼠血清样本。分析前取150ìl采集到的血浆与350ìl磷酸盐缓冲液(50mmol K2HPO4/NaH2PO4,pH 7.4,0.9%NaCl,99.9%D2O)混和,在4℃下以10000r/min的转速离心10min,离心后取450ìl上清液存储于5mm NMR管中并于4℃温度下保存直至实验。1D 1H NMR实验在Bruker Advance 850MHZ谱仪上进行。预饱和序列来抑制水峰信号并结合CPMG(Carr-Purcell-Meiboom-Gill)自旋回波序列来采集血清谱。试验温度设置为25℃,每一个谱共累加64次,谱宽17KHz。将所采集到的fid数据末尾补零至2的整次幂个数据点,进行快速傅里叶变换、相位校正,得到样品的核磁共振波谱数据,并将每个样本的核磁共振波谱数据保存为16382维向量。
2、数据预处理
选择所述波谱数据的无谱峰的2000个数据点计算噪声水平ó,并对所述波谱数据进行高斯平滑,得到高斯平滑后的波谱数据。采用公式
Figure BDA0002869917870000041
对所述高斯平滑后的波谱数据进行尺度变换,其中,ysmooth表示所述高斯平滑后的波谱数据,16382表示所述高斯平滑后的波谱数据的维度,436906874表示波谱数据中谱峰强度的最大值。
3、谱峰识别和参数设定
对所述预处理后的波谱数据进行谱峰识别,并计算每个谱峰的半高宽。将半径R设定为谱峰半高宽的最大值。
4、计算基线模型
以所述预处理后的波谱数据维度为横坐标,谱峰强度为纵坐标绘制波谱。将基线初始化为16382维零向量。对于所述预处理后的波谱数据的每一个点,以R为半径作圆,使所作的圆与所述预处理后的波谱数据曲线相切,且圆心在所述预处理后的波谱数据曲线下方,将该点处的圆弧纵坐标值作为该点处的基线值。
5、尺度恢复
采用公式
Figure BDA0002869917870000042
对所述预处理后的波谱数据的基线进行尺度变换,得到波谱数据的基线。
6、用原始波谱数据减去最终的基线,实现核磁共振波谱的基线校正。
图1为本发明实施例基线校正效果示意图。其中,(a)为样品的核磁共振波谱数据;
(b)为应用本发明方法构建的基线模型;虚线曲线(·)表示应用本方法构建的基线模型;
(c)为应用本发明方法进行基线校正后的样品核磁共振波谱数据。本发明通过构建圆弧曲线方程,让圆弧曲线与核磁共振波谱信号相切,利用圆弧滑动,得到一条切线点轨迹作为波谱的候选基线点,通过平滑处理,最后得到样品的基线。计算速度快,可以自适应地求出波谱不同位置上的基线,且避免基线校正引入的谱峰畸变,具有重要的意义。
以上所述,仅为本发明较佳实施例而已,故不能依此限定本发明实施的范围,即依本发明专利范围及说明书内容所作的等效变化与修饰,皆应仍属本发明涵盖的范围内。

Claims (5)

1.一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于包括以下步骤:
1)数据收集和预处理:将所采集到的fid数据末尾补零至2的整次幂个数据点,进行快速傅里叶变换、相位校正,得到样品的核磁共振波谱数据,并进行数据预处理;
2)圆弧方程:对样品核磁共振波谱进行谱峰识别,计算每个谱峰的半高宽,设置圆弧半径R为谱峰半高宽的中值,构建圆弧方程;
3)圆弧切线模型:滑动圆弧的圆心位置,计算样品核磁共振波谱与圆弧的切线点,得到圆弧滑动轨迹,根据圆弧滑动轨迹计算得到波谱基线点;
4)基线校正:对样品核磁共振波谱基线点进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正。
2.如权利要求1所述一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于在步骤1)中,数据预处理的方法包括:
(1)高斯平滑:选择所述波谱数据的一段无谱峰区间,计算该区间数据点的标准差σ作为噪声水平,并对所述波谱数据进行窗宽为2σ的高斯平滑,得到高斯平滑后的波谱数据;
(2)尺度归一化:采用公式
Figure FDA0003254644930000011
对高斯平滑后的波谱数据进行尺度变换,其中,y表示尺度变换后的波谱数据,ysmooth表示高斯平滑后的波谱数据,n表示高斯平滑后的波谱数据点数,ymax表示高斯平滑后的波谱数据的峰强度最大值。
3.如权利要求1所述一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于在步骤2)中,所述构建圆弧方程的具体方法为:
(1)以预处理后的波谱数据维度为横坐标,谱峰强度为纵坐标构建坐标系并绘制波谱;
(2)以参数xo、yo分别表示圆心横坐标、纵坐标,以R为半径构建含参数的圆弧方程:
Figure FDA0003254644930000012
x∈(xo-R,xo+R)。
4.如权利要求1所述一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于在步骤3)中,所述得到波谱基线点的具体方法为:
(1)滑动圆弧中心横坐标,滑动范围从1~n,滑动步长为1,滑动后纵向平移圆弧,使得所述圆弧与所述波谱相切;
(2)计算样品核磁共振波谱与圆弧的切线点,确定yo的值,得到确定的圆弧方程,并记录圆弧滑动轨迹;
(3)对于波谱的每一数据点,选择该数据点所对应的圆弧滑动轨迹的最大值,作为该数据点的基线值。
5.如权利要求1所述一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于在步骤4)中,所述对样品核磁共振波谱基线点进行平滑和尺度逆变换,采用公式
Figure FDA0003254644930000021
对预处理后的波谱数据的基线进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正,其中,b为波谱数据的基线,btemp为预处理后的波谱数据的基线,n为高斯平滑后的波谱数据点数,ymax为高斯平滑后波谱的波谱数据的峰强度最大值。
CN202011594412.4A 2020-12-29 2020-12-29 一种基于圆弧切线模型的核磁共振波谱基线校正方法 Expired - Fee Related CN112731234B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011594412.4A CN112731234B (zh) 2020-12-29 2020-12-29 一种基于圆弧切线模型的核磁共振波谱基线校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011594412.4A CN112731234B (zh) 2020-12-29 2020-12-29 一种基于圆弧切线模型的核磁共振波谱基线校正方法

Publications (2)

Publication Number Publication Date
CN112731234A CN112731234A (zh) 2021-04-30
CN112731234B true CN112731234B (zh) 2021-10-29

Family

ID=75607870

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011594412.4A Expired - Fee Related CN112731234B (zh) 2020-12-29 2020-12-29 一种基于圆弧切线模型的核磁共振波谱基线校正方法

Country Status (1)

Country Link
CN (1) CN112731234B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102360502A (zh) * 2011-09-07 2012-02-22 中国科学院武汉物理与数学研究所 一种自动基线校正方法
DE102013005455A1 (de) * 2012-03-29 2013-10-02 Jasco Corporation Zirkulardichroismus-Spektrometer mit einem Ausrichtungsmechanismus
CN103635136A (zh) * 2011-04-29 2014-03-12 磁共振创新有限公司 组织相似度图
CN104458785A (zh) * 2014-12-12 2015-03-25 中国科学院武汉物理与数学研究所 一种核磁共振波谱谱峰对齐及谱峰提取方法
US9702810B1 (en) * 2014-09-22 2017-07-11 University Of Puerto Rico Dual cell holder system
CN108287137A (zh) * 2017-12-22 2018-07-17 必欧瀚生物技术(合肥)有限公司 一种基于分段多项式拟合的基线校正方法
CN108844939A (zh) * 2018-03-14 2018-11-20 西安电子科技大学 基于非对称加权最小二乘的拉曼光谱检测基线校正方法
CN111161171A (zh) * 2019-12-18 2020-05-15 三明学院 爆破振动信号基线漂零校正和消噪方法、装置、设备及系统
CN111337606A (zh) * 2020-03-19 2020-06-26 南通乐尔环保科技有限公司 一种应用于色谱分析的重叠峰处理办法及装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103635136A (zh) * 2011-04-29 2014-03-12 磁共振创新有限公司 组织相似度图
CN102360502A (zh) * 2011-09-07 2012-02-22 中国科学院武汉物理与数学研究所 一种自动基线校正方法
DE102013005455A1 (de) * 2012-03-29 2013-10-02 Jasco Corporation Zirkulardichroismus-Spektrometer mit einem Ausrichtungsmechanismus
US9702810B1 (en) * 2014-09-22 2017-07-11 University Of Puerto Rico Dual cell holder system
CN104458785A (zh) * 2014-12-12 2015-03-25 中国科学院武汉物理与数学研究所 一种核磁共振波谱谱峰对齐及谱峰提取方法
CN108287137A (zh) * 2017-12-22 2018-07-17 必欧瀚生物技术(合肥)有限公司 一种基于分段多项式拟合的基线校正方法
CN108844939A (zh) * 2018-03-14 2018-11-20 西安电子科技大学 基于非对称加权最小二乘的拉曼光谱检测基线校正方法
CN111161171A (zh) * 2019-12-18 2020-05-15 三明学院 爆破振动信号基线漂零校正和消噪方法、装置、设备及系统
CN111337606A (zh) * 2020-03-19 2020-06-26 南通乐尔环保科技有限公司 一种应用于色谱分析的重叠峰处理办法及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Iterative baseline correction algorithm for dead time truncated one-dimensional solid-state MAS NMR spectra;Maxime Yon *;《Solid State Nuclear Magnetic Resonance》;20201204;第1-8页 *
Tomasz P. Zielin'ski.Fast MinMax energy-based phase correction method for NMR spectra with linear phase distortion.《Journal of Magnetic Resonance》.2017, *
核磁共振代谢轮廓的基线校正和去卷积新方法;龙怡驰;《厦门大学研究生院官网》;20210520;第1-2页 *
核磁共振谱自动基线校正新方法;钟彩娇;《波谱学杂志》;20141231;第523-533页 *
采用曲线拟合方法进行EPR波谱基线校正;邹曦露;《波谱学杂志》;19961231;第602-609页 *

Also Published As

Publication number Publication date
CN112731234A (zh) 2021-04-30

Similar Documents

Publication Publication Date Title
Cobas et al. A new general-purpose fully automatic baseline-correction procedure for 1D and 2D NMR data
CN109145727B (zh) 一种基于vmd参数优化的轴承故障特征提取方法
CN105279379B (zh) 基于凸组合核函数主成分分析的太赫兹光谱特征提取方法
CN104899906B (zh) 基于自适应正交基的磁共振图像重建方法
CN110619296B (zh) 一种基于奇异分解的信号降噪方法
Hoang Wavelet-based spectral analysis
Bao et al. A new automatic baseline correction method based on iterative method
CN103941091A (zh) 基于改进emd端点效应的电力系统hht谐波检测方法
CN110503060B (zh) 一种光谱信号去噪方法及其系统
CN110501664A (zh) 磁共振波谱相位校正方法及应用该方法的磁共振成像系统
CN102008306B (zh) 磁共振并行成像方法
Cobas Applications of the Whittaker smoother in NMR spectroscopy
CN112731234B (zh) 一种基于圆弧切线模型的核磁共振波谱基线校正方法
CN109102818B (zh) 一种基于信号频率概率密度函数分布的去噪音频采样算法
CN113095353B (zh) 一种基于ap聚类的欠定盲源分离方法
CN117290788B (zh) 一种基于改进小波变换算法的配电网故障辨识方法及系统
Antoniadis et al. Peaks detection and alignment for mass spectrometry data
CN104360295B (zh) 基于字典学习的磁共振并行成像方法和装置
CN105954223A (zh) 一种提高汽油性质预测精度的方法
CN112816779A (zh) 一种解析信号生成的谐波实信号参数估计方法
CN108537738A (zh) 一种矩阵补全方法
CN109460614B (zh) 基于瞬时带宽的信号时间-频率分解方法
CN112054803B (zh) 一种基于压缩感知的通信信号的分选方法
CN106934783B (zh) 一种高频次遥感时间序列数据的时域重建方法
KR20120115729A (ko) 스펙트럼 데이터의 이상값 영향을 최소화하기 위한 전처리 방법

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211029

CF01 Termination of patent right due to non-payment of annual fee