CN112731234B - 一种基于圆弧切线模型的核磁共振波谱基线校正方法 - Google Patents
一种基于圆弧切线模型的核磁共振波谱基线校正方法 Download PDFInfo
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 title claims abstract description 40
- 238000001228 spectrum Methods 0.000 claims abstract description 68
- 230000003595 spectral effect Effects 0.000 claims abstract description 44
- 238000009499 grossing Methods 0.000 claims abstract description 26
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 230000001131 transforming effect Effects 0.000 claims abstract description 5
- 230000009466 transformation Effects 0.000 claims description 11
- 238000013480 data collection Methods 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 210000002966 serum Anatomy 0.000 description 2
- 238000005160 1H NMR spectroscopy Methods 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000008363 phosphate buffer Substances 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000000079 presaturation Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 206010039073 rheumatoid arthritis Diseases 0.000 description 1
- AJPJDKMHJJGVTQ-UHFFFAOYSA-M sodium dihydrogen phosphate Chemical compound [Na+].OP(O)([O-])=O AJPJDKMHJJGVTQ-UHFFFAOYSA-M 0.000 description 1
- 229910000162 sodium phosphate Inorganic materials 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000000264 spin echo pulse sequence Methods 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/46—NMR spectroscopy
-
- 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
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)尺度归一化:采用公式对高斯平滑后的波谱数据进行尺度变换,其中,y表示尺度变换后的波谱数据,ysmooth表示高斯平滑后的波谱数据,n表示高斯平滑后的波谱数据点数,ymax表示高斯平滑后的波谱数据的峰强度最大值。
在步骤2)中,所述构建圆弧方程的具体方法为:
(1)以预处理后的波谱数据维度为横坐标,谱峰强度为纵坐标构建坐标系并绘制波谱;
在步骤3)中,得到波谱基线点的具体方法为:
(1)滑动圆弧中心横坐标,滑动范围从1~n,滑动步长为1,滑动后纵向平移圆弧,使得所述圆弧与所述波谱相切;
(2)计算样品核磁共振波谱与圆弧的切线点,确定yo的值,得到确定的圆弧方程,并记录圆弧滑动轨迹;
(3)对于波谱的每一数据点,选择该数据点所对应的圆弧滑动轨迹的最大值,作为该数据点的基线值。
在步骤4)中,所述对基线进行尺度逆变换,采用公式对预处理后的波谱数据的基线进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正,其中,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个数据点计算噪声水平ó,并对所述波谱数据进行高斯平滑,得到高斯平滑后的波谱数据。采用公式对所述高斯平滑后的波谱数据进行尺度变换,其中,ysmooth表示所述高斯平滑后的波谱数据,16382表示所述高斯平滑后的波谱数据的维度,436906874表示波谱数据中谱峰强度的最大值。
3、谱峰识别和参数设定
对所述预处理后的波谱数据进行谱峰识别,并计算每个谱峰的半高宽。将半径R设定为谱峰半高宽的最大值。
4、计算基线模型
以所述预处理后的波谱数据维度为横坐标,谱峰强度为纵坐标绘制波谱。将基线初始化为16382维零向量。对于所述预处理后的波谱数据的每一个点,以R为半径作圆,使所作的圆与所述预处理后的波谱数据曲线相切,且圆心在所述预处理后的波谱数据曲线下方,将该点处的圆弧纵坐标值作为该点处的基线值。
5、尺度恢复
6、用原始波谱数据减去最终的基线,实现核磁共振波谱的基线校正。
图1为本发明实施例基线校正效果示意图。其中,(a)为样品的核磁共振波谱数据;
(b)为应用本发明方法构建的基线模型;虚线曲线(·)表示应用本方法构建的基线模型;
(c)为应用本发明方法进行基线校正后的样品核磁共振波谱数据。本发明通过构建圆弧曲线方程,让圆弧曲线与核磁共振波谱信号相切,利用圆弧滑动,得到一条切线点轨迹作为波谱的候选基线点,通过平滑处理,最后得到样品的基线。计算速度快,可以自适应地求出波谱不同位置上的基线,且避免基线校正引入的谱峰畸变,具有重要的意义。
以上所述,仅为本发明较佳实施例而已,故不能依此限定本发明实施的范围,即依本发明专利范围及说明书内容所作的等效变化与修饰,皆应仍属本发明涵盖的范围内。
Claims (5)
1.一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于包括以下步骤:
1)数据收集和预处理:将所采集到的fid数据末尾补零至2的整次幂个数据点,进行快速傅里叶变换、相位校正,得到样品的核磁共振波谱数据,并进行数据预处理;
2)圆弧方程:对样品核磁共振波谱进行谱峰识别,计算每个谱峰的半高宽,设置圆弧半径R为谱峰半高宽的中值,构建圆弧方程;
3)圆弧切线模型:滑动圆弧的圆心位置,计算样品核磁共振波谱与圆弧的切线点,得到圆弧滑动轨迹,根据圆弧滑动轨迹计算得到波谱基线点;
4)基线校正:对样品核磁共振波谱基线点进行平滑和尺度逆变换,得到波谱基线,实现样品核磁共振波谱的基线校正。
4.如权利要求1所述一种基于圆弧切线模型的核磁共振波谱基线校正方法,其特征在于在步骤3)中,所述得到波谱基线点的具体方法为:
(1)滑动圆弧中心横坐标,滑动范围从1~n,滑动步长为1,滑动后纵向平移圆弧,使得所述圆弧与所述波谱相切;
(2)计算样品核磁共振波谱与圆弧的切线点,确定yo的值,得到确定的圆弧方程,并记录圆弧滑动轨迹;
(3)对于波谱的每一数据点,选择该数据点所对应的圆弧滑动轨迹的最大值,作为该数据点的基线值。
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)
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 | 南通乐尔环保科技有限公司 | 一种应用于色谱分析的重叠峰处理办法及装置 |
-
2020
- 2020-12-29 CN CN202011594412.4A patent/CN112731234B/zh not_active Expired - Fee Related
Patent Citations (9)
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)
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 |