CN108918499A - 拉曼图谱中去除拉曼基线漂移的方法 - Google Patents
拉曼图谱中去除拉曼基线漂移的方法 Download PDFInfo
- Publication number
- CN108918499A CN108918499A CN201810683051.7A CN201810683051A CN108918499A CN 108918499 A CN108918499 A CN 108918499A CN 201810683051 A CN201810683051 A CN 201810683051A CN 108918499 A CN108918499 A CN 108918499A
- Authority
- CN
- China
- Prior art keywords
- raman
- trough
- abscissa
- baseline
- map
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/65—Raman scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2201/00—Features of devices classified in G01N21/00
- G01N2201/12—Circuits of general importance; Signal processing
- G01N2201/127—Calibration; base line adjustment; drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明公开了拉曼图谱中去除拉曼基线漂移的方法,包括:第一步:对原始拉曼图谱进行平滑,通过中间值比两边值大,求得平滑后的拉曼波峰,通过中间值比两边值小,求得平滑后的波谷,得到平滑后的拉曼图谱;第二步:对平滑后的拉曼区域进行滑动窗口式的线性拟合,获得拟合后的拉曼基线;第三步:去除拟合后的拉曼基线中的荧光包;第四步:原始拉曼峰图谱减去步骤三获得的去除荧光包的拉曼基线,即得到最终的去除拉曼基线漂移的拉曼图谱。该方法能够保证校正过程中的稳定性,加快基线校正的速度。
Description
技术领域
本发明涉及光学领域中的拉曼基线漂移的处理方法,具体是指拉曼图谱中去除拉曼基线漂移的方法。
背景技术
拉曼光谱在生物分析和检测领域有重要作用。近年来,拉曼光谱被运用于细胞的分类和分析中。这往往需要更加准确的数据来进行算法分析。在拉曼检测中,背景荧光是一种最普遍的背景干扰光,背景荧光外观通常比拉曼峰要宽得多,看起来就像是拉曼光谱缓慢变化的基线。背景荧光会导致拉曼基线发生漂移,从而导致拉曼数据的不精确。
常见的基线校正方法有求导、多项式拟合、小波变换等。Savitzky—Golay导数光谱是比较常用的导数光谱,它是根据一系列给定的光谱点拟合一个多项式,计算数据段的中点处拟合多项式的导数。然而导数光谱用于基线校正会降低信噪比或导致原始光谱变形。多项式拟合的基线校正方法是从光谱中选取拟合数据点,再根据光谱的特点定义合适的拟合阶数,然后根据最小二乘原理计算拟合数据点的最佳逼近函数。然而拟合阶数的难以确定,将会增加运算时间,同时拟合基线会因为拟合阶数的改变而发生振荡。小波变换在荧光包存在问题上,会导致荧光背景的残留。
目前小波变化在拉曼中的基线校正是普遍的,通过小波变化进行高低频的分离,同时采用最小二乘法原理拟合数据点,去除荧光背景对光谱的干扰。这种方法的使用在专利201310094703.0的一种拉曼光谱预处理方法中,也有所提及和使用。
目前的小波变换去除拉曼基线虽然在荧光造成的基线偏移上有一定的作用,但是在荧光包的处理上存在问题,不能很好的去除荧光包。
发明内容
本发明的目的是提供拉曼图谱中去除拉曼基线漂移的方法,该方法能够保证校正过程中的稳定性,加快基线校正的速度。
本发明的上述目的通过如下技术方案来实现:拉曼图谱中去除拉曼基线漂移的方法,其特征在于,该方法包括如下步骤:
第一步:对原始拉曼图谱进行平滑,通过中间值比两边值大,求得平滑后的拉曼波峰,通过中间值比两边值小,求得平滑后的波谷,得到平滑后的拉曼图谱;
第二步:对平滑后的拉曼区域进行滑动窗口式的线性拟合,获得拟合后的拉曼基线;
第三步:去除拟合后的拉曼基线中的荧光包;
第四步:原始拉曼峰图谱减去步骤三获得的去除荧光包的拉曼基线,即得到最终的去除拉曼基线漂移的拉曼图谱。
本发明中,所述第二步的具体过程如下:选择平滑后的拉曼图谱中横坐标为X的波谷,下一个波谷的横坐标为X1,下一个波谷的横坐标为X2,下一个波谷的横坐标为X3,求得横坐标为X的波谷与横坐标为X1的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a1,求得横坐标为X的波谷与横坐标为X2的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a2,比较a1与a2的大小,若a1>a2,则选择横坐标为X的波谷与横坐标为X1的波谷进行线性拟合;若a1<a2,继续向下比较,求得横坐标为X的波谷与横坐标为X3的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a3,若a2>a3,不再继续比较,则选择横坐标为X的波谷与横坐标为X2的波谷进行线性拟合;若a2<a3,则继续向下比较,直至出现a(N)>a(N+1),N为横坐标为X的波谷下的第N个波谷,则选择横坐标为X的波谷与横坐标为X(N)的波谷进行线性拟合,得到拟合后的拉曼基线,继续以横坐标X(N)为波谷起点,重复上述步骤,直到所有波谷都被遍历,获得最终拟合的拉曼基线,线性拟合均采用现有最小二乘法的拟合方法。
所述第三步的具体过程如下:拟合后的拉曼基线内的波谷可能存在荧光包,需要去除,设置一个荧光系数的阖值,如果波谷的荧光系数超过这个阖值,则该波谷会被视为荧光包,
当荧光系数大于阖值时,拟合后的拉曼基线进行上拉来去除荧光包,获得去除荧光包的拉曼基线。
本发明的方法在校正拉曼基线漂移的过程中,不会出现拟合过度或者欠拟合的情况,保证了校正过程中的稳定性,因为是线性拟合,所以不需要确定拟合的阶数,因此大大的加快了基线校正的速度。
附图说明
下面结合附图和具体实施例对本发明做进一步详细说明。
图1是本发明拉曼图谱中去除拉曼基线漂移的方法中对平滑后的拉曼区域进行滑动窗口式的线性拟合,获得拟合后的拉曼基线的图示,为平滑后的拉曼图谱与其拟合基线的复合图,其中,连续曲线为平滑后的拉曼图谱,点划曲线为平滑后的拉曼图谱的拟合基线,通过比较两个波谷与y负半轴连线的余弦值,找出适合进行线性拟合的两个波谷的位置;
图1A是图1中平滑后的拉曼图谱;
图1B是图1中平滑后的拉曼图谱的拟合基线;
图2是本发明拉曼图谱中去除拉曼基线漂移的方法中去除拟合后的拉曼基线中的荧光包的图示;为有荧光包的平滑后的拉曼图谱的拟合基线与其拟合基线的复合图,其中,连续曲线为有荧光包的平滑后的拉曼图谱,点划曲线为有荧光包的平滑后的拉曼图谱的拟合基线,该图为有荧光的平滑后的拉曼图谱的拟合基线效果图,从图中可以看出,平滑后拉曼图谱的荧光包被很好的拟合,而拉曼峰没有出现过拟合和欠拟合的情况;
图2A是图2中有荧光包的平滑后的拉曼图谱;
图2B是图2中有荧光包的平滑后的拉曼图谱的拟合基线;
图3是本发明拉曼图谱中去除拉曼基线漂移的方法中获得去除荧光包的拉曼基线的图示;为有荧光包的拉曼图谱与其拟合基线的复合图,其中,连续曲线为有荧光包的拉曼图谱,点划曲线为有荧光包的拉曼图谱的拟合基线,将获得的基线与拉曼图谱进行作图,从图中可以看出,拉曼图谱的荧光包被很好的拟合,而拉曼峰也没有出现过拟合和欠拟合的情况;
图3A是图3中有荧光包的拉曼图谱;
图3B是图3中有荧光包的拉曼图谱的拟合基线;
图4是本发明拉曼图谱中去除拉曼基线漂移的方法中基线校正以后的拉曼图谱,通过拉曼图谱减去平滑后拉曼图谱的拟合基线从图中可以看出,荧光包被去除,而拉曼峰强度并未出现实质上的损失。
具体实施方式
本发明拉曼图谱中去除拉曼基线漂移的方法,包括如下步骤:
第一步:对原始拉曼图谱进行平滑,通过中间值比两边值大,求得平滑后的拉曼波峰,通过中间值比两边值小,求得平滑后的波谷,得到平滑后的拉曼图谱;
第二步:对平滑后的拉曼区域进行滑动窗口式的线性拟合,获得拟合后的拉曼基线;
选择平滑后的拉曼图谱中横坐标为X的波谷,下一个波谷的横坐标为X1,下一个波谷的横坐标为X2,下一个波谷的横坐标为X3,求得横坐标为X的波谷与横坐标为X1的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a1,求得横坐标为X的波谷与横坐标为X2的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a2,比较a1与a2的大小,若a1>a2,则选择横坐标为X的波谷与横坐标为X1的波谷进行线性拟合;若a1<a2,继续向下比较,求得横坐标为X的波谷与横坐标为X3的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a3,若a2>a3,不再继续比较,则选择横坐标为X的波谷与横坐标为X2的波谷进行线性拟合;若a2<a3,则继续向下比较,直至出现a(N)>a(N+1),N为横坐标为X的波谷下的第N个波谷,则选择横坐标为X的波谷与横坐标为X(N)的波谷进行线性拟合,得到拟合后的拉曼基线,继续以横坐标X(N)为波谷起点,重复上述步骤,直到所有波谷都被遍历,获得最终拟合的拉曼基线,线性拟合均采用现有最小二乘法的拟合方法;
如图1、图1A、图1B所示,以(x=1651)处的波谷为例,我们通过两个波谷的线性拟合作为拉曼基线,该拉曼的下一个波谷为(x=1607),但是这个波谷并不合适作为线性拟合的波谷,因为这时拉的基线会出现欠拟合的情况,而下一个波谷(x=1554)适合作为线性拟合的波谷,下一个波谷(x=1522)并不适合作为线性拟合的波谷,因为拉曼基线会拉到拉曼图谱上方,而出现拉曼基线的过拟合。
判断(x=1651)与哪个波谷进行线性拟合,求得(x=1651)与(x=1607)的直线与y负半轴的Cos记a1(a1为-0.906),求得(x=1651)与(x=1554)的直线与y负半轴的Cos记a2(a2为-0.293),比较a1与a2的大小,a1<a2,继续向下比较,求得(x=1651)与(x=1522)的直线与y负半轴的Cos记为a3(a3为-0.515),因为a2>a3,那么此时不再继续比较,我们选择(x=1651)和(x=1544)处的波谷进行线性拟合。
这样做的好处:可以防止拉曼出现欠拟合,同时又不会因为基线上漂,出现过拟合的情况。能够很好的实现和去除拉曼的基线漂移。
第三步:去除拟合后的拉曼基线中的荧光包;
拟合后的拉曼基线内的波谷可能存在荧光包,需要去除。
如图2、图2A、图2B所示,在拉曼分区后得到的(x=1242)波谷与(x=1406)波谷的区域存在明显的拉曼荧光包,这样的荧光包实际并不存在拉曼峰,但是如果不能通过基线去除的话,会影响到后续的数据分析。
设置一个荧光系数的阖值,数值为1.5,如果波谷的荧光系数超过这个阖值,则该波谷会被视为荧光包,不对称系数的公式中基础高度为0.03。
当荧光系数大于阖值时,拟合后的拉曼基线进行上拉来去除荧光包,获得去除荧光包的拉曼基线,如图3、图3A、图3B所示;
第四步:原始拉曼峰图谱减去步骤三获得的去除荧光包的拉曼基线,即得到最终的去除拉曼基线漂移的拉曼图谱,去除后的拉曼图谱如图4所示。
本发明的上述实施例并不是对本发明保护范围的限定,本发明的实施方式不限于此,凡此种种根据本发明的上述内容,按照本领域的普通技术知识和惯用手段,在不脱离本发明上述基本技术思想前提下,对本发明上述结构做出的其它多种形式的修改、替换或变更,均应落在本发明的保护范围之内。
Claims (3)
1.拉曼图谱中去除拉曼基线漂移的方法,其特征在于,该方法包括如下步骤:
第一步:对原始拉曼图谱进行平滑,通过中间值比两边值大,求得平滑后的拉曼波峰,通过中间值比两边值小,求得平滑后的波谷,得到平滑后的拉曼图谱;
第二步:对平滑后的拉曼区域进行滑动窗口式的线性拟合,获得拟合后的拉曼基线;
第三步:去除拟合后的拉曼基线中的荧光包;
第四步:原始拉曼峰图谱减去步骤三获得的去除荧光包的拉曼基线,即得到最终的去除拉曼基线漂移的拉曼图谱。
2.根据权利要求1所述的拉曼图谱中去除拉曼基线漂移的方法,其特征在于,所述第二步的具体过程如下:选择平滑后的拉曼图谱中横坐标为X的波谷,下一个波谷的横坐标为X1,下一个波谷的横坐标为X2,下一个波谷的横坐标为X3,求得横坐标为X的波谷与横坐标为X1的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a1,求得横坐标为X的波谷与横坐标为X2的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a2,比较a1与a2的大小,若a1>a2,则选择横坐标为X的波谷与横坐标为X1的波谷进行线性拟合;若a1<a2,继续向下比较,求得横坐标为X的波谷与横坐标为X3的波谷相连的直线与y轴负半轴之间的夹角的余弦值记为a3,若a2>a3,不再继续比较,则选择横坐标为X的波谷与横坐标为X2的波谷进行线性拟合;若a2<a3,则继续向下比较,直至出现a(N)>a(N+1),N为横坐标为X的波谷下的第N个波谷,则选择横坐标为X的波谷与横坐标为X(N)的波谷进行线性拟合,得到拟合后的拉曼基线,继续以横坐标X(N)为波谷起点,重复上述步骤,直到所有波谷都被遍历,获得最终拟合的拉曼基线,线性拟合均采用现有最小二乘法的拟合方法。
3.根据权利要求2所述的拉曼图谱中去除拉曼基线漂移的方法,其特征在于,所述第三步的具体过程如下:拟合后的拉曼基线内的波谷可能存在荧光包,需要去除,设置一个荧光系数的阖值,如果波谷的荧光系数超过这个阖值,则该波谷会被视为荧光包,
当荧光系数大于阖值时,拟合后的拉曼基线进行上拉来去除荧光包,获得去除荧光包的拉曼基线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810683051.7A CN108918499B (zh) | 2018-06-28 | 2018-06-28 | 拉曼图谱中去除拉曼基线漂移的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810683051.7A CN108918499B (zh) | 2018-06-28 | 2018-06-28 | 拉曼图谱中去除拉曼基线漂移的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108918499A true CN108918499A (zh) | 2018-11-30 |
CN108918499B CN108918499B (zh) | 2019-08-02 |
Family
ID=64422893
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810683051.7A Active CN108918499B (zh) | 2018-06-28 | 2018-06-28 | 拉曼图谱中去除拉曼基线漂移的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108918499B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109993155A (zh) * | 2019-04-23 | 2019-07-09 | 北京理工大学 | 针对低信噪比紫外拉曼光谱的特征峰提取方法 |
CN110108694A (zh) * | 2019-05-17 | 2019-08-09 | 广西科技大学 | 拉曼光谱仪波数迁移测量精度的提升方法 |
CN113109318A (zh) * | 2021-03-26 | 2021-07-13 | 中国科学院西安光学精密机械研究所 | 基于光谱峰高直接提取的拉曼光谱定量分析方法及*** |
CN114993891A (zh) * | 2022-07-14 | 2022-09-02 | 长光辰英(杭州)科学仪器有限公司 | 基于余弦相似度的颗粒物拉曼检测方法 |
CN116030032A (zh) * | 2023-02-16 | 2023-04-28 | 四川省肿瘤医院 | 基于拉曼光谱数据的乳腺癌分析设备、***及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002003857A2 (en) * | 2000-07-11 | 2002-01-17 | Lightouch Medical, Inc. | Method of tissue modulation for noninvasive measurement of an analyte |
CN103529012A (zh) * | 2013-09-18 | 2014-01-22 | 武汉钢铁(集团)公司 | 一种适用于高炉瓦斯灰中碳来源的拉曼光谱定量检测方法 |
CN105806824A (zh) * | 2016-03-22 | 2016-07-27 | 武汉轻工大学 | 基于多源光谱特征融合的食用油种类快速鉴别方法 |
CN107345911A (zh) * | 2017-06-06 | 2017-11-14 | 浙江大学 | 一种表面增强拉曼光谱定性定量分析米鱼肌肉中组胺的方法 |
-
2018
- 2018-06-28 CN CN201810683051.7A patent/CN108918499B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002003857A2 (en) * | 2000-07-11 | 2002-01-17 | Lightouch Medical, Inc. | Method of tissue modulation for noninvasive measurement of an analyte |
CN103529012A (zh) * | 2013-09-18 | 2014-01-22 | 武汉钢铁(集团)公司 | 一种适用于高炉瓦斯灰中碳来源的拉曼光谱定量检测方法 |
CN105806824A (zh) * | 2016-03-22 | 2016-07-27 | 武汉轻工大学 | 基于多源光谱特征融合的食用油种类快速鉴别方法 |
CN107345911A (zh) * | 2017-06-06 | 2017-11-14 | 浙江大学 | 一种表面增强拉曼光谱定性定量分析米鱼肌肉中组胺的方法 |
Non-Patent Citations (1)
Title |
---|
覃赵军等: "《分段式线性拟合校正拉曼光谱基线漂移》", 《光谱学与光谱分析》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109993155A (zh) * | 2019-04-23 | 2019-07-09 | 北京理工大学 | 针对低信噪比紫外拉曼光谱的特征峰提取方法 |
CN109993155B (zh) * | 2019-04-23 | 2021-01-01 | 北京理工大学 | 针对低信噪比紫外拉曼光谱的特征峰提取方法 |
CN110108694A (zh) * | 2019-05-17 | 2019-08-09 | 广西科技大学 | 拉曼光谱仪波数迁移测量精度的提升方法 |
CN110108694B (zh) * | 2019-05-17 | 2021-11-02 | 广西科技大学 | 拉曼光谱仪波数迁移测量精度的提升方法 |
CN113109318A (zh) * | 2021-03-26 | 2021-07-13 | 中国科学院西安光学精密机械研究所 | 基于光谱峰高直接提取的拉曼光谱定量分析方法及*** |
CN114993891A (zh) * | 2022-07-14 | 2022-09-02 | 长光辰英(杭州)科学仪器有限公司 | 基于余弦相似度的颗粒物拉曼检测方法 |
CN114993891B (zh) * | 2022-07-14 | 2024-04-19 | 长光辰英(杭州)科学仪器有限公司 | 基于余弦相似度的颗粒物拉曼检测方法 |
CN116030032A (zh) * | 2023-02-16 | 2023-04-28 | 四川省肿瘤医院 | 基于拉曼光谱数据的乳腺癌分析设备、***及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108918499B (zh) | 2019-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108918499B (zh) | 拉曼图谱中去除拉曼基线漂移的方法 | |
CN109696788B (zh) | 一种基于显示面板的快速自动对焦方法 | |
CN105117725B (zh) | 一种基于车牌可信度的车牌自纠正识别方法 | |
CN103886553B (zh) | 图像非局部均值去噪方法和*** | |
CN108896278A (zh) | 一种滤光片丝印缺陷检测方法、装置及终端设备 | |
CN110111301A (zh) | 基于频域变换的金属氧化表面缺陷视觉检测方法 | |
CN107369176B (zh) | 一种柔性ic基板氧化面积检测***及方法 | |
WO2023070593A1 (zh) | 线宽测量方法、装置、计算处理设备、计算机程序及计算机可读介质 | |
CN105139391B (zh) | 一种雾霾天气交通图像边缘检测方法 | |
CN110033439B (zh) | 一种制丝过程中皮带输送机堵料的视觉检测方法 | |
CN102663384A (zh) | 基于贝塞尔控制点搜索的曲线识别方法及装置 | |
CN110533616A (zh) | 一种图像天空区域分割的方法 | |
CN107038444A (zh) | 一种指针式表盘的图像识别方法 | |
CN107977989A (zh) | 一种确定图像偏移量的方法及设备 | |
CN108918093A (zh) | 一种滤光片镜面缺陷检测方法、装置及终端设备 | |
CN117541582A (zh) | 一种高频变频器igbt绝缘质量检测方法 | |
CN112750089B (zh) | 基于局部块最大和最小像素先验的光学遥感影像去雾方法 | |
CN114140416A (zh) | 一种基于机器视觉的玻璃边缘检测方法和*** | |
CN117689574A (zh) | 一种用于肿瘤射频消融诊疗定位的医疗影像处理方法 | |
CN110489810B (zh) | 一种基于数据块的趋势自动提取方法 | |
CN108460735A (zh) | 基于单幅图像的改进暗通道去雾方法 | |
CN105957061A (zh) | 一种用于旋冠盖的阈值图制作和图片检测方法 | |
CN109146900A (zh) | 应用于全息投影显示的Blob区域改进提取方法 | |
CN105989351B (zh) | 一种指纹图像背景分割的方法 | |
CN108734180A (zh) | 一种基于计算方式优化的sift特征点梯度生成方法 |
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 |