CN101539402B - 实时测量数据的处理方法 - Google Patents

实时测量数据的处理方法 Download PDF

Info

Publication number
CN101539402B
CN101539402B CN2009100686693A CN200910068669A CN101539402B CN 101539402 B CN101539402 B CN 101539402B CN 2009100686693 A CN2009100686693 A CN 2009100686693A CN 200910068669 A CN200910068669 A CN 200910068669A CN 101539402 B CN101539402 B CN 101539402B
Authority
CN
China
Prior art keywords
data
measured object
standard measured
value
sensor
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
Application number
CN2009100686693A
Other languages
English (en)
Other versions
CN101539402A (zh
Inventor
陈炳生
刘宝瑛
周涛
戴清华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Eastloong Opto-Electronic Measuring & Control Technology Co Ltd
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN2009100686693A priority Critical patent/CN101539402B/zh
Publication of CN101539402A publication Critical patent/CN101539402A/zh
Application granted granted Critical
Publication of CN101539402B publication Critical patent/CN101539402B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)

Abstract

一种实时测量数据的处理方法,它是采用系统误差和随机误差分类处理,并逐一修正的四步数据修正方法。本发明采用高阶多项式的数学模型来对被测物处在径向不同位置上的测量数据进行非线性修正;采用一阶多项式的数学模型将被测物处在相关方向不同位置上的测量数据进行偏离修正;为消除系统误差中固定不变的部分而采用常数修正;为了消除随机误差而利用过往数据对实时数据进行滚动滤波修正。本发明的处理方法保证了对数据的实时响应,又去除了以往数据处理方法对随机误差不能有效消除的弱点,以较低的处理工作量,实现了高效、高稳定的对实时测量数据的误差修正。

Description

实时测量数据的处理方法
技术领域
本发明属于数据处理领域,涉及一种针对实时测量数据的误差和干扰的处理方法。
背景技术
工业现场的多种测量仪器,如光电测量、超声测量、激光测量等,在线测量具有数据量大、易受干扰、变化规律难以预测等特点。对于这类数据的一般处理方法比较简单,仅采用大数平均法,即对所获得的数据依照一定的时间间隔进行平均处理;或者是使用统计学中的标准差等概念进行衡量和评价。经过平均处理的数据虽可以减小干扰的影响,减少一部分随机误差,但实时性差,不能反映出测量数据的真实情况。而标准差等统计概念对数据的评价也是不全面的,更不能实现对有用数据中包含的噪声和各种干扰因素的处理。
发明内容
本发明的目的是为了克服上述的缺陷,提出一种实时测量数据的处理方法,它是将系统误差和随机误差进行分类后逐一修正的数据处理方法。
本发明的解决方案是:
一种实时测量的数据处理方法,其特征在于:它包括以下步骤:
一、校准:
1、使用至少两路有重合区域的传感器对同一标准被测物进行测量;
2、获取修正系数:
使用导轨将标准被测物置于上述的重合区域测量区,并在被修正传感器径向上能够使得标准被测物移动固定的步长,采集每次移动后静止时被修正传感器测得的标准被测物两侧边缘值x1、x2,相减获得每次移动后静止时的标准被测物外径x,由于标准被测物的实际外形尺寸y是已知的,根据移动步长指定标准被测物左侧边缘值p1,并指定右侧边缘值p2为左侧边缘值p1与y的和,将这两组数据和修正前的两组位置值数据x1、x2,代入公式(1),并根据公式(1),列出超定方程组,求解得到非线性修正系数a1i、a2i
p 1 = Σ i = 0 n a 1 i x 1 i p 2 = Σ i = 0 n a 2 i x 2 i - - - ( 1 )
二、修正:
1、非线性修正:
采用高阶多项式的数学模型来对标准被测物处在径向不同位置上的测量数据x进行修正;使用上述步骤求得的非线性系数,将测得的标准被测物两侧边缘值x1、x2代入式(2),获得修正后的位置值p1’、p2’,将修正后的位置值相减得到标准被测物直径的测量数据x’;
p 1 ' = Σ i = 0 n a 1 i x 1 i p 2 ' = Σ i = 0 n a 2 i x 2 i - - - ( 2 )
2、动态偏离:
采用一阶多项式的数学模型,将标准被测物处在法向平行位置上的测量数据进行修正;
通过计算在位置关系上与被修正传感器相关方向放置的其它组传感器的测量数据xv,和标准被测物在被修正传感器法向测量区中间位置的测量值xm,还有被修正传感器相关方向不同位置的测量误差平均值Ee及被修正传感器测量值x’,获得偏离修正后的测量值x”为:
x″=(xv-xm)·Ee+x′        (2)
3、校零:
校零是通过数据零点的移动来进行固定的系统误差的补偿,目的是为了降低系统误差中固定不变的部分,将标准被测物放置于测量区中心位置,采集被修正传感器的测量数据,由被修正传感器经过动态偏离修正后的测量数据x”和标准被测物的已知实际外形尺寸y,根据公式(4),计算得到校零系数c,在测量时,便将校零系数作为固定的系统误差,加入测量数据;
c=y-x″        (4)
4、滚动滤波:
滚动滤波是为了消除随机误差所引起的野值而设计的,原理是通过在时间上连续的前n个时刻数据的平均值来设定滤波门限,对当前时刻采集的数据进行判断筛选,如果在门限之内,则将数据送入后续处理过程,否则重新采集数据进行判断。
本发明具有以下优点:它解决了工业现场高温或温度变化剧烈情况下的对在线产品宽度、厚度的高精密测量装置的形变问题。通过采用微晶石为原料的测量装置结构,极大地减小了装置形变所引起的测量误差。
附图说明
图1是本发明的示意图。
图2是本发明的数据处理流程图。
图3是本发明中滚动滤波的计算流程图。
图4是采样时序图。
图5是本发明一种实施方式的测量示意图。
具体实施方式
实施例:本发明的一种实时测量系统是依照图1固定好两路光电传感器1、2的发射端和接收端,确保两传感器光路垂直;将标准被测物安装于移动导轨上,置于两传感器测量区的重合部分内,即可对于传感器1开始进行非线性修正、动态偏离、校零和滚动滤波四个步骤的数据处理参数的计算和获取;
本发明的一种针对实时测量领域的数据处理方法,以下由一种测量标准被测物外径的光电装置为例进行说明,由传感器直接获得的测量数据为标准被测物外径的两侧边缘值,如图5所示。一种实时测量数据的处理方法,包括校准和修正,校准是为了获得数据处理所需要的上述四步骤系数,然后使用该系数对测量数据进行修正处理。本实施例中的两传感器互相垂直放置,
一、校准的具体步骤如下:
1、标准被测物放置在移动导轨上,将移动导轨移动方向调整到与传感器1的径向方向平行,之后,标准被测物由导轨带动,在两传感器重合区域内运动,每次运动固定的步长1mm,对标准被测物进行测量,采集每次移动后静止时传感器1测量的标准被测物两侧边缘值x1、x2
2、将上述获得的标准被测物两侧边缘值相减,得到每次移动后静止时的标准被测物外径x,如图2所示;
3、获取修正系数:
标准被测物的标称外径y是已知的,根据移动步长指定标准被测物左侧边缘值p1,并指定右侧边缘值p2为左侧边缘值p1与y的和,将这两组数据和修正前的两组位置值数据x1、x2,代入公式(1),并根据公式(1),列出超定方程组,求解得到非线性修正系数a1i、a2i
p 1 = Σ i = 0 n a 1 i x 1 i p 2 = Σ i = 0 n a 2 i x 2 i - - - ( 1 )
二、修正的具体步骤如下:
1、非线性修正:使用上述步骤求得的非线性修正系数,将测得的两侧边缘值x1、x2代入式(2),获得修正后的位置值p1’、p2’,将修正后的位置值相减得到标准被测物直径的测量数据x’;
p 1 ' = Σ i = 0 n a 1 i x 1 i p 2 ' = Σ i = 0 n a 2 i x 2 i - - - ( 2 )
修正前后的测量数据如表1所示:
表1修正前后的测量数据
Figure GSB00000048600300043
  9322.941   11550.35   17.998   22.999
  9805.023   12038.57   18.999   24.002
  10286.49   12527.3   19.998   25.003
  10769.49   13016.67   21.001   26.003
  11251.62   13506.87   22.000   27.001
  11733.53   13996.24   22.998   27.995
  12218.29   14490.22   24.000   28.996
  12702.92   14987   25.001   30.000
  13187.64   15483.2   26.001   31.001
  13671.22   15981.16   26.999   32.003
  14154.91   16476.71   28.000   32.998
2、动态偏离:
(1)将导轨移动方向调整到与传感器1的法向方向平行;
(2)之后,标准被测物由导轨带动,在两传感器重合区域内运动,以每次运动固定步长1mm为例,采集每次移动后静止时传感器2测量的标准被测物两侧边缘值;
(3)将传感器2直接读取的标准被测物两侧边缘值相减,获得测量数据xv
(4)当标准被测物移动到传感器1法向测量区中间位置时,将传感器1直接读取的标准被测物两侧边缘值相加除以2,获得标准被测物在法向测量区中间位置的测量值xm
(5)经过上述非线性修正后获得传感器1输出的两侧边缘值之差为x’;
(6)采集传感器1在法向测量区中所有移动位置的两侧边缘值,将其分别相减,获得在这些位置上的标准被测物外径,再与标准被测物真实值相减之后求平均,获得测量误差平均值Ee
根据公式(3),计算得到经动态偏离修正后的测量值x”:
x″=(xv-xm)·Ee+x′        (3)
3、校零:将标准被测物放置于测量区中心位置,采集传感器1的测量数据,由传感器1经过动态偏离修正后的测量数据x”和标准被测物的已知实际外形尺寸y,根据公式(4),计算得到校零系数c,在测量时,便将校零系数作为固定的系统误差,加入测量数据;
c=y-x″        (4)
4、滚动滤波:将经过以上修正步骤的非线性修正、动态偏离和校零三步处理的数据,再根据图3的流程进行滚动滤波处理:
(1)首先通过如采样时序图4中所示,采集在时间上连续的前n个时刻(t1至tn时刻)的测量数据,
(2)再将上述所获得的测量数据进行平均,获得平均值,
(3)之后,由过往数据的平均值的(1±f)%设定滤波门限,其中滤波门限中的f为手动输入,
(4)随后,对当前时刻采集的数据进行判断筛选,
A、如果在门限之内,
a、将当前时刻数据送入平均数据队列,
b、再将前第n时刻数据送入后续其他数据处理或直接显示并存储,
B、如不在门限之内,重新采集数据进行判断。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (1)

1.一种实时测量数据的处理方法,其特征在于:它包括以下步骤:校准和修正,
(1)校准的具体步骤如下:
A、使用至少两路有重合区域相互垂直放置的传感器对同一标准被测物进行测量;
B、获取修正系数:
使用导轨将标准被测物置于上述重合区域测量区,并在被修正传感器径向上能够使得标准被测物移动固定的步长,采集每次移动后静止时被修正传感器测得的标准被测物两侧边缘值x1、x2,相减获得每次移动后静止时的标准被测物外径x,由于标准被测物的实际外形尺寸y是已知的,根据移动步长指定标准被测物左侧边缘值p1,并指定右侧边缘值p2为左侧边缘值p1与y的和,将这两组数据和修正前的两组位置值数据x1、x2,代入公式(1),并根据公式(1),列出超定方程组,求解得到非线性修正系数a1i、a2i
p 1 = Σ i = 0 n a 1 i x 1 i p 2 = Σ i = 0 n a 2 i x 2 i - - - ( 1 )
(2)、修正的具体步骤如下:
A、非线性修正:使用上述步骤求得的非线性修正系数,将测得的标准被测物两侧边缘值X1、X2代入式(2),获得修正后的位置值p1’、p2’,将修正后的位置值相减得到标准被测物直径的测量数据x’;
p 1 ′ = Σ i = 0 n a 1 i x 1 i p 2 ′ = Σ i = 0 n a 2 i x 2 i - - - ( 2 )
B、动态偏离:
(a)将导轨移动方向调整到与被修正传感器的法向方向平行或0-90度内的其他固定角度的不同位置上的测量数据进行修正;
(b)之后,标准被测物由导轨带动,在各传感器重合区域内运动,每次运动步长固定,采集每次移动后静止时其它路传感器测量的标准被测物两侧边缘值;
(c)将其它路传感器直接读取的标准被测物两侧边缘值相减,获得测量数据xv
(d)当标准被测物移动到被修正传感器法向测量区中间位置时,将被修正传感器直接读取的标准被测物两侧边缘值相加除以2,获得标准被测物在法向测量区中间位置的测量值xm
(e)经过上述非线性修正后获得被修正传感器输出的两侧边缘值之差为x’;
(f)采集被修正传感器在法向测量区中所有移动位置的两侧边缘值,将其分别相减,获得在这些位置上的标准被测物外径,再与标准被测物真实值相减之后求平均,获得测量误差平均值Ee
根据公式(3),计算得到经动态偏离修正后的测量值x”:
x″=(xv-xm)·Ee+x′        (3)
C、校零:将标准被测物放置于测量区中心位置,采集被修正传感器的测量数据,由被修正传感器经过动态偏离修正后的测量数据x”和标准被测物的已知实际外形尺寸y,根据公式(4),计算得到校零系数c,在测量时,便将校零系数作为固定的系统误差,加入测量数据;
c=y-x″      (4)
D、滚动滤波:将经过以上修正步骤的非线性修正、动态偏离和校零三步处理的数据,再进行滚动滤波处理:
(a)首先通过采集在时间上连续的前n个时刻t1至tn时刻的测量数据,
(b)再将上述所获得的测量数据进行平均,获得平均值,
(c)之后,由过往数据的平均值的(1±f)%设定滤波门限,其中滤波门限中的f为手动输入,
(d)随后,对当前时刻采集的数据进行判断筛选,
1)、如果在门限之内,
I、将当前时刻数据送入平均数据队列,
II、再将前第n时刻数据送入后续其他数据处理或直接显示并存储,
2)、如不在门限之内,重新采集数据进行判断。
CN2009100686693A 2009-04-29 2009-04-29 实时测量数据的处理方法 Active CN101539402B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100686693A CN101539402B (zh) 2009-04-29 2009-04-29 实时测量数据的处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100686693A CN101539402B (zh) 2009-04-29 2009-04-29 实时测量数据的处理方法

Publications (2)

Publication Number Publication Date
CN101539402A CN101539402A (zh) 2009-09-23
CN101539402B true CN101539402B (zh) 2011-08-17

Family

ID=41122714

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100686693A Active CN101539402B (zh) 2009-04-29 2009-04-29 实时测量数据的处理方法

Country Status (1)

Country Link
CN (1) CN101539402B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813494B (zh) * 2010-02-25 2012-03-07 上海乐春重工机电设备有限公司 消除传感器自身和外界干扰的方法
CN102493240B (zh) * 2011-11-15 2014-07-16 江苏赛福天钢索股份有限公司 电梯钢丝绳直径在线控制系统及方法
CN103134615B (zh) * 2011-11-29 2014-10-01 上海开通数控有限公司 空压机系统传感器的自动校准方法
TWI451336B (zh) * 2011-12-20 2014-09-01 Univ Nat Cheng Kung 預測模型之建模樣本的篩選方法及其電腦程式產品
CN105630158B (zh) * 2015-12-16 2018-06-29 广东欧珀移动通信有限公司 传感器数据处理方法、装置及终端设备
CN106770861B (zh) * 2016-11-21 2018-01-16 浙江大学 油浸式变压器在线监控数据可用性的评价方法
CN109100125B (zh) * 2018-06-28 2020-06-23 中国直升机设计研究所 一种旋转放大器输出信号修正方法
CN109696665B (zh) * 2018-12-28 2021-02-19 百度在线网络技术(北京)有限公司 超声波传感器测量数据的处理方法、装置及设备
CN110986751B (zh) * 2019-12-24 2020-10-23 深圳市交通工程试验检测中心有限公司 一种北斗、gnss形变监测方法
CN114113153A (zh) * 2020-08-26 2022-03-01 郑高山 对射式测量仪器在线零点校准装置及校准方法
CN113970297A (zh) * 2021-11-19 2022-01-25 江苏智库智能科技有限公司 堆垛机货位列数据的标定测量方法
CN114910565B (zh) * 2022-07-19 2022-09-27 天津市特种设备监督检验技术研究院(天津市特种设备事故应急调查处理中心) 一种非线性超声检测中相对非线性系数的修正方法

Also Published As

Publication number Publication date
CN101539402A (zh) 2009-09-23

Similar Documents

Publication Publication Date Title
CN101539402B (zh) 实时测量数据的处理方法
Shokravi et al. Vehicle-assisted techniques for health monitoring of bridges
Fu et al. A non-contact measuring system for in-situ surface characterization based on laser confocal microscopy
Almeida et al. On the acoustic filtering of the pipe and sensor in a buried plastic water pipe and its effect on leak detection: an experimental investigation
Gui et al. A component decomposition model for 3D laser scanning pavement data based on high-pass filtering and sparse analysis
Mullapudi et al. Shaping streamflow using a real-time stormwater control network
Nasrin et al. Impact of short duration intense rainfall events on sanitary sewer network performance
Selbig et al. The effect of particle size distribution on the design of urban stormwater control measures
Muñoz et al. Effect of the resolution of tipping-bucket rain gauge and calculation method on rainfall intensities in an Andean mountain gradient
Bezak et al. Estimation of suspended sediment loads using copula functions
Rys Investigation of weigh-in-motion measurement accuracy on the basis of steering axle load spectra
Kırbaş IRI sensitivity to the influence of surface distress on flexible pavements
Dong et al. Reconstruction of 3D pavement texture on handling dropouts and spikes using multiple data processing methods
Fereidoon et al. Predicting rainfall and runoff through satellite soil moisture data and SWAT modelling for a poorly gauged basin in Iran
Jeon et al. Traffic measurement on multiple drive lanes with wireless ultrasonic sensors
Pedret Rodés et al. GPR spectra for monitoring asphalt pavements
De Blasiis et al. Assessing of the road pavement roughness by means of LiDAR technology
Huang et al. Ground-based radar interferometry for monitoring the dynamic performance of a multitrack steel truss high-speed railway bridge
Kaloop et al. Adjustment and assessment of the measurements of low and high sampling frequencies of GPS real-time monitoring of structural movement
Chang et al. A case study for the application of an operational two-dimensional real-time flooding forecasting system and smart water level gauges on roads in Tainan City, Taiwan
Bahrani et al. Alternate method of pavement assessment using geophones and accelerometers for measuring the pavement response
Matsuoka et al. Experimental validation of Non-Marker simple image displacement measurements for railway bridges
Chen et al. Comparative analysis of four baseflow separation methods in the south Atlantic-Gulf Region of the US
Laucka et al. Research of the equipment self-calibration methods for different shape fertilizers particles distribution by size using image processing measurement method
Hou et al. Variational mode decomposition based time-varying force identification of stay cables

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Assignee: Tianjin EastLoong Opto-electronic Measuring & Control Technology Co., Ltd.

Assignor: Chen Bingsheng

Contract record no.: 2011120000171

Denomination of invention: Method for processing real-time measuring data

Granted publication date: 20110817

License type: Exclusive License

Open date: 20090923

Record date: 20110830

ASS Succession or assignment of patent right

Owner name: TIANJIN EAST LOONG OPTO-ELECTRONIC MEASURING AND C

Free format text: FORMER OWNER: CHEN BINGSHENG

Effective date: 20130114

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 300190 NANKAI, TIANJIN TO: 300384 NANKAI, TIANJIN

TR01 Transfer of patent right

Effective date of registration: 20130114

Address after: The 300384 Tianjin City Huayuan Industrial Zone Development Road six No. 6 Haitai green industry base K1-4-204

Patentee after: Tianjin EastLoong Opto-electronic Measuring & Control Technology Co., Ltd.

Address before: 300190, Tianjin, Nankai District, Tian Jiang, Nanchang, South Ningxia, 8-3-31

Patentee before: Chen Bingsheng