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
- focal plane
- matrix
- overbar
- infrared focal
- uniformity
- 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 42
- 230000001419 dependent effect Effects 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 68
- 238000010606 normalization Methods 0.000 claims description 11
- 239000000203 mixture Substances 0.000 claims description 9
- 230000001915 proofreading effect Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 6
- 230000008901 benefit Effects 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 238000009432 framing Methods 0.000 claims description 2
- 238000012937 correction Methods 0.000 abstract description 35
- 238000003331 infrared imaging Methods 0.000 abstract description 7
- 238000001514 detection method Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 abstract description 3
- 238000007796 conventional method Methods 0.000 abstract 2
- 230000004044 response Effects 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000013500 data storage Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 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
- 239000004615 ingredient Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004297 night vision Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 239000000523 sample 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 (23)
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 | 广东正业科技股份有限公司 | 一种视觉定位装置的位置补偿***与方法 |
CN108225571A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108226205A (zh) * | 2017-12-12 | 2018-06-29 | 中国科学院地质与地球物理研究所 | 一种离子或电子扫描成像图像漂移的自适应校正方法 |
CN108225570A (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 | 深圳市诺金***集成有限公司 | 一体机触控操作的实现方法及*** |
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 (37)
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 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225570A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225570B (zh) * | 2017-12-31 | 2020-01-14 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
CN108225571A (zh) * | 2017-12-31 | 2018-06-29 | 北京华科德科技有限公司 | 一种短波红外焦平面自适应非均匀校正算法 |
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 | 深圳市诺金***集成有限公司 | 一体机触控操作的实现方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN101776486B (zh) | 2011-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101776486A (zh) | 一种基于红外焦平面非均匀性指纹模式的校正方法 | |
CN101666682B (zh) | 基于场景统计的神经网络非均匀性校正方法 | |
CN101776487B (zh) | 一种红外焦平面非均匀性指纹提取及图像校正方法 | |
CN106017695B (zh) | 基于运动状态估计的自适应红外非均匀性校正方法 | |
Zuo et al. | Improved interframe registration based nonuniformity correction for focal plane arrays | |
US10931901B2 (en) | Method and apparatus for selectively correcting fixed pattern noise based on pixel difference values of infrared images | |
CN106197690B (zh) | 一种宽温范围条件下的图像校准方法及*** | |
CN111968066A (zh) | 红外图像的校正方法、装置、设备及制冷红外成像*** | |
CN104008543A (zh) | 一种图像融合质量评价方法 | |
Jara et al. | Joint de-blurring and nonuniformity correction method for infrared microscopy imaging | |
CN105987758A (zh) | 一种无挡片红外热像仪的非均匀性校正方法 | |
CN113049118B (zh) | 一种红外焦平面探测器闪元测试装置及测试方法 | |
Gladysz | Absolute and differential G-tilt in turbulence: theory and applications | |
CN112710397B (zh) | 一种基于温度替代的两点校正方法及*** | |
CN104599261A (zh) | 短波红外焦平面的均值直方图神经网络非均匀性校正算法 | |
CN105092043A (zh) | 一种基于场景的变积分时间的非均匀性校正方法 | |
CN104065853B (zh) | 一种红外相机串扰消除方法 | |
US20200410645A1 (en) | Image processing apparatus, image processing method, and storage medium | |
Zhao et al. | Fast iterative adaptive nonuniformity correction with gradient minimization for infrared focal plane arrays | |
CN103868601B (zh) | Irfpa探测器非均匀响应的双边全变分正则化校正方法 | |
Yang et al. | Real‐time illumination adjustment for video deflectometers | |
Chrzanowski et al. | Virtual MRTD–an indirect method to measure MRTD of thermal imagers using computer simulation | |
JPWO2018235777A1 (ja) | 画像データ加工装置、植物栽培システム及び画像データ加工方法 | |
CN115115653A (zh) | 一种冷热冲击试验箱精细化温度校准方法 | |
US9152880B1 (en) | Method for modeling human visual discrimination task performance of dynamic scenes |
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 |