CN105911046A - 一种基于频数分布原理的滤波方法 - Google Patents
一种基于频数分布原理的滤波方法 Download PDFInfo
- Publication number
- CN105911046A CN105911046A CN201610521973.9A CN201610521973A CN105911046A CN 105911046 A CN105911046 A CN 105911046A CN 201610521973 A CN201610521973 A CN 201610521973A CN 105911046 A CN105911046 A CN 105911046A
- Authority
- CN
- China
- Prior art keywords
- raman spectrum
- group
- data
- raman
- array
- 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
Abstract
本发明涉及一种基于频数分布原理的滤波方法,其包括以下步骤:采集N组拉曼光谱并求其峰值强度存入数组Ix;将Ix分为M个区间;按照分区对数组Ix进行分组,统计各区间频数;求出频数最大的区间K,并求该区间的中值Ixm(K);设定拉曼光谱强度波动最大偏差A,若第N+1组拉曼光谱峰值强度在区间[Ixm(K)‑A,Ixm(K)+A,]内,该拉曼光谱有效,进行预处理后存入寄存器,否则该拉曼光谱无效,舍弃,数组Ix长度固定,每新采集一组拉曼光谱求其峰值强度放入队列尾,扔掉原队列首的一个数据,保持N个数据最新,寄存器内有效拉曼光谱达到L组后进行叠加。本发明滤波效果明显,实时性好,可有效排除瞬态噪声干扰,提高拉曼散射海水温度激光遥感测量系统的抗干扰性和水温测量精度。
Description
技术领域
本发明属于海洋监测领域,具体地说是一种适用于拉曼散射海水温度激光遥感测量系统的基于频数分布原理的滤波方法。
背景技术
海水水下温度测量对研究海洋环境和气候监测及自然灾害的早期预报等十分重要。由于蓝绿激光在水中良好的透射性,用于海水水下温度的遥感测量可以快速获取其水平和垂直分布,具有浮标等定点站位测量和微波辐射或红外辐射等表面温度测量方式不可比拟的优点。目前国内外很多研究人员对拉曼散射海水温度激光遥感测量技术进行了广泛深入的研究,但是海水拉曼散射受波浪、外界光线、水下悬浮物质、CDOM荧光散射等因素影响信号强度波动大,导致水温测量误差较大。虽然现有根据海水拉曼光谱反演海水温度的方法中已经采用光谱强度归一化的方式降低拉曼光谱信号强度波动对水温测量结果的影响,如2015年“Applied Optics”第54卷第10期发表的“Ultimate sensing resolution of watertemperature by remote Raman spectroscopy”和2001年“Journal of the KoreanPhysical Society”第38卷第6期发表的“LIDAR measurement of water temperature byusing frequency-shifted Raman scattering”。但是该方法尚不能很好地消除所述外界瞬态噪声干扰信号对水温测量精度的影响,严重限制了拉曼散射海水温度激光遥感测量技术的实际应用。迄今为止,尚未见针对拉曼光谱受瞬态噪声干扰信号强度波动的相关滤波算法报道。
发明内容
本发明所要解决的技术问题是针对水的拉曼散射光谱受噪声干扰信号强度波动大的问题,提供一种方法简单、合理,实现方便,滤波效果明显,能够有效排除外界噪声信号干扰,提高拉曼散射海水温度激光遥感测量系统的实时测温精度的基于频数分布原理的滤波方法。
本发明解决上述技术问题采用的技术方案是:一种基于频数分布原理的滤波方法,其特征是:其包括以下步骤:
步骤(1):连续采集多组拉曼光谱数据I1(λ)、I2(λ)……,分别求得各组光谱的峰值强度Iimax,把连续N个Iimax值看成一个队列,存入数组Ix;
Ix(i)=[I1max I2max … INmax],i=1,2,…N
步骤(2):将数组Ix划分为M个区间,根据数组Ix中数据的最大值Ixmax和最小值Ixmin来确定分区间距ΔIx;
第j个分区区间的中间值表示为:
Ixm(j)=Ixmin-0.5ΔIx+j·ΔIx,j=1,2,…M
步骤(3):按照分区对数组Ix中的元素进行分组,若Ix(i)属于指定第j个区间,则函数γj[Ix(i)]值为1,否则函数γj[Ix(i)]值为0,即:
统计各分区区间的频数Fj;
步骤(4):根据上述频数分布统计结果求出频数最大值Fmax所处的区间编号K,并求得该区间对应的中间值Ixm(K),作为第N+1组拉曼光谱筛选的参照值;
步骤(5):设定相邻两组拉曼光谱强度波动允许的最大偏差值为A,若第N+1组拉曼光谱数据的峰值强度满足Ixm(K)-A<I(N+1)max<Ixm(K)+A,则认为这一组拉曼光谱有效,进行光谱预处理后存入寄存器用于拉曼光谱叠加,否则认为这一组拉曼光谱无效,光谱数据舍弃,数组Ix的长度固定为N,每采集一组新的拉曼光谱数据并求得其对应的峰值强度Iimax后都放在队列尾部,并扔掉原来队列首部的一个峰值强度数据,保持这N个数据始终是最近的数据,当寄存器内有效拉曼光谱数据达到L组以后,将L组预处理后的拉曼光谱进行叠加。
根据本发明,将L组拉曼光谱叠加后的结果用于光谱特征值提取,进而将光谱特征值用于物质结构及成分鉴别、材料特性分析等,例如可根据所述滤波叠加后的拉曼光谱提取光谱特征值,进而根据光谱特征值计算海水盐度或温度。本发明所述一种基于频数分布原理的滤波方法具有实现方便,滤波效果明显,实时性好的优点,可有效排除瞬态噪声信号的干扰,提高拉曼散射海水温度激光遥感测量系统的抗干扰性和水温测量精度。
附图说明
下面结合附图对本发明作进一步的描述。
图1是拉曼光谱峰值强度的频数分布原理示意图。
图2是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时,灯光长时间干扰基于频数分布原理的滤波方法测试结果。
图3是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时,灯光短时间干扰基于频数分布原理的滤波方法测试结果。
图4是设定相邻两组拉曼光谱强度波动允许的最大偏差值为200时,灯光长时间干扰基于频数分布原理的滤波方法测试结果。
图5是设定相邻两组拉曼光谱强度波动允许的最大偏差值为200时,灯光短时间干扰基于频数分布原理的滤波方法测试结果。
具体实施方式
本发明一种基于频数分布原理的滤波方法,其包括以下具体步骤:
步骤(1):连续采集多组拉曼光谱数据I1(λ)、I2(λ)……,分别求得各组光谱的峰值强度Iimax,把连续N个Iimax值看成一个队列,存入数组Ix;
Ix(i)=[I1max I2max … INmax],i=1,2,…N (1)
式(1)中Ix(i)是数组Ix的第i个元素。
步骤(2):将数组Ix划分为M个区间,根据数组Ix中数据的最大值Ixmax和最小值Ixmin来确定分区间距ΔIx;
式(2)中Ixmax是输入数组Ix中的最大强度值,Ixmin是数组Ix中的最小强度值,每个区间的中间值可以表示为:
Ixm(j)=Ixmin-0.5ΔIx+j·ΔIx,j=1,2,…M (3)
式(3)中Ixm(j)是第j个区间的中间值。
步骤(3):按照分区对数组Ix中的每个元素进行分组,若Ix(i)属于指定第j个区间,则函数γj[Ix(i)]值为1,否则函数γj[Ix(i)]值为0,即
式(4)中γj[Ix(i)]反映数组Ix中第i个元素是否属于第j个分区的情况。统计各分区区间的频数Fj;
式(5)中Fj是数组Ix中每个元素在第j个分区区间的频数。频数分布统计结果如图1所示,直接反映拉曼光谱峰值强度的变化情况。
步骤(4)根据所述频数分布统计结果求出频数最大值Fmax所处的区间编号K,并求得该区间对应的中间值Ixm(K),作为第N+1组拉曼光谱筛选的参照值。
步骤(5)设定相邻两组拉曼光谱强度波动允许的最大偏差值为A,若第N+1组拉曼光谱数据的峰值强度满足Ixm(K)-A<I(N+1)max<Ixm(K)+A,则认为这一组拉曼光谱有效,进行光谱预处理后存入寄存器用于拉曼光谱叠加,否则认为这一组拉曼光谱无效,光谱数据舍弃,数组Ix的长度固定为N,每采集一组新的拉曼光谱数据并求得其对应的峰值强度Iimax后都放在队列尾部,并扔掉原来队列首部的一个峰值强度数据,保持这N个数据始终是最近的数据,当寄存器内有效拉曼光谱数据达到L组以后,将L组预处理后的拉曼光谱进行叠加。
根据本发明,将L组拉曼光谱叠加后的结果用于光谱特征值提取,进而将光谱特征值用于物质结构及成分鉴别、材料特性分析等,例如可根据所述滤波叠加后的拉曼光谱提取光谱特征值,进而根据光谱特征值计算海水盐度或温度。本发明所述一种基于频数分布原理的滤波方法具有实现方便,滤波效果明显,实时性好的优点,可有效排除瞬态噪声信号的干扰,提高拉曼散射海水温度激光遥感测量系统的抗干扰性和水温测量精度。
为了验证所述一种基于频数分布原理的滤波方法的抗干扰能力,在实验室采用灯光扰动的方式进行了测试。测试中首先实时采集20组拉曼光谱,求得其对应的光强极大值Iimax,存入数组Ix。然后将Ix中20个元素分为4个区间,统计拉曼光强在各个区间的分布,以频数最高的一个区间K的中间值Ixm(K)作为下一组拉曼光谱筛选的参照值。当设定相邻两组拉曼光谱强度波动允许的最大偏差值A=500时,若第21组拉曼光谱的峰值强度满足Ixm(K)-500<I21max<Ixm(K)+500,则认为这一组拉曼光谱有效,进行预处理后存入寄存器用于拉曼光谱叠加;否则认为这一组拉曼光谱无效,舍弃该组拉曼光谱。数组Ix的长度固定为20,每次采集到一组新的拉曼光谱并求得其对应的峰值强度Iimax后都放在队列尾部,并扔掉原来队列首部的一个峰值强度数据,保持这20个数据始终是最新的数据。最后,当寄存器内有效拉曼光谱数据达到5组以后,将5组预处理后的拉曼光谱进行叠加,计算光谱特征值,进而计算水温并实时显示。图2是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时,灯光长时间干扰基于频数分布原理的滤波方法测试结果,其中虚线是高精度温度传感器的测温结果,实线是拉曼散射海水温度激光遥感测量系统测温结果,t1处产生的尖峰是由于系统运行过程中打开室内照明灯8s,然后关闭,这一扰动导致拉曼光谱强度增大,从而改变了测温结果,t2处的产生的尖峰对应开灯时间为7s。图3是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时,灯光短时间干扰基于频数分布原理的滤波方法测试结果,其中虚线是高精度温度传感器的测温结果,实线是拉曼散射海水温度激光遥感测量系统测温结果,测量过程中多次开关照明灯,每次照明时间缩短到5s以内,由于频数分布滤波算法的滤波作用,干扰信号对温度测量结果没有影响。
图4是当设定相邻两组拉曼光谱强度波动允许的最大偏差值A=200时灯光长时间干扰的测试结果,其中虚线是温度传感器同步测量的结果,实线是拉曼散射海水温度激光遥感测量系统测温结果,t3处产生的尖峰对应照明时间为7s,而t4处产生的尖峰对应开灯时间为10s。由于灯光扰动时间长导致拉曼光谱强度增大,从而改变了测温结果。图5是当设定相邻两组拉曼光谱强度波动允许的最大偏差值A=200时灯光短时间干扰的测试结果,其中虚线是温度传感器同步测量的结果,实线是拉曼散射海水温度激光遥感测量系统测温结果,测量过程中多次开关照明灯,每次照明时间缩短到5s以内,由于频数分布滤波算法的滤波作用,干扰信号对温度测量结果没有影响。
在上述说明中参照着附图,通过具体实施方式适当且充分地对本发明进行了说明,但是应该认识到,对本领域技术人员来说,能够容易的对上述的实时方式进行变更和/或改良。因此对于本领域技术人员所实施的变更方式或改良方式,只要是不脱离权利要求书中记载的权利要求的权力范围的层次,则都应被解释为被包括在该权利要求的权力范围内。
Claims (1)
1.一种基于频数分布原理的滤波方法,其特征是:其包括以下步骤:
步骤(1):连续采集多组拉曼光谱数据I1(λ)、I2(λ)......,分别求得各组光谱的峰值强度Iimax,把连续N个Iimax值看成一个队列,存入数组Ix;
Ix(i)=[I1max I2max … INmax],i=1,2,…N
步骤(2):将数组Ix划分为M个区间,根据数组Ix中数据的最大值Ixmax和最小值Ixmin来确定分区间距ΔIx;
第j个分区区间的中间值表示为:
Ixm(j)=Ixmin-0.5ΔIx+j·ΔIx,j=1,2,…M
步骤(3):按照分区对数组Ix中的元素进行分组,若Ix(i)属于指定第j个区间,则函数γj[Ix(i)]值为1,否则函数γj[Ix(i)]值为0,即:
统计各分区区间的频数Fj;
步骤(4):根据上述频数分布统计结果求出频数最大值Fmax所处的区间编号K,并求得该区间对应的中间值Ixm(K),作为第N+1组拉曼光谱筛选的参照值;
步骤(5):设定相邻两组拉曼光谱强度波动允许的最大偏差值为A,若第N+1组拉曼光谱数据的峰值强度满足Ixm(K)-A<I(N+1)max<Ixm(K)+A,则认为这一组拉曼光谱有效,进行光谱预处理后存入寄存器用于拉曼光谱叠加,否则认为这一组拉曼光谱无效,光谱数据舍弃,数组Ix的长度固定为N,每采集一组新的拉曼光谱数据并求得其对应的峰值强度Iimax后都放在队列尾部,并扔掉原来队列首部的一个峰值强度数据,保持这N个数据始终是最近的数据,当寄存器内有效拉曼光谱数据达到L组以后,将L组预处理后的拉曼光谱进行叠加。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610521973.9A CN105911046B (zh) | 2016-07-04 | 2016-07-04 | 一种基于频数分布原理的滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610521973.9A CN105911046B (zh) | 2016-07-04 | 2016-07-04 | 一种基于频数分布原理的滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105911046A true CN105911046A (zh) | 2016-08-31 |
CN105911046B CN105911046B (zh) | 2018-09-04 |
Family
ID=56753916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610521973.9A Expired - Fee Related CN105911046B (zh) | 2016-07-04 | 2016-07-04 | 一种基于频数分布原理的滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105911046B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106886684A (zh) * | 2017-02-24 | 2017-06-23 | 科博达技术有限公司 | 数字量的等级确定方法 |
CN110530841A (zh) * | 2019-07-16 | 2019-12-03 | 河北伊诺光学科技股份有限公司 | 一种谱图有效性验证方法、系统、电子设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5765948A (en) * | 1995-03-07 | 1998-06-16 | Kabushiki Kaisha Toshiba | Light-temperature distribution sensor using back scattering light produced by incident light pulse and temperature distribution measuring method |
US20090324498A1 (en) * | 2008-06-26 | 2009-12-31 | Sony Corporation | Method of measuring quantity of in-vivo substance by use of coherent anti-stokes raman scattered light |
CN101825567A (zh) * | 2010-04-02 | 2010-09-08 | 南开大学 | 一种近红外光谱和拉曼光谱波长的筛选方法 |
CN103217409A (zh) * | 2013-03-22 | 2013-07-24 | 重庆绿色智能技术研究院 | 一种拉曼光谱预处理方法 |
WO2013183271A1 (ja) * | 2012-06-05 | 2013-12-12 | パナソニック株式会社 | 信号処理装置 |
CN105700019A (zh) * | 2016-02-01 | 2016-06-22 | 电子科技大学 | 一种基于Born-Jordan时频分布的地震信号时频峰值滤波方法 |
-
2016
- 2016-07-04 CN CN201610521973.9A patent/CN105911046B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5765948A (en) * | 1995-03-07 | 1998-06-16 | Kabushiki Kaisha Toshiba | Light-temperature distribution sensor using back scattering light produced by incident light pulse and temperature distribution measuring method |
US20090324498A1 (en) * | 2008-06-26 | 2009-12-31 | Sony Corporation | Method of measuring quantity of in-vivo substance by use of coherent anti-stokes raman scattered light |
CN101825567A (zh) * | 2010-04-02 | 2010-09-08 | 南开大学 | 一种近红外光谱和拉曼光谱波长的筛选方法 |
WO2013183271A1 (ja) * | 2012-06-05 | 2013-12-12 | パナソニック株式会社 | 信号処理装置 |
CN103217409A (zh) * | 2013-03-22 | 2013-07-24 | 重庆绿色智能技术研究院 | 一种拉曼光谱预处理方法 |
CN105700019A (zh) * | 2016-02-01 | 2016-06-22 | 电子科技大学 | 一种基于Born-Jordan时频分布的地震信号时频峰值滤波方法 |
Non-Patent Citations (2)
Title |
---|
XIUYUN REN ET AL.: "Theoretical and experimental investigations on measuring underwater temperature by the coherent Brillouin scattering method", 《APPLIED OPTICS》 * |
王翥等: "基于 WSN的温室大棚温湿度监测系统的设计", 《仪表技术与传感器》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106886684A (zh) * | 2017-02-24 | 2017-06-23 | 科博达技术有限公司 | 数字量的等级确定方法 |
CN106886684B (zh) * | 2017-02-24 | 2019-06-04 | 科博达技术股份有限公司 | 数字量的等级确定方法 |
CN110530841A (zh) * | 2019-07-16 | 2019-12-03 | 河北伊诺光学科技股份有限公司 | 一种谱图有效性验证方法、系统、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105911046B (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kerns et al. | Large‐scale precipitation tracking and the MJO over the Maritime Continent and Indo‐Pacific warm pool | |
Quinn | Analysis of temporal changes in fish assemblages in Serpentine Creek, Queensland | |
Mazzocchi et al. | Stability and resilience in coastal copepod assemblages: The case of the Mediterranean long-term ecological research at Station MC (LTER-MC) | |
CN102252973B (zh) | 土壤含水量遥感监测方法 | |
Ghose et al. | Rain-fed rice yield fluctuation to climatic anomalies in Bangladesh | |
CN103411846B (zh) | 基于高光谱技术的叶面降尘量测定方法 | |
Crowley et al. | Validation of coral temperature calibrations | |
CN205910119U (zh) | 基于太赫兹波的植物叶片含水量检测系统 | |
CN103868860A (zh) | 一种基于高光谱植被指数监测湿地植被冠层氮浓度的方法 | |
Skok et al. | Analysis of tropical cyclone precipitation using an object-based algorithm | |
CN106053370A (zh) | 一种基于hico模拟的反演近海岸水体透明度的方法 | |
CN108876172B (zh) | 一种基于改进型modis植被供水指数的地表土壤含水量评估方法 | |
CN105353371B (zh) | 基于ar谱扩展分形的海面雷达目标检测方法 | |
CN108286962A (zh) | 一种水环境光谱库的建立方法及系统 | |
CN113553907A (zh) | 一种基于遥感技术的森林生态环境状况评价方法 | |
CN105911046A (zh) | 一种基于频数分布原理的滤波方法 | |
CN109060771A (zh) | 一种基于光谱不同特征集的共识模型构建方法 | |
Mandal et al. | Hilbert-Huang transform analysis of surface wavefield under tropical cyclone Hudhud | |
Lee et al. | The structure and dynamics of the stratospheric northern annular mode in CMIP5 simulations | |
Ling et al. | Remote sensing estimation of colored dissolved organic matter (CDOM) from GOCI measurements in the Bohai Sea and Yellow Sea | |
Krupica et al. | The utility of body size indices derived from optical plankton counter data for the characterization of marine zooplankton assemblages | |
Xu et al. | Algorithmic characterization of lake stratification and deep chlorophyll layers from depth profiling water quality data | |
Lamouroux et al. | Objective assessment of the contribution of the RECOPESCA network to the monitoring of 3D coastal ocean variables in the Bay of Biscay and the English Channel | |
Lee et al. | Mixing depth estimation from operational JMA and KMA wind-profiler data and its preliminary applications: Examples from four selected sites | |
Blondel et al. | Acoustic soundscapes and biodiversity–comparing metrics, seasons and depths with data from the Neptune Ocean observatory offshore British Columbia |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180904 Termination date: 20190704 |
|
CF01 | Termination of patent right due to non-payment of annual fee |