CN112730712B - 提高lc-ms数据信噪比的方法 - Google Patents
提高lc-ms数据信噪比的方法 Download PDFInfo
- Publication number
- CN112730712B CN112730712B CN202011614196.5A CN202011614196A CN112730712B CN 112730712 B CN112730712 B CN 112730712B CN 202011614196 A CN202011614196 A CN 202011614196A CN 112730712 B CN112730712 B CN 112730712B
- Authority
- CN
- China
- Prior art keywords
- peak
- chromatographic peak
- obtaining
- coefficient
- noise ratio
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8603—Signal analysis with integration or differentiation
- G01N30/8617—Filtering, e.g. Fourier filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
- G01N30/72—Mass spectrometers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8631—Peaks
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Optical Recording Or Reproduction (AREA)
- Treatment Of Liquids With Adsorbents In General (AREA)
Abstract
Description
技术领域
本发明涉及色谱,特别涉及提高LC-MS数据信噪比的方法。
背景技术
液相色谱-质谱联用仪(LC-MS)是对复杂混合物样品进行定性分析的重要工具。它融合了液相色谱和质谱两种仪器单独使用时的优点,具有更大的优势,广泛应用于环境、食品、地质检测等各个领域。原始LC-MS数据是LC-MS进行定性分析的基础,但它会包含来自仪器噪声、中性粒子、其他化合物等的干扰,若直接应用原始LC-MS数据对混合物样品进行定性分析,产生的分析结果是不可靠的。此外LC-MS***的复杂性常常使分析人员难以满足方法检测限制。因此需要一种能提高LC-MS信噪比,同时保证信号不失真的算法。
为了提高LC-MS数据的信噪比,常用一些去噪算法进行处理,包括:
1.傅里叶滤波法,傅里叶滤波法对于非平稳信号的处理能力不足且无法获得各频率出现的时刻。
2.小波去噪方法,可以处理非平稳信号,但是小波变换的阈值参数选取困难。
3.多项式平滑算法也称为SG滤波法,对窗口内的数据进行加权,能够平滑噪声,但是它不具有自适应性。
此外,上述算法都是同时对整条谱线进行处理,在滤除噪声的同时可能会引起信号的失真。
LC-MS数据的信号具有局部突变性的特点,在用传统的去噪算法进行处理时,去除噪声的同时往往会把信号的尖端也去掉,从而使信号峰的响应降低,这非常不利于后续的分析。
发明内容
为解决上述现有技术方案中的不足,本发明提供了一种去噪效果好、有用信号不失真的提高LC-MS数据信噪比的方法。
本发明的目的是通过以下技术方案实现的:
提高LC-MS数据信噪比的方法,所述提高LC-MS数据信噪比的方法包括以下步骤:
(A1)检测色谱峰s(t)的峰值位置p和宽度w,获得所述色谱峰s(t)的开始位置w1和结束位置w2;
(A2)将所述色谱峰s(t)划分为峰部分s1(t)和非峰部分s2(t),分别进入步骤(B1)和步骤(C1);
(B1)根据所述色谱峰s(t)设计基函数,所述基函数离散化表示为F,根据F构成四个长度为L的正交滤波器;L·D、H·D和L·R、H·R,L·R是F的归一化表示,H·R是L·R的正交反向滤波器,L·D是L·R反向,H·D是H·R的反向;
(B2)得到峰部分s1(t)的第一层近似系数a1和第一层细节系数d1;
(B3)利用步骤(B2)的方式得到峰部分s1(t)的多层近似系数an和各层细节系数d1,d2···dn,n为大于2的整数;
(B4)将近似系数an和各层细节系数d1,d2···dn分别进行重构,得到峰部分s1(t)的近似部分s1a(t)和细节部分s1b(t);
(C1)对于非峰部分s2(t),按照步骤(B1)-(B4)的方式获得非峰部分s2(t)的近似部分s2a(t)和细节部分;
(C2)整合所述峰部分信息s′1(t)和近似部分s2a(t),得到去噪的色谱峰s′(t)。
与现有技术相比,本发明具有的有益效果为:
去噪效果好、有用信号不失真;
采用峰检测算法,识别出信号峰的位置,根据信号峰的特点,为峰部分和非峰部分构造不同的基函数进行处理,对信号波采取多层次提取有用信息的去噪算法,对非波位置采取强制消噪处理,从而在有效消除噪声的同时不会引起信号的失真,有效地提高信噪比;
对于信号峰的处理,将噪声所在尺度的细节部分进行了imf分量分解并重构,能够避免阈值系数选取的困难,同时保留噪声所在尺度的有用信号,可以使信号不失真。
附图说明
参照附图,本发明的公开内容将变得更易理解。本领域技术人员容易理解的是:这些附图仅仅用于举例说明本发明的技术方案,而并非意在对本发明的保护范围构成限制。图中:
图1是根据本发明实施例的原始信号示意图;
图2是根据小波变换后的信号示意图;
图3是根据本发明实施例的方法处理后的信号示意图。
具体实施方式
图1-3和以下说明描述了本发明的可选实施方式以教导本领域技术人员如何实施和再现本发明。为了教导本发明技术方案,已简化或省略了一些常规方面。本领域技术人员应该理解源自这些实施方式的变型或替换将在本发明的范围内。本领域技术人员应该理解下述特征能够以各种方式组合以形成本发明的多个变型。由此,本发明并不局限于下述可选实施方式,而仅由权利要求和它们的等同物限定。
实施例1:
本发明实施例的提高LC-MS数据信噪比的方法的结构示意图,所述提高LC-MS数据信噪比的方法包括步骤:
(A1)检测色谱峰s(t)的峰值位置p和宽度w,获得所述色谱峰s(t)的开始位置w1和结束位置w2;
(A2)将所述色谱峰s(t)划分为峰部分s1(t)和非峰部分s2(t),分别进入步骤(B1)和步骤(C1);
(B1)根据所述色谱峰s(t)设计基函数,所述基函数离散化表示为F,根据F构成四个长度为L的正交滤波器;L·D、H·D和L·R、H·R,L·R是F的归一化表示,H·R是L·R的正交反向滤波器,L·D是L·R反向,H·D是H·R的反向;
(B2)得到峰部分s1(t)的第一层近似系数a1和第一层细节系数d1;
(B3)利用步骤(B2)的方式得到峰部分s1(t)的多层近似系数an和各层细节系数d1,d2···dn,n为大于2的整数;
(B4)将近似系数an和各层细节系数d1,d2···dn分别进行重构,得到峰部分s1(t)的近似部分s1a(t)和细节部分s1b(t);
(C1)对于非峰部分s2(t),按照步骤(B1)-(B4)的方式获得非峰部分s2(t)的近似部分s2a(t)和细节部分;
(C2)整合所述峰部分信息s′1(t)和近似部分s2a(t),得到去噪的色谱峰s′(t)。
为了降低噪声,并不失真,进一步地,在步骤(B1)中,
F=[0.014 -0.015 -0.124 0.012 0.561 0.640 0.141 -0.028 0.021 0.010],L=10。
为了降低噪声,并不失真,进一步地,在步骤(B2)中,峰部分s1(t)与低通滤波器L·D卷积,得到a1;
峰部分s1(t)与高通滤波器H·D卷积,得到d1。
为了降低噪声,并不失真,进一步地,在步骤(C1)中,基函数离散化表示为:
F=[0.014 -0.015 -0.124 0.012 0.283 0.321 0.141 -0.028 0.021 0.010],L=10。
为了准确寻峰,进一步地,在步骤(A1)中,检测色谱峰s(t)的峰值位置p的方式为:
为了准确地获得峰的宽度,进一步地,在步骤(A1)中,检测色谱峰s(t)的宽度w的方式为:
实施例2:
根据本发明实施例1的提高LC-MS数据信噪比的方法在农残检测中的应用例。
在该应用例中,提高LC-MS数据信噪比的方法包括步骤:
(A1)获得色谱峰s(t),如图1所示,检测色谱峰s(t)的峰值位置p和宽度w,获得所述色谱峰s(t)的开始位置w1和结束位置w2;
峰值位置p的获得方式为:
色谱峰s(t)的宽度w的获得方式为:
(A2)将所述色谱峰s(t)划分为峰部分s1(t)和非峰部分s2(t),分别进入步骤(B1)和步骤(C1);
(B1)根据所述色谱峰s(t)设计基函数,所述基函数离散化表示为F=[0.014 -0.015 -0.124 0.012 0.561 0.640 0.141 -0.028 0.021 0.010],根据F构成四个长度为L的正交滤波器;L·D、H·D和L·R、H·R,L·R是F的归一化表示,H·R是L·R的正交反向滤波器,L·D是L·R反向,H·D是H·R的反向;L=10;
(B3)利用步骤(B2)的方式得到峰部分s1(t)的多层近似系数an和各层细节系数d1,d2···dn,n为大于2的整数;本实施例为三层;
(B4)将近似系数an和各层细节系数d1,d2···dn分别进行重构,得到峰部分s1(t)的近似部分s1a(t)和细节部分s1b(t);
(C1)对于非峰部分s2(t),按照步骤(B1)-(B4)的方式获得非峰部分s2(t)的近似部分s2a(t)和细节部分;
针对非峰部分s2(t)的基函数的离散化为;
F=[0.014 -0.015 -0.124 0.012 0.283 0.321 0.141 -0.028 0.021 0.010],L=10。
(C2)整合所述峰部分信息s′1(t)和近似部分s2a(t),得到去噪的色谱峰s′(t),如图3所示。
图2示意性地给出了利用小波变化处理原始信号后的信号,与本实施例的处理相比,信噪比低,且在非峰部分信号失真。
Claims (6)
1.提高LC-MS数据信噪比的方法,所述提高LC-MS数据信噪比的方法包括以下步骤:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011614196.5A CN112730712B (zh) | 2020-12-31 | 2020-12-31 | 提高lc-ms数据信噪比的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011614196.5A CN112730712B (zh) | 2020-12-31 | 2020-12-31 | 提高lc-ms数据信噪比的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112730712A CN112730712A (zh) | 2021-04-30 |
CN112730712B true CN112730712B (zh) | 2022-06-07 |
Family
ID=75611190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011614196.5A Active CN112730712B (zh) | 2020-12-31 | 2020-12-31 | 提高lc-ms数据信噪比的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112730712B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116136518B (zh) * | 2023-04-20 | 2023-08-01 | 杭州泽天春来科技有限公司 | 色谱仪 |
CN117907511B (zh) * | 2024-03-20 | 2024-06-14 | 浙江灵析精仪科技发展有限公司 | 一种多组分重叠峰的自动化解析方法、装置及电子设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017665A (zh) * | 2012-12-04 | 2013-04-03 | 上海交通大学 | 数字散斑干涉条纹快速滤波*** |
CN103020907A (zh) * | 2012-12-04 | 2013-04-03 | 上海交通大学 | 基于二维集合经验模态分解的dspi条纹滤波*** |
CN103903233A (zh) * | 2014-04-10 | 2014-07-02 | 北京工业大学 | 基于双树离散小波包和信噪比估计的图像去噪方法 |
US9569843B1 (en) * | 2015-09-09 | 2017-02-14 | Siemens Healthcare Gmbh | Parameter-free denoising of complex MR images by iterative multi-wavelet thresholding |
CN106951835A (zh) * | 2017-03-03 | 2017-07-14 | 东华大学 | 一种脑电信号噪声去除方法 |
CN107782706A (zh) * | 2017-09-05 | 2018-03-09 | 燕山大学 | 一种基于集合经验模态分解的多环芳烃荧光光谱去噪方法 |
CN109697320A (zh) * | 2018-12-25 | 2019-04-30 | 北京万维盈创科技发展有限公司 | 拖尾峰处理方法及装置 |
CN111965499A (zh) * | 2020-07-27 | 2020-11-20 | 广东电网有限责任公司广州供电局 | 基于经验模态分解的局部放电数据去噪方法和*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9503054B2 (en) * | 2013-03-14 | 2016-11-22 | Avago Technolgies General IP (Singapore) Pte. Ltd. | Linear phase FIR biorthogonal wavelet filters with complementarity for image noise reduction |
AU2017350807B2 (en) * | 2016-10-25 | 2022-07-07 | Regeneron Pharmaceuticals, Inc. | Methods and systems for chromatography data analysis |
-
2020
- 2020-12-31 CN CN202011614196.5A patent/CN112730712B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103017665A (zh) * | 2012-12-04 | 2013-04-03 | 上海交通大学 | 数字散斑干涉条纹快速滤波*** |
CN103020907A (zh) * | 2012-12-04 | 2013-04-03 | 上海交通大学 | 基于二维集合经验模态分解的dspi条纹滤波*** |
CN103903233A (zh) * | 2014-04-10 | 2014-07-02 | 北京工业大学 | 基于双树离散小波包和信噪比估计的图像去噪方法 |
US9569843B1 (en) * | 2015-09-09 | 2017-02-14 | Siemens Healthcare Gmbh | Parameter-free denoising of complex MR images by iterative multi-wavelet thresholding |
CN106951835A (zh) * | 2017-03-03 | 2017-07-14 | 东华大学 | 一种脑电信号噪声去除方法 |
CN107782706A (zh) * | 2017-09-05 | 2018-03-09 | 燕山大学 | 一种基于集合经验模态分解的多环芳烃荧光光谱去噪方法 |
CN109697320A (zh) * | 2018-12-25 | 2019-04-30 | 北京万维盈创科技发展有限公司 | 拖尾峰处理方法及装置 |
CN111965499A (zh) * | 2020-07-27 | 2020-11-20 | 广东电网有限责任公司广州供电局 | 基于经验模态分解的局部放电数据去噪方法和*** |
Non-Patent Citations (2)
Title |
---|
ENSEMBLE EMPIRICAL MODE DECOMPOSITION:A NOISE-ASSISTED DATA ANALYSIS METHOD;ZHAOHUA WU et al.;《Advances in Adaptive Data Analysis》;20091231;第1卷(第1期);第1-41页 * |
平稳小波去噪算法中的参数选择;张恒 等;《国防科技大学学报》;20190831;第41卷(第4期);第165-170页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112730712A (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112730712B (zh) | 提高lc-ms数据信噪比的方法 | |
CN109883547B (zh) | 一种基于小波阈值差分的宽波段光谱信号去噪方法 | |
CN101647215B (zh) | 用于信号处理***中分析噪声的方法和装置 | |
Neelamani et al. | Wavelet-based deconvolution for ill-conditioned systems | |
CN105279379B (zh) | 基于凸组合核函数主成分分析的太赫兹光谱特征提取方法 | |
Galiana-Merino et al. | De-noising of short-period seismograms by wavelet packet transform | |
CN110763913B (zh) | 一种基于信号分段分类的导数谱平滑处理方法 | |
CN105182417B (zh) | 一种基于形态成分分析的面波分离方法及*** | |
de Oliveira Lyrio et al. | Efficient automatic denoising of gravity gradiometry data | |
CN110646841B (zh) | 时变稀疏反褶积方法及*** | |
CN109724693B (zh) | 一种基于平稳小波的融合光谱去噪方法 | |
CN111505709B (zh) | 一种基于稀疏谱分解的衰减定性分析的方法 | |
CN109670531A (zh) | 一种基于Hodrick-Prescott滤波器的近红外光谱信号的去噪方法 | |
CN112213773B (zh) | 一种地震分辨率提高方法及电子设备 | |
CN112764108B (zh) | 一种基于改进经验小波变换的新型地震资料噪声压制算法 | |
CN107907542A (zh) | 一种ivmd及能量估计相结合的dspi相位滤波方法 | |
CN116165713A (zh) | 一种井中das地震数据去噪方法和装置 | |
Kuang et al. | Joint empirical mode decomposition and optimal frequency band estimation for adaptive low-frequency noise suppression | |
CN113419276B (zh) | 一种自适应相位估算的时变子波提取方法 | |
CN113324918B (zh) | 一种岩石光谱去噪方法 | |
CN109917458B (zh) | 油气分析方法及装置 | |
Lees | Reshaping spectrum estimates by removing periodic noise: Application to seismic spectral ratios | |
CN109884706B (zh) | 非平稳微分加权叠加地震资料处理方法 | |
CN112505779B (zh) | 基于特征分解的采集脚印去除方法 | |
CN112241024B (zh) | 地震资料信噪比提高方法、计算机存储介质及*** |
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 |