一种数字TDI红外探测器的实时盲元检测方法技术

本发明专利技术公开了一种数字TDI红外探测器的实时盲元检测方法。该方法对数字TDI红外探测器工作时的图像进行实时采集，将当前帧的图像数据按象元顺序与前一帧缓存器的图像数据、特征缓存器中的特征数据进行比较累加，计算新的象元统计特征。当统计到一定程度的象元统计特征后，在扫描方向以行为单位读出象元特征并进行判决，每一次的判决结果累加到盲元概率表中。利用盲元概率表计算TDI累加权重，通过TDI累加权重实时去除盲元在TDI累加过程中的影响，完成盲元补偿。本发明专利技术在不打断TDI成像工作状态，无需额外辅助设备的情况下，准确实时的完成盲元的检测和补偿。

【技术实现步骤摘要】

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

【技术保护点】
一种数字TDI红外探测器的实时盲元检测与补偿方法，其特征在于：所述方法步骤如下：(1)实时采集数字TDI红外探测器工作时的图像数据，其采集到的图像数据为未进行TDI累加前的探测器原始输出信号；对一个M·N的数字TDI红外探测器，每次采集到的图像帧大小为M·N，其中其中M为行数，N为列数，工作时的扫描方向垂直于列方向；(2)将当前帧的图像数据按象元顺序与前一帧缓存器的图像数据、特征缓存器中的特征数据进行比较累加，计算新的象元统计特征并将结果存储回特征缓存器中；其中通过缓存器存储前一帧的图像数据，其存储大小与图像大小相同，特征缓存器存储每一个象元对应的特征数据，其大小由图像大小与所选择的特征数共同决定；常用的象元特征数据特征包括，象元的累积变化长度d，计算公式如下：其中n为帧计数，In(x,y)为位置为(x,y)的象元的在当前图像的响应值，In‑1(x,y)为位置为象元的在前一帧的响应值，dn‑1(x,y)为特征缓存器存储的累积变化长度，dn(x,y)为新的象元累积变化长度；象元的累积正向(或负向)变化次数,计算公式如下：ptn(x，y)＝ptn‑1(x，y)+((In(x，y)‑In‑1(x，y))＞Tp)      (2) 其中ptn(x,y)为特征缓存器存储的象元累积正向变化次数，ptn‑1(x,y)为新的象元累积正向变化次数，Tp为人工设置的变化有效阈值；象元在统计周期的变化幅度，计算公式如下：dvn(x，y)＝max(Imax(x，y)，In(x，y))‑min(Imin(x，y)，In(x，y))      (3) 其中Imax(x,y)，Imin(x,y)为特征缓存器存储的统计周期内象元响应最大值和最小值，dvn(x,y)为新的象元在统计周期的变化幅度；(3)当累加的特征次数达到统计周期Nf时，从特征缓存器中按行读出象元的统计特征送入判决器中，并将特征缓存器中的对应行象元特征清零，重新开始统计；(4)判决器通过对同一行象元的统计特征进行计算判别，决定出该行的盲元与正常象元；判决器对输入的同一行象元特征按下式进行判决处理：其中a为象元统计特征权重，Td为人工设置的盲元判别阈值，D(x,y)为象元的盲元判决结果，其中计算的是第x行象元第k个特征的均值，也可以用第x行象元第k个特征的中值进行替代；(5)将当前象元判决结果与该象元在盲元概率表中存储的盲元概率计算更新盲元概率，并将结果存储回盲元概率表中。盲元概率计算更新公式如下：pm(x,y)＝β1·pm‑1(x,y)+β2·D(x,y)     (5) 其中β1、β2为概率收敛系数，由人工进行设置，且β1+β2＝1，p(x,y)为象元的盲元概率，下标为盲元概率的累加次数；(6)从盲元概率表中读出待补偿象元所在行象元的盲元概率，利用盲元概率生成累加权重对该象元进行补偿。累加权重由下式计算得到：其中Pb为人工设置的盲元概率阈值，W(x,y)为对应象元的累加权重。利用累加权重对该象元进行补偿，计算公式如下：Inew(x,y)＝W(x,y)·I(x,y)        (7) 其中Inew(x,y)为补偿后的象元信号；利用补偿后的信号进行TDI累加，去除盲元在图像中的影响。...

【技术特征摘要】
1.一种数字TDI红外探测器的实时盲元检测与补偿方法，其特征在于：所述方法步骤如下：
(1)实时采集数字TDI红外探测器工作时的图像数据，其采集到的图像数据为未进行TDI累加前的探测器原始输出信号；对一个M·N的数字TDI红外探测器，每次采集到的图像帧大小为M·N，其中其中M为行数，N为列数，工作时的扫描方向垂直于列方向；
(2)将当前帧的图像数据按象元顺序与前一帧缓存器的图像数据、特征缓存器中的特征数据进行比较累加，计算新的象元统计特征并将结果存储回特征缓存器中；其中通过缓存器存储前一帧的图像数据，其存储大小与图像大小相同，特征缓存器存储每一个象元对应的特征数据，其大小由图像大小与所选择的特征数共同决定；
常用的象元特征数据特征包括，象元的累积变化长度d，计算公式如下：
其中n为帧计数，In(x,y)为位置为(x,y)的象元的在当前图像的响应值，In-1(x,y)为位置为象元的在前一帧的响应值，dn-1(x,y)为特征缓存器存储的累积变化长度，dn(x,y)为新的象元累积变化长度；象元的累积正向(或负向)变化次数,计算公式如下：
ptn(x，y)＝ptn-1(x，y)+((In(x，y)-In-1(x，y))＞Tp)(2)
其中ptn(x,y)为特征缓存器存储的象元累积正向变化次数，ptn-1(x,y)为新的象元累积正向变化次数，Tp为人工设置的变化有效阈值；象元在统计周期的变化幅度，计算公式如下：
dvn(x，y)＝max(Imax(x，y)，In(x，y))-min(Imin(x，y)，In(x，y))(3)
其中Imax(x,y)，Imin(x,y)为特征缓存器存储的统计周期内象元响应最大值和最小值，dvn(x,y)为新的象元在统计周期的变化幅度；
(3)当累加的特征次数达到统...

【专利技术属性】
技术研发人员：陈忻，韩冰，苏晓锋，夏晖，李夜金，黄茂潼，饶鹏，孙胜利，
申请(专利权)人：中国科学院上海技术物理研究所，
类型：发明
国别省市：上海;31