一种数字tdi红外探测器的实时盲元检测方法
【技术领域】
[0001] 本发明涉及到红外探测器信号处理技术，具体指一种应用于数字TDI红外探测器 的盲元检测及补偿方法。
【背景技术】
[0002] 红外遥感成像技术作为当今迅速发展的高新技术之一，越来越广泛地应用于气候 预测、自然灾害检测、地球环境检测、导航、农业、天文学、军事等领域。为了满足不断提高的 分辨率要求和与之相对应的光通量的要求，在不降低轨道高度的情况下采用时间延迟积分 (Time Delay Intergration，TDI)技术提高分辨率和信噪比。由于TDI技术不仅具有在不牺 牲空间分辨率和成像系统工作速度的情况下获得高灵敏度，而且它能够不同程度地增加系 统的信噪比和提高探测器响应均匀性的特性使其在空间红外遥感应用中具有广泛的应用 前景。TDI的成像过程如图1所示，TDI探测器通过对同一景物在不同时间多次曝光后进行延 迟累加，N次曝光得到的图像按照一定的顺序进行累加即可得到目标景物的N级TDI结果 [1]。 按照实现形式的不同，TDI红外探测器分为数字TDI与模拟TDI两类，模拟TDI将探测器获得 模拟电信号在电荷域进行累加，实现TDI功能，而数字TDI指的是将探测器获得模拟电信号 先转换成数字信号，之后通过在数字信号域中进行累加以实现TDI的功能，其处理过程如图 2所示。由于数字化得信号可以通过处理进行信号的补偿与校正，因此数字TDI红外探测器 与模拟TDI红外探测器比较具有更高性噪比、更高动态范围并且级数可调使用灵活等优点， 所以数字TDI红外探测器更广泛的被使用。
[0003] 由于制造材料、工艺等因素的影响（如材料的不均匀性、掩模误差、缺陷等），红外 焦平面阵列器件存在不可避免的盲元(包括死像元和过热像元）问题。因此在TDI延迟累加 前需要对探测器的盲元进行定位和补偿，否则会大大降低成像信号的输出信噪比和输出动 态范围。目前常用的盲元检测与补偿算法主要分为两大类 [2]:-类是基于黑体定标的检测 方法，另一类为基于场景的检测方法。基于黑体定标类的方法目前在各类红外成像系统中 得到了最广泛的应用，基于黑体定标的盲元检测算法原理简单，但是需要黑体的配合，无法 处理实际应用中因环境温度变化而随机出现的新盲元，其检测效率较低 [3][4]。基于场景检 测法不依赖额外的设备，能够有效地校正红外探测器本身所固有的盲元和因环境温度变化 而产生的随机盲元。但目前基于于场景的检测方法多在图像空间进行处理和检测，容易受 到临域非均匀性和连续盲元的干扰，且运算量大，消耗存储资源多，难以做到硬件实现和实 时检测处理 [5]。
[0004] [1]谢恺，孙振江等.新型空间红外扫描传感器[J].红外技术，2006，V〇1.28(4): 230-233
[0005] [2]赖睿，刘上乾等.红外焦平面阵列盲元检测技术研究[J].《半导体光电》2005年 6月第26卷第3期
[0006] [3]张晓琳.红外焦平面阵列的多判据盲元检测方法.中国发明专利授权号 CN103076156 B
[0007] [4]隋修宝，陈钱等.红外焦平面阵列盲元检测算法[J].光电工程，2008.8， Vol.35,N0.8:107-111
[0008] [5]顾国华.基于滑动窗口与多帧补偿的自适应盲元检测与补偿算法[J].红外技 术，2010 · 6, Vol · 32,N0 · 7:420:423
【发明内容】
[0009] 为了克服上述现有技术的不足，本发明提供了一种数字TDI红外探测器的实时盲 元检测方法。该方法利用数字TDI红外探测器成像过程中的数据，独立统计象元在成像过程 中的变化差异，利用正常象元与盲元的累积变化差异检测出盲元。具有无需额外辅助设备， 运算量低，资源消耗少等优点，能够在不打断数字TDI红外探测器工作状态的情况下，实现 实时动态的盲元检测和补偿。
[0010] 为了解决上述技术问题，本发明是这样实现的：
[0011] 步骤一实时采集数字TDI红外探测器工作时的图像数据。采集到的图像数据为未 进行TDI累加前的探测器原始输出。即对一个Μ · N的数字TDI红外探测器，每次采集到的图 像帧大小为Μ · Ν，其中其中Μ为行数，Ν为列数，工作时的扫描方向垂直于列方向；
[0012] 步骤二将当前帧的图像数据按象元顺序与前一帧缓存器的图像数据、特征缓存器 中的特征数据进行比较累加，计算新的象元统计特征并将结果存储回特征缓存器中。象元 统计特征主要有:象元的累积变化长度，计算公式如下：
[0013] djx,y) =：dnA(.V, γ) + λ/ι +(/"(-ν,.)·)-/" ,(.ν, r))2 -1 ⑴：
[0014] 其中η为帧计数，In(x,y)为位置为(x,y)的象元的在当前图像的响应值，In-i(x,y) 为位置为象元的在前一帧的响应值，d n-i (X, y)为特征缓存器存储的累积变化长度，dn (X, y) 为新的象元累积变化长度。象元的累积正向（或负向）变化次数，计算公式如下：
[0015] ptn(x,y)=ptn-i(x,y) + ((In(x,y)-In-i(x,y))>TP) (2)
[001 ?]其中ptn (X，y)为特征缓存器存储的象元累积正向变化次数，ptn-1 (X，y)为新的象 元累积正向变化次数，TP为人工设置的变化有效阈值。象元在统计周期的变化幅度，计算公 式如下：
[0017] dvn(x，y) =max(Imax(x，y)，In(x，y) )_min(Imin(x，y)，In(x，y)) (3)
[0018] 其中Imax(x，y)，Imin(x，y)为特征缓存器存储的统计周期内象元响应最大值和最小 值，dv n(x，y)为新的象元在统计周期的变化幅度。在实际应用中采用上述统计特征中的一 种或是多种组合。将象元的统计特征以向量的形式表达如下所示：
[0019] E(X，y，k) = [d(X，y)，pt(X，y)，...，dv(X，y)] (4)
[0020] 其中k为特征索引，其值域为[1，K]，K为所采用的特征数；
[0021] 步骤三当累加的特征次数达到统计周期Nf时，从特征缓存器中按行读出象元的统 计特征送入判决器中，并将特征缓存器中的对应行象元特征清零，重新开始统计；
[0022] 步骤四判决器通过对同一行象元的统计特征进行计算判别，决定出该行的盲元与 正常象元。判决器的判决公式如下：
[0023]
(5)
[0024] 其中a为象元统计特征权重，Td为人工设置的盲元判别阈值，D(x，y)为象元的盲元 判决结果；
[0025] 步骤五将当前象元判决结果与该象元在盲元概率表中存储的盲元概率计算更新 盲元概率，并将结果存储回盲元概率表中。盲元概率计算公式如下：
[0026] pm(x,y)=0i · pm-i(x,y)+^2 · D(x,y) (6)
[0027] 其中仇、此为概率收敛系数，由人工进行设置，且fo+i32 = l，p(X，y)为象元的盲元概 η 率，下标w=为盲元概率的累加次数；
[0028] 步骤六从盲元概率表中读出待补偿象元所在行象元的盲元概率，利用盲元概率生 成累加权重对该象元进行补偿。累加权重由下式计算得到： r(x，v)=产。，W ⑴
[0029] Σ[(Α?(·ν,/)<^] /-1
[0030] 其中Pb为人工设置的盲元概率阈值，W( X, y)为对应象元的累加权重。利用累加权 重对该象元进行补