CN107255521B - 一种红外图像非均匀性校正方法及系统 - Google Patents
一种红外图像非均匀性校正方法及系统 Download PDFInfo
- Publication number
- CN107255521B CN107255521B CN201710507734.2A CN201710507734A CN107255521B CN 107255521 B CN107255521 B CN 107255521B CN 201710507734 A CN201710507734 A CN 201710507734A CN 107255521 B CN107255521 B CN 107255521B
- Authority
- CN
- China
- Prior art keywords
- image
- correction
- assistant images
- radiation intensity
- detector
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 173
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000005855 radiation Effects 0.000 claims abstract description 62
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 6
- 238000003702 image correction Methods 0.000 claims 1
- 238000003331 infrared imaging Methods 0.000 abstract description 10
- 238000004422 calculation algorithm Methods 0.000 description 17
- 238000002474 experimental method Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005057 refrigeration 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- 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
- G01J2005/0077—Imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
Abstract
本发明公开了一种红外图像非均匀性校正方法及系统,涉及红外图像处理领域,包括步骤:对待校正图像进行一次校正,得到第一图像;采集待校正图像在多个不同辐射强度时的均匀背景图像,并对每幅不同辐射强度时的均匀背景图像进行校正,得到多幅不同辐射强度时的辅助图像;计算辅助图像的偏置本底,根据加权法建立当前场景下的补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解补偿校正模型中辅助图像的加权系数,得到当前场景下的补偿校正公式,并对第一图像进行补偿校正,得到第二图像。本发明提供的红外图像非均匀性校正方法,可以消除探测器温度和场景变化对红外成像的影响,提高红外成像质量。
Description
技术领域
本发明涉及红外图像处理领域,具体涉及一种红外图像非均匀性校正方法及系统。
背景技术
随着红外成像技术的不断发展,已广泛应用于民用、军事等多种领域。在红外成像过程中,红外图像的非均匀性容易受探测器温度、图像场景变化的影响,严重降低红外探测器的成像质量,因此,需要对红外图像进行非均匀性校正,消除外界因素对成像质量的影响。
红外图像非均匀性校正通常采用无快门定标算法,现有的无快门定标算法较多,参见专利号为GB2445254A提出的基于相邻像素差平方和最小的无快门校正算法,该算法首先需要在探测器工作温度范围内采集几十幅均匀本底,然后计算其对应的偏置本底,在实时校正时,通过计算所有偏置本底的加权和来获得实时的偏置校正参数,加权系数可根据图像场景的变化不断地进行自适应更新,但该算法只考虑了探测器温度对图像质量的影响,并未考虑外界场景变化对图像非均匀性的影响,仅实现了基于探测器温度的红外图像的非均匀性校正,该算法校正后的成像质量较差。
另参见专利号为CN201210232737.7提出的一种基于模板法的无快门非均匀性校正算法,该算法选择一幅固定的偏置本底来计算校正参数,随着探测器实际温度和该偏置本底对应的探测器温度差异越来越大,计算误差也会越来越大,导致校正后的图像质量越差。
现有技术中的无快门校正算法得到的红外图像质量较差,需要提出一种新的算法,来提高红外成像质量。
发明内容
针对现有技术中存在的缺陷,本发明的目的在于提供一种红外图像非均匀性校正方法,可以消除探测器温度和场景变化对红外成像的影响,提高红外成像质量。
为达到以上目的,本发明采取的技术方案是:
本发明提供了一种红外图像非均匀性校正方法,包括步骤:
采集探测器在不同温度下的均匀背景图像,采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到一幅第一图像;
采集待校正图像在多个不同辐射强度时的均匀背景图像,并对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
采用单点校正法计算辅助图像的偏置本底,根据加权法建立当前场景下的补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解所述补偿校正模型中辅助图像的加权系数,得到当前场景下的补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
在上述技术方案的基础上,所述基于探测器温度的非均匀性校正方法包括步骤:
计算探测器的图像增益系数及多个温度点下的偏置本底,采用插值方式估算探测器在当前温度下的偏置本底,得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
在上述技术方案的基础上,所述基于探测器温度的非均匀性校正公式具体为:
Yij=Gainij(Xij+Offsetij(t))
式中,Gainij是图像坐标(i,j)处的增益系数,Xij是待校正图像坐标(i,j)处的像素灰度值,Yij是一次校正后的图像坐标(i,j)处的像素灰度值,Offsetij(t)是探测器在当前工作温度下图像坐标(i,j)处的偏置本底。
在上述技术方案的基础上,采用单点校正法计算辅助图像的偏置本底的计算公式为:
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数。
在上述技术方案的基础上,所述当前场景下的补偿校正公式为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数。
在上述技术方案的基础上,求解所述补偿校正模型中辅助图像的加权系数具体过程为:
根据补偿校正模型,得到校正后图像和期望图像的像素差平方和S,计算公式为:
式中,i=1,2…M,j=1,2…N,M为图像坐标的行数,N为图像坐标的列数,Zij是第二图像在(i,j)处的像素灰度值,为Zij的期望值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φ1),oij(Φq)分别为辐射强度为Φ1,Φq的辅助图像的偏置本底,分别为oij(Φ1),oij(Φq)的期望值,a1,aq为辐射强度为Φ1,Φq的辅助图像的加权系数,q为辅助图像的个数;
使平方和S的导数为0,求解所有辅助图像的加权系数a1,a2…aq。
本发明还提供了一种红外图像非均匀性校正系统,包括:
一次校正模块,其用于采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到第一图像;
采集模块,其用于采集待校正图像在多个不同辐射强度时的均匀背景图像;
辅助处理模块,其用于对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
补偿校正模块,其用于采用单点校正法计算辅助图像的偏置本底,根据加权法建立当前场景下的补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解所述补偿校正模型中辅助图像的加权系数,得到当前场景下的补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
在上述技术方案的基础上,所述一次校正模块包括:
计算模块,用于计算探测器的图像增益系数及多个温度点下的偏置本底,并采用插值方式估算探测器在当前温度下的偏置本底;
校正处理模块,用于根据探测器的图像增益系数和探测器在当前温度下的偏置本底得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
在上述技术方案的基础上,采用单点校正法计算辅助图像的偏置本底的计算公式为:
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数。
在上述技术方案的基础上,所述补偿校正公式为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数。
与现有技术相比,本发明的优点在于:
(1)本发明的红外图像非均匀性校正方法,先对红外图像进行一次校正,消除探测器温度变化对红外图像的影响,再对其进行补偿校正,消除图像场景变化对红外成像质量的影响,从而使得校正后的红外成像质量好。
(2)本发明的红外图像非均匀性校正方法,计算简单,运算量小,易于实现。
附图说明
图1为本发明实施例中红外图像非均匀性校正方法的流程图;
图2为本发明实施例中第一组实验的计算结果图;
图3为本发明实施例中第二组实验的计算结果图。
具体实施方式
以下结合附图及实施例对本发明作进一步详细说明。
本发明实施例提供一种红外图像非均匀性校正方法,包括步骤:
采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到一幅第一图像;
采集待校正图像在多个不同辐射强度时的均匀背景图像,并对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
采用单点校正法计算辅助图像的偏置本底,根据加权法建立补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解补偿校正模型中辅助图像的加权系数,得到补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
在本发明实施例中,基于探测器温度的非均匀性校正方法为:计算探测器的图像增益系数及多个温度点下的偏置本底,采用插值方式估算探测器在当前温度下的偏置本底,得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
参见图1所示,本发明实施例的红外图像非均匀性校正方法,详细步骤包括:
S1:计算增益系数:探测器在同一温度下采集待校正图像在两个不同辐射强度时的均匀背景图像,并将这两幅均匀背景图像分别记为BlackH和BlackL,计算探测器的图像增益系数;
式中,(i,j)表示图像坐标位置,M、N分别表示图像坐标的行数和列数,Gainij是图像坐标(i,j)处的增益系数,BlackHij是均匀背景图像BlackH在坐标(i,j)处的像素灰度值,BlackLij是均匀背景图像BlackL在图像坐标(i,j)处的像素灰度值;
S2:计算探测器在各个温度点下的偏置本底:探测器在不同温度下采集待校正图像在相同辐射强度时的均匀背景图像,并记录相应的若干探测器温度点,探测器温度点的个数为L,探测器温度点记为tn,探测器温度点tn的均匀背景图像记为F(tn),计算探测器在各个温度点下的偏置本底公式如下:
式中,Fij(tn)是温度点tn的均匀背景图像F(tn)在坐标(i,j)处的像素灰度值,为温度点tn的均匀背景图像F(tn)所有坐标点的像素灰度均值,M、N分别表示图像的行数和列数,Gainij是图像坐标(i,j)处的增益系数,offsetij(tn)为温度点tn下图像坐标(i,j)处的偏置本底;
S3:估算探测器在当前工作温度下的偏置本底:测量探测器当前温度,从若干探测器温度点中选取与当前温度最近的P个探测器温度点,得到相应探测器温度点的偏置本底,采用插值方式计算探测器在当前温度下的偏置本底Offsetij(t),计算公式如下;
式中,offset(T1),…offset(TP)是与温度点T1,T2…TP对应的偏置本底,offsetij(tn)为温度点tn下图像坐标(i,j)处的偏置本底,Kn(t)为探测器温度点Tn对应的插值基函数,
式中,t为探测器的当前工作温度;
S4:一次校正:根据图像增益系数和探测器在当前温度下的偏置本底,得到基于探测器温度的非均匀性校正公式,并使用基于探测器温度的非均匀性校正公式对待校正图像进行一次校正,得到第一图像;
Yij=Gainij(Xij+Offsetij(t))
式中,Gainij是图像坐标(i,j)处的增益系数,Xij是待校正图像坐标(i,j)处的像素灰度值,Yij是一次校正后的图像坐标(i,j)处的像素灰度值,Offsetij(t)是探测器在当前工作温度下图像坐标(i,j)处的偏置本底;
S5:校正不同辐射强度时的均匀背景图像:探测器在同一温度下采集待校正图像在q个不同辐射强度时的均匀背景图像,记录该探测器温度点tt和相应的q个辐射强度,将q个不同的辐射强度分别记为Φ1,Φ2,…Φq,将与之对应的均匀背景图像分别记为f(Φ1),f(Φ2),…f(Φq),
根据基于探测器温度的非均匀性校正公式,对每幅不同辐射强度时的均匀背景图像分别进行校正,相应得到q幅不同辐射强度的辅助图像,得到辅助图像的计算公式为:
ffij(Φn)=Gainij(fij(Φn)+Offsetij(tt))
式中,fij(Φn)是均匀背景图像f(Φn)在坐标(i,j)处的像素灰度值,tt为探测器温度点,Offsetij(tt)为探测器在当前温度tt下的偏置本底,Gainij是图像坐标(i,j)处的增益系数,ffij(Φn)是辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值;
S6:采用单点校正法计算辅助图像的偏置本底;
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数;
S7:得到当前场景下的补偿校正公式:根据加权法计算当前场景下实时的偏置本底,计算公式为:
Oij=a1×oij(Φ1)+a2×oij(Φ2)+…+ap×oij(Φp)
式中,a1,a2…aq为加权系数,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,Oij为当前场景下实时的偏置本底;
建立补偿校正模型为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数;
通过使补偿校正后图像和期望图像的像素差平方和最小求解补偿校正模型中辅助图像的加权系数a1,a2…aq,得到补偿校正公式;
S8:补偿校正:根据补偿校正公式,对第一图像进行补偿校正,得到第二图像。
上述步骤S3中,插值温度点个数P的选取需要根据采集均匀本底时的温度间隔和探测器本身的特性而定。由经验得知,当温度间隔大约为5℃时,在保证图像质量的条件下,P最大取值为5。
步骤S7的求解补偿校正模型中辅助图像的加权系数具体过程为:
S71:根据补偿校正模型,得到校正后图像和期望图像的像素差平方和S,计算公式为:
式中,i=1,2…M,j=1,2…N,M为图像坐标的行数,N为图像坐标的列数,Zij是第二图像在(i,j)处的像素灰度值,为Zij的期望值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φ1),oij(Φq)分别为辐射强度为Φ1,Φq的辅助图像的偏置本底,分别为oij(Φ1),oij(Φq)的期望值,a1,aq为辐射强度为Φ1,Φq的辅助图像的加权系数,q为辅助图像的个数;
S72:分别对每个加权系数均进行求导,并使之为0,求导公式如下:
为了方便化简,令:
其中的简写;
化简求导公式,得到:
其中
把上式代入计算当前场景下实时偏置本底的公式中,计算出实时的偏置本底。
本发明实施例还提供了一种红外图像非均匀性校正系统,包括:
一次校正模块,其用于采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到第一图像;
采集模块,其用于采集待校正图像在多个不同辐射强度时的均匀背景图像;
辅助处理模块,其用于对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
补偿校正模块,其用于采用单点校正法计算辅助图像的偏置本底,根据加权法建立补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解补偿校正模型中辅助图像的加权系数,得到补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
一次校正模块包括:
计算模块,用于计算探测器的图像增益系数及多个温度点下的偏置本底,并采用插值方式估算探测器在当前温度下的偏置本底;
校正处理模块,用于根据探测器的图像增益系数和探测器在当前温度下的偏置本底得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
采用单点校正法计算辅助图像的偏置本底的计算公式为:
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数。
补偿校正公式为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数。
通常剩余非均匀性NUES是用来衡量红外探测器焦平面阵列响应非均匀性的常用指标,被看作是焦平面阵列响应非均匀性的数字定义,当作为红外图像质量的一种评价标准时,只适用于背景辐射均匀的红外图像,其定义如下:
其中,Vi,j是探测器焦平面第i行,第j列的响应输出,d为焦平面阵列中的死像元数,h为焦平面阵列中的过热象元数,Vavg是有效像素的响应平均值,M和N分别代表红外图像的行数和列数。
采用红外非制冷型探测器型号是ULIS PICO384P,探测器的分辨率是384×288,工作温度范围是-40到60℃,分别对专利号为GB2445254A提出的基于相邻像素差平方和最小的无快门校正算法(简称为MSSE算法)、专利号为CN201210232737.7提出的一种基于模板法的无快门非均匀性校正算法(简称为TBS算法)、本发明实施例中的红外图像非均匀性校正方法进行多组实验对比。
第一组实验:场景温度为20℃,探测器工作温度由低温到高温连续变化时采集的31幅均匀背景图像,分别通过上述三种算法对同一红外图像进行非均匀性校正,计算结果如图2所示。
经本发明实施例中的红外图像非均匀性校正方法校正后,红外图像的剩余非均匀性NUES是最小的,并且随着探测器工作温度(探测器焦平面温度)的变化,校正效果也比较稳定,波动较小。
由于TBS算法在每个温度段只选用一幅固定的均匀本底作为参考,很难精确计算出探测器在当前工作温度下的偏置本底,当探测器当前工作温度和选用的均匀本底对应的探测器温度相差较小,计算偏置本底的误差就越小,剩余非均匀性NUES的值越小,图像质量也就越好,反之,两者的温差越大,计算误差就越大,剩余非均匀性NUES的值越大,图像质量也就越差,导致TBS算法校正后的红外图像的剩余非均匀性NUES呈周期性变化,稳定性差。
经MSSE校正后的图像剩余非均匀性曲线基本上没有太大的波动,也比较稳定,但与本发明实施例中的红外图像非均匀性校正方法校正后红外图像的剩余非均匀性NUES相比,仍然偏大,校正效果偏差。
可见,采用本发明实施例中的红外图像非均匀性校正方法校正后,红外图像的剩余非均匀性NUES最小,而且稳定性好,图像质量更好,较好地抑制了探测器温度对图像非均匀性的影响。
第二组实验:探测器工作温度在20℃附近,场景温度在5℃到40℃之间连续变化,分别通过上述三种算法对同一红外图像进行非均匀性校正,计算结果如图3所示。
经本发明实施例中的红外图像非均匀性校正方法校正后,红外图像的剩余非均匀性NUES是最小的,校正后图像质量好,并且随着场景温度的变化,校正效果也比较稳定,波动较小。
经TBS校正后图像的剩余非均匀性NUES曲线不再呈现周期性的波动,但稳定性仍然较差,且剩余非均匀性NUES偏大,图像质量较差。
经MSSE校正后的图像剩余非均匀性NUES曲线波动较大,稳定性差,随着场景温度的升高,剩余非均匀性NUES的值逐渐增大,在常温以下,MSSE校正后的图像剩余非均匀性NUES的值小于BS校正后图像的剩余非均匀性NUES,在常温以上,MSSE校正后的图像剩余非均匀性NUES的值大于TBS校正后图像的剩余非均匀性NUES,MSSE算法在校正的过程中没有考虑图像场景温度变化对图像非均匀性的影响。
可见,采用本发明实施例中的红外图像非均匀性校正方法校正后,红外图像的剩余非均匀性NUES最小,而且稳定性好,图像质量更好,较好地抑制了场景温度对图像非均匀性的影响。
由上可知,本发明实施例中的红外图像非均匀性校正方法,先对红外图像进行一次校正,消除探测器温度变化对红外图像的影响,再对其进行补偿校正,消除图像场景变化对红外成像质量的影响,从而使得校正后的红外成像质量好;而且该红外图像非均匀性校正方法,计算简单,运算量小,易于实现。
本发明不局限于上述实施方式,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围之内。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (10)
1.一种红外图像非均匀性校正方法,其特征在于,包括步骤:
采集探测器在不同温度下的均匀背景图像,采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到一幅第一图像;
采集待校正图像在多个不同辐射强度时的均匀背景图像,并对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
采用单点校正法计算辅助图像的偏置本底,根据加权法建立当前场景下的补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解所述补偿校正模型中辅助图像的加权系数,得到当前场景下的补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
2.如权利要求1所述的红外图像非均匀性校正方法,其特征在于,所述基于探测器温度的非均匀性校正方法包括步骤:
计算探测器的图像增益系数及多个温度点下的偏置本底,采用插值方式估算探测器在当前温度下的偏置本底,得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
3.如权利要求2所述的红外图像非均匀性校正方法,其特征在于,所述基于探测器温度的非均匀性校正公式具体为:
Yij=Gainij(Xij+Offsetij(t))
式中,Gainij是图像坐标(i,j)处的增益系数,Xij是待校正图像坐标(i,j)处的像素灰度值,Yij是一次校正后的图像坐标(i,j)处的像素灰度值,Offsetij(t)是探测器在当前工作温度下图像坐标(i,j)处的偏置本底。
4.如权利要求1所述的红外图像非均匀性校正方法,其特征在于,采用单点校正法计算辅助图像的偏置本底的计算公式为:
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数。
5.如权利要求4所述的红外图像非均匀性校正方法,其特征在于,所述当前场景下的补偿校正公式为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数。
6.如权利要求5所述的红外图像非均匀性校正方法,其特征在于,求解所述补偿校正模型中辅助图像的加权系数具体过程为:
根据补偿校正模型,得到校正后图像和期望图像的像素差平方和S,计算公式为:
式中,i=1,2…M,j=1,2…N,M为图像坐标的行数,N为图像坐标的列数,Zij是第二图像在(i,j)处的像素灰度值,为Zij的期望值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φ1),oij(Φq)分别为辐射强度为Φ1,Φq的辅助图像的偏置本底,分别为oij(Φ1),oij(Φq)的期望值,a1,aq为辐射强度为Φ1,Φq的辅助图像的加权系数,q为辅助图像的个数;
使平方和S的导数为0,求解所有辅助图像的加权系数a1,a2…aq。
7.一种红外图像非均匀性校正系统,其特征在于,包括:
一次校正模块,其用于采用基于探测器温度的非均匀性校正方法,对待校正图像进行一次校正,得到第一图像;
采集模块,其用于采集待校正图像在多个不同辐射强度时的均匀背景图像;
辅助处理模块,其用于对每幅不同辐射强度时的均匀背景图像进行基于探测器温度的非均匀性校正,得到多幅不同辐射强度时的辅助图像;
补偿校正模块,其用于采用单点校正法计算辅助图像的偏置本底,根据加权法建立当前场景下的补偿校正模型,通过使补偿校正后图像和期望图像的像素差平方和最小求解所述补偿校正模型中辅助图像的加权系数,得到当前场景下的补偿校正公式,并对第一图像进行补偿校正,得到第二图像。
8.如权利要求7所述的红外图像非均匀性校正系统,其特征在于,所述一次校正模块包括:
计算模块,用于计算探测器的图像增益系数及多个温度点下的偏置本底,并采用插值方式估算探测器在当前温度下的偏置本底;
校正处理模块,用于根据探测器的图像增益系数和探测器在当前温度下的偏置本底得到基于探测器温度的非均匀性校正公式,采用两点校正法对待校正图像进行一次校正。
9.如权利要求7所述的红外图像非均匀性校正系统,其特征在于,采用单点校正法计算辅助图像的偏置本底的计算公式为:
式中,oij(Φn)为辐射强度为Φn的辅助图像的偏置本底,ffij(Φn)为辐射强度为Φn的辅助图像在坐标(i,j)处的像素灰度值,是ffij(Φn)的像素灰度均值,q为辅助图像的个数。
10.如权利要求9所述的红外图像非均匀性校正系统,其特征在于,所述补偿校正公式为:
式中,Zij是第二图像在(i,j)处的像素灰度值,Yij是第一图像在(i,j)处的像素灰度值,oij(Φm)为辐射强度为Φm的辅助图像的偏置本底,am为辐射强度为Φm的辅助图像的加权系数,q为辅助图像的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710507734.2A CN107255521B (zh) | 2017-06-28 | 2017-06-28 | 一种红外图像非均匀性校正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710507734.2A CN107255521B (zh) | 2017-06-28 | 2017-06-28 | 一种红外图像非均匀性校正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107255521A CN107255521A (zh) | 2017-10-17 |
CN107255521B true CN107255521B (zh) | 2019-03-26 |
Family
ID=60024380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710507734.2A Active CN107255521B (zh) | 2017-06-28 | 2017-06-28 | 一种红外图像非均匀性校正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107255521B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110595630B (zh) * | 2019-08-28 | 2020-08-28 | 武汉华中数控股份有限公司 | 基于探测器温度的多点非均匀校正的方法 |
CN110852976B (zh) * | 2019-11-22 | 2023-04-18 | 昆明物理研究所 | 一种红外图像明暗不均匀校正方法及计算机程序产品 |
CN110850500B (zh) * | 2019-11-24 | 2022-02-08 | 北京长峰科威光电技术有限公司 | 一种红外图像多段单点校正参数修正方法 |
CN111507915B (zh) * | 2020-04-13 | 2024-03-29 | 烟台艾睿光电科技有限公司 | 基于模糊配准的实时红外非均匀性校正方法、设备及介质 |
CN112150377B (zh) * | 2020-08-27 | 2024-03-29 | 北京空间机电研究所 | 一种红外图像非均匀性校正系数解析解交替迭代优化方法 |
CN113188664B (zh) * | 2021-04-02 | 2022-09-06 | 烟台艾睿光电科技有限公司 | 红外快速打快门方法、装置以及红外热像仪 |
CN113436088B (zh) * | 2021-06-09 | 2022-07-26 | 浙江兆晟科技股份有限公司 | 一种红外图像热窗效应实时抑制方法及系统 |
CN113724171A (zh) * | 2021-09-14 | 2021-11-30 | 浙江大华技术股份有限公司 | 图像的校正方法、装置、存储介质及电子装置 |
CN114140345A (zh) * | 2021-11-05 | 2022-03-04 | 中国安全生产科学研究院 | 一种基于红外图像的非均匀矫正方法、装置 |
CN115861088B (zh) * | 2022-10-20 | 2023-06-20 | 国科天成科技股份有限公司 | 红外摄像头的非均匀性漂移实时校正方法和系统 |
CN115375590B (zh) * | 2022-10-27 | 2023-04-07 | 之江实验室 | 一种亮度非均匀性在线补偿的图像处理方法 |
CN115760653B (zh) * | 2023-01-09 | 2023-05-09 | 武汉中导光电设备有限公司 | 图像校正方法、装置、设备及可读存储介质 |
CN117314774B (zh) * | 2023-09-26 | 2024-04-09 | 广州市星飞达电子科技有限公司 | 一种基于探测器温度的红外图像校正系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102768071A (zh) * | 2012-07-06 | 2012-11-07 | 武汉高德红外股份有限公司 | 基于模板法的非制冷热像仪无快门非均匀性校正方法 |
CN105737990A (zh) * | 2016-02-24 | 2016-07-06 | 华中科技大学 | 一种基于探测器温度的红外图像非均匀性校正方法及系统 |
WO2016182961A1 (en) * | 2015-05-08 | 2016-11-17 | Flir Systems, Inc. | Isothermal image enhancement systems and methods |
CN106595870A (zh) * | 2016-12-08 | 2017-04-26 | 苏州长风航空电子有限公司 | 一种红外成像装置及其非均匀性处理方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10018504B2 (en) * | 2015-11-12 | 2018-07-10 | Sensors Unlimited, Inc. | Gain normalization |
-
2017
- 2017-06-28 CN CN201710507734.2A patent/CN107255521B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102768071A (zh) * | 2012-07-06 | 2012-11-07 | 武汉高德红外股份有限公司 | 基于模板法的非制冷热像仪无快门非均匀性校正方法 |
WO2016182961A1 (en) * | 2015-05-08 | 2016-11-17 | Flir Systems, Inc. | Isothermal image enhancement systems and methods |
CN105737990A (zh) * | 2016-02-24 | 2016-07-06 | 华中科技大学 | 一种基于探测器温度的红外图像非均匀性校正方法及系统 |
CN106595870A (zh) * | 2016-12-08 | 2017-04-26 | 苏州长风航空电子有限公司 | 一种红外成像装置及其非均匀性处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107255521A (zh) | 2017-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107255521B (zh) | 一种红外图像非均匀性校正方法及系统 | |
CN108846805B (zh) | 一种基于场景自适应的红外热图像两点非均匀校正方法 | |
CN100535618C (zh) | 一种红外焦平面阵列非均匀性自适应校正方法 | |
CN109272520B (zh) | 一种联合运动指导与边缘检测的自适应红外焦平面非均匀校正方法 | |
CN110211056B (zh) | 基于局部中值直方图的自适应红外图像去条纹算法 | |
CN106197673B (zh) | 一种自适应宽温度范围非均匀校正方法及系统 | |
CN104515599B (zh) | 基于fpga的红外焦平面阵列盲元检测系统及方法 | |
CN109685877A (zh) | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 | |
CN102778296A (zh) | 基于总变分的自适应红外焦平面非均匀性校正方法 | |
CN103335724A (zh) | 基于定标的场景自适应irfpa非均匀性校正方法 | |
CN106017695B (zh) | 基于运动状态估计的自适应红外非均匀性校正方法 | |
CN102589707B (zh) | 红外焦平面阵列探测器非均匀性校正残差的实时补偿方法 | |
CN104159001A (zh) | 基于高精度大幅面扫描仪系统的图像多点分区域校正方法 | |
CN110782403A (zh) | 一种红外图像非均匀性校正方法 | |
CN106855435A (zh) | 长波线阵红外相机星上非均匀性实时校正方法 | |
CN105869129B (zh) | 针对非均匀校正后的热红外图像剩余非均匀噪声去除方法 | |
CN111076821A (zh) | 一种积分时间自适应切换的红外图像非均匀校正方法 | |
CN109813442A (zh) | 一种基于多帧处理的内部杂散辐射非均匀性校正方法 | |
CN109741376A (zh) | 一种基于改进ransac算法的中、长波红外图像配准方法 | |
CN110363714A (zh) | 自适应学习速率的基于场景帧间配准的非均匀性校正方法 | |
CN104268870A (zh) | 基于小波变换直方图的短波红外焦平面非均匀性校正算法 | |
CN108055487B (zh) | 像感器阵列不均匀性的一致校正方法及系统 | |
CN103868601B (zh) | Irfpa探测器非均匀响应的双边全变分正则化校正方法 | |
WO2019183843A1 (zh) | 基于帧间配准和自适应步长的红外图像非均匀性校正方法 | |
CN104050661A (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 |