CN109297604B - 一种获取红外热像仪两点校正定标参数的方法及系统 - Google Patents
一种获取红外热像仪两点校正定标参数的方法及系统 Download PDFInfo
- Publication number
- CN109297604B CN109297604B CN201811137971.5A CN201811137971A CN109297604B CN 109297604 B CN109297604 B CN 109297604B CN 201811137971 A CN201811137971 A CN 201811137971A CN 109297604 B CN109297604 B CN 109297604B
- Authority
- CN
- China
- Prior art keywords
- average value
- xsor
- zsor
- obtaining
- time average
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000001931 thermography Methods 0.000 claims 5
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000006467 substitution reaction 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
-
- 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/80—Calibration
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明涉及一种获取红外热像仪两点校正定标参数的方法及系统,分别采集T1、T2温度下的黑体数据:使用待处理红外热像仪连续采集M帧黑体数据获取该M帧黑体数据的时间平均值Xsori,j、Zsori,j,根据时间平均值得到列噪声值Xsorj、Zsorj;获取时间平均值得到一帧图像,计算该帧图像所有点的平均值;根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j。本发明通过对采集的数据分离列噪声;使用分离列噪声后的数据计算两点校正参数值,得到一个更符合红外热像仪探测元的特性的校正参数值。
Description
技术领域
本发明涉及一种定标参数的方法及系统,尤其涉及一种获取红外热像仪两点校正定标参数的方法及系统。
背景技术
红外热像仪的主要部件红外焦平面阵列,由于生产工艺的原因,存在固有非均匀性,需要进行非均匀性校正。目前普遍应用的两点定标的校正方法,在定标计算的过程中,都包含了红外焦平面阵列读出电路的列噪声,影响了红外热像仪的非均匀性校正效果。
发明内容
本发明解决的技术问题是:构建一种获取红外热像仪两点校正定标参数的方法及系统,克服现有技术非均匀性校正效果不佳的技术问题。
本发明的技术方案是:构建一种获取红外热像仪两点校正定标参数的方法及系统,包括如下步骤:
分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j,根据时间平均值得到列噪声值Xsorj;将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,根据时间平均值得到列噪声值Zsorj;
获取时间平均值Xsori,j得到一帧图像,计算该帧图像所有点的平均值Xsor;获取时间平均值Zsori,j得到一帧图像,计算该帧图像所有点的平均值Zsor;
根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:所述时间平均值Xsori,j、时间平均值Zsori,j,计算过程如下:
其中:XNi,j为M帧中采集的第N帧图像数据,ZNi,j为P帧中采集的第N帧图像数据,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:计算时间平均值数据每一列的平均值,得到每一列的列噪声值Xsorj、列噪声值Zsorj,计算过程如下:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:计算所采集图像的时间平均值Xsori,j得到的一帧图像上的所有点的平均值Xsor,计算所采集图像的时间平均值Zsori,j得到的一帧图像上的所有点的平均值Zsor:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:两点校正定标参数获取过程如下:
Xi,j=Xsori,j-Xsorj+Xsor
Zi,j=Zsori,j-Zsorj+Zsor
其中:Xi,j为第一定标参数,Zi,j为第二定标参数,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的技术方案是:构建一种获取红外热像仪两点校正定标参数的系统,包括采集模块、时间平均值获取模块、列噪声值获取模块、平均值获取模块、定标参数获取模块,所述采集模块分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数;所述时间平均值获取模块获取该M帧黑体数据的时间平均值Xsori,j,所述时间平均值获取模块获取该P帧黑体数据的时间平均值Zsori,j,所述列噪声值获取模块根据时间平均值Xsori,j得到列噪声值Xsorj,所述列噪声值获取模块根据时间平均值Zsori,j得到列噪声值Zsorj,所述平均值获取模块根据时间平均值Xsori,j得到的一帧图像计算该帧图像所有点的平均值Xsor;所述平均值获取模块根据时间平均值Zsori,j得到的一帧图像计算该帧图像所有点的平均值Zsor;所述定标参数获取模块根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,所述定标参数获取模块根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:所述时间平均值获取模块获取过程如下:
其中:XNi,j为M帧中采集的第N帧图像数据,ZNi,j为P帧中采集的第N帧图像数据,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:计算时间平均值数据每一列的平均值,得到每一列的列噪声值Xsorj、列噪声值Zsorj,所述列噪声值获取模块获取过程如下:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:所述平均值获取模块计算所采集图像的时间平均值Xsori,j得到的一帧图像上的所有点的平均值Xsor,计算所采集图像的时间平均值Zsori,j得到的一帧图像上的所有点的平均值Zsor:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的进一步技术方案是:所述定标参数获取模块获取过程如下:
Xi,j=Xsori,j-Xsorj+Xsor
Zi,j=Zsori,j-Zsorj+Zsor
其中:Xi,j为第一定标参数,Zi,j为第二定标参数,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的技术效果是:构建一种获取红外热像仪两点校正定标参数的方法及系统,分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j,根据时间平均值得到列噪声值Xsorj;将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,根据时间平均值得到列噪声值Zsorj;获取时间平均值Xsori,j得到一帧图像,计算该帧图像所有点的平均值Xsor;获取时间平均值Zsori,j得到一帧图像,计算该帧图像所有点的平均值Zsor;根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。本发明通过对采集的数据分离列噪声;使用分离列噪声后的数据计算两点校正参数值,得到一个更符合红外热像仪探测元的特性的校正参数值。
附图说明
图1为本发明的结构示意图。
具体实施方式
下面结合具体实施例,对本发明技术方案进一步说明。
如图1所示,本发明的具体实施方式是:构建一种获取红外热像仪两点校正定标参数的方法及系统,构建一种获取红外热像仪两点校正定标参数的系统,包括采集模块、时间平均值获取模块、列噪声值获取模块、平均值获取模块、定标参数获取模块,包括如下步骤:
分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j,根据时间平均值得到列噪声值Xsorj;将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,根据时间平均值得到列噪声值Zsorj;
所述采集模块1采集两种温度下的黑体数据:即将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j。时间平均值是一种定义概念,即:在一段时间内采集了多帧图像,该多帧图像的采集时间间隔相同,将该多帧图像对应像素点位置的像素灰度平均值即为该多帧图像的时间平均值。所述时间平均值获取模块2获取过程如下,具体来说:
时间平均值Xsori,j
将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,所述时间平均值获取模块2获取过程如下,具体来说:
其中:XNi,j为M帧中采集的第N帧图像数据,ZNi,j为P帧中采集的第N帧图像数据,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
获取时间平均值Xsori,j得到一帧图像,计算该帧图像所有点的平均值Xsor;获取时间平均值Zsori,j得到一帧图像,计算该帧图像所有点的平均值Zsor;
根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
根据时间平均值得到列噪声值Xsorj,计算时间平均值数据每一列的平均值,得到每一列的列噪声值Xsorj,所述列噪声值获取模块3获取过程如下:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
根据时间平均值得到列噪声值Zsorj,计算时间平均值数据每一列的平均值,得到每一列的列噪声值Zsorj,所述列噪声值获取模块获取过程如下:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
所述平均值获取模块4计算所采集图像的时间平均值Xsori,j得到的一帧图像上的所有点的平均值Xsor:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
所述平均值获取模块4计算所采集图像的时间平均值Zsori,j得到的一帧图像上的所有点的平均值Zsor:
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
定标参数获取模块5获取两点校正定标参数获取过程如下:
Xi,j=Xsori,j-Xsorj+Xsor
Zi,j=Zsori,j-Zsorj+Zsor
其中:Xi,j为第一定标参数,Zi,j为第二定标参数,i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
本发明的技术效果是:构建一种获取红外热像仪两点校正定标参数的方法及系统,分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j,根据时间平均值得到列噪声值Xsorj;将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,根据时间平均值得到列噪声值Zsorj;获取时间平均值Xsori,j得到一帧图像,计算该帧图像所有点的平均值Xsor;获取时间平均值Zsori,j得到一帧图像,计算该帧图像所有点的平均值Zsor;根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。本发明通过对采集的数据分离列噪声;使用分离列噪声后的数据计算两点校正参数值,得到一个更符合红外热像仪探测元的特性的校正参数值。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (8)
1.一种获取红外热像仪两点校正定标参数的方法,其特征在于,包括如下步骤:
分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,获取该M帧黑体数据的时间平均值Xsori,j,根据时间平均值得到列噪声值Xsorj;将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数,获取该P帧黑体数据的时间平均值Zsori,j,根据时间平均值得到列噪声值Zsorj;计算时间平均值数据每一列的平均值,得到每一列的列噪声值Xsorj、列噪声值Zsorj;
获取时间平均值Xsori,j得到一帧图像,计算该帧图像所有点的平均值Xsor;获取时间平均值Zsori,j得到一帧图像,计算该帧图像所有点的平均值Zsor;
根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,
Xi,j=Xsori,j-Xsorj+Xsor
Zi,j=Zsori,j-Zsorj+Zsor
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
5.一种获取红外热像仪两点校正定标参数的系统,其特征在于,包括采集模块、时间平均值获取模块、列噪声值获取模块、平均值获取模块、定标参数获取模块,所述采集模块分别采集T1、T2温度下的黑体数据:将黑体温度设置为T1,使用待处理红外热像仪连续采集M帧黑体数据,其中,M=2N,N为大于或等于6的自然数,将黑体温度设置为T2,使用待处理红外热像仪连续采集P帧黑体数据,其中,P=2N,N为大于或等于6的自然数;所述时间平均值获取模块获取该M帧黑体数据的时间平均值Xsori,j,所述时间平均值获取模块获取该P帧黑体数据的时间平均值Zsori,j,所述列噪声值获取模块根据时间平均值Xsori,j得到列噪声值Xsorj,所述列噪声值获取模块根据时间平均值Zsori,j得到列噪声值Zsorj,计算时间平均值数据每一列的平均值,得到每一列的列噪声值Xsorj、列噪声值Zsorj,所述平均值获取模块根据时间平均值Xsori,j得到的一帧图像计算该帧图像所有点的平均值Xsor;所述平均值获取模块根据时间平均值Zsori,j得到的一帧图像计算该帧图像所有点的平均值Zsor;所述定标参数获取模块根据时间平均值Xsori,j、列噪声值Xsorj、平均值Xsor得到第一校正定标参数Xi,j,所述定标参数获取模块根据时间平均值Zsori,j、列噪声值Zsorj、平均值Zsor得到第二校正定标参数Zi,j,
Xi,j=Xsori,j-Xsorj+Xsor
Zi,j=Zsori,j-Zsorj+Zsor
其中:i为1到n的自然数,j为1到m的自然数,n×m为红外热像仪探测器阵列大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811137971.5A CN109297604B (zh) | 2018-09-28 | 2018-09-28 | 一种获取红外热像仪两点校正定标参数的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811137971.5A CN109297604B (zh) | 2018-09-28 | 2018-09-28 | 一种获取红外热像仪两点校正定标参数的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109297604A CN109297604A (zh) | 2019-02-01 |
CN109297604B true CN109297604B (zh) | 2020-06-09 |
Family
ID=65164543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811137971.5A Active CN109297604B (zh) | 2018-09-28 | 2018-09-28 | 一种获取红外热像仪两点校正定标参数的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109297604B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113639880B (zh) * | 2021-08-13 | 2023-09-12 | 中国科学院上海技术物理研究所 | 一种真空红外定标辐射响应一致性方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076096A (zh) * | 2013-01-07 | 2013-05-01 | 南京理工大学 | 基于中值直方图均衡的红外非均匀性校正算法 |
WO2014173464A1 (en) * | 2013-04-26 | 2014-10-30 | Mtech Imaging Pte Ltd | Noise correction in an infra-red imaging system |
CN104406699A (zh) * | 2014-11-26 | 2015-03-11 | 浙江红相科技股份有限公司 | 基于自适应红外图像校正算法红外热像仪 |
CN107341780A (zh) * | 2017-07-12 | 2017-11-10 | 成都中昊英孚科技有限公司 | 一种红外图像预处理校正方法 |
-
2018
- 2018-09-28 CN CN201811137971.5A patent/CN109297604B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103076096A (zh) * | 2013-01-07 | 2013-05-01 | 南京理工大学 | 基于中值直方图均衡的红外非均匀性校正算法 |
WO2014173464A1 (en) * | 2013-04-26 | 2014-10-30 | Mtech Imaging Pte Ltd | Noise correction in an infra-red imaging system |
CN104406699A (zh) * | 2014-11-26 | 2015-03-11 | 浙江红相科技股份有限公司 | 基于自适应红外图像校正算法红外热像仪 |
CN107341780A (zh) * | 2017-07-12 | 2017-11-10 | 成都中昊英孚科技有限公司 | 一种红外图像预处理校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109297604A (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9947086B2 (en) | Image adjustment based on locally flat scenes | |
EP3289759B1 (en) | Compact row column noise filter for an imaging system | |
US10605668B2 (en) | Thermography process for converting signal to temperature in a thermal imaging system | |
US10128808B2 (en) | Gain calibration for an imaging system | |
US10742913B2 (en) | Shutterless calibration | |
US20160065848A1 (en) | Thermography for a thermal imaging camera | |
CN103528694B (zh) | 一种用红外热像仪测量目标物体的温度的方法 | |
JP2020522937A5 (zh) | ||
CN110595630B (zh) | 基于探测器温度的多点非均匀校正的方法 | |
CN104677501B (zh) | 非制冷红外焦平面阵列非均匀性校正的方法和装置 | |
US10467736B2 (en) | Image adjustment based on locally flat scenes | |
US20170243326A1 (en) | Pixel decimation for an imaging system | |
CN110631706B (zh) | 红外图像的校正方法、装置及存储介质 | |
US20190102869A1 (en) | Apparatus for monitoring surroundings of vehicle and method of calibrating the same | |
WO2017177392A1 (zh) | 一种红外测温方法及装置 | |
US10890490B2 (en) | Thermography process for converting signal to temperature in a thermal imaging system | |
CN107084797A (zh) | 国产红外探测器模拟输出自适应调整装置及其调整方法 | |
CN109297604B (zh) | 一种获取红外热像仪两点校正定标参数的方法及系统 | |
CN113138026A (zh) | 一种基于测温系统的测温方法、装置、计算机设备 | |
CN106500855A (zh) | 一种红外探测器盲元检测方法 | |
US20180315171A1 (en) | Image adjustment based on locally flat scenes | |
CN107566763B (zh) | 大尺寸图像传感器及其图像校正方法 | |
CN105704406A (zh) | 一种图像处理中坏点去除的方法 | |
CN104166960B (zh) | 一种基于场景的自适应非均匀固定噪声去除方法 | |
KR102015087B1 (ko) | 테라헤르츠 이미지 센서용 화소 어레이 구조, 테라헤르츠 이미지 센서, 그 제조방법 및 이를 포함하는 영상장치 |
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 |