CN101776486A - 一种基于红外焦平面非均匀性指纹模式的校正方法 - Google Patents
一种基于红外焦平面非均匀性指纹模式的校正方法 Download PDFInfo
- Publication number
- CN101776486A CN101776486A CN200910273496A CN200910273496A CN101776486A CN 101776486 A CN101776486 A CN 101776486A CN 200910273496 A CN200910273496 A CN 200910273496A CN 200910273496 A CN200910273496 A CN 200910273496A CN 101776486 A CN101776486 A CN 101776486A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mover
- focal plane
- msubsup
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012937 correction Methods 0.000 claims abstract description 66
- 230000001419 dependent effect Effects 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 52
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000010606 normalization Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 10
- 230000007613 environmental effect Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 238000003702 image correction Methods 0.000 claims description 3
- 238000003331 infrared imaging Methods 0.000 abstract description 9
- 238000001514 detection method Methods 0.000 abstract description 8
- 238000007796 conventional method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 8
- 230000003044 adaptive effect Effects 0.000 description 5
- 230000036760 body temperature Effects 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004297 night vision Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Landscapes
- Radiation Pyrometers (AREA)
- Image Processing (AREA)
Abstract
一种基于红外焦平面非均匀性指纹模式的校正方法,属于红外成像技术领域。目的是在无背景帧的情况下,对非制冷红外焦平面探测器获取的红外图像进行非均匀性校正。本发明包括原始数据采集步骤、非均匀性指纹提取步骤、校正处理步骤。所谓非均匀性指纹,是指每个红外焦平面探测器都有相对稳定的非均匀性模式及其随温度变化的规律,这两者统称为非均匀性指纹。利用这些可以事先估算出来的非均匀性模式和规律,便可以在红外焦平面探测器获取实际红外图像后,对其进行非均匀性校正。与常规的红外焦平面非均匀性校正方法比较本发明能够有效地减小了探测装置的体积,也不需要像常规方法每次校正时都需要利用均匀挡板获取背景帧,从而大大简化了校正过程。
Description
技术领域
本发明属于红外成像技术领域,具体涉及一种基于红外焦平面非均匀性指纹模式的校正方法。
背景技术
随着红外探测器的发展,红外成像系统也相继发展。在第一代红外成像系统中,采用线列探测器,通过一维光机扫描成像。到了20世纪70年代中期,IRFPA(Infrared FocalPlane Array,红外焦平面阵列)探测器的出现标志着第二代红外成像系统——凝视红外成像系统的诞生。与线列探测器相比,焦平面探测器成像具有空间分辨率高、探测能力强、帧频高等优点,正迅速成为红外成像技术的主流器件。目前凝视红外成像系统已开始广泛应用于夜视、海上营救搜索、天文、工业热探测和医学等民用领域,是红外成像系统的发展方向。然而,红外焦平面阵列存在的非均匀性与无效像元严重影响了系统的成像质量,降低了系统的空间分辨率、温度分辨率、探测距离以及辐射量的正确度量,直接制约着系统的最终性能。尽管随着器件制作工艺的改进,焦平面的非均匀性和无效像元问题有了较大改善,但离完全解决问题还有很大距离,仍是当前红外焦平面阵列成像系统必须解决的首要问题。
针对红外焦平面阵列的非均匀性问题所提出的相应的校正方法,主要分为两大类:一类是基于定标的校正方法,如一点法,两点法等。该类方法原理简洁,硬件易于实现和集成;校正精度高,可用于场景温度的度量;对目标没有任何要求,是实际IRFPA组件产品中主要采用的方法。但这类方法受限于IRFPA响应漂移带来的校正误差;实际校正时需要参考源进行标定,使得设备装置相对复杂;同时在实际应用中需要进行周期性的定标,定标频率取决于系统的稳定性,对于实际探测器不易做到快速反应。另一类是基于场景类的自适应校正方法,如时域高通滤波校正法,神经网络校正法和恒定统计约束校正法等。这类方法可以在一定程度上克服IRFPA响应漂移带来的校正误差,不要求或只需要简单的定标,根据场景信息适应性的更新校正系数,但这类算法应用时计算量大,往往需要特殊并行计算机结构来实现,不利于系统硬件的实现、集成以及对场景的实时处理。
华中科技大学易新建等人在《红外与激光工程》2004年第33卷第1期《红外焦平面阵列非均匀性的两点校正及依据》一文中以普朗克(Plank)辐射定律和红外探测元的线性响应模型为基础,在理论上完整地推导了红外焦平面非均匀性的两点校正方法。文章从理论上论证了两点法的物理依据,表明如果IRFPA的响应是稳定的、线性的,则两点校正的算法没有误差,但实际上IRFPA探测元响应都是非线性的,而且存在响应漂移的问题,因此,用两点校正方法存在较大的剩余误差。
华中科技大学图像识别与人工智能研究所张天序等人在《红外与毫米波学报》2005年第4期《红外焦平面非均匀性噪声的空间频率特性及空间自适应非均匀性校正方法改进》一文中分析了红外焦平面阵列非均匀性噪声的空间频率特性,指出空间低频噪声为其中的主要成分。针对传统空域自适应校正方法去除低频空间噪声存在的不足,提出采用一点校正与神经网络自适应校正相结合的方法。该方法在空间低频噪声占优时能获得较好的校正效果,但存在目标退化问题。
华中科技大学图像识别与人工智能研究所石岩等人在《红外与毫米波学报》2005年第5期《红外焦平面阵列非均匀性自适应校正算法中目标退化与伪像的消除方法》一文中,从场景中边缘信息获取的角度出发,分析到了上述问题出现的原因,并提出了采用边缘指导的神经网络自适应智能型校正算法(ED_NN_NUC)来消除目标退化(fade-out,表明目标图像变模糊,即是目标的信噪比降低、与背景的反差变小、融于背景,不易被识别)。该算法在自适应非均匀校正过程中,自适应地提取当前帧校正后图像边缘信息,以此指导校正参数的更新环节。该算法在能准确获取场景边缘信息的前提下,能较好地抑制目标边缘的退化和伪像,特别有良好的保留弱小亮目标的作用,只是在运算速度方面较慢。
中山大学信息科学与技术学院汪民等人在《红外技术》2007年第6期《一种非制冷焦平面阵列图像漂移的双温度补偿新方法》一文以温度漂移为研究对象,分析环境温度与机芯温度的规律和相互之间的关系,提出了一种双变量线性回归模型进行图像温度漂移补偿的新方法,较好地解决非制冷IRFPA的图像输出因漂移造成图像测温偏差大的工业应用难题。
南京理工大学电子工程与光电技术学院白俊奇等人在《基于环境温度补偿的红外焦平面探测器非均匀性校正模型》一文中通过分析环境变化对红外焦平面阵列探测器输出的影响,建立了一种基于环境温度和目标温度非线性的非均匀性理论模型,实验表明能够提高非均匀性校正的精度。但利用该模型对红外焦平面探测器进行非均匀性校正需要记录大量的数据,对红外焦平面探测器存储硬件要求较高。
分析以往的两大类非均匀性校正方法,基于定标的校正方法,如两点法,由于IRFPA存在响应漂移,实际校正时需要利用均匀参考源作为挡板周期性采集面源黑体图像作为背景帧用于定标,否则校正后会存在较大的剩余误差;基于场景类的自适应校正方法,如神经网络校正法,应用时计算量大,不利于系统的实时处理,而且校正后往往还会出现目标退化和伪像等问题。本发明从研究IRFPA相对稳定的非均匀性模式及其随温度变化的规律(本发明称之为IRFPA的非均匀性指纹)出发,挖掘出相关规律信息,并精炼成代表非均匀性指纹的少量数据,存储在红外焦平面探测器的存储单元中。利用本发明对实际红外图像进行校正时,红外焦平面探测器不需要置于恒温箱中,也无需利用背景帧用于定标,只需要利用存储在红外焦平面探测器存储单元中的少量IRFPA非均匀性指纹数据,经过环境温度的即时测量和简单的运算即可有效地校正IRFPA的非均匀性。
发明内容
本发明提出一种基于红外焦平面非均匀性指纹模式的校正方法,目的是无需利用背景帧用于定标,运算简单,利于硬件实现和实时处理。
一种基于红外焦平面非均匀性指纹模式的图像校正方法,包括(本发明中温度单位均为摄氏度,所用IRFPA的规格为R×C):
(1)原始数据采集步骤:将红外焦平面置于恒温箱中,在恒温箱内的环境温度从设定的下限温度起,每间隔一恒定温度增量采集一组帧数固定的图像序列,直至持续到设定的上限温度,从而得到多组不同环境温度的图像序列;
(2)非均匀性指纹提取步骤:以上述多组图像序列中红外焦平面工作稳定以后的图像为基础数据,利用归一化和最小二乘法计算出红外焦平面非均匀性指纹;
(3)校正处理步骤:首先,以环境温度作为输入参数,根据上述提取出来的非均匀性指纹计算出非均匀性特征图imgNU,再利用如下公式对红外焦平面实际所成图像进行校正:
AftimgX(i,j)=imgX(i,j)-imgNU(i,j)+mean(imgNU),
其中,AftimgX为校正后的图像,(i,j)表示图像矩阵上第i行第j列的像素位置,mean()表示对矩阵求均值。
进一步地,所述的原始数据采集步骤包括:把红外焦平面探测器置于恒温箱内,恒温箱的温度能在一定的温度范围内进行调节,把红外焦平面探测器工作时周围的环境温度记为TS,红外焦平面探测器在恒温箱设定的环境温度下工作,对面源黑体进行成像,在面源黑体温度TB恒定的情况下,恒温箱从设定的环境温度下限TSL变化到环境温度上限TSH,红外焦平面探测器每间隔环境恒定的温度增量ΔTS采集一组面源黑体图像序列作为实验数据,红外焦平面探测器在恒温箱设定的每个环境温度下的采集时间相同,均为t秒;采集到的面源黑体图像帧数也一样,用FrmNumber表示,FrmNumber为自然数,红外焦平面探测器共采集到M组图像序列,每组图像序列为FrmNumber帧,M为自然数,计算公式为:
M=(TSH-TSL)/ΔTS+1,
取采集面源黑体图像序列时的M个环境温度组成TSA,TSA为M×1的矩阵,TSA的各个元素取值为:
TSA(n)=TSL+(n-1)ΔTS,n=1,2,…M,
其中,M组面源黑体图像序列每组都具有如下规律:随着采集时间的增大,面源黑体图像每个像素的灰度值都呈上升趋势,到达一定的时间ts后,面源黑体图像每个像素的灰度值都趋于稳定,ts是一个和具体红外焦平面阵列(IRFPA)有关的常量,记M组面源黑体图像序列每组趋于稳定的黑体图像帧数均为N。
进一步地,,所述的非均匀性指纹包括非均匀性指纹图f1、f2,指纹参数a,b,c,d和p1,p2,p3,非均匀性指纹图f1、f2用于记录红外焦平面探测器相对稳定的非均匀性模式,非均匀性参数a,b,c,d和p1,p2,p3分别为两条曲线方程的系数,这两条曲线记录红外焦平面探测器非均匀性随环境温度变化的规律,非均匀性指纹的提取步骤包括:
(1)非均匀性指纹图f1的计算:把M组面源黑体图像序列每组中趋于稳定后的N帧图像分别取平均:
(n=1,2…M)表示(n=1,2…M)归一化后的结果,是和规格一样的矩阵,./表示点除,该符号两边操作数必须为两个规格一致的矩阵,表示该两矩阵对应位置的元素相除,结果为一个矩阵,规格和操作数矩阵一致,
从矩阵(n=1,2…M)中选取一个参考位置(ist,jst),(ist,jst)的选取原则如下:为矩阵上最小的元素,当n从1到M取值时,组成规格为M×1的矩阵,记为Dst,Dst各个元素的取值如下式所示:
选择以下的数学模型,对M组数据(TSA(n),Dst(n))(n=1,2…M)按最小二乘原理进行曲线拟合,即可求得非均匀性指纹参数a,b,c,d:
其中,T表示环境温度,为自变量,f为因变量,a,b,c,d为待求的非均匀性指纹参数,*表示数乘,如果该符号两边的操作数都为标量,则表示两个标量相乘,结果仍为标量;如果该符号两边的操作数为标量和矩阵,表示该标量与矩阵每个元素相乘,结果是一个矩阵;
(3)非均匀性指纹参数p1,p2,p3的计算:记(n=1,2…M)中位置分别为(ira,jra)和(ist,jst)两个元素的差值为(ist,jst)为步骤(2)中确定的参考位置,(ira,jra)为矩阵中任意一个位置,的计算公式为:
y=p1*T2+p2*T+p3,
其中T表示环境温度,为自变量,y为因变量,p1,p2,p3为待求的非均匀性指纹参数;
进一步地,所述的校正处理步骤中,
校正处理步骤是等待红外焦平面探测器输出的红外图像趋于稳定后,对红外焦平面探测器当前输出的实际红外图像imgX进行校正。此步骤中,红外焦平面探测器不再置于恒温箱中,而是置于实际的工作环境中,红外焦平面探测器的环境温度TS随着周围工作环境的温度而变化。红外焦平面探测器装有温度敏感元,用于测量即时环境温度。记TSW为红外焦平面探测器工作时,作为输入参数用于计算非均匀性特征图imgNU的温度,红外焦平面探测器温度敏感元输出的即时环境温度为TSC。非均匀性特征图记录了红外焦平面探测器当前输出的红外图像的非均匀性,具体包括:
(1)记红外焦平面探测器趋于稳定后获取第一帧实际红外图像时温度敏感元输出的环境温度为TSG,此时TSW=TSG;
(2)以TSW为输入参数,根据步骤(2)中存储的非均匀性指纹图和非均匀性指纹参数,估算出该环境温度下红外焦平面探测器相对稳定不变的非均匀性特征图imgNU;
(3)通过如下公式对红外焦平面探测器当前输出的实际红外图像imgX进行校正:AftimgX(i,j)=imgX(i,j)-imgNU(i,j)+mean(imgNU),0<i<R,0<j<C,
其中,AftimgX为校正后的图像,R、C表示红外图像的规格和IRFPA的规格一致,为R行,C列,(i,j)表示红外图像第i行第j列的像元,mean(imgNU)表示对imgNU求均值:
(4)红外焦平面探测器获取下一帧实际图像。如果下式成立:
|TSC-TSW|>ΔTSW,
则
TSW=TSC,
并转步骤(2),否则TSW不变,并转步骤(3),其中,ΔTSW为设定的判断TSW是否需要更新的温度阈值。
所述的一种基于红外焦平面非均匀性指纹模式的校正方法,其特征在于所述校正处理步骤中,非均匀性特征图imgNU通过以下公式计算得出:
imgNU的计算公式为:
imgNU=f1·*Temp,
其中,T表示作为所述输入参数的环境温度,a,b,c,d,p1,p2,p3为存储的非均匀性指纹数据,·*表示点乘,该符号两边操作数必须为两个规格一致的矩阵,表示该两矩阵对应位置的元素相乘,结果为一个矩阵,规格和操作数矩阵一致。
本发明利用曲线拟合的方法来研究IRFPA响应随环境温度的变化规律,把相关规律用少量数据进行表示,并把这些数据存储在红外焦平面探测器存储单元中。校正时,以红外焦平面探测器的温度敏感元的输出为输入参数,利用存储的数据经过简单的运算即可估算出当前环境温度下红外焦平面探测器的非均匀性,用红外焦平面探测器采集到的实际红外图像和估算出的红外焦平面探测器的非均匀性进行简单运算,即可对实际的红外图像中的非均匀性进行有效的校正。与常规的方法相比,本发明具有以下优点:无需将红外焦平面探测器置于恒温箱中,有效减小探测装置的体积和复杂度;无需利用均匀挡板获取背景帧用于定标,运算简单,利于硬件实现和实时处理,而且有效地克服了IRFPA响应漂移带来的校正误差。
附图说明
图1为本发明的流程图;
图2中的曲线为环境温度TS=20度,面源黑体温度TB=40度时采集的面源黑体图像序列中位置为(129,99)的像元的灰度值随时间的变化。
图7中为环境温度TS=20度,面源黑体温度TB=40度时采集的面源黑体图像序列第36帧的原始图和校正图,图7(a)为原始图,图7(b)为校正图;
图8中为环境温度TS=-10度,面源黑体温度TB=40度时采集的面源黑体图像序列第36帧的原始图和校正图,图8(a)为原始图,图8(b)为校正图;
图9中为环境温度TS=20度,面源黑体温度TB=40度时采集的面源黑体图像序列的第36帧原始图和校正图的直方图,图9(a)为原始图的直方图,图9(b)为校正图的直方图。
图10中为环境温度TS=-10度,面源黑体温度TB=40度时采集的面源黑体图像序列的第36帧原始图和校正图的直方图,图10(a)为原始图的直方图,图10(b)为校正图的直方图。
具体实施方式
下面结合附图和具体实施方式对本发明作详细说明,其中选用了实际规格为R×C=288×384的IRFPA(本发明中温度单位均为摄氏度)。
(1)原始数据采集步骤。采集方法:把红外焦平面探测器置于恒温箱内,恒温箱的温度能在一定的温度范围内进行调节。面源黑体温度恒定于TB=40度,恒温箱环境温度下限TSL=-40度到环境温度上限TSH=20度,红外焦平面探测器以每间隔ΔTS=10度的环境温度增量采集面源黑体图像作为实验数据。红外焦平面探测器在每个环境温度下的采集时间t相同,均为4800秒;采集到的面源黑体图像帧数FrmNumber也一样,均为36帧。红外焦平面探测器共采集到M=7组图像序列。
取采集面源黑体图像序列时的7个环境温度组成TSA。TSA为7×1的矩阵,TSA的各个元素取值为:
TSA(n)=TSL+(n-1)ΔTS,n=1,2,…7,
其中,7组面源黑体图像序列每组都具有如下规律:随着采集时间的增大,面源黑体图像每个像素的灰度值都呈上升趋势,到达一定的时间ts后,面源黑体图像每个像素的灰度值都趋向稳定。7组面源黑体图像序列每组趋于稳定的黑体图像帧数为6。
图2所示曲线为环境温度TS=20度时,面源黑体图像序列中位置为(129,99)的像元的灰度值随时间的变化。
(2)非均匀性指纹提取步骤。
(2.1)f1的计算:把7组面源黑体图像序列每组中趋于稳定后的6帧图像分别取平均:
(n=1,2…7)表示(n=1,2…7)归一化后的结果,是和规格一样的矩阵。./表示点除,该符号两边操作数必须为两个规格一致的矩阵,表示该两矩阵对应位置的元素相除,结果为一个矩阵,规格和操作数矩阵一致。
从矩阵(n=1,2…7)中选取一个参考位置(ist,jst)=(20,20)。(ist,jst)的选取原则如下:为矩阵上最小的元素。取(n=1,2…7)矩阵上的同一位置(20,20)的元素组成规格为7×1的矩阵,记为Dst,Dst各个元素的取值如下式所示:
选择以下的数学模型,对7组数据(TSA(n),Dst(n))(n=1,2…7)按最小二乘原理进行曲线拟合:
即可求得非均匀性指纹参数a=0.5191,b=0.5839,c=-0.02007,d=0.4609。其中,T表示环境温度,为自变量。f为因变量。*表示数乘,如果该符号两边的操作数都为标量,则表示两个标量相乘,结果仍为标量;如果该符号两边的操作数为标量和矩阵,表示该标量与矩阵每个元素相乘,结果是一个矩阵。
(2.3)非均匀性指纹参数p1,p2,p3的计算:记(n=1,2…7)中位置分别为(ira,jra)和(ist,jst)两个元素的差值为(ist,jst)为(2.2)中确定的参考位置(20,20),(ira,jra)为矩阵中任意一个位置。的计算公式为:
图4中所示两条曲线分别表示(n=1,2…7)矩阵中,位置分别为(129,99)和(20,20)两个元素的值随环境温度TS变化的趋势。图5中曲线表示矩阵中位置为(129,99)和(20,20)的两个元素的差值随环境温度TS变化的趋势。
y=p1*T2+p2*T+p3,
即可求得非均匀性指纹参数p1=-0.000167,p2=-0.02007,p3=0.4609。其中T表示环境温度,为自变量。y为因变量。
得到第二个非均匀性指纹图f2。
(3)校正处理步骤。校正处理步骤是等待红外焦平面探测器输出的红外图像趋于稳定后,对红外焦平面探测器当前输出的实际红外图像imgX进行校正。此步骤中,红外焦平面探测器不再置于恒温箱中,而是置于实际的工作环境中,红外焦平面探测器的环境温度TS随着周围工作环境的温度而变化。红外焦平面探测器装有温度敏感元,用于测量即时环境温度。记TSW为红外焦平面探测器工作时,作为输入参数用于计算非均匀性特征图imgNU的温度,红外焦平面探测器温度敏感元输出的即时环境温度为TSC。非均匀性特征图记录了红外焦平面探测器当前输出的红外图像的非均匀性。
(3.1)记红外焦平面探测器趋于稳定后获取第一帧实际红外图像时温度敏感元输出的环境温度为TSG,此时TSW=TSG。
(3.2)以TSW为输入参数,根据步骤(2)中存储的非均匀性指纹图和非均匀性指纹参数,估算出该环境温度下红外焦平面探测器相对稳定不变的非均匀性特征图imgNU。
imgNU的计算公式为:
imgNU=f1.*Temp,
其中,T表示作为所述输入参数的环境温度。a,b,c,d,p1,p2,p3为存储的非均匀性指纹数据。·*表示点乘,该符号两边操作数必须为两个规格一致的矩阵,表示该两矩阵对应位置的元素相乘,结果为一个矩阵,规格和操作数矩阵一致。
(3.3)通过如下公式对红外焦平面探测器当前输出的实际红外图像imgX进行校正:
AftimgX(i,j)=imgX(i,j)-imgNU(i,j)+mean(imgNU),
,
0<i<288,0<j<384
其中,AftimgX为校正后的图像.(i,j)表示红外图像第i行第j列的像元,0<i<288,0<j<384,mean(imgNU)表示对imgNU求均值:
(4)红外焦平面探测器获取下一帧实际图像。如果下式成立:
|TSC-TSW|>ΔTSW,
则
TSW=TSC,
并转步骤(2),否则TSW不变,并转步骤(3)。其中,ΔTSW为设定的判断TSW是否需要更新的温度阈值,ΔTSW=5度。
校正前后的图像见图7和图8。由图9和图10可以看出,校正前的面源黑体图像的灰度值分布在较大的范围内,即非均匀性比较严重;校正后的面源黑体图像的灰度值分布比较集中,非均匀性得到明显改善。为定量表示本发明的校正效果,应用以下公式计算校正前后面源黑体图像的非均匀性:
其中m为IRFPA探测单元数量,d为焦平面上无效像元的个数,Vs为IRFPA的响应输出电压,means(Vs)为整个IRFPA的平均响应输出电压。
计算结果见表1,由结果可以看出,面源黑体图像的非均匀性得到较大改善。
表1校正前后图像的非均匀性UR
环境温度 | -40度 | -30度 | -20度 | -10度 | 0度 | 10度 | 20度 |
校正前 | 0.1685 | 0.1519 | 0.1365 | 0.1223 | 0.0979 | 0.0846 | 0.0761 |
校正后 | 0.0013 | 0.0011 | 0.0015 | 0.0023 | 0.0032 | 0.0031 | 0.0003 |
Claims (5)
1.一种基于红外焦平面非均匀性指纹模式的图像校正方法,包括:
一、原始数据采集步骤:将红外焦平面置于恒温箱中,在恒温箱内的环境温度从设定的下限温度起,每增加一恒定温度增量采集一组帧数固定的图像序列,直至设定的上限温度,从而得到多组不同环境温度的图像序列;
二、非均匀性指纹提取步骤:以上述多组图像序列中红外焦平面工作稳定以后的图像为基础数据,利用归一化和最小二乘算法计算出红外焦平面非均匀性指纹;
三、校正处理步骤:首先,以环境温度作为输入参数,根据上述提取出来的非均匀性指纹计算出非均匀性特征图imgNU,再利用校正公式对红外焦平面实际所成图像进行校正:
AftimgX(i,j)=imgX(i,j)-imgNU(i,j)+mean(imgNU),
其中,AftimgX为校正后的图像,(i,j)表示图像矩阵上第i行第j列的像素位置,mean()表示对矩阵求均值。
2.如权利要求1所述的一种基于红外焦平面非均匀性指纹模式的图像校正方法,其特征在于,所述的原始数据采集步骤包括:
把红外焦平面探测器置于恒温箱内,恒温箱的温度能在一定的温度范围内进行调节,把红外焦平面探测器工作时周围的环境温度记为TS,红外焦平面探测器在恒温箱设定的环境温度下工作,对面源黑体进行成像,在面源黑体温度TB恒定的情况下,恒温箱从设定的环境温度下限TSL变化到环境温度上限TSH,红外焦平面探测器每间隔恒定的环境温度增量ΔTS采集一组面源黑体图像序列作为实验数据,红外焦平面探测器在恒温箱设定的每个环境温度下的采集时间相同,均为t秒,采集到的面源黑体图像帧数也一样,均为FrmNumber帧,FrmNumber为自然数,红外焦平面探测器共采集到M组图像序列,M为自然数,计算公式为:
M=(TSH-TSL)/ΔTS+1,
取采集面源黑体图像序列时的M个环境温度组成TSA,TSA为M×1的矩阵,TSA的各个元素取值为:
TSA(n)=TSL+(n-1)ΔTS,n=1,2,…M,
3.如权利要求1或2所述的一种基于红外焦平面非均匀性指纹模式的图像校正方法,其特征在于,所述的非均匀性指纹包括非均匀性指纹图f1、f2,非均匀性指纹参数a,b,c,d和p1,p2,p3,非均匀性指纹的提取步骤包括:
(1)非均匀性指纹图f1的计算:
记M组面源黑体图像序列每组趋于稳定的黑体图像帧数均为N,把M组面源黑体图像序列每组的所述N帧图像分别取平均:
(2)非均匀性指纹参数a,b,c,d的计算:
从矩阵中选取一个参考位置(ist,jst),(ist,jst)的选取原则如下:为矩阵上最小的元素,取矩阵上的所有同一位置(ist,jst)的元素组成规格为M×1的矩阵,记为Dst,Dst各个元素的取值如下式所示:
选择以下的数学模型,对M组数据(TSA(n),Dst(n))按最小二乘原理进行曲线拟合,即可求得非均匀性指纹参数a,b,c,d:
其中,T表示环境温度,为自变量,f为因变量,a,b,c,d为待求的非均匀性指纹参数,*表示数乘,如果该符号两边的操作数都为标量,则表示两个标量相乘,结果仍为标量;如果该符号两边的操作数为标量和矩阵,表示该标量与矩阵每个元素相乘,结果是一个矩阵;
(3)非均匀性指纹参数p1,p2,p3的计算:
y=p1*T2+p2*T+p3,
其中T表示环境温度,为自变量,y为因变量,p1,p2,p3为待求的非均匀性指纹参数;
(4)非均匀性指纹图f2的计算:
4.如权利要求1或2或3所述的一种基于红外焦平面非均匀性指纹模式的图像校正方法,其特征在于,所述的校正处理步骤中,
红外焦平面探测器置于实际的工作环境中,探测器上装有温度敏感元,用于测量即时环境温度,记TSW为红外焦平面探测器工作时,作为输入参数用于计算非均匀性特征图imgNU的温度,红外焦平面探测器温度敏感元输出的即时环境温度为TSC,具体步骤为:
(A)记红外焦平面探测器趋于稳定后获取第一帧实际红外图像时温度敏感元输出的环境温度为TSG,此时TSW=TSG;
(B)以TSW为输入参数,根据所述的非均匀性指纹和该输入参数,计算出该即时环境温度下红外焦平面探测器的非均匀性特征图imgNU;
(C)通过所述的校正公式进行校正;
(D)红外焦平面探测器获取下一帧实际图像,如果下式成立:
|TSC-TSW|>ΔTSW,
则
TSW=TSC,
并转步骤(B),否则TSW不变,并转步骤(C),其中,ΔTSW为设定的判断所述输入参数TSW是否需要更新的温度阈值。
5.如权利要求1-4任一项所述的一种基于红外焦平面非均匀性指纹模式的校正方法,其特征在于所述校正处理步骤中,所述的非均匀性特征图imgNU通过以下公式计算得出:
imgNU=f1.*Temp,
其中,T表示作为所述输入参数的环境温度,a,b,c,d,p1,p2,p3为存储的非均匀性指纹数据,.*表示点乘,该符号两边操作数必须为两个规格一致的矩阵,表示该两矩阵对应位置的元素相乘,结果为一个矩阵,规格和操作数矩阵一致。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102734969A CN101776486B (zh) | 2009-12-31 | 2009-12-31 | 一种基于红外焦平面探测器非均匀性指纹模式的校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102734969A CN101776486B (zh) | 2009-12-31 | 2009-12-31 | 一种基于红外焦平面探测器非均匀性指纹模式的校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101776486A true CN101776486A (zh) | 2010-07-14 |
CN101776486B CN101776486B (zh) | 2011-08-10 |
Family
ID=42513004
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102734969A Expired - Fee Related CN101776486B (zh) | 2009-12-31 | 2009-12-31 | 一种基于红外焦平面探测器非均匀性指纹模式的校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101776486B (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012122788A1 (zh) * | 2011-03-14 | 2012-09-20 | 浙江兆晟科技有限公司 | 图像校正系数的获取方法、非均匀图像校正方法及系统 |
CN103219997A (zh) * | 2013-02-06 | 2013-07-24 | 北京空间机电研究所 | 多谱段多焦面拼接红外探测器控制与信号采样电路 |
CN103279947A (zh) * | 2013-05-06 | 2013-09-04 | 西安交通大学 | 一种基于单幅图像灰度特征的红外图像校正方法 |
CN103528690A (zh) * | 2013-09-24 | 2014-01-22 | 电子科技大学 | 一种红外热像仪的非均匀校正方法 |
CN105466566A (zh) * | 2015-12-05 | 2016-04-06 | 中国航空工业集团公司洛阳电光设备研究所 | 一种红外非均匀性校正实时补偿方法 |
CN105987758A (zh) * | 2015-02-05 | 2016-10-05 | 南京理工大学 | 一种无挡片红外热像仪的非均匀性校正方法 |
CN106846292A (zh) * | 2016-12-14 | 2017-06-13 | 中国资源卫星应用中心 | 一种无效像元在轨检测方法 |
CN107036715A (zh) * | 2017-03-30 | 2017-08-11 | 智来光电科技(苏州)有限公司 | 一种红外图像无挡片非均匀性校正装置及其校正方法 |
CN107121202A (zh) * | 2017-05-27 | 2017-09-01 | 珠海格力电器股份有限公司 | 红外传感器温度校验方法、装置、红外传感器及设备 |
CN107784648A (zh) * | 2017-10-18 | 2018-03-09 | 广东正业科技股份有限公司 | 一种视觉定位装置的位置补偿系统与方法 |
CN108226205A (zh) * | 2017-12-12 | 2018-06-29 | 中国科学院地质与地球物理研究所 | 一种离子或电子扫描成像图像漂移的自适应校正方法 |
CN108225570A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225571A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108240863A (zh) * | 2016-12-23 | 2018-07-03 | 南京理工大学 | 用于非制冷红外热像仪的实时两点非均匀性校正方法 |
CN109813439A (zh) * | 2019-02-02 | 2019-05-28 | 湖南大学 | 一种主动式红外热成像热像图序列处理方法 |
WO2019127059A1 (zh) * | 2017-12-26 | 2019-07-04 | 西安电子科技大学 | 基于引导滤波和高通滤波的红外图像非均匀性校正方法 |
CN105376499B (zh) * | 2015-12-11 | 2019-07-05 | 上海兴芯微电子科技有限公司 | 红外摄像装置的死点校正方法、系统及校正系统 |
CN110298274A (zh) * | 2019-06-18 | 2019-10-01 | Oppo广东移动通信有限公司 | 光学指纹参数升级方法及相关产品 |
CN110553739A (zh) * | 2018-12-19 | 2019-12-10 | 中科和光(天津)应用激光技术研究所有限公司 | 一种红外热成像的无挡片非均匀性校正方法 |
CN111595462A (zh) * | 2020-07-23 | 2020-08-28 | 南京智谱科技有限公司 | 红外成像测温系统标定方法、装置、计算设备及存储介质 |
CN112862902A (zh) * | 2021-02-24 | 2021-05-28 | 中国资源卫星应用中心 | 一种空间线阵相机的相对辐射校正方法 |
CN114140345A (zh) * | 2021-11-05 | 2022-03-04 | 中国安全生产科学研究院 | 一种基于红外图像的非均匀矫正方法、装置 |
CN115237257A (zh) * | 2022-08-04 | 2022-10-25 | 深圳市诺金系统集成有限公司 | 一体机触控操作的实现方法及系统 |
CN118392324A (zh) * | 2024-07-01 | 2024-07-26 | 中国科学院光电技术研究所 | 一种基于背景板扣除内部杂散辐射的红外辐射测量方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108200436B (zh) * | 2017-12-29 | 2019-06-18 | 华中科技大学 | 一种面向线列扫描红外成像非均匀性残差的无损压缩方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL173418A (en) * | 2006-01-29 | 2013-10-31 | Rafael Advanced Defense Sys | Correction of unevenness of characters produced by staring detectors |
CN100535618C (zh) * | 2007-04-19 | 2009-09-02 | 华中科技大学 | 一种红外焦平面阵列非均匀性自适应校正方法 |
CN101776487B (zh) * | 2009-12-31 | 2011-05-18 | 华中科技大学 | 一种红外焦平面非均匀性指纹提取及图像校正方法 |
-
2009
- 2009-12-31 CN CN2009102734969A patent/CN101776486B/zh not_active Expired - Fee Related
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012122788A1 (zh) * | 2011-03-14 | 2012-09-20 | 浙江兆晟科技有限公司 | 图像校正系数的获取方法、非均匀图像校正方法及系统 |
CN103219997A (zh) * | 2013-02-06 | 2013-07-24 | 北京空间机电研究所 | 多谱段多焦面拼接红外探测器控制与信号采样电路 |
CN103219997B (zh) * | 2013-02-06 | 2016-08-10 | 北京空间机电研究所 | 多谱段多焦面拼接红外探测器控制与信号采样电路 |
CN103279947A (zh) * | 2013-05-06 | 2013-09-04 | 西安交通大学 | 一种基于单幅图像灰度特征的红外图像校正方法 |
CN103279947B (zh) * | 2013-05-06 | 2015-12-02 | 西安交通大学 | 一种基于单幅图像灰度特征的红外图像校正方法 |
CN103528690A (zh) * | 2013-09-24 | 2014-01-22 | 电子科技大学 | 一种红外热像仪的非均匀校正方法 |
CN103528690B (zh) * | 2013-09-24 | 2016-03-02 | 电子科技大学 | 一种红外热像仪的非均匀校正方法 |
CN105987758B (zh) * | 2015-02-05 | 2019-07-12 | 南京理工大学 | 一种无挡片红外热像仪的非均匀性校正方法 |
CN105987758A (zh) * | 2015-02-05 | 2016-10-05 | 南京理工大学 | 一种无挡片红外热像仪的非均匀性校正方法 |
CN105466566A (zh) * | 2015-12-05 | 2016-04-06 | 中国航空工业集团公司洛阳电光设备研究所 | 一种红外非均匀性校正实时补偿方法 |
CN105376499B (zh) * | 2015-12-11 | 2019-07-05 | 上海兴芯微电子科技有限公司 | 红外摄像装置的死点校正方法、系统及校正系统 |
CN106846292B (zh) * | 2016-12-14 | 2019-08-27 | 中国资源卫星应用中心 | 一种无效像元在轨检测方法 |
CN106846292A (zh) * | 2016-12-14 | 2017-06-13 | 中国资源卫星应用中心 | 一种无效像元在轨检测方法 |
CN108240863A (zh) * | 2016-12-23 | 2018-07-03 | 南京理工大学 | 用于非制冷红外热像仪的实时两点非均匀性校正方法 |
CN107036715B (zh) * | 2017-03-30 | 2019-06-14 | 智来光电科技(苏州)有限公司 | 一种红外图像无挡片非均匀性校正装置及其校正方法 |
CN107036715A (zh) * | 2017-03-30 | 2017-08-11 | 智来光电科技(苏州)有限公司 | 一种红外图像无挡片非均匀性校正装置及其校正方法 |
CN107121202A (zh) * | 2017-05-27 | 2017-09-01 | 珠海格力电器股份有限公司 | 红外传感器温度校验方法、装置、红外传感器及设备 |
CN107784648B (zh) * | 2017-10-18 | 2021-08-24 | 广东正业科技股份有限公司 | 一种视觉定位装置的位置补偿系统与方法 |
CN107784648A (zh) * | 2017-10-18 | 2018-03-09 | 广东正业科技股份有限公司 | 一种视觉定位装置的位置补偿系统与方法 |
CN108226205A (zh) * | 2017-12-12 | 2018-06-29 | 中国科学院地质与地球物理研究所 | 一种离子或电子扫描成像图像漂移的自适应校正方法 |
WO2019127059A1 (zh) * | 2017-12-26 | 2019-07-04 | 西安电子科技大学 | 基于引导滤波和高通滤波的红外图像非均匀性校正方法 |
CN108225571B (zh) * | 2017-12-31 | 2020-01-10 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225571A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225570A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225570B (zh) * | 2017-12-31 | 2020-01-14 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN110553739B (zh) * | 2018-12-19 | 2021-03-12 | 中科和光(天津)应用激光技术研究所有限公司 | 一种红外热成像的无挡片非均匀性校正方法 |
CN110553739A (zh) * | 2018-12-19 | 2019-12-10 | 中科和光(天津)应用激光技术研究所有限公司 | 一种红外热成像的无挡片非均匀性校正方法 |
CN109813439B (zh) * | 2019-02-02 | 2020-04-03 | 湖南大学 | 一种主动式红外热成像热像图序列处理方法 |
CN109813439A (zh) * | 2019-02-02 | 2019-05-28 | 湖南大学 | 一种主动式红外热成像热像图序列处理方法 |
CN110298274A (zh) * | 2019-06-18 | 2019-10-01 | Oppo广东移动通信有限公司 | 光学指纹参数升级方法及相关产品 |
CN111595462A (zh) * | 2020-07-23 | 2020-08-28 | 南京智谱科技有限公司 | 红外成像测温系统标定方法、装置、计算设备及存储介质 |
CN111595462B (zh) * | 2020-07-23 | 2020-11-13 | 南京智谱科技有限公司 | 红外成像测温系统标定方法、装置、计算设备及存储介质 |
CN112862902A (zh) * | 2021-02-24 | 2021-05-28 | 中国资源卫星应用中心 | 一种空间线阵相机的相对辐射校正方法 |
CN112862902B (zh) * | 2021-02-24 | 2024-05-07 | 中国资源卫星应用中心 | 一种空间线阵相机的相对辐射校正方法 |
CN114140345A (zh) * | 2021-11-05 | 2022-03-04 | 中国安全生产科学研究院 | 一种基于红外图像的非均匀矫正方法、装置 |
CN115237257A (zh) * | 2022-08-04 | 2022-10-25 | 深圳市诺金系统集成有限公司 | 一体机触控操作的实现方法及系统 |
CN115237257B (zh) * | 2022-08-04 | 2023-07-04 | 深圳市诺金系统集成有限公司 | 一体机触控操作的实现方法及系统 |
CN118392324A (zh) * | 2024-07-01 | 2024-07-26 | 中国科学院光电技术研究所 | 一种基于背景板扣除内部杂散辐射的红外辐射测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101776486B (zh) | 2011-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101776486B (zh) | 一种基于红外焦平面探测器非均匀性指纹模式的校正方法 | |
CN101776487B (zh) | 一种红外焦平面非均匀性指纹提取及图像校正方法 | |
Zuo et al. | Scene-based nonuniformity correction algorithm based on interframe registration | |
CN104580894B (zh) | 一种红外焦平面的多点校正方法及系统 | |
CN107255521A (zh) | 一种红外图像非均匀性校正方法及系统 | |
Wang et al. | Projected Barzilai–Borwein method for large-scale nonnegative image restoration | |
CN108230249B (zh) | 基于各向异性的l1范数全变分正则化非均匀性校正方法 | |
CN110503137B (zh) | 基于交叉融合的遥感影像时空融合基础图像对的确定方法 | |
CN105737990A (zh) | 一种基于探测器温度的红外图像非均匀性校正方法及系统 | |
CN112393807B (zh) | 红外图像处理方法、装置、系统以及计算机可读存储介质 | |
CN113705788A (zh) | 基于全卷积神经网络的红外图像温度估计方法及系统 | |
Li et al. | A method for the destriping of an orbita hyperspectral image with adaptive moment matching and unidirectional total variation | |
Shi et al. | Multi-scale thermal radiation effects correction via a fast surface fitting with Chebyshev polynomials | |
Li et al. | A novel stripe noise removal model for infrared images | |
Zhang et al. | Cmos fixed pattern noise removal based on low rank sparse variational method | |
CN112710397B (zh) | 一种基于温度替代的两点校正方法及系统 | |
Lee et al. | Infrared Image Deconvolution Considering Fixed Pattern Noise | |
Zhang et al. | Fuzzy c-means clustering algorithm with deformable spatial information for image segmentation | |
Lv et al. | Statistical scene-based non-uniformity correction method with interframe registration | |
Charbal et al. | Regularised digital-level corrections for infrared image correlation | |
CN110111277B (zh) | 一种平面热像图修补方法及装置 | |
CN108921116B (zh) | 遥感影像时变信息提取方法 | |
Zhou et al. | Study of Surface Temperature Distribution for High-Temperature U75V Rail Steel Plates in Rolling Process by Colorimetry Thermometry | |
Wen et al. | Image sensor correction algorithm for photon transfer curve based on neural network | |
Yang et al. | Novel real-time nonuniformity correction solution for infrared focal plane arrays based on S-curve model |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110810 |