CN113899457B - 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 - Google Patents
一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 Download PDFInfo
- Publication number
- CN113899457B CN113899457B CN202110988247.9A CN202110988247A CN113899457B CN 113899457 B CN113899457 B CN 113899457B CN 202110988247 A CN202110988247 A CN 202110988247A CN 113899457 B CN113899457 B CN 113899457B
- Authority
- CN
- China
- Prior art keywords
- flame
- image
- temperature
- emissivity
- calibration
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 27
- 238000001228 spectrum Methods 0.000 title claims abstract description 13
- 230000005855 radiation Effects 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000005457 Black-body radiation Effects 0.000 claims description 3
- 239000004215 Carbon black (E152) Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 3
- 229930195733 hydrocarbon Natural products 0.000 claims description 3
- 150000002430 hydrocarbons Chemical class 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims 1
- 238000009529 body temperature measurement Methods 0.000 abstract description 2
- 238000002485 combustion reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
- 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/10—Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,包括以下步骤:(1)设定火焰温度图像上有M×N个像素点,任意一点表示为(i,j),i的取值范围为从1到M的整数,j的取值范围为从1到N的整数,将针对火焰中心方向对应的像素点设为(i0,j0)。本发明中,利用光谱仪获得火焰中心点的辐射强度,再利用彩色CCD相机获得二维温度彩色图像,建立R、G、B值与辐射强度之间的关系,求得火焰中心点的标定系数,利用标定系数标定二维图像上的每一个像素点,采用迭代的方法求得每一个像素点的温度与发射率,完成二维火焰彩色图像标定,本方法设定辐射率是波长的一阶函数,能够实时检测出每个波长下的辐射率,提高测温结果的准确性。
Description
技术领域
本发明涉及热工测试技术领域,尤其涉及一种采用光谱检测的火焰彩色图像标定及温度图像检测方法。
背景技术
火焰温度和发射率是反映燃烧进程的两个重要参数,准确测量这两个参数对深入燃烧机理研宄具有指导作用,同时也对监控燃烧过程具有重要意义。现有技术中,通过双色法温度检测方法进行测定,这种方法操作简单,是目前利用火焰图像进行温度场重建的常用测量方法,这种方法通过忽略火焰黑度随波长变化,假设火焰满足灰性来测量火焰温度,根据辐射定律,利用两个波长下的单色辐射强度可以计算得到被测对象的温度,但是利用这种方法检测温度的前提是假设待测火焰满足灰性条件,事实上火焰的辐射率是随波长变化的,因此,得到的检测结果并不准确,目前,许多学者使用黑体炉对温度场进行标定,也是近几十年来在温度检测领域对温度进行标定的常用方法,但是黑体炉体积大,控温时间长,升温元件容易受到损伤,大多数情况下都是在实验室中进行应用,很难实现现场的、实时的标定。而工程中的探测器不可避免的会积灰,导致信号偏差,为了保证长时间的测量精度就需要周期性的标定,工程应用给传统的黑体炉标定带来了困难。所以,使用黑体炉对彩色火焰图像传感器进行标定的方法具有一定的局限性,发展一种采用光谱检测的火焰彩色图像标定及温度图像检测方法具有现实意义。
发明内容
本发明的目的在于:为了解决上述问题,而提出的一种采用光谱检测的火焰彩色图像标定及温度图像检测方法。
为了实现上述目的,本发明采用了如下技术方案:
一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,包括以下步骤:
(1)设定火焰温度图像上有M×N个像素点,任意一点表示为(i,j),i的取值范围为从1到M的整数,j的取值范围为从1到N的整数,将针对火焰中心方向对应的像素点设为(i0,j0);
(2)利用光谱仪对准火焰中心方向,获得火焰光谱辐射强度Iλ(i0,j0);
(3)利用彩色CCD相机在与光谱仪相同方向上获得的红(R)、绿(G)、蓝(B)三基色值分别为R(i0,j0)、G(i0,j0)、B(i0,j0);设CCD传感器R、G、B三基色代表性波长分别为λR、λG、λB,则火焰中心方向图像像素R(i0,j0)、G(i0,j0)、B(i0,j0)对应的在其代表性波长λR、λG、λB下的火焰辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)根据光谱仪检测结果直接获得;
(4)引入修正系数,cR为对R值的修正系数,cG为对G值的修正系数,cB为对B值的修正系数,建立图像中心像素点三基色信息R(i0,j0)、G(i0,j0)、B(i0,j0)与相应的辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)之间的线性关系:
IR(i0,j0)=cR·R(i0,j0) (1)
IG(i0,j0)=cG·G(i0,j0) (2)
IB(i0,j0)=cB·B(i0,j0) (3)
获得修正系数cR、cG、cB:
cR=IR(i0,j0)/R(i0,j0) (4)
cG=IG(i0,j0)/G(i0,j0) (5)
cB=IB(i0,j0)/B(i0,j0) (6);
(5)根据维恩定律,将发射率表示为波长的一阶函数辐,射强度的计算公式为:
式中,I为辐射强度,λ为波长,T为温度,ε为发射率,Ib为黑体辐射强度,a0、a1表示未知的多项式系数,c1和c2为普朗克常数;在工程实际中,碳氢火焰温度校正只需将发射率函数表达为一阶函数即可满足测量、标定精度要求;
(6)由式(10),得到;
(7)三个方程中含有三个未知数,所以将式(11)和(12)联立求得:
(8)将式(14)、(15)带入式(16)中
采用迭代算法求出温度T(i,j)和发射率ε(i,j)=a0(i,j)+a1(i,j)λ,完成对二维彩色温度图像的标定以及温度和辐射率图像的检测,同时得到二维温度图像T(i,j)和辐射率图像ε(i,j)。
优选地,所述步骤(4)中用修正系数cR、cG、cB将火焰图像全部像素的三基色信息转化为三个代表性波长下单色辐射强度分布图像IR(i,j)、IG(i,j)、IB(i,j);
IR(i,j)=cR·R(i,j) (7)
IG(i,j)=cG·G(i,j) (8)
IB(i,j)=cB·B(i,j) (9)。
综上所述,由于采用了上述技术方案,本发明的有益效果是:
本申请中,首先利用光谱仪对准火焰中心点,获得火焰中心点的辐射强度,再利用彩色CCD相机获得二维温度彩色图像,建立R、G、B值与辐射强度之间的关系,求得火焰中心点的标定系数,由于火焰表面并不严格满足灰体条件,发射率随火焰温度变化而变化,利用标定系数标定二维图像上的每一个像素点,采用迭代的方法求得每一个像素点的温度与发射率,从而完成二维火焰彩色图像标定,本方法设定辐射率是波长的一阶函数,能够实时检测出每个波长下的辐射率,提高测温结果的准确性。
附图说明
图1示出了根据本发明实施例提供的光谱仪标定结构示意图;
图2示出了根据本发明实施例提供的方法总体流程图;
图3示出了根据本发明实施例提供的火焰原始图像示意图;
图4示出了根据本发明实施例提供的CCD的R、G、B图像示意图;
图5示出了根据本发明实施例提供的CCD的IR、IG、IB图像示意图;
图6示出了根据本发明实施例提供的温度分布图;
图7示出了根据本发明实施例提供的发射率分布图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1-7,本发明提供一种技术方案:
1.一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,包括以下步骤:
(1)设定火焰温度图像上有M×N个像素点,任意一点表示为(i,j),i的取值范围为从1到M的整数,j的取值范围为从1到N的整数,将针对火焰中心方向对应的像素点设为(i0,j0);
(2)利用光谱仪对准火焰中心方向,获得火焰光谱辐射强度Iλ(i0,j0);
(3)利用彩色CCD相机在与光谱仪相同方向上获得的红(R)、绿(G)、蓝(B)三基色值分别为R(i0,j0)、G(i0,j0)、B(i0,j0);设CCD传感器R、G、B三基色代表性波长分别为λR、λG、λB,则火焰中心方向图像像素R(i0,j0)、G(i0,j0)、B(i0,j0)对应的在其代表性波长λR、λG、λB下的火焰辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)根据光谱仪检测结果直接获得;
(4)引入修正系数,cR为对R值的修正系数,cG为对G值的修正系数,cB为对B值的修正系数,建立图像中心像素点三基色信息R(i0,j0)、G(i0,j0)、B(i0,j0)与相应的辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)之间的线性关系:
IR(i0,j0)=cR·R(i0,j0) (1)
IG(i0,j0)=cG·G(i0,j0) (2)
IB(i0,j0)=cB·B(i0,j0) (3)
获得修正系数cR、cG、cB:
cR=IR(i0,j0)/R(i0,j0) (4)
cG=IG(i0,j0)/G(i0,j0) (5)
cB=IB(i0,j0)/B(i0,j0) (6);
(5)根据维恩定律,将发射率表示为波长的一阶函数辐,射强度的计算公式为:
式中,I为辐射强度,λ为波长,T为温度,ε为发射率,Ib为黑体辐射强度,a0、a1表示未知的多项式系数,c1和c2为普朗克常数;在工程实际中,碳氢火焰温度校正只需将发射率函数表达为一阶函数即可满足测量、标定精度要求;
(6)由式(10),得到;
(7)三个方程中含有三个未知数,所以将式(11)和(12)联立求得:
(8)将式(14)、(15)带入式(16)中
采用迭代算法求出温度T(i,j)和发射率ε(i,j)=a0(i,j)+a1(i,j)λ,完成对二维彩色温度图像的标定以及温度和辐射率图像的检测,同时得到二维温度图像T(i,j)和辐射率图像ε(i,j)。
具体的,如图3至图5所示,步骤(4)中用修正系数cR、cG、cB将火焰图像全部像素的三基色信息转化为三个代表性波长下单色辐射强度分布图像IR(i,j)、IG(i,j)、IB(i,j);
IR(i,j)=cR·R(i,j) (7)
IG(i,j)=cG·G(i,j) (8)
IB(i,j)=cB·B(i,j) (9)。
采用光谱仪测量某电厂90%负荷下的火焰辐射能量,如表1;采用CCD相机获取火焰原始图像,如图3;火焰R、G、B图像如图4;根据中心点计算标定系数为:CR=IR(i0,j0)/R(i0,j0)=6.8936e+05,CG=IG(i0,j0)/G(i0,j0)=1.3350e+05,CB=IB(i0,j0)/B(i0,j0)=2.3761e+04;将中心点标定系数应用于所有像素,得到IR、IG、IB图像如图5,其中,三个代表性波长:λR=0.6108um,λG=0.5108um
λB=0.4358um;温度分布如图6;发射率分布如图7。
表1实施例数据
综上所述,本实施例所提供的一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,首先利用光谱仪对准火焰中心点,获得火焰中心点的辐射强度,再利用彩色CCD相机获得二维温度彩色图像,建立R、G、B值与辐射强度之间的关系,求得火焰中心点的标定系数,由于火焰表面并不严格满足灰体条件,发射率随火焰温度变化而变化,利用标定系数标定二维图像上的每一个像素点,采用迭代的方法求得每一个像素点的温度与发射率,从而完成二维火焰彩色图像标定。
实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (2)
1.一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,其特征在于,包括以下步骤:
(1)设定火焰温度图像上有M×N个像素点,任意一点表示为(i,j),i的取值范围为从1到M的整数,j的取值范围为从1到N的整数,将针对火焰中心方向对应的像素点设为(i0,j0);
(2)利用光谱仪对准火焰中心方向,获得火焰光谱辐射强度Iλ(i0,j0);
(3)利用彩色CCD相机在与光谱仪相同方向上获得的红(R)、绿(G)、蓝(B)三基色值分别为R(i0,j0)、G(i0,j0)、B(i0,j0);设CCD传感器R、G、B三基色代表性波长分别为λR、λG、λB,则火焰中心方向图像像素R(i0,j0)、G(i0,j0)、B(i0,j0)对应的在其代表性波长λR、λG、λB下的火焰辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)根据光谱仪检测结果直接获得;
(4)引入修正系数,cR为对R值的修正系数,cG为对G值的修正系数,cB为对B值的修正系数,建立图像中心像素点三基色信息R(i0,j0)、G(i0,j0)、B(i0,j0)与相应的辐射强度IR(i0,j0)、IG(i0,j0)、IB(i0,j0)之间的线性关系:
IR(i0,j0)=cR·R(i0,j0) (1)
IG(i0,j0)=cG·G(i0,j0) (2)
IB(i0,j0)=cB·B(i0,j0) (3)
获得修正系数cR、cG、cB:
cR=IR(i0,j0)/R(i0,j0) (4)
cG=IG(i0,j0)/G(i0,j0) (5)
cB=IB(i0,j0)/B(i0,j0) (6);
(5)根据维恩定律,将发射率表示为波长的一阶函数辐,射强度的计算公式为:
式中,I为辐射强度,λ为波长,T为温度,ε为发射率,Ib为黑体辐射强度,a0、a1表示未知的多项式系数,c1和c2为普朗克常数;在工程实际中,碳氢火焰温度校正只需将发射率函数表达为一阶函数即可满足测量、标定精度要求;
(6)由式(10),得到;
(7)三个方程中含有三个未知数,所以将式(11)和(12)联立求得:
(8)将式(14)、(15)带入式(16)中
采用迭代算法求出温度T(i,j)和发射率ε(i,j)=a0(i,j)+a1(i,j)λ,完成对二维彩色温度图像的标定以及温度和辐射率图像的检测,同时得到二维温度图像T(i,j)和辐射率图像ε(i,j)。
2.根据权利要求1所述的一种采用光谱检测的火焰彩色图像标定及温度图像检测方法,其特征在于,所述步骤(4)中用修正系数cR、cG、cB将火焰图像全部像素的三基色信息转化为三个代表性波长下单色辐射强度分布图像IR(i,j)、IG(i,j)、IB(i,j);
IR(i,j)=cR·R(i,j) (7)
IG(i,j)=cG·G(i,j) (8)
IB(i,j)=cB·B(i,j) (9)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110988247.9A CN113899457B (zh) | 2021-08-26 | 2021-08-26 | 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110988247.9A CN113899457B (zh) | 2021-08-26 | 2021-08-26 | 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113899457A CN113899457A (zh) | 2022-01-07 |
CN113899457B true CN113899457B (zh) | 2024-02-02 |
Family
ID=79188018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110988247.9A Active CN113899457B (zh) | 2021-08-26 | 2021-08-26 | 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113899457B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403639A (zh) * | 2008-11-17 | 2009-04-08 | 华中科技大学 | 碳氢火焰的温度图像及黑度图像检测方法 |
CN103557965A (zh) * | 2013-11-25 | 2014-02-05 | 北京汉能清源科技有限公司 | 水泥回转窑温度测定及其温度场在线检测方法、装置 |
JP2017058370A (ja) * | 2015-09-16 | 2017-03-23 | 三菱マテリアル株式会社 | 被測定物の温度、粉塵の温度及び粉塵の濃度を計測する方法 |
CN110160657A (zh) * | 2019-05-31 | 2019-08-23 | 华中科技大学 | 一种基于可见光热成像的高温分布检测方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104849260B (zh) * | 2015-05-18 | 2017-08-01 | 华中科技大学 | 一种锅炉燃烧火焰中的气相碱金属浓度的在线检测方法 |
-
2021
- 2021-08-26 CN CN202110988247.9A patent/CN113899457B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101403639A (zh) * | 2008-11-17 | 2009-04-08 | 华中科技大学 | 碳氢火焰的温度图像及黑度图像检测方法 |
CN103557965A (zh) * | 2013-11-25 | 2014-02-05 | 北京汉能清源科技有限公司 | 水泥回转窑温度测定及其温度场在线检测方法、装置 |
JP2017058370A (ja) * | 2015-09-16 | 2017-03-23 | 三菱マテリアル株式会社 | 被測定物の温度、粉塵の温度及び粉塵の濃度を計測する方法 |
CN110160657A (zh) * | 2019-05-31 | 2019-08-23 | 华中科技大学 | 一种基于可见光热成像的高温分布检测方法及装置 |
Non-Patent Citations (1)
Title |
---|
基于图像处理的火焰温度及辐射率图像检测法;姜志伟, 周怀春, 娄春, 於正前;华中科技大学学报(自然科学版)(第09期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113899457A (zh) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101943604B (zh) | 测温成像系统及其测量方法 | |
CN103913237B (zh) | 三波段红外辐射精确测温方法 | |
Fu et al. | Measurement performance of an optical CCD-based pyrometer system | |
CN110207830A (zh) | 一种非黑体辐射源的图像传感器标定装置和标定方法 | |
CN111044153B (zh) | 一种图谱关联系统红外光谱的非线性定标方法及装置 | |
CN110017904A (zh) | 一种基于ccd相机的多光谱辐射测温方法 | |
CN106644102A (zh) | 一种基于彩色ccd相机的碳氢火焰温度测量方法 | |
CN105354859B (zh) | 一种火焰可见光辐射标定方法 | |
CN104062010B (zh) | 一种优化定标算法的分光光源颜色照度测量仪器 | |
CN112102271B (zh) | 一种基于普通数字相机的实时在线火焰温度测量方法 | |
CN113375815A (zh) | Ccd和红外热像仪相结合的物体表面温度测量方法及系统 | |
CN115265803A (zh) | 一种远距离点目标红外光谱辐射特性测量方法 | |
Yue et al. | Color crosstalk correction for synchronous measurement of full-field temperature and deformation | |
CN114509165B (zh) | 一种光谱发射率测量装置及表面温度测量方法 | |
CN113899457B (zh) | 一种采用光谱检测的火焰彩色图像标定及温度图像检测方法 | |
CN114449079A (zh) | 一种基于手机摄像头的高温测量装置和方法 | |
CN110160660A (zh) | 基于光场相机的高温部件多光谱测温方法及系统 | |
CN212133888U (zh) | 基于辐射光谱的彩色ccd自标定测温装置 | |
CN105606220B (zh) | 一种优化波长矫正方法及采用该方法的分光测色仪 | |
CN112525951A (zh) | 一种加热成像装置及关联辐射图像与积灰温度的方法 | |
CN111649830A (zh) | 基于辐射光谱的彩色ccd自标定测温装置和方法 | |
CN116256069A (zh) | 一种双波长温度测量系统定标方法 | |
CN112013970A (zh) | 一种级联式红外辐射测量系统标定方法 | |
CN106768384B (zh) | 基于辅助调频光源的彩色成像温度场测量装置及方法 | |
CN113405672A (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 |