CN106197690B - 一种宽温范围条件下的图像校准方法及系统 - Google Patents

一种宽温范围条件下的图像校准方法及系统 Download PDF

Info

Publication number
CN106197690B
CN106197690B CN201610482300.7A CN201610482300A CN106197690B CN 106197690 B CN106197690 B CN 106197690B CN 201610482300 A CN201610482300 A CN 201610482300A CN 106197690 B CN106197690 B CN 106197690B
Authority
CN
China
Prior art keywords
temperature
correction coefficient
multinomial
equation
correction
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
CN201610482300.7A
Other languages
English (en)
Other versions
CN106197690A (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.)
Hubei Jiuzhiyang Infrared System Co Ltd
Original Assignee
Hubei Jiuzhiyang Infrared System Co Ltd
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 Hubei Jiuzhiyang Infrared System Co Ltd filed Critical Hubei Jiuzhiyang Infrared System Co Ltd
Priority to CN201610482300.7A priority Critical patent/CN106197690B/zh
Publication of CN106197690A publication Critical patent/CN106197690A/zh
Application granted granted Critical
Publication of CN106197690B publication Critical patent/CN106197690B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/52Radiation pyrometry, e.g. infrared or optical thermometry using comparison with reference sources, e.g. disappearing-filament pyrometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/80Calibration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077Imaging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

本发明公开了一种宽温范围条件下的图像校准方法及系统,该方法包括以下步骤:S1、使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;S2、在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;S3、在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而完成对红外图像的非均匀校正。本发明能够满足红外图像在大动态范围的非均匀性校正的要求,消除由非均匀性引起的图像噪声问题。

Description

一种宽温范围条件下的图像校准方法及系统
技术领域
本发明涉及红外图像的非均匀性校正领域,尤其涉及一种宽温范围条件下的图像校准方法及系统。
背景技术
红外成像由于具有比较强的环境适应性,对伪装目标具有较强的识别能力,其设备具有体积小、功耗低的等特点。红外成像技术由于这些特点在军事和民用领域都起到了重要的作用。然而红外探测器的自身特性,使得红外图像中存在非均匀性,并且严重影响其成像质量,所以为了得到清晰的红外图像,研究非均匀性校正技术成为必不可少的内容。
非均匀性校正一般可以分为参考辐射源的定标和基于场景的校正方法。常见的温度定标方法有一点温度校正和两点温度校正,该类方法假设在一定时间的间隔内,探测单元的响应特性保持不变;基于场景的校正方法,如时域高通滤波,神经网络等,该类方法通过场景序列的迭代估计,自适应更新校正参数。
对于宽温范围条件下的图像非均匀性校正,由于常用的一点/两点校正,其在标定点外的非均匀性较差,采用线性模型,温度动态范围较窄,并不适合宽温范围;基于场景的计算量大,容易存在“鬼影”,工程中较少应用。目前对宽温范围条件下的图像校准技术较少有应用,相关研究也比较少,因此有必要提出一种宽温范围条件下的图像校准方法。
发明内容
本发明要解决的技术问题在于针对现有技术中宽温度范围图像校正的非均匀性较差,温度动态范围较窄的缺陷,提供一种能够满足红外图像在大动态范围的非均匀性校正的要求,消除由非均匀性引起的图像噪声问题的宽温范围条件下的图像校准方法及系统。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种宽温范围条件下的图像校准方法,包括以下步骤:
S1、使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
S2、在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
S3、引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
进一步地,本发明的步骤S1中建立多项级次方程的方法具体为:
从低温到饱和的高温,设置三个以上的温度段,在每个温度段内通过黑体辐射源设置n个温度点,红外探测器某一像元的输出为x,校正后的输出为y,它们之间关系表示为多项级数方程:y=K0+K1x+K2x2+....+Knxn
进一步地,本发明的步骤S2中求解一个温度段内的校正系数的方法具体为:
在一个温度段内,从曲线拟合的精度考虑n取5,即将多项级数方程转化为y=K5x5+K4x4+K3x3+A(x+B)(x+C),通过计算得出校正系数K5、K4、K3、A、B、C。
进一步地,本发明的步骤S3中对校正系数进行修正的方法具体为:
经过试验和计算得出多项级数方程的四阶以上系数较小,仅对多项级数方程得到的校正系数A、B、C进行修正;
将得到的红外图像为多个子块,分子块得到的校正系数用下标加数字来表示,对多项级数方程中的二次参数采用几何分解的方法,对加入漂移补偿因子M1、M2、M3、m1、m2、m3的多项级数方程进行分解,最终得到校正后的探测器两温度段之间的某像元的输出为:
y=x3+A′(x-M1B1-M2B2-M3)(x-m1C1-m2C2-m3)
通过最小二乘法求解式中的M1、M2、M3、m1、m2、m3
本发明提供一种宽温范围条件下的图像校准系统,包括:
多项级次方程建立单元,用于使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
校正系数计算单元,用于在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
校正系数修正单元,用于引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
本发明产生的有益效果是:本发明的宽温范围条件下的图像校准方法,采用了多段的多项级数方程求解校正系数,较一元线性方程的两点温度校正,提高了拟合曲线的精度,同时增加了动态范围,满足宽温范围条件下的图像校准要求;增加了对环境温度的校正,进一步提高了红外图像非均匀性的处理效果,更好的适配宽温范围条件下的图像校准技术要求;能够满足红外图像在大动态范围的非均匀性校正的要求,消除由非均匀性引起的图像噪声问题。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的宽温范围条件下的图像校准方法的校正方法流程图;
图2是本发明实施例的宽温范围条件下的图像校准方法的具体实施例的校正方法流程图;
图3是本发明实施例的宽温范围条件下的图像校准方法的红外焦平面辐射图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明实施例的宽温范围条件下的图像校准方法,包括以下步骤:
S1、使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
建立多项级次方程的方法具体为:
从低温到饱和的高温,设置三个以上的温度段,在每个温度段内通过黑体辐射源设置n个温度点,红外探测器某一像元的输出为x,校正后的输出为y,它们之间关系表示为多项级数方程:y=K0+K1x+K2x2+....+Knxn
S2、在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
求解一个温度段内的校正系数的方法具体为:
在一个温度段内,从曲线拟合的精度考虑n取5,即将多项级数方程转化为y=K5x5+K4x4+K3x3+A(x+B)(x+C),通过计算得出校正系数K5、K4、K3、A、B、C。
S3、引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
对校正系数进行修正的方法具体为:
经过试验和计算得出多项级数方程的四阶以上系数较小,仅对多项级数方程得到的校正系数A、B、C进行修正;
将得到的红外图像为多个子块,分子块得到的校正系数用下标加数字来表示,对多项级数方程中的二次参数采用几何分解的方法,对加入漂移补偿因子M1、M2、M3、m1、m2、m3的多项级数方程进行分解,最终得到校正后的探测器两温度段之间的某像元的输出为:
y=x3+A′(x-M1B1-M2B2-M3)(x-m1C1-m2C2-m3)
通过最小二乘法求解式中的M1、M2、M3、m1、m2、m3
如图2所示,在本发明的另一个具体实施例中,该方法包括以下步骤:
首先将参考辐射源(黑体)从低温到高温分三个以上温度段,在每个温度段内取5个温度点,利用公式(1),计算出校正系数K5、K4、K3、A、B、C:
y=K5x5+K4x4+K3x3+A(x+B)(x+C) (1)
如图3所示,可知红外焦平面成像除了与目标源的温度还和环境背景温度有关,也就是认为,环境温度同样也影响着图像的非均匀性:
取两段多项级数方程,设一段环境温度为Ts1时,通过求解计算得到其校正系数为A1、B1、C1,另一段环境温度为Ts2时,同样得到此时其校正系数为A2、B2、C2。当温度由Ts1到Ts2,探测器整体受环境温度的影响发生漂移,非均匀性变差。可得校正模型如下式:
式中M1、M2、M3、m1、m2、m3为漂移补偿因子。
然后,将红外图像按照一定等分,分多个像元集块,(这里建议按照像面大小的十等分为宜),在各个像元集块内,计算像元非均匀性按照公式(1)可表示为:
Y-(K5x5+K4x4+K3x3)=A′(X-M1B1-M2B2-M3)(X-m1C1-m2C2-m3)(3)
在每个像元集块校正后的Y数据中,搜索小于一定阈值的像素对,则相邻像素之间的差异可以表示为:
考虑到简化计算,这里计算漂移补偿因子时设:
Q1=A′(a)X(a)-A′(b)X(b)
Q2=A′(a)B1(a)-A′(b)B1(b)
Q3=A′(a)B2(a)-A′(b)B2(b)
Q4=A′(a)-A′(b) (5)
对公式(3)简化计算,如下:
Y-(K5x5+K4x4+K3x3)=A′(X-M1B1-M2B2-M3)2 (6)
依照y=a2(x-b)2采用最小二乘法简化计算漂移补偿因子,对所有像素求取平方和可得:
在满足等式左边最小的情况下,求取相关系数:
对M1、M2、M3分别求偏导可得:
求解上述方程(8)可得M1、M2、M3,采用公式(5)-(8)同理可得m1、m2、m3,从而获得整个校正系数。根据调整后的校正系数完成非均匀性校正,然后输出最终效果图,完成了宽温度范围条件下的图像校准。
本发明实施例的宽温范围条件下的图像校准系统,包括:
多项级次方程建立单元,用于使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
校正系数计算单元,用于在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
校正系数修正单元,用于引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (5)

1.一种宽温范围条件下的图像校准方法,其特征在于,包括以下步骤:
S1、使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
S2、在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
S3、引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
2.根据权利要求1所述的宽温范围条件下的图像校准方法,其特征在于,步骤S1中建立多项级次方程的方法具体为:
从低温到饱和的高温,设置三个以上的温度段,在每个温度段内通过黑体辐射源设置n个温度点,红外探测器某一像元的输出为x,校正后的输出为y,它们之间关系表示为多项级次方程:y=K0+K1x+K2x2+....+Knxn
3.根据权利要求2所述的宽温范围条件下的图像校准方法,其特征在于,步骤S2中求解一个温度段内的校正系数的方法具体为:
在一个温度段内,从曲线拟合的精度考虑n取5,即将多项级次方程转化为y=K5x5+K4x4+K3x3+A(x+B)(x+C),通过计算得出校正系数K5、K4、K3、A、B、C。
4.根据权利要求3所述的宽温范围条件下的图像校准方法,其特征在于,步骤S3中对校正系数进行修正的方法具体为:
经过试验和计算得出多项级次方程的四阶以上系数较小,仅对多项级次方程得到的校正系数A、B、C进行修正;
将得到的红外图像为多个子块,分子块得到的校正系数用下标加数字来表示,对多项级次方程中的二次参数采用几何分解的方法,对加入漂移补偿因子M1、M2、M3、m1、m2、m3的多项级次方程进行分解,最终得到校正后的探测器两温度段之间的某像元的输出为:
y=x3+A′(x-M1B1-M2B2-M3)(x-m1C1-m2C2-m3)
通过最小二乘法求解式中的M1、M2、M3、m1、m2、m3
5.一种宽温范围条件下的图像校准系统,其特征在于,包括:
多项级次方程建立单元,用于使用不同温度的均匀黑体辐射作为参考源,将其从低温到高温分为多个温度段,在每个温度段内设置多个温度点,根据红外探测器某一像元校正前后的输出值,建立多项级次方程;
校正系数计算单元,用于在一个温度段内,根据曲线拟合的精度确定选取的温度点数量,设置多项级次方程的级数为温度点数量,进而求解出该温度段内的校正系数;并求解各个温度段的校正系数;
校正系数修正单元,用于引入环境温度对校正系数的影响,在多项级次方程中加入漂移补偿因子并对其进行分解,根据最小二乘法求解修正后的校正系数,进而根据修正后的校正系数完成对红外图像的非均匀校正。
CN201610482300.7A 2016-06-27 2016-06-27 一种宽温范围条件下的图像校准方法及系统 Active CN106197690B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610482300.7A CN106197690B (zh) 2016-06-27 2016-06-27 一种宽温范围条件下的图像校准方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610482300.7A CN106197690B (zh) 2016-06-27 2016-06-27 一种宽温范围条件下的图像校准方法及系统

Publications (2)

Publication Number Publication Date
CN106197690A CN106197690A (zh) 2016-12-07
CN106197690B true CN106197690B (zh) 2019-01-22

Family

ID=57461316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610482300.7A Active CN106197690B (zh) 2016-06-27 2016-06-27 一种宽温范围条件下的图像校准方法及系统

Country Status (1)

Country Link
CN (1) CN106197690B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106679817B (zh) * 2016-12-23 2019-08-09 武汉星网光测科技有限公司 一种用于标定红外热像仪的方法
CN109711517B (zh) * 2018-12-13 2022-04-15 湖南岛千岛物联网科技有限责任公司 集成自校准温度传感器的rfid标签
CN111964786B (zh) * 2020-07-20 2021-08-31 四川虹美智能科技有限公司 温度检测方法及装置
CN111964789B (zh) * 2020-08-14 2021-08-10 深圳市视美泰技术股份有限公司 温度测量方法、装置、计算机设备和存储介质
CN112304436B (zh) * 2020-10-23 2022-07-19 北京百度网讯科技有限公司 用于创建温度校正模型的方法、装置、电子设备和介质
CN113514155B (zh) * 2021-04-13 2023-04-28 武汉华中数控股份有限公司 一种无快门的非均匀校正的方法
CN113658054B (zh) * 2021-07-06 2024-03-29 北京空间机电研究所 一种基于温漂特征线逼近的红外图像拼接校正方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564598A (zh) * 2012-01-06 2012-07-11 电子科技大学 一种红外探测器测温的定标和校正方法及相应的测温方法
CN103076097A (zh) * 2013-01-06 2013-05-01 河北汉光重工有限责任公司 基于参照源的分段线性非均匀矫正方法
CN103674259A (zh) * 2012-09-25 2014-03-26 中国航天科工集团第二研究院二〇七所 一种适用于红外热像定量处理的多点非均匀性校正方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR0133637B1 (ko) * 1994-06-21 1998-04-23 정명세 새로운 내삽공식을 이용하여 복사온도계를 교정하는 방법
US8282273B2 (en) * 2009-03-24 2012-10-09 K-Space Associates, Inc. Blackbody fitting for temperature determination
WO2016033297A1 (en) * 2014-08-28 2016-03-03 Seek Thermal, Inc. Thermography for a thermal imaging camera

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102564598A (zh) * 2012-01-06 2012-07-11 电子科技大学 一种红外探测器测温的定标和校正方法及相应的测温方法
CN103674259A (zh) * 2012-09-25 2014-03-26 中国航天科工集团第二研究院二〇七所 一种适用于红外热像定量处理的多点非均匀性校正方法
CN103076097A (zh) * 2013-01-06 2013-05-01 河北汉光重工有限责任公司 基于参照源的分段线性非均匀矫正方法

Also Published As

Publication number Publication date
CN106197690A (zh) 2016-12-07

Similar Documents

Publication Publication Date Title
CN106197690B (zh) 一种宽温范围条件下的图像校准方法及系统
Ancuti et al. I-HAZE: A dehazing benchmark with real hazy and haze-free indoor images
US7899271B1 (en) System and method of moving target based calibration of non-uniformity compensation for optical imagers
Cao et al. Effective strip noise removal for low-textured infrared images based on 1-D guided filtering
CN108846805B (zh) 一种基于场景自适应的红外热图像两点非均匀校正方法
Hardie et al. Scene-based nonuniformity correction with reduced ghosting using a gated LMS algorithm
CN106197673B (zh) 一种自适应宽温度范围非均匀校正方法及系统
CN103218778A (zh) 一种图像和视频的处理方法及装置
CN109813442B (zh) 一种基于多帧处理的内部杂散辐射非均匀性校正方法
Lai et al. Spatiotemporal adaptive nonuniformity correction based on BTV regularization
CN112348754B (zh) 一种低照度彩色图像增强方法及装置
CN105928627A (zh) 一种消除制冷型热像仪画面复现非均匀性的方法
CN111080561A (zh) 一种时域高通滤波方法
CN112710397B (zh) 一种基于温度替代的两点校正方法及系统
CN106023108A (zh) 基于边界约束和上下文正则化的图像去雾算法
CN103491318A (zh) 一种红外焦平面探测器图像校正方法及系统
CN110782403A (zh) 一种红外图像非均匀性校正方法
Li et al. Single-frame-based column fixed-pattern noise correction in an uncooled infrared imaging system based on weighted least squares
Li et al. Interior radiation noise reduction method based on multiframe processing in infrared focal plane arrays imaging system
CN105092043B (zh) 一种基于场景的变积分时间的非均匀性校正方法
Yu et al. An improved retina-like nonuniformity correction for infrared focal-plane array
Bekhtin et al. Algorithmic supply of IR sensors with FPN using texture homogeneity levels
Jara Chavez et al. Acceleration algorithm for constant-statistics method applied to the nonuniformity correction of infrared sequences
CN111047521A (zh) 一种基于图像熵的红外图像非均匀性参数化校正优选方法
CN115002360A (zh) 一种基于稳健估计的红外视频非均匀性校正方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant