红外非均匀性校正的算法研究

1、IRFPA响应的非均匀性校正的算法研究CONTENTS一 IRFPA响应的非均匀性产生原因二 常见的非均匀性校正方法的原理三 两点标定法及改进思路四 自适应算法的介绍和改进五 盲元替代的实现附录一：用于实现两点标定法和盲元处理的MFC源程序参考文献一 IRFPA响应的非均匀性产生原因为了对IRFPA响应的非均匀性进行校正，我们有必要先对它的定义和成因进行介绍。 红外成像过程可以描述为目标和背景的红外辐射通过大气和光学系统传输后到达红外焦平面阵列（infrared focal plane arrays, IRFPA），红外探测器把辐射信号转换为电信号，然后经过读出电路输出显示的过程。因此，IRF

2、PA响应输出是目标辐射特性、大气传输特性、光学系统特性、器件响应特性等诸多因素共同作用的结果。依此，我们可以把影响红外成像质量的因素分为四种：响应的非均匀性，响应的漂移特性，盲元，目标辐射的对比度。对于目标辐射的对比度，由观测对象决定，但可以通过图像增强的方法进行改善。本文主要讨论前三种。对于理想的IRFPA，假定输入辐射均匀且相同的话，那么每个探测器的输出信号应该完全相同。但事实上，由于制作工艺，材料质量等因素的影响，每个探测器在阻抗，容抗，热敏面积，电阻温度系数等参数方面均有差别，因此输出信号幅度并不相同，即产生固定图案噪声( Fixed Pattern Noise , FPN )，这种不

3、一致就是IRFPA响应的非均匀性。1999年中华人民共和国国家标准红外焦平面阵列参数测试技术规范中的关于非均匀性(NU)的定义如下：(1-1)其中表示红外焦平面阵列上所有有效像元的输出信号平均值(在计算输出信号的和以及非均匀性时，均不包括无效像元的信号值)，M和N分别为焦平面阵列的行数和列数，d为焦平面阵列中的死像元数，h为焦平面阵列中的过热像元数。盲元包括死像素和过热像素。死像素指响应率小于平均响应率的110的像素；而过热像素就是响应率大于平均响应率10倍的像素。它严重影响了红外系统的成像质量.因此，在实际应用中，必须要先对其进行非均匀校正(Non-uniformity Correction

4、, NUC) 。不同的IRFPA响应的非均匀性并不相同。常用的HgCdTe和InSb型IRFPA的非均匀性为10%左右，而测微辐射热计IRFPA的非均匀性可达20%。若不校正的话，我们就无法获得清晰的图像。与此同时，非均匀性还会随着时间和工作条件（如环境温度、偏置电压）的变化而产生漂移。影响漂移的因素间也是相互关联，很难对其中一个单独校正。工作条件的改变是对整个器件施加的影响，因此具有一种整体性变化趋势，这种趋势表现为直流成分，可以通过图像增强的方法消除；对于每一个探测元，工作条件的改变引起的变化略有不同，其响应的漂移也略有差别，这种不同探测器的差别性变化可以通过校正系数的更新来消除。盲元是指

5、响应过高或过低的探测单元，在图像中表现为过亮点和过暗点。如果某个探测器的非均匀性无法通过校正消除的话，就认为它是盲元。盲元数量和分布直接影响探测器的成像质量，若盲元多而集中，图像就会出现无响应区域，称为盲区，严重限制了红外成像器件的应用。因此，只有开发有效的盲元检测算法对其进行补偿，才能更好发挥探测器的性能。 对于盲元的定义，主要是从器件对黑体辐射的响应程度作为量化指标的。国家标准从器件对黑体辐射的响应程度作为盲元的量化指标，包括死像元和过热像元。死像元是指像元响应率小于110平均响应率的像元；过热像元是指像元噪声电压大于平均噪声电压10倍的像元。像元响应率和平均响应率的定义如下。像元响应率：

6、像元响应率R(i ,j)是单位辐射功率产生的输出信号电压，表达式如下： (1-2)上式中，V(i,j)是第i行第j列像元对于辐射功率P的响应电压，P表示红外焦平面阵列接收的辐射功率。平均响应率：平均响应率表示红外焦平面阵列各个有效像元响应率的平均值，表达式如下： (1-3)上式中：M、N分别表示红外焦平面阵列的行数和列数；d、h分别表示死像元数和过热像元数。实际测量中，d和h是经多次迭代计算得到。盲元率：红外焦平面阵列的盲元数占总像元数的百分比，由下式表示： (1-4)二常见的非均匀性校正方法的原理非均匀性校正的方法有很多，大致可以分为两类：基于标定技术的算法和基于场景技术的算法。基于标定技术

7、的算法是指在实验室里利用均匀的高温和低温黑体对红外焦平面进行标定，对指定焦平面器件进行参数的提取，从而计算出增益和偏移系数，对探测器输出进行校正的方法，此类算法常见的有的两点标定法（TPC）、多点温度校正法（ETPC）等，其中的两点标定法最为常见，也是本文着重介绍的一种算法。定标校正法具有较高的校正精度、技术成熟、硬件实现容易、校正效果明显、图像较为均匀。其缺点是校正时，需要许多辅助器件(如：黑体源、光学设备等) 定标，因此大大增加了探测器的体积和成本；在标定过程中成像系统还必须暂停工作；此外它还没有补偿各探测元差别性漂移的能力，对探测器长时间工作的漂移需要周期性采样定标。这严重制约了该类算法

8、的使用范围； 基于场景技术的算法,它直接利用每帧图像的场景信息及统计特性，获得每个像素的校正系数，进行非均匀校正。基于场景的非均匀性校正算法有：神经网络法(ANNC)、时域高通滤波法(THPFC)、恒定统计法(CSC)、代数校正算法、Kalman滤波法等。这类算法克服了标定法的不足，能够随着IRFPA非均匀性的变化自适应地对图像进行校正，其校正参数的提取不需要对系统进行干预，系统正常使用过程中自动对当前探测器输出的非均匀性进行参数的统计、计算和提取。但其技术还不成熟，好的算法实现起来繁琐，硬件要求资源多、速度快，能实现的算法需要经过近似、简化，校正效果也不尽相同但往往具有较大的计算量和存储量,

9、因此很难实时处理；收敛速度慢，收敛时间长，算法的实时性较差。三 两点标定法两点标定法建立在两个假设条件下：(1)每个探测单元的响应是线性的。这一假设在探测单元响应的大部分区域是成立的，只有在辐射输入很小和接近饱和时才有比较大的非线性，在忽略两端误差的情况下可以认为探测单元的响应是线性的。根据线性模型的数学公式，只需对直线上两点进行定标测量，便可求出直线，从而对非均匀性进行校正。在此条件下，红外焦平面阵列在均匀辐射背景下任一像元的响应输出可以表示为： (3-1)其中为辐射通量，为像元(i,j)的响应，为其增益，为其偏移系数。在假设(1)的条件下，IRFPA的响应非均匀性体现为每个像元的增益和偏移

10、系数的不同。(2)探测单元的响应分布相同，并具有时间的稳定性，否则，定标数据在使用时就会失去意义。非均匀校正的基本思想就是是相同辐射条件下的不同像素的响应曲线重合于同一条曲线。为此可先设定一条基准曲线，将各探测元的响应曲线分别作旋转和平移使之与基准曲线重合，这个旋转和平移系数即为该像元的校正系数。下图给出了两点校正法的原理图，两点法校正将所有像元在高温和低温下的响应分别规格化为和。图3.1 两点校正法示意图（直线a和b分别代表一个像元的响应特性，直线c代表校正后所有像元的响应特性，既基准曲线）假定基准响应曲线在输入辐射分别是和时的响应为和，像元（i,j）对应的响应分别为和，则可以求出像元（i,

11、j）的增益校正系数和偏移校正系数分别为 (3-1)而像元（i，j）经两点校正后的输出表示为 (3-2)式中为像元（i,j）未校正前的输出，为两点非均匀校正后的输出。由该式我们可以看出获得任意像元的校正系数并不需要知道实际的增益和偏移量，而只需要知道和。基准响应和可以任意选取，通常为了简便，我们取所有像素在输入辐射分别为和时的响应平均值。利用两点标定法获取校正系数的过程大致为：(1). 调节热成像仪系统，使焦平面通过光学系统与平面黑体源对准，并使辐射充满焦平面的整个视场。(2). 调节黑体辐射源的温度到低温；测量N组焦平面每个像元的响应值，求出每个像素在下的平均值以及所有像元的响应平均值。(3)

12、. 调节黑体辐射源的温度到高温；测量N组焦平面每个像元的响应值，求出每个像素在下的平均值以及所有像元的响应平均值。(4). 利用算法计算每个像元的增益和偏移系数，并将其存储在查找表内，以供校正时使用。根据上述过程，对所给定的4251文件进行处理后得到的图像如下：图3.2 使用MFC经两点标定法校正后图像两点校正法的缺点及改进思路：根据假定(1),两点校正法假定像元的响应特性在所感兴趣的温度范围内是线性的，实际情况并非如此。当目标温度和定标温度差距较大时，将会有较大的残留空间噪声,而这会影响成像质量，因此适合于非线性不严重的焦平面探测器或小动态范围成像系统。为弥补两点校正法的这一不足，必须把温度

13、定标点设置在目标温度附近，也可以可进一步采用选择多温度点进行多点校正。根据假定(2),探测单元的响应分布相同，并具有时间的稳定性，但长时间的使用后，将会产生相应的漂移。对于各像元的差别性漂移，它表现为增益和偏移系数的细微变化，因此利用标定系数获得的永久性校正系数将不再使用，此时我们就必须重新进行标定。这时我们可以通过电机波动放置在IRFPA前面的黑体挡板来重新标定。四 其他算法的介绍基于标定技术的算法无法永久适用于各像元的差别性漂移，从而使得我们长时间使用时必须进行标定。这在实际使用中带来很大不便，自适应非均匀性校正方法是非均匀性校正发展的必然趋势。基于场景的非均匀性校正算法有：神经网络法(A

14、NNC)、时域高通滤波法(THPFC)、恒定统计法(CSC)、代数校正算法、Kalman滤波法等。本文仅介绍时域高通滤波法(THPFC)、神经网络法(ANNC)这两种算法。1990年以来，美国海军研究实验室的D.A.Scribner等人基于人眼视觉神经机理的现有研究成果，提出了两种自适应非均匀校正算法。一种是时域高通滤波校正算法(Temporal High-pass Filter, THPF-NUC)，另一种是一神经网络计算为核心的人工神经网络算法(artificial neural networks, ANNC）。 下面分别介绍其原理和优缺点。4.1 时域高通滤波校正算法可以认为空中环境尤其

15、是温度对于探测器而言是个慢变化量，而图像本身包括目标、背景噪声是个快变化量，因而图像可以看作是高频部分，而探测器漂移的非均匀性是探测器的固有的响应，分布在低频(变化较慢)的部分，它能用时域低通表现出来，从而可以利用增益补偿和高通滤波算法来完成红外焦平面阵列的非均匀性校正。图4.1 THPF-NUC原理图假设探测器的响应输出为，经过增益补偿后变为，低通滤波器输出，即为固定图案噪声(FPN)，则校正后图像像元输出为 (4-1) 其中的低通采样输出为： (4-2)这样,在不显著增加计算复杂性的前提下,通过简单的迭代处理完成FPN 参数初值的计算, 避免了使用黑体所导致成像系统光机结构复杂、 易受周围

16、环境温度影响的缺陷。其中N为预先设置好的帧数，G为增益补偿系数。可以先假设某像素在第n帧经增益校正后得到一个滤波器输入量，经过无限冲击响应滤波器低通处理后得到了低通输出，再通过减法运算得到一个高通输出。经z变换得到传递函数为： (4-3)令(为采样时间)，则可以求出系统的频率响应为： (4-4)可求出截止频率为： (4-5)其中为焦平面阵列图像的帧速率。积累帧数N越大，截止频率越小，在实际应用中可以根据具体要求选取合适的N值，过小的积累帧采用时域高通滤波将会破坏到运动的目标。时域高通滤波法的优点：时域高通滤波法的算法简单可以实时实施；不需要定标；合理选择时间常数，空间噪声可以降低到时域噪声水平

17、；对低频率的非均匀性效果显著。时域高通滤波法的缺点：只进行了偏移的校正，需要极好的增益均匀性，或是先进行增益的定标校正；时域高通滤波法要求视场中的景物做随机运动，否则图像会退化，这在有些应用场合是一个非常苛刻的条件；这种方法的难点之一是设置的帧数N的确定。较小的N值会较快达到稳定，但是对景物的随机性要求也更强。反之，较大的N值会使系统需要较长时间才能稳定，但景物的短时非随机的影响也小一些。过大的N值有可能使系统的稳定时间超过探测器的漂移速度，从而无法彻底补偿系统的漂移。时域高通滤波法可能会减弱信号。对于静止场景而言，基于时域高通的算法可能会将场景与固定噪声一起通过低通滤波器去除，产生图像的消隐

18、现象。另外，因为低通滤波器是通过时域上的加权累加进行更新，新加入的帧的权重低于以往的累加，因此这样可以避免异常的噪声对滤波器性能的影响。但是，当存在长时间静止的亮目标突然产生运动，新加入的场景亮度较原有亮度差距较大时，由于系数更新缓慢，就会在局部出现鬼影现象，鬼影现象是自适应更新系数的校正算法中共有的问题。4.2 人工神经网络校正算法可以完全不对FPA进行标定(或自动标定)，是红外成像系统的理想境界。依赖于神经网络方法自动实时地进行校正系数的更新是目前实验室研究的热点之一。对人眼视觉研究表明，人眼视网膜的各个光接收元之间总存在一定的差异，但人可根据具体的景像通过视网膜上水平细胞元进行调节，每一

19、个水平细胞元都与邻近的几个光接收元连在一起。受此启发，国外有学者采用周期性神经网络结构进行自适应的非均匀性校正，该法采用一个隐含层计算某像素邻域输出的平均值，并以此作为该像素的输出，回馈给线性校正神经元计算NUC系数。具体办法是：让每一个神经元连接一个阵元，再设计一个隐含层，它的每一个神经元就像水平细胞元那样与邻近的几个阵元连接起来，得到它们的平均输出值反馈到它的上层神经元去计算非均匀性，采用最陡下降法依据实际景像逐帧迭代，直至达到最佳校正状态。 D.A.Scribner 提出的 BP 人工神经网络算法用于非均匀性校正的实现模型，其原理如下图所示。图4.2 BP 神经网络校正算法原理图可见，网

20、络由输入层、NUC 层、一个隐含层和输出层四层组成，是 BP 人工神经网络算法应用于非均匀性校正的具体实现形式。它是经典的BP人工神经网络和二点定标校正算法的结合，即在选用最简单的三层BP人工神经网络算法结构的基础上，结合红外非均匀性校正的实际情况，引入了二点定标NUC（非均匀性校正）层。因此这种实现形式既具备了人工神经网络的自适应性（即：参数不断向着最优化调整） ，又结合了红外非均匀性校正的具体问题，使得算法更具针对性。(1)输入层的输入为连续帧的二维图像灰度值，(i,j)表示像素所在的行值、列值。(2)隐含层的作用是对图像进行线性平滑，平滑的结果作为期望的输出图像。因为红外图像具有低对比度

21、的特性，因此相邻的阵列元所接收的红外辐射应当比较接近，由于各阵列元的响应特性参数近似不相关，因此可以把某一像素的邻域平均作为该像素的理想输出。对于第 n 帧图像的第i行，第 j列的像素，取其四邻域像素的灰度值进行平均，作为该像素的期望输出值。即 (4-6)并将期望输出值反馈至NUC层。(3)NUC层除了要对输入灰度值进行二点非均匀性校正，如果实际输出与期望输出间的误差大于设定的误差阈值，还要利用隐含层反馈的期望输出值，对增益校正系数和偏置校正系数进行修正。NUC层行二点非均匀性校正的公式如下： (4-7)式中n为输入图像的帧序号，为增益校正系数，为偏置校正系数。修正过程如下：定义关于、的误差函

22、数为（为简便起见，以下公式均略去下标）： (4-8)根据最小梯度法，此函数的梯度为: (4-9)则趋进于最小误差的最陡下降路径为（即沿着负梯度方向）: (4-10) 其中是迭代步长。从的表达式可以看出误差函数的形状为抛物谷面，谷面底部满足方程，即 (4-11) 由于人工神经网络具有鲁棒性，所以当足够小时，就可以保证迭代过程稳定收敛，即误差值在误差函数抛物谷面上逐渐向着谷面底部靠近。(4)输出层输出经过校正的二维灰度图像，即 (4-12)理论分析和试验结果表明，神经网络算法是稳定和收敛的，即能连续修正每一像元的增益和偏置系数，直至稳定为止。与定标方式的校正算法相比，自适应校正算法更能满足实际需要

23、。即只要一次定标测量，就能连续修正光敏元的增益和偏置系数，因此对探测器的变化具有更强的适应性；减小了空间噪声电平，降低了红外焦平面阵列的非均匀性。当背景变化时可以跟踪探测器单元响应的非线性；空间噪声可以降低到低于时域噪声。神经网络算法的缺点：算法的计算量比较大，实时实现有困难；结构比较复杂，硬件实现有困难；与时域高通滤波算法相类似，这种方法同样要求视场内的景物作随机运动，否则会引起图像退化；在寻求最优的过程中，可能陷入局部最小值而达不到全局最优解，以及收敛较慢、需要较多训练次数；对于空间频率低的非均匀性效果不明显。4.4 对人工神经网络算法的进一步分析和改进 对上述BP人工神经网络算法进行分析

24、，可以得出，算法实际上由三部分组成：（1）二点定标线性校正算法；（2）线性平滑算法；（3）系数修正算法。其算法原理框图如图4.4.1所示。图4.4.1 BP算法的原理框图下面分析这三部分算法。(1) 二点定标线性校正算法这部分算法本质上是和二点定标法是完全相同的，只是其增益校正系数和偏置校正系数是处于不断的自适应调整过程，使得输出在外界条件变化时仍可接近期望值。(2) 线性平滑算法 BP人工神经网络算法中的线性平滑算法与普通的线性平滑算法并不相同，它没有采取直接代替的方法，而是将线性平滑的结果作为期望输出值，来校正二点定标线性校正算法中的增益校正系数和偏置校正系数。这相当于对图像进行了二维低通

25、滤波，所以当目标长时间固定在某一个像素点上时，将会使目标模糊。邻域大小可任意选择，但将直接影响算法的复杂度、计算量等性能。(3) 系数修正算法BP 人工神经网络的系数修正算法实质上就是最小均方误差算法（LMS）。最小均方误差算法（LMS）的工作原理就是把增益和偏移校正系数作为权系数，把期望响应和输出信号的差值作为误信号，根据误信号的大小自动调节权系数，从而使得滤波输出信号的均方误差最小。在LMS滤波器中，期望响应对性能具有重要影响，但它又是不可获知的，因此必须选用某值代替。常用方法是把目标在运动轨迹上均值作为理想期望。由于运动轨迹上的无规则性，LMS校正后的剩余非均匀性近似等于过去轨迹运动后所

26、有像元平均后的非均匀性。这要求我们保持系统与场景间的相对运动，否则将会退化为像元自身的响应，滤波器将失效。在分析了经典的BP 人工神经网络算法之后，我们可以依次进行算法思路的改进。改进思路一：在算法试验仿真中，注意到经BP校正算法中的增益校正系数（简称“增益”）的调整量是偏置校正系数（简称“偏置”）调整量的倍，这就意味着在LMS算法中增益的调整量要远大于偏置的调整量，而通常增益处于以 1 为中心的一个区间内，这种调整量数量级上的差别导致算法实现中步长值难以选取。按照经典BP算法进行校正，当步长取很小时，会造成收敛缓慢；而当步长取过大时，又造成迭代不收敛。这使算法在应用于实际时遇到了障碍。针对这

27、一缺点，人们提出了归一化BP人工神经网络校正算法，既利用图像的灰度平均值对图像的灰度数据进行归一化， 使其成为1附近的数值，对归一化的灰度数据，利用经典BP校正算法中的增益和偏置调节公式进行非均匀性校正，在处理完成后，对数据进行还原，使其恢复到原始的数量级。可以看出此时对增益和偏置进行调节，我们可以获得相同数量级的调整量，这样就避免了迭代步长a 难以选取的问题，克服了BP算法实际应用的障碍。可见，归一化BP人工神经网络算法，处理的目的就是为了消除增益和偏置的调整量的数量级的差别。归一化BP校正算法具体实现如下。在进行BP校正前，对需要处理的每一帧图像的灰度值进行归一化： (4-13)为当前帧的

28、灰度均值。对归一化的图像灰度值利用经典BP算法进行非均匀性校正。而后，对归一化输出结果进行还原，实际校正输出为： (4-14)至此归一化BP校正算法对一帧图像的处理过程结束。此外，为了保证算法尽可能地迭代收敛到最优值，我们可以设定增益和偏置量的动态范围，当调整后的增益和偏置量超出其动态范围，用下式进行循环回调，直至其处于动态范围之内。 (4-15)改进思路二：BP人工神经网络算法中的线性平滑算法采用邻域平均法估计真实的图形值，但是在图像不同灰度的边界处，这种方法会产生很大的偏差。因此我们可以根据邻域像元的灰度为其分配不同的权系数，然后通过加权求和得到估计的真实值。即根据两像元灰度级差值的倒数确

29、定权系数，二者灰度级越接近，差值就越小，权值就越大，这就是所谓的“梯度倒数加权法”。此外，如果当前像元位于图像的边界处时，周围的像元可能与当前像元处于不同区域，我们可以将周围像元值与当前像元值进行比较，如果其差值大于选定的阈值Th,就认为它与当前像元处于不同区域，将其权系数置为0。上述改进是基于像元的空间相关性，我们还可以利用图像序列的时间相关性估算图像的真值，确定其权系数。采用相邻的三帧图像，在每帧图像内采用最近邻的四个像素和当前像素的数据对图像的真值进行估计，用表示第K帧图像的地(i,j)像元的原始值，表示第K帧图像的地(i,j)像元的校正值，为像元的权系数，可以得到以下的计算公式： (4

30、-16)其中l=k-2,k-1,k, (4-17) (4-18) (4-19)其中sigh（x）为符号函数，Th是选定阈值。当前像素和它最近的4个像素的校正值均送入隐含层，通过比较估计的真值和实际输出值的误差，来修正上一层的增益系数和偏置。图 采用单帧图像估计真值图像的神经网络结构图针对上述这些算法的实现步骤，分析ANNC算法中步长和误差阈值大小对算法的影响。(1) 在远离算法收敛的终点（最优值）时，不论误差平面是陡是缓，较大的步长值可使算法快速收敛，而且较大的步长值甚至有助于跳过崎岖的局部极小值，如图所示图 较大步长跨越局部最小的情况(2) 在接近收敛的终点时，较大的步长会形成振荡，导致算法

31、不收敛。如图所示图 较大步长形成振荡的情况(3) 小的步长值固然能保证算法的收敛，但是过小的步长使算法收敛缓慢，在期望的时间内得不到较好的校正效果(4) 误差阈值的选择也很重要，大的误差阈值可使算法很快的达到稳定的输出，但是校正输出的效果较差，就是说，在没有达到最优值时，就停止了修正过程，而不能得到最优值。过大的误差阈值还可使算法完全失效，起不到任何校正作用。小的误差阈值情况下，收敛较为缓慢，但是校正输出的非均匀性较好，而过小的误差阈值会引起算法的振荡，而得不到稳定的输出。总而言之，步长和误差阈值的选择，对于算法来说就是稳定性和效率二者的选择。在实践中，常常根据具体情况和要求对二者进行折衷。五

32、. 盲元替代的实现红外焦平面阵列像元响应曲线如下图所示,图5.1 像元响应曲线曲线A和C分别表示过热像元和死像元温度响应，它们的响应输出不随红外焦平面阵列接收的辐射功率而改变，曲线B表示正常像元的温度响应，随着接收的辐射功率增加，像元的响应随之改变，并且在较窄的动态范围内，像元响应和辐射功率呈线性关系。目前，常用的盲元检测算法有定义法、直方图法、小波滤波法、定标法等。盲元包括死像素和过热像素。IRFPA中的盲元与正常像素在响应特性上有很大差异，在一定的动态范围内，正常像素的温度响应特性曲线是线性的，而盲元的动态范围远离正常的像素动态范围，盲元的温度响应特性曲线多为非线性的，并且变化斜率偏高或者

33、偏低，对两个不同温度的响应，其两点差值偏离正常的像素的两点差值。针对盲元响应特性的这一特点，周慧鑫等提出了一种基于两点参考辐射源的盲元自动检测技术。这种算法的缺陷在于比较门限值很难确定，对于不同的两幅图像，门限值也可能不同。因此，这给硬件实现带来困难。最直接的盲元测试技术就是根据盲元的定义进行测试。盲元的定义是按IRFPA器件对黑体辐射的响应程度作为量化指标，它是基于实验室测量结果给出的关系式。因此，采用均匀辐射的黑体源照射IRFFA，得到一组响应数据，由于死像素和过热像素是随机的，其统计平均响应值应接近为零，总的响应平均值与剔除盲元后的响应平均值基本相同。因此，可采用总的响应平均值与各个像素

34、的响应值相比较。确定盲元的位置。在盲元补偿方面，现有的方法有相邻像元替代法、中值滤波法和插值法等。这里主要介绍插值法和中值滤波法。5.1 插值法盲元补偿技术是根据相邻像素或前后帧图像的响应相关性对盲元位置的信息进行预测和替代的过程。红外成像系统主要是对景物的实时动态成像，根据图像信息理论，其相邻两点或两帧之间的图像数据具有极高的相关性。根据插值理论，对于连续变化的函数，任何一点都可以通过邻域进行插值预测。插值分为线性插值和非线性插值。非线性插值计算复杂，对于实时性要求很高的红外热成像系统显然并不合适。目前均采用线性插值技术对盲元进行替代。假设像素（i,j）为盲元，则替代方程为 (5-1) 式中

35、为时空邻域或空间邻域，为邻域权系数，。权系数要根据邻域情况做出选择，如果邻域内某个像素也同为盲元，则其对应的权系数应该为0，其权值应按比例分配给其他邻域像素。替代方程计算量与邻域大小成正比，实际应用中，应选择合适的邻域。常见的方式为帧内4邻域替代方式，其邻域为盲元的上下左右4个像素。如果邻域内不存在盲元，则其替代方程为 (5-2)如果邻域内存在盲元，则替代方程的权系数必须按上述原则调整。盲元替代的实现方式对红外成像系统的性能具有一定影响。从检索到的文献来看，目前均采用中断方式实现盲元替代。盲元的位置已事先做成查找表存储于热成像系统内部。当红外图像被校正后，首先被保存下来。在两帧图像的间隙，由中

36、断触发盲元替代程序，逐一查表替代盲元。中断程序结束以后，红外图像才能输出。传统的盲元替代方式不仅增加了存储空间的要求，也导致图像输出存在12帧的延迟。在设计的热成像系统内部，将采用更为简化的盲元替代方式，利用盲元的左右像素的平均值替代盲元，其方程为 (5-3)利用该方程，我们发展了一种新型的替代结构，它处于整个信号处理流水线的一环，由FPGA来实现。系统为每个像素分配一个标志位，表明其是否为盲元。每个像素的标识位和它的增益校正系数存储在一个字内（因为所用AD为14位，所以增益校正系数只需15位，剩余的一位作为盲元标识位）。工作时，标识位一直跟随流水线移动。当校正结束后，标识位判断是否需要替代。

37、这种方式结构简单，不需要把校正后的图像存储起来，降低了存储空间的要求；不存在延迟，提高了成像系统的实时性。5.2 中值滤波法 中值滤波是一种非线性滤波，其输出是观测样本分布中心的最小绝对值均值估计，是一种优化的保持边界与滤除高频噪声干扰的滤波方法，尤其对处理突变尖峰之类的脉冲噪声非常有效。中值滤波器于1971年由Tukey提出并应用于一维信号分析中，后来被二维图像处理技术引用。 中值滤波就是用一个奇数点的移动窗口(中值窗)，将窗口中心点的值用窗口内各点按其数值大小排序后中间的那个数(即中值)代替取中值窗长度L为奇数，对采样序列X(1),X(2),X(N)进行中值滤波，是从输入序列中相继抽出L个

38、数X(i-k), ,X(i-1),X(i),X(i+1), ,X(i+k),其中X(i)为窗口中心点值，k=(L-1)/2；再将这L个点按数值大小排列，取中值作为滤波输出，用公式表示为 (5-4)二维中值滤波可 表示为 (5-5)中值滤波对于某些类型的随机噪声具有非常理想的降噪能力，对于线性平滑滤波而言，在处理像素邻域之内的噪声点时，噪声的存在总会或多或少影响该点的像素计算，但在中值滤波中噪声点则常常直接被忽略掉，而且与线性平滑滤波器相比，中值滤波在降噪同时引起的模糊效应降低。对4251文件，使用MFC进行编程处理盲元后得到的图像如下：图5.2 使用中值滤波进行盲元处理后得到的图像附录一：用于

39、实现两点标定法和盲元处理的MFC源程序热像仪每个像素校正后的系数表void CIR_ImageView:OnAdjust() / TODO: Add your command handler code here / 读入低温数据文件，并存储到数组中 FILE *plow=fopen("B42520.txt","r"); int i, j; for (i = 0; i < 240; i+) for (j = 0; j < 320; j+) fscanf(plow, "%d", &mlowij); fclose(plo

40、w); / 读入高温数据文件，并存储到数组中FILE *phigh=fopen("B42540.txt","r"); for (i = 0; i < 240; i+) for (j = 0; j < 320; j+) fscanf(phigh, "%d", &mhighij);fclose(phigh); /求两组数据的平均值int y1,y2;int temp1=0;int temp2=0;for (i=0;i<240;i+) for (j=0;j<320;j+)temp1=temp1+mlowij;

41、temp2=temp2+mhighij; y1=temp1/(240*320); y2=temp2/(240*320);/求校正系数for (i=0;i<240;i+) for (j=0;j<320;j+) Gij=(float)(y2-y1)/(mhighij-mlowij+0.001); Oij=(float)(y2*mlowij-y1*mhighij)/(mlowij-mhighij+0.001);/将校正系数读入文件中FILE *pG=fopen("G.txt","w");FILE *pO=fopen("O.txt"

42、;,"w"); for (i = 0; i < 240; i+) for (j = 0; j < 320; j+) fprintf(pG, "%2.3f ", Gij);fprintf(pO, "%6.1f ", Oij); fprintf(pG, "%n");fprintf(pO, "%n");fclose(pG);fclose(pO);利用校正系数表对此文件进行图像处理void CIR_ImageView:OnProcess() / TODO: Add your command

43、handler code here/载入数据FILE *pIR=fopen("4251.txt","r"); int i, j; for (i = 0; i < 240; i+)for (j = 0; j < 320; j+)fscanf(pIR, "%d", &mIRij);fclose(pIR); /两点校正for (i=0;i<240;i+) for (j=0;j<320;j+) mIRnewij=(int)(Gij*mIRij+Oij); m_Image.LoadFromFile("t

44、est.bmp");/打开图像 flag=TRUE;/图像打开后，将标记变量置为真m_Image.RgbToGrade();LPBYTElpDst;/ 指向缓存图像的指针LPBYTElpNewDIBBits;/ 指向缓存DIB图像的指针HLOCALhNewDIBBits;LPBYTE lpDIBBits=m_Image.GetData();/找到原图像的起始位置/ 暂时分配内存，以保存新图像hNewDIBBits = LocalAlloc(LHND, 240 * 320);lpNewDIBBits = (LPBYTE )LocalLock(hNewDIBBits);/ 初始化新分配的内存，设定初始值为0lpDst = (LPBYTE)lpNewDIBBits;memset(lpDst, (BYTE)0, 240 * 320);/将数组转化为图像for(i = 0; i<240; i+)for(j = 0;j<320; j+)/ 指向目标图像倒数第i行，第j个像素的指针lpDst