CN113405672A - 一种基于对数多项式的实时高温温度场测量方法 - Google Patents
一种基于对数多项式的实时高温温度场测量方法 Download PDFInfo
- Publication number
- CN113405672A CN113405672A CN202110689424.3A CN202110689424A CN113405672A CN 113405672 A CN113405672 A CN 113405672A CN 202110689424 A CN202110689424 A CN 202110689424A CN 113405672 A CN113405672 A CN 113405672A
- Authority
- CN
- China
- Prior art keywords
- temperature
- intensity
- logarithm
- polynomial
- primary colors
- 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.)
- Pending
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 239000003086 colorant Substances 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 9
- 206010020843 Hyperthermia Diseases 0.000 claims 1
- 230000036031 hyperthermia Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000005259 measurement Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 abstract description 4
- 230000005855 radiation Effects 0.000 abstract description 4
- 238000001514 detection method Methods 0.000 abstract description 3
- 230000003287 optical effect Effects 0.000 abstract description 3
- 238000013500 data storage Methods 0.000 abstract description 2
- 230000009897 systematic effect Effects 0.000 abstract description 2
- 238000009529 body temperature measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 101100248200 Arabidopsis thaliana RGGB gene Proteins 0.000 description 1
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000003723 Smelting Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000002309 gasification Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000010310 metallurgical process Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011819 refractory material Substances 0.000 description 1
- 239000004071 soot Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0014—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
- G01J5/0018—Flames, plasma or welding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Plasma & Fusion (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
Description
技术领域
本发明涉及高温温度场的测量技术领域,具体涉及一种基于对数多项式的实时高温温度场测量方法。
背景技术
工业领域,有非常多的行业涉及到高温过程,诸如火力发电、煤气化、冶金等。以冶金过程为例,可靠且连续的金属表面温度测量对有效操作控制至关重要。有效的温度控制,可以提高熔炼率,降低燃料消耗,并延长耐火材料寿命。其他工业过程也是类似情况,都只有在可靠且连续的测量出高温物体表面温度场的基础上,才能够针对性对过程进行调整,降本增效。
通常高温物体通常会自发辐射可见光。而可见光的辐射强度与光源的温度直接相关。国内外已开发了诸多基于光学方法测量高温物体表面温度的技术。专利公布号CN112556859 A的专利公开了一种碳烟火焰温度测量方法,该方法需要对同一高温表面透过两个波长相近的滤光片拍摄两次,通过比色测温法计算问题。但由于更换滤光片需要一定时间,必须要求火焰稳定。因此该法只适用于实验室等场景,对环境要求较高,且不能连续获取温度。
专利公布号CN 101403639 A的专利公开了一种碳氢火焰的温度图像及黑度图像检测方法,是一种利用CCD相机的测温方法,通过事先经黑体炉标定CCD探测器,拟合两基色比值的多项式与温度之间的函数关系。在已知该函数关系的基础上测量火焰温度。但由于该方法利用的数据为RGB24位图,基色值以8位字节存储。该数据已经经过图像压缩,对于原始数据而言存在一些非线性损失,并不能直接反映出火焰的辐射强度,因此,将会对温度测量产生较大数值误差。且在其应用黑体炉进行标定时,将两基色比值与温度直接回归,将过多的非线性性交给数据进行拟合,这对数据采集的质量要求非常高,且可能会导致模型过于复杂,泛用性不强。综合考虑该方法数值误差较大,在理论上有所欠缺。
发明内容
本发明基于现有技术中火焰温度测量存在的不足,提供一种基于对数变换和倒数变换的实时高温温度场测量方法,该方法不仅能够在工业场景中实时连续测量,还考虑了图像压缩过程中的损失,降低数值计算问题所带来的误差,测量精度高,检测误差可低于2.2%。
为实现上述目的,本发明采用的技术方案是:
一种基于对数多项式的实时高温温度场测量方法,包括如下步骤:
(1)拍摄原始图像:采用普通相机对待测高温物体进行拍摄,获得原始图像;
(2)获取RGB空间图像:将步骤(1)拍摄获取的原始图像导出颜色滤镜矩阵,进行Bayer插值,获得待测高温物体的RGB空间图像;
(3)计算基色强度对数r:计算待测高温物体的RGB空间图像中任意两基色的比例并取对数,获得任意两基色的强度比值对数,称为基色强度对数r;
步骤(3)中,基色强度对数r的计算公式如下:
r=ln(R/G)
或,r=ln(R/B)
或,r=ln(B/G)
其中,R、G、B分别代表三种基色的强度数据。
优选地,步骤(3)中,不同的待测高温物体颜色选用相同的两种基色计算基色强度对数r。可以提高后续温度计算的稳定性和精度,比如红黄色火焰选用R和G。
其中,n为多项式最大的次数,i表示多项式次数,ki表示i次多项式的系数。
其中,ki系数的获取方法为:采用普通相机拍摄不同温度的温度场,获得一系列已知温度的图像,从而确定基色强度对数r,和r所对应的温度T,根据1/T与r之间的非线性程度设定n值;在设定n之后,根据最小二乘原理,确定所有的ki。
优选地,根据1/T与r之间的非线性程度,设定n,非线性程度越大n越大,通常n可以取2。
进一步地,获取ki系数的具体公式为:
其中,K=[k0 k1...ki...kN],表示i=0,1,…,n的ki构成的向量;
与现有技术相比,本发明具有以下有益效果:
本发明的测量方法为一种非接触式的实时高温温度场测量方法,相较于其他方法,充分考虑光学传感器所得数据与辐射强度的关系,减小由于数据存储问题引入的系统性误差。同时,也考虑了求解计算过程中的数值误差问题。对温度的检测误差可以小于2.2%。
附图说明
图1为本发明的实验装置结构示意图,其中1为待测高温物体,2为拍摄相机,3为计算机。
图2为实施例中相机拍摄获取的原始图像。
图3为实施例中获得的RGB空间图像。
图4为实施例中基色强度对数r与黑体炉温度T的拟合曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。本领域技术人员在理解本发明的技术方案基础上进行修改或等同替换,而未脱离本发明技术方案的精神和范围,均应涵盖在本发明的保护范围内。
本发明的实验装置结构示意图如图1所示,其中1为待测的高温物体,采用拍摄相机2对其进行拍摄,通过计算机3对拍摄相机2的图像进行获取和计算处理。
实施例
在本实施例中待测高温物体为黑体炉,采用数字相机对黑体炉进行拍摄,其中计算机为与相机集成在一起的具有一定计算功能的计算设备。所使用的数字相机为佳能MARKIII G7X,该相机的颜色滤镜矩阵排布方式为RGGB的拜尔滤镜分布。首先,设定数字相机的相机参数,白平衡模式选择为晴天,曝光补偿设为0,设定感光度和快门速度为手动模式,感光度设为1/125,快门时间为1/2000。在使用数字相机拍摄黑体炉后,将相机产生RAW格式数据传输至计算机中,进行计算处理。
设定黑体炉温度从850℃到1200℃以50℃为间隔升温,在每一个温度点都拍摄一次图像,针对每一次拍摄都进行以下步骤:
首先将相机拍摄得到的原始图像,如图2所示,将图像的RAW格式数据进行Bayer插值,获得如图3所示的RGB空间图像,这里的RGB空间为uint16的数据,将该数据改为用double格式数据存储,用以计算出基色强度对数r;本实施例中选择r=ln(R/G),R、G分别为红色和绿色基色的强度数据;
即,
根据该函数关系式,得到黑体炉根据图像计算的温度,表1为黑体炉计算测定的结果与黑体炉设定值比较,从表1中可以看出,通过本发明的方法计算的黑体炉温度误差小于1.8%。
表1黑体炉设定温度、计算温度及误差比例
Claims (9)
2.根据权利要求1所述的基于对数多项式的实时高温温度场测量方法,其特征在于,步骤(3)中,基色强度对数r的计算公式如下:
r=ln(R/G)
或,r=ln(R/B)
或,r=ln(B/G)
其中,R、G、B分别代表三种基色的强度数据。
3.根据权利要求1或2所述的基于对数多项式的实时高温温度场测量方法,其特征在于,步骤(3)中,不同的高温物体颜色选用相同的两种基色计算基色强度对数r。
5.根据权利要求4所述的基于对数多项式的实时高温温度场测量方法,其特征在于,其中ki系数的获取方法为:采用普通相机拍摄不同温度的温度场,获得一系列已知温度的图像,从而确定基色强度对数r,和r所对应的温度T,根据1/T与r之间的非线性程度设定n值;在设定n之后,根据最小二乘原理,确定所有的ki。
6.根据权利要求5所述的基于对数多项式的实时高温温度场测量方法,其特征在于,1/T与r之间的非线性程度越大,n取值越大。
7.根据权利要求5所述的基于对数多项式的实时高温温度场测量方法,其特征在于,n取值2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110689424.3A CN113405672A (zh) | 2021-06-22 | 2021-06-22 | 一种基于对数多项式的实时高温温度场测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110689424.3A CN113405672A (zh) | 2021-06-22 | 2021-06-22 | 一种基于对数多项式的实时高温温度场测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113405672A true CN113405672A (zh) | 2021-09-17 |
Family
ID=77682150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110689424.3A Pending CN113405672A (zh) | 2021-06-22 | 2021-06-22 | 一种基于对数多项式的实时高温温度场测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113405672A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114862795A (zh) * | 2022-04-29 | 2022-08-05 | 吉林大学 | 物体表面温度测量系统、测量方法、智能终端及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4885633A (en) * | 1988-06-13 | 1989-12-05 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Quantitative surface temperature measurement using two-color thermographic phosphors and video equipment |
CN101403639A (zh) * | 2008-11-17 | 2009-04-08 | 华中科技大学 | 碳氢火焰的温度图像及黑度图像检测方法 |
CN104180908A (zh) * | 2014-08-07 | 2014-12-03 | 华能国际电力股份有限公司 | 一种raw图像辐射测温装置及方法 |
CN106644102A (zh) * | 2017-01-17 | 2017-05-10 | 浙江大学 | 一种基于彩色ccd相机的碳氢火焰温度测量方法 |
CN112102271A (zh) * | 2020-09-02 | 2020-12-18 | 浙江大学 | 一种基于普通数字相机的实时在线火焰温度测量方法 |
CN112556859A (zh) * | 2020-12-11 | 2021-03-26 | 中国科学院广州能源研究所 | 一种炭烟火焰温度测量方法 |
-
2021
- 2021-06-22 CN CN202110689424.3A patent/CN113405672A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4885633A (en) * | 1988-06-13 | 1989-12-05 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Quantitative surface temperature measurement using two-color thermographic phosphors and video equipment |
CN101403639A (zh) * | 2008-11-17 | 2009-04-08 | 华中科技大学 | 碳氢火焰的温度图像及黑度图像检测方法 |
CN104180908A (zh) * | 2014-08-07 | 2014-12-03 | 华能国际电力股份有限公司 | 一种raw图像辐射测温装置及方法 |
CN106644102A (zh) * | 2017-01-17 | 2017-05-10 | 浙江大学 | 一种基于彩色ccd相机的碳氢火焰温度测量方法 |
CN112102271A (zh) * | 2020-09-02 | 2020-12-18 | 浙江大学 | 一种基于普通数字相机的实时在线火焰温度测量方法 |
CN112556859A (zh) * | 2020-12-11 | 2021-03-26 | 中国科学院广州能源研究所 | 一种炭烟火焰温度测量方法 |
Non-Patent Citations (2)
Title |
---|
李汉舟等: "基于数字图像处理的温度检测算法研究", 《中国电机工程学报》 * |
格伦.兰德等: "《美国布鲁克斯摄影学院数码摄影教程》", 31 January 2014, 中国摄影出版社 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114862795A (zh) * | 2022-04-29 | 2022-08-05 | 吉林大学 | 物体表面温度测量系统、测量方法、智能终端及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101282489B (zh) | 一种光源检测设备和方法及图像处理方法 | |
CN106644102B (zh) | 一种基于彩色ccd相机的碳氢火焰温度测量方法 | |
CN112102271B (zh) | 一种基于普通数字相机的实时在线火焰温度测量方法 | |
CN106679817A (zh) | 一种用于标定红外热像仪的方法 | |
CN113375815B (zh) | Ccd和红外热像仪相结合的物体表面温度测量方法及系统 | |
CN110160657B (zh) | 一种基于可见光热成像的高温分布检测方法及装置 | |
CN111721225A (zh) | 一种高温环境下温度变形动态测量方法及装置 | |
CN111707382B (zh) | 一种温度变形同步测量的动态光学补偿方法及装置 | |
CN113405672A (zh) | 一种基于对数多项式的实时高温温度场测量方法 | |
CN111337132B (zh) | 一种温度测量方法、装置以及数字图像采集设备 | |
TWI442032B (zh) | 非接觸式溫度量測方法 | |
CN113865717A (zh) | 一种基于高速相机的瞬态高温比色测温装置 | |
CN111207833B (zh) | 一种基于图像数据归一化技术的测温方法 | |
TWI465702B (zh) | 非接觸式溫度量測方法 | |
CN109738068B (zh) | 一种校正多光谱相机响应值非线性的方法 | |
Yu et al. | A method for color calibration based on simulated annealing optimization | |
CN113758683B (zh) | 一种基于平均色彩饱和度的相机系统辐照后光谱退化评估方法 | |
CN115937093A (zh) | 融合hsl空间和改进暗通道技术的烟雾浓度检测方法 | |
CN106595868B (zh) | 一种基于改进三色法的高炉燃烧带温度场检测方法 | |
TW201122309A (en) | Method of calibrating a light source | |
CN105744267B (zh) | 基于齐次多项式数码相机参数可变的获取三刺激值方法 | |
JP2014137780A (ja) | 画像処理方法及び画像処理装置 | |
CN113899457B (zh) | 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 | |
US20140340547A1 (en) | Imaging apparatus | |
Czech et al. | The linearization of the relationship between scene luminance and digital camera output levels |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210917 |