CN106908169B - 一种分布式光纤温度传感器数据处理方法 - Google Patents
一种分布式光纤温度传感器数据处理方法 Download PDFInfo
- Publication number
- CN106908169B CN106908169B CN201710114234.2A CN201710114234A CN106908169B CN 106908169 B CN106908169 B CN 106908169B CN 201710114234 A CN201710114234 A CN 201710114234A CN 106908169 B CN106908169 B CN 106908169B
- Authority
- CN
- China
- Prior art keywords
- temperature
- edge
- continuous
- data
- edge detection
- 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
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
Abstract
本发明公开了一种分布式光纤温度传感器数据处理方法,属于传感器技术领域,包括:在利用滤波器模板对分布式温度传感器测量的连续温度数据进行滤波处理之前,对连续温度数据进行边缘检测;在边缘检测结果为温度突变处时,滤波器模板在进行平滑滤波时跳过此边缘。本发明通过在对分布式光纤温度传感器测量的数据进行平滑滤波处理之前,先对温度数据中的边缘进行检测,然后在平滑滤波时跳过每个边缘。这就避免了对边缘位置使用平滑时,使得边缘位置变得平缓,从而降低了系统的空间定位精准度和测温的精准度。
Description
技术领域
本发明涉及传感器技术领域,特别涉及一种分布式光纤温度传感器数据处理方法。
背景技术
分布式光纤温度传感器是一种实时连续测量沿光纤方向温度场分布的技术,该技术大多采用光纤中的拉曼背向散射光强对温度的敏感来实现温度的调节。但是由于拉曼散射光非常微弱,目前,普遍采用累加平均的方法来提供高信噪比,理论上可以认为,只要进行足够多的平均次数,总能够使噪声水平控制在可接受的范围内。但是在实际应用中,由于温度测量的实时性对于分布式光纤传感来说也是一个重要的参数,在平均次数过多的情况下,温度测量的实时性也会变差。因此,在实际应用中,一般通过信号累加平均和信号滤波相结合的方式来提高整个系统的信噪比。
一般对于数字信号来说,经典的滤波方法包括:傅立叶变换滤波、短时傅立叶变换滤波、小波变换滤波和平滑滤波等。对于分布式传感,目前广泛使用的只有频域中的小波变换滤波而不采用空域中的平滑滤波。因为空域中的平滑滤波虽然能够有效的去除噪声的影响,但是具有降低分布式温度数据中包含的空间分辨率、增加部分温度数据的测量误差这些较为严重的缺陷。
但是,长期的实践表明,虽然数字信号在频域中的小波变换滤波方法能有效的改善系统的信噪比,但是该方法也存在一些明显的缺陷:一是,该方法需要针对不同性能的传感器进行单独的调试,实现及运算过程相对复杂。二是,该方法难以在嵌入式系统中进行数据处理的相关硬件中实现,其应用受到了极大的限制。
发明内容
本发明的目的在于提供一种分布式光纤温度传感器数据处理方法,以保证在不降低系统定位精确度和测温精确度的条件下简单的实现平滑滤波。
为实现以上目的,本发明采用的技术方案为:提供一种分布式光纤温度传感器数据处理方法,包括:
在利用滤波器模板对分布式温度传感器测量的连续温度数据进行滤波处理之前,对连续温度数据进行边缘检测;
在边缘检测结果为温度突变处时,滤波器模板在进行平滑滤波时跳过此边缘。
与现有技术相比,本发明存在以下技术效果:本发明通过在对分布式光纤温度传感器测量的数据进行平滑滤波处理之前,先对温度数据中的边缘进行检测,然后在平滑滤波时跳过每个边缘。这就避免了对边缘位置使用平滑时,使得边缘位置变得平缓,从而降低了系统的空间定位精准度和测温的精准度。
附图说明
图1是本发明一实施例提供的一种分布式光纤温度传感器数据处理方法的流程示意图;
图2是本发明一实施例中利用滤波器模板进行空间滤波的示意图;
图3是本发明一实施例中的滤波器模板示意图;
图4是本发明一实施例中利用现有平滑滤波方法对温度测量数据进行滤波的实际效果图;
图5是本发明一实施例中利用改进后的平滑滤波方法对温度测量数据进行滤波的实际效果图。
具体实施方式
下面结合图1至图5所示,对本发明做进一步详细叙述。
如图1所示,本实施例公开了一种分布式光纤温度传感器数据处理方法,该方法包括如下步骤S1至S2:
S1、在利用滤波器模板对分布式温度传感器测量的连续温度数据进行滤波处理之前,对连续温度数据进行边缘检测;
S2、在边缘检测结果为温度突变处时,滤波器模板在进行平滑滤波时跳过此边缘。
进一步地,在步骤S1中,对连续温度数据进行边缘检测,具体包括如下步骤:
对连续温度数据中的相邻不同温度点进行做差处理;
当存在至少两个连续相邻不同温度点之间的差值的绝对值均大于预设阈值时,则判断该至少两个连续相邻不同温度点之间为温度突变处。
其中,在实际应用中,滤波器模板是采用自由游走的方式,也就是说模板的游走方式是不受限制的。如果分布式光纤温度传感器连续测量的两个不同温度区域时,必然会在温度分布图中有一个温度上升沿或者温度下降沿,温度上升沿或者温度下降沿称为温度突变处。而本实施例中在采用滤波器模板进行平滑滤波时,跳过这个温度上升沿或者下降沿,因此,能够保留温度信号的上升沿和下降沿的大小不变,提高了平滑滤波后的温度数据的真实性。
进一步地,本实施例中的边缘检测方法包括但不限于作差法。比如本实施例中进行边缘检测的方法还可以为对两个连续相邻不同温度点间的曲线做微分求导,然后根据边缘斜率的大小判断此处是否为温度突变处。其中,当斜率趋于零时,该两个连续相邻不同温度点之间不存在上升沿或者下降沿。在斜率大于零或者小于零时,该两个连续相邻不同温度点之间存在上升沿或者下降沿。因此,本实施例只要保护的是在平滑滤之前对温度数据中的边缘进行检测的思想,本领域技术人员可以根据实际情况采用不同的边缘检测方法。
进一步地,在步骤S2之前,还包括:
判断边缘检测结果是否为温度突变处,如果是,则执行步骤S2;
如果否,则利用滤波器模板对此边缘进行平滑滤波。
需要说明的是,在实际应用中,温度数据中的噪声信号本身就是无规律变化的,噪声信号之间也可能存在信号突变。因此,需要对边缘检测结果进行判断识别,以区分当前检测到的信号突变是温度信号突变还是噪声信号突变。
进一步地,判断边缘检测的结果是否为温度突变处,具体包括:
将边缘检测结果与预设的精度阈值进行比较;
在边缘检测结果大于精度阈值时,则此边缘为温度突变处;
在边缘检测结果小于或等于精度阈值时,则此边缘为噪声突变处。
其中,本实施例中预设的精度阈值是一个用来与边缘检测结果进行比较的值,本实施例中预设的精度阈值的具体指为0.8。
进一步地,上述的步骤S2中,还包括,在两个连续相邻温度点之间不存在上升沿或者下降沿时,采用滤波器模板对连续的温度数据进行平滑滤波处理。
具体地,在分布式温度数据中模板就像一个小标尺,标尺的刻度可以有3个、5个或者更大,刻度的间距对应分布式温度数据中的两点间距,平滑的过程是将分布式温度数据中的某个值用其相邻的值进行改变,具体步骤为:
(1)将模板在分布式温度数据中漫游,并将尺子的中心与分布式温度中的某个数据重合;
(2)将模板中的各个系数与对应的数据值相乘,系数均在[0,1]范围内且系数的和为1;
(3)将所有的乘积累加;
(4)将上述的累加结果赋值给模板中心位置的数据。
如图2所示,图a是分布式光纤温度传感器测量的一部分连续数据,将图b是一个刻度为5的模板,模板系数在模板内标出,如果将k0所在的位置与分布式温度数据T5重合,如图c所示:那么模板中对数据处理后输出的结果为:
R=k0T5+k1T6+k2T4+k3T3+k4T7;
将R取代T5的数据就完成了该温度数据的滤波,同样对原始温度数据中的每个数据都按照这种方式进行处理,就可以实现对所有温度数据的平滑滤波。
进一步地,将采用现行平滑滤波方法对分布式温度数据进行滤波的效果与采用本实施例公开的平滑滤波方法对分布式温度数据进行滤波的效果进行比较如下:
图4为采用现行平滑滤波方法,利用图3所示的滤波器模板对一次恒温箱温度测量数据进行滤波的效果图,从图4中可以看出,这种方法对温度数据起到了很好的平滑作用,但是同时使得温度数据的上升沿和下降沿变得平缓,间接的影响了分布式温度传感系统的空间分辨率,同时使得上升沿和下降沿附近的温度测量结果误差变大。
图5为采用本实施例公开的方法对分布式温度数据进行滤波处理后的效果图,从图中可以看出,该方法不仅能对温度数据起到很好的平滑作用,还能保留温度信号在上升沿和下降沿的大小不变,保证了分布式温度传感系统的空间分辨率和温度测量结果的准确性。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种分布式光纤温度传感器数据处理方法,其特征在于,包括如下步骤S1至S2:
S1、在利用滤波器模板对分布式温度传感器测量的连续温度数据进行滤波处理之前,对连续温度数据进行边缘检测;
S2、在根据边缘检测结果判断有至少两个连续相邻不同温度点之间为温度突变处时,滤波器模板在进行平滑滤波时跳过此边缘。
2.如权利要求1所述的方法,其特征在于,步骤S1中的对连续温度数据进行边缘检测,具体包括:
对连续温度数据中的相邻不同温度点进行做差处理;
当存在至少两个连续相邻不同温度点之间的差值的绝对值大于预设阈值时,则判断该至少两个连续相邻不同温度点之间为疑似温度突变处。
3.如权利要求1或2所述的方法,其特征在于,在所述步骤S2之前,还包括:
判断边缘是否为温度突变处,如果是,则执行步骤S2;
如果否,则利用滤波器模板对此边缘进行平滑滤波。
4.如权利要求3所述的方法,其特征在于,所述的判断边缘是否为温度突变处,具体包括:
将边缘检测结果与预设的精度阈值进行比较;
在边缘检测结果大于精度阈值时,则此边缘为温度突变处;
在边缘检测结果小于或等于精度阈值时,则此边缘为噪声突变处。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710114234.2A CN106908169B (zh) | 2017-02-28 | 2017-02-28 | 一种分布式光纤温度传感器数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710114234.2A CN106908169B (zh) | 2017-02-28 | 2017-02-28 | 一种分布式光纤温度传感器数据处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106908169A CN106908169A (zh) | 2017-06-30 |
CN106908169B true CN106908169B (zh) | 2020-03-17 |
Family
ID=59208421
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710114234.2A Active CN106908169B (zh) | 2017-02-28 | 2017-02-28 | 一种分布式光纤温度传感器数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106908169B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111879389B (zh) * | 2020-06-30 | 2022-02-08 | 厦门汉印电子技术有限公司 | 称重方法、介质及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10318853A (ja) * | 1997-05-20 | 1998-12-04 | Pub Works Res Inst Ministry Of Constr | 路面温度測定方法 |
WO2011049630A2 (en) * | 2009-10-23 | 2011-04-28 | SensorTran, Inc | Stimulated brillouin system with multiple fbg's |
CN202736125U (zh) * | 2012-06-15 | 2013-02-13 | 黄淮学院 | 一种医学图像的边缘检测系统 |
CN105953941A (zh) * | 2016-04-29 | 2016-09-21 | 深圳艾瑞斯通技术有限公司 | 一种基于拉曼散射的分布式光纤测温方法及装置 |
CN106441386A (zh) * | 2016-09-29 | 2017-02-22 | 西南交通大学 | 基于分布式光纤传感系统的数据处理方法及装置 |
-
2017
- 2017-02-28 CN CN201710114234.2A patent/CN106908169B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10318853A (ja) * | 1997-05-20 | 1998-12-04 | Pub Works Res Inst Ministry Of Constr | 路面温度測定方法 |
WO2011049630A2 (en) * | 2009-10-23 | 2011-04-28 | SensorTran, Inc | Stimulated brillouin system with multiple fbg's |
CN202736125U (zh) * | 2012-06-15 | 2013-02-13 | 黄淮学院 | 一种医学图像的边缘检测系统 |
CN105953941A (zh) * | 2016-04-29 | 2016-09-21 | 深圳艾瑞斯通技术有限公司 | 一种基于拉曼散射的分布式光纤测温方法及装置 |
CN106441386A (zh) * | 2016-09-29 | 2017-02-22 | 西南交通大学 | 基于分布式光纤传感系统的数据处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106908169A (zh) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109239360B (zh) | 一种反应曲线异常检测方法及装置 | |
JP6336982B2 (ja) | その場測定のための光計測 | |
Son et al. | Vibration displacement measurement technology for cylindrical structures using camera images | |
US8594975B2 (en) | Systems and methods for wafer edge feature detection and quantification | |
CN112985639B (zh) | 基于去啁啾和时域子啁啾脉冲提取的分布式光纤传感方法 | |
KR20170134309A (ko) | 기준 광섬유 간섭계 데이터와 함께 측정된 광섬유 간섭계 데이터의 등록 | |
CN108896274B (zh) | 一种基于子集窗长度优化算法的分布式光纤应变解调方法 | |
EP2985566B1 (en) | Data generation method and data generation apparatus | |
CN113838139B (zh) | 图像传感器的参数检测方法、装置、电子设备及存储介质 | |
CN112345759A (zh) | 一种用于荧光强度峰检测的方法 | |
CN106908169B (zh) | 一种分布式光纤温度传感器数据处理方法 | |
CN109696702B (zh) | 一种海水放射性核素k40检测的重叠峰判断方法 | |
CN1713744B (zh) | 用于测量图像信号中的噪声的装置及其方法 | |
CN117168337B (zh) | Ofdr应变边缘优化方法及测量方法 | |
CN109669205B (zh) | 一种海水放射性核素k40元素的寻峰方法 | |
US6738511B1 (en) | Reduced noise sensitivity method and apparatus for converting an interferogram phase map to a surface profile map | |
JP7209960B2 (ja) | 振動特性検出装置、および振動特性検出方法 | |
CN109870500B (zh) | 一种基于交流电磁场检测的缺陷实时判别方法及系统 | |
CN112034268A (zh) | 一种用于消除电荷累积效应测量空间电场的方法及系统 | |
CN105548935A (zh) | 磁场测量仪分辨率的检测方法和装置 | |
CN111751609A (zh) | 一种用于50Hz轨道信号检测及状态鉴别的方法及系统 | |
US6560547B1 (en) | Real time sampling system and method for measuring an interrupted surface | |
KR102305190B1 (ko) | 샘플 표면의 3d 단차 측정 장치 | |
CN115077728B (zh) | 一种多波长检测方法、装置、系统 | |
CN111272060B (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 |