CN104867122A - 一种红外自适应非均匀性校正及细节增强级联处理方法 - Google Patents

Abstract

本发明公开的一种红外自适应非均匀性校正及细节增强级联处理方法，涉及一种基于FPGA的红外自适应非均匀性校正及细节增强级联处理方法，属于红外成像技术领域。本发明在FPGA处理平台上实现自适应非均匀性校正及细节增强级联处理，所述的自适应非均匀性校正通过基于灰度相关的时域高通非均匀性校正算法实现，用于减小非均匀性漂移导致的噪声。所述的细节增强通过基于双平台直方图的双边滤波细节增强算法实现，用于提高红外图像细节清晰度。本发明基于FPGA硬件实现处理延迟小于20行，资源消耗少，并可减小或消除非均匀性漂移导致的噪声，提高红外图像清晰度和成像质量。本发明可应用于带有FPGA处理平台的红外成像系统中。

Description

一种红外自适应非均匀性校正及细节增强级联处理方法

技术领域

本发明涉及一种非均匀性校正及细节增强级联处理方法，尤其涉及一种基于FPGA的红外自适应非均匀性校正及细节增强级联处理方法，属于红外成像技术领域。

技术背景

红外热成像系统广泛应用于军事、工业、农业、医疗等各领域。红外辐射经光学系统进入红外焦平面探测器，探测器输出模拟信号经A\D转换为数字信号进入红外数字图像处理系统处理，生成视频输出成像。目前红外数字图像处理系统中普遍有盲元补偿、非均匀性校正和高动态范围压缩三种处理技术，用非均匀性校正技术去除图像中的固定图案噪声，分为基于场景的自适应校正和基于参考源的定标校正两种，基于场景的自适应校正技术由于不需要在视场中插入参考源，成为非均匀性校正技术的研究重点，使用高动态压缩技术将14bit的红外数据压缩至8bit由显示器输出，分为线性压缩、直方图均衡化等，基于双边滤波的细节增强技术在压缩图像数位的同时保留住场景的细节成为高动态压缩技术的研究重点。常用现场可编辑门阵列FPGA+数字信号处理器DSP的系统架构形式实现上述功能，目前常用FPGA实现基于参考源的非均匀性校正技术和线性压缩技术，DSP实现直方图变化技术，这些级联处理技术难以应对探测器的非均匀性漂移导致的噪声和红外图像本身细节不清晰的问题，同时算法的处理延迟往往不可控，在FPGA处理平台上实现基于场景的自适应校正技术和基于双边滤波的细节增强技术的级联处理是未来改进的目标。

发明内容

本发明的目的是为了解决红外图像因非均匀性漂移导致的噪声和红外图像本身细节不清晰的问题，同时实现硬件处理延迟小于20行。本发明公开的一种红外自适应非均匀性校正及细节增强级联处理方法，基于FPGA硬件实现处理延迟小于20行，资源消耗少，并可减小或消除非均匀性漂移导致的噪声，提高红外图像清晰度和成像质量。

本发明的目的是通过下述技术方案实现的。

本发明公开的一种红外自适应非均匀性校正及细节增强级联处理方法，在FPGA处理平台上实现自适应非均匀性校正及细节增强级联处理，所述的自适应非均匀性校正通过基于灰度相关的时域高通非均匀性校正算法实现，用于减小非均匀性漂移导致的噪声。所述的细节增强通过基于双平台直方图的双边滤波细节增强算法实现，用于提高红外图像细节清晰度。由于FPGA处理延迟可控，且所述的自适应非均匀性校正及细节增强级联处理方法减少FPGA硬件资源消耗，在实现处理延迟小于20行条件下可提高红外图像质量。

本发明公开的一种红外自适应非均匀性校正及细节增强级联处理方法，包括如下步骤：

步骤一：在FPGA处理平台实现基于灰度相关的时域高通非均匀性校正模块和实现基于双平台直方图的双边滤波细节增强模块。

步骤二：通过基于灰度相关的时域高通非均匀性校正模块减小非均匀性漂移导致的噪声，具体实现方法包括步骤2.1、2.2、2.3：

步骤2.1：红外焦平面探测器单元的入射辐射经A\D转换为数字信号进入FPGA数字处理模块中，从SRAM中读取预校正偏置值，根据时域阈值条件对输入数字信号预校正，满足时域阈值条件时，用SRAM输出的偏置矩阵进行非均匀性预校正，不满足时域阈值条件时，预校正偏置为0。

所述的预校正方法为：

$y_k(i,j,t)=\left\{\begin{array}{cc}{z}_k(i,j,t)+b_{k-1}(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ z_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(1\right)$

其中，z_k(i,j,t)为温度t时第k帧探测元(i,j)入射辐射值，b_k-1(i,j,t)为温度t时第k-1帧得到探测元(i,j)的加性非均匀，y_k(i,j,t)为温度t时第k帧探测元(i,j)预校正结果，所述的时域阈值条件T_th为：

$T_{\mathrm{th}}=d_2\&CenterDot;\frac{Z_{\max }-Z_{\min }}{100}---\left(2\right)$

其中，Z_max和Z_min探测器整幅图像灰度值的最大值和最小值，d₂是空域阈值稀释系数，其数值范围为0.5～5。

步骤2.2：根据预校正结果，由FPGA内部控制单元在行有效期间计算单帧的偏置估计值。用FPGA实现自适应选择滤波器，计算入射辐射的空域估计值。由FPGA内2×S级FIFO首尾级联形成的(2×S+1)×(2×S+1)空域均值滤波器的滤波窗口，由FPGA内部控制单元，在场消隐期间得到空域阈值，不满足空域阈值条件的像素不作为空域均值滤波的样本。

所述的偏置估计方法为：

${\hat{x}}_k(i,j,t)=\frac{\underset{m=-s}{\overset{s}{\mathrm{\&Sigma;}}}\underset{n=-s}{\overset{s}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)\&CenterDot;y_k(i+m,j+n)}{\underset{m=-s}{\overset{s}{\mathrm{\&Sigma;}}}\underset{n=-s}{\overset{s}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)}---\left(3\right)$

${\hat{b}}_k(i,j,t)=y_k(i,j,t)-{\hat{x}}_k(i,j,t)---\left(4\right)$

其中，为温度t时第k帧探测元(i,j)入射辐射值估计值，(m,n)是窗口内像素距离中心像素(i,j)的距离，为温度t时第k帧探测元(i,j)偏置估计值，选择因子δ为：

$\mathrm{\&delta;}(i+m,j+n)=\left\{\begin{array}{cc}1& |I(i+m,j+n)-I(i,j)|<T_{\mathrm{sp}}\\ 0& |I(i+m,j+n)-I(i,j)|>T_{\mathrm{sp}}\end{array}\right.---\left(5\right)$

其中，空域阈值条件T_sp为：

$T_{\mathrm{sp}}=d_1(\mathrm{\&alpha;}\frac{\underset{i=2}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i-1,j,t)|}{(m-1)N}+\mathrm{\&beta;}\frac{\underset{i=1}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=2}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i,j-1,t)|}{M(N-1)})---\left(6\right)$

其中，d₁是空域阈值稀释系数，其在0.5至3间变化，α和β的数值与非均匀性的形式有关，如果非均匀性形状为“横向”条纹，α＝0，β＝1；如果非均匀性形状为“纵向”条纹，α＝1，β＝0；如果非均匀性形状为“网格”或者“水纹”，α＝0.5，β＝0.5。

步骤2.3：根据单帧的偏置值和时域阈值条件，由FPGA内部控制单元在行有效期间计算非均匀性偏置，存入SRAM中，根据该偏置完成基于灰度相关的时域高通非均匀性校正。满足时域阈值条件时，非均匀性偏置值为单帧的偏置估计值的时域平均，不满足时域阈值条件时，非均匀性偏置值为0。步骤2.1和2.3中采用乒乓结构将一块SRAM按地址分为存储上一帧和当前帧的非均匀性偏置值的两块区域，由FPGA内部控制单元生成SRAM控制器。

所述自适应非均匀性校正方法为：

$b(i,j,t)=\left\{\begin{array}{cc}\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\hat{b}}_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ 0& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(7\right)$

$x_k(i,j,t)=y_k(i,j,t)-\hat{b}(i,j,t)---\left(8\right)$

其中，x_k(i,j,t)为温度t时第k帧探测元(i,j)自适应非均匀性校正后的输出，为非均匀性校正的偏置值。

步骤三：非均匀性校正后，通过基于双平台直方图的双边滤波细节增强模块解决红外图像不清晰的问题，具体实现方法包括步骤3.1、3.2：

步骤3.1：根据非均匀性校正后的图像，从SRAM中读取双平台直方图映射表，压缩该图像。由FPGA内部控制单元在行有效期间，统计直方图信息，在场消隐期间，计算更新直方图映射表。采用乒乓结构按地址将一块SRAM划分为存储当前帧直方图、上一帧和当前帧的直方图映射表的三块区域，由FPGA内部控制单元生成SRAM控制器。

步骤3.2：由FPGA内部控制单元在行有效期间进行双边滤波分离出细节图像，自适应线性映射放大细节图像，与步骤3.1中的压缩图像形成合并形成输出图像，完成基于双平台直方图的双边滤波细节增强。由FPGA内FIFO首尾级联形成的双边滤波器的滤波窗口。

所述的细节图像分离方法为：

$B_k(i,j)=\frac{1}{w(i,j)}\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)x_k(i,j)---\left(9\right)$

D_k(i,j)＝x_k(i,j)-B_k(i,j)  (10)

$w(i,j)=\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)---\left(11\right)$

其中，φ为双边滤波器，w(i,j)为双边滤波权重，B_k(i,j)为第k帧探测元(i,j)的双边滤波器输出，D_k(i,j)为分离出的第k帧探测元(i,j)细节图像。

重复上述步骤一至三，处理每个探测元，直到完成整个红外图像的处理。

有益效果：

1、本发明公开的一种红外自适应非均匀性校正及细节增强级联处理方法，在FPGA处理平台上实现自适应非均匀性校正及细节增强级联处理，通过基于灰度相关的时域高通非均匀性校正算法减小非均匀性漂移导致的噪声，通过基于双平台直方图的双边滤波细节增强算法提高红外图像细节清晰度。

2、由于FPGA处理延迟可控，且自适应非均匀性校正及细节增强级联处理方法可减少FPGA硬件资源消耗，便于硬件实现。本发明可采用单片FPGA和两片SRAM来实现基于灰度映射的时域高通非均匀性校正和基于双平台直方图的双边滤波细节增强级联处理，资源消耗情况少：10486 SliceRegisters、12047 Slice LUTs和36 DSP48Es；处理延迟小：小于20行；本发明可应用于带有FPGA处理平台的红外成像系统中。

附图说明

图1是本发明的一种红外自适应非均匀性校正及细节增强级联处理方法的流程图；

图2是本发明的一种红外自适应非均匀性校正及细节增强级联处理方法模块方框图。

具体实施方式

下面结合附图和实施例对本发明做进一步详细的说明。

实施例1：

如图1所示，本实施例的一种红外自适应非均匀性校正及细节增强级联处理方法，红外探测器入射辐射经A\D转换为数字信号进入FPGA数字图像处理系统中，经基于灰度相关的时域高通非均匀性校正模块和实现基于双平台直方图的双边滤波细节增强模块先后处理，校正因非均匀性漂移导致的噪声，解决红外图像本身细节不清晰的问题。

本实施例的一种红外自适应非均匀性校正及细节增强级联处理方法，包括如下步骤：

步骤一：在FPGA处理平台实现基于灰度相关的时域高通非均匀性校正模块和实现基于双平台直方图的双边滤波细节增强模块。

步骤二：通过基于灰度相关的时域高通非均匀性校正模块减小非均匀性漂移导致的噪声，具体实现步骤包括：

步骤2.1：对于第k帧探测元(i,j)，从SRAM中读取预校正偏置值，根据时域阈值条件对输入数字信号预校正，满足时域阈值条件时，用SRAM输出的偏置矩阵进行非均匀性预校正，不满足时域阈值条件时，预校正偏置为0。

所述的预校正方法为：

$y_k(i,j,t)=\left\{\begin{array}{cc}{z}_k(i,j,t)+b_{k-1}(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ z_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(12\right)$

其中，z_k(i,j,t)为温度t时第k帧探测元(i,j)入射辐射值，b_k-1(i,j,t)为温度t时第k-1帧得到探测元(i,j)的加性非均匀，y_k(i,j,t)为温度t时第k帧探测元(i,j)预校正结果，所述的时域阈值条件T_th为：

$T_{\mathrm{th}}=d_2\&CenterDot;\frac{Z_{\max }-Z_{\min }}{100}---\left(13\right)$

其中，Z_max和Z_min探测器整幅图像灰度值的最大值和最小值，d₂是空域阈值稀释系数，其数值范围为0.5～5。

步骤2.2：根据预校正结果，由FPGA内部控制单元在行有效期间计算单帧的偏置估计值。用FPGA实现自适应选择滤波器，计算入射辐射的空域估计值。由FPGA内6级FIFO首尾级联形成的7×7空域均值滤波器的滤波窗口，由FPGA内部控制单元，在场消隐期间得到空域阈值，不满足空域阈值条件的像素不作为空域均值滤波的样本。

所述的偏置估计方法为：

${\hat{x}}_k(i,j,t)=\frac{\underset{m=-3}{\overset{3}{\mathrm{\&Sigma;}}}\underset{n=-3}{\overset{3}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)\&CenterDot;y_k(i+m,j+n)}{\underset{m=-3}{\overset{3}{\mathrm{\&Sigma;}}}\underset{n=-3}{\overset{3}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)}---\left(14\right)$

${\hat{b}}_k(i,j,t)=y_k(i,j,t)-{\hat{x}}_k(i,j,t)---\left(15\right)$

其中，为温度t时第k帧探测元(i,j)入射辐射值估计值，(m,n)是窗口内像素距离中心像素(i,j)的距离，为温度t时第k帧探测元(i,j)偏置估计值，选择因子δ为：

$\mathrm{\&delta;}(i+m,j+n)=\left\{\begin{array}{cc}1& |I(i+m,j+n)-I(i,j)|<T_{\mathrm{sp}}\\ 0& |I(i+m,j+n)-I(i,j)|>T_{\mathrm{sp}}\end{array}\right.---\left(16\right)$

其中，空域阈值条件T_sp为：

$T_{\mathrm{sp}}=d_1(\mathrm{\&alpha;}\frac{\underset{i=2}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i-1,j,t)|}{(m-1)N}+\mathrm{\&beta;}\frac{\underset{i=1}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=2}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i,j-1,t)|}{M(N-1)})---\left(17\right)$

其中，d₁是空域阈值稀释系数，其在0.5至3间变化，α和β的数值与非均匀性的形式有关，如果非均匀性形状为“横向”条纹，α＝0，β＝1；如果非均匀性形状为“纵向”条纹，α＝1，β＝0；如果非均匀性形状为“网格”或者“水纹”，α＝0.5，β＝0.5。

步骤2.3：根据单帧的偏置值和时域阈值条件，由FPGA内部控制单元在行有效期间计算非均匀性偏置，存入SRAM中，根据该偏置完成基于灰度相关的时域高通非均匀性校正。满足时域阈值条件时，非均匀性偏置值为单帧的偏置估计值的时域平均，不满足时域阈值条件时，非均匀性偏置值为0。步骤2.1和2.3中采用乒乓结构将一块SRAM按地址分为存储上一帧和当前帧的非均匀性偏置值的两块区域，由FPGA内部控制单元生成SRAM控制器。

所述自适应非均匀性校正方法为：

$b(i,j,t)=\left\{\begin{array}{cc}\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\hat{b}}_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ 0& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(18\right)$

$x_k(i,j,t)=y_k(i,j,t)-\hat{b}(i,j,t)---\left(19\right)$

其中，x_k(i,j,t)为温度t时第k帧探测元(i,j)自适应非均匀性校正后的输出，为非均匀性校正的偏置值。

步骤三：非均匀性校正后，通过基于双平台直方图的双边滤波细节增强模块解决红外图像不清晰的问题，具体实现步骤包括：

步骤3.1：根据非均匀性校正后的图像，从SRAM中读取双平台直方图映射表，压缩该图像。由FPGA内部控制单元在行有效期间，统计直方图信息，在场消隐期间，计算更新直方图映射表。采用乒乓结构按地址将一块SRAM划分为存储当前帧直方图、上一帧和当前帧的直方图映射表的三块区域，由FPGA内部控制单元生成SRAM控制器。

步骤3.2：由FPGA内部控制单元在行有效期间进行双边滤波分离出细节图像，自适应线性映射放大细节图像，与步骤3.1中的压缩图像形成合并形成输出图像，完成基于双平台直方图的双边滤波细节增强。由FPGA内FIFO首尾级联形成的双边滤波器的滤波窗口。

所述的细节图像分离方法为：

$B_k(i,j)=\frac{1}{w(i,j)}\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)x_k(i,j)---\left(20\right)$

D_k(i,j)＝x_k(i,j)-B_k(i,j)  (21)

$w(i,j)=\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)---\left(22\right)$

其中，φ为双边滤波器，w(i,j)为双边滤波权重，B_k(i,j)为第k帧探测元(i,j)的双边滤波器输出，D_k(i,j)为分离出的第k帧探测元(i,j)细节图像。

重复上述步骤一至三，处理每个探测元，直到完成整个红外图像的处理。

如图2所示，本实施例的一种红外自适应非均匀性校正及细节增强级联处理方法，具体实现装置包括基于FPGA实现的基于灰度相关的时域高通非均匀性校正模块、基于双平台直方图的双边滤波细节增强模块和外部存储器SRAM。基于FPGA实现的基于灰度相关的时域高通非均匀性校正模块和基于双平台直方图的双边滤波细节增强模块用于实现非均匀性校正和双边滤波细节增强级联处理，用于减小或消除非均匀性漂移导致的噪声，提高红外图像清晰度和成像质量。2片外部存储器SRAM用来存储直方图映射表和偏置校正矩阵。本实施例基于FPGA硬件实现处理延迟小于20行，资源消耗少，便于硬件实现。

本实施例可采用Xilinx公司的Virtex5 XQ5VLX50T芯片硬件实现，包括时钟管理模块、SRAM控制器模块、基于灰度值映射的时域高通非均匀性校正模块、双平台直方图压缩模块、基于双边滤波的细节增强模块。其外围电路包括有源晶体振荡器、电源模块、9Mbits静态存储器GS880Z18T-150。

本发明保护范围不仅局限于实施例，本实施例用于解释本发明，凡与本发明在相同原理和构思条件下的变更或修改均在本发明公开的保护范围之内。

Claims (5)

1.一种红外自适应非均匀性校正及细节增强级联处理方法，其特征在于：在FPGA处理平台上实现自适应非均匀性校正及细节增强级联处理，所述的自适应非均匀性校正通过基于灰度相关的时域高通非均匀性校正算法实现，用于减小非均匀性漂移导致的噪声；所述的细节增强通过基于双平台直方图的双边滤波细节增强算法实现，用于提高红外图像细节清晰度。

2.如权利要求1所述的一种红外自适应非均匀性校正及细节增强级联处理方法，其特征在于：具体实现包括如下步骤，

步骤一：在FPGA处理平台实现基于灰度相关的时域高通非均匀性校正模块和实现基于双平台直方图的双边滤波细节增强模块；

步骤二：通过基于灰度相关的时域高通非均匀性校正模块减小非均匀性漂移导致的噪声；

步骤三：非均匀性校正后，通过基于双平台直方图的双边滤波细节增强模块提高红外图像清晰度；

重复上述步骤一至三，处理每个探测元，直到完成整个红外图像的处理。

3.如权利要求2所述的一种红外自适应非均匀性校正及细节增强级联处理方法，其特征在于：

所述的步骤二具体实现方法包括步骤2.1、2.2、2.3：

步骤2.1：红外焦平面探测器单元的入射辐射经A\D转换为数字信号进入FPGA数字处理模块中，从SRAM中读取预校正偏置值，根据时域阈值条件对输入数字信号预校正，满足时域阈值条件时，用SRAM输出的偏置矩阵进行非均匀性预校正，不满足时域阈值条件时，预校正偏置为0；

所述的预校正方法为：

$y_k(i,j,t)=\left\{\begin{array}{cc}{z}_k(i,j,t)+b_{k-1}(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ z_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(1\right)$

其中，z_k(i,j,t)为温度t时第k帧探测元(i,j)入射辐射值，b_k-1(i,j,t)为温度t时第k-1帧得到探测元(i,j)的加性非均匀，y_k(i,j,t)为温度t时第k帧探测元(i,j)预校正结果，所述的时域阈值条件T_th为：

$T_{\mathrm{th}}=d_2\&CenterDot;\frac{Z_{\max }-Z_{\min }}{100}---\left(2\right)$

其中，Z_max和Z_min探测器整幅图像灰度值的最大值和最小值，d₂是空域阈值稀释系数，其数值范围为0.5～5；

步骤2.2：根据预校正结果，由FPGA内部控制单元在行有效期间计算单帧的偏置估计值；用FPGA实现自适应选择滤波器，计算入射辐射的空域估计值；由FPGA内2×S级FIFO首尾级联形成的(2×S+1)×(2×S+1)空域均值滤波器的滤波窗口，由FPGA内部控制单元，在场消隐期间得到空域阈值，不满足空域阈值条件的像素不作为空域均值滤波的样本；

所述的偏置估计方法为：

${\hat{x}}_k(i,j,t)=\frac{\underset{m=-s}{\overset{s}{\mathrm{\&Sigma;}}}\underset{n=-s}{\overset{s}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)\&CenterDot;y_k(i+m,j+n)}{\underset{m=-s}{\overset{s}{\mathrm{\&Sigma;}}}\underset{n=-s}{\overset{s}{\mathrm{\&Sigma;}}}\mathrm{\&delta;}(i+m,j+n)}---\left(3\right)$

${\hat{b}}_k(i,j,t)=y_k(i,j,t)-{\hat{x}}_k(i,j,t)---\left(4\right)$

其中，为温度t时第k帧探测元(i,j)入射辐射值估计值，(m,n)是窗口内像素距离中心像素(i,j)的距离，为温度t时第k帧探测元(i,j)偏置估计值，选择因子δ为：

$\mathrm{\&delta;}(i+m,j+n)=\left\{\begin{array}{cc}1& |I(i+m,j+n)-I(i,j)|<T_{\mathrm{sp}}\\ 0& |I(i+m,j+n)-I(i,j)|>T_{\mathrm{sp}}\end{array}\right.---\left(5\right)$

其中，空域阈值条件T_sp为：

$T_{\mathrm{sp}}=d_1(\mathrm{\&alpha;}\frac{\underset{i=2}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=1}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i-1,j,t)|}{(M-1)N}+\mathrm{\&beta;}\frac{\underset{i=1}{\overset{M}{\mathrm{\&Sigma;}}}\underset{j=2}{\overset{N}{\mathrm{\&Sigma;}}}|y_k(i,j,t)-y_k(i,j-1,t)|}{M(N-1)})---\left(6\right)$

其中，d₁是空域阈值稀释系数，其在0.5至3间变化，α和β的数值与非均匀性的形式有关，如果非均匀性形状为“横向”条纹，α＝0，β＝1；如果非均匀性形状为“纵向”条纹，α＝1，β＝0；如果非均匀性形状为“网格”或者“水纹”，α＝0.5，β＝0.5；

步骤2.3：根据单帧的偏置值和时域阈值条件，由FPGA内部控制单元在行有效期间计算非均匀性偏置，存入SRAM中，根据该偏置完成基于灰度相关的时域高通非均匀性校正；满足时域阈值条件时，非均匀性偏置值为单帧的偏置估计值的时域平均，不满足时域阈值条件时，非均匀性偏置值为0；步骤2.1和2.3中采用乒乓结构将一块SRAM按地址分为存储上一帧和当前帧的非均匀性偏置值的两块区域，由FPGA内部控制单元生成SRAM控制器；

所述自适应非均匀性校正方法为：

$b(i,j,t)=\left\{\begin{array}{cc}\frac{1}{K}\underset{k=1}{\overset{K}{\mathrm{\&Sigma;}}}{\hat{b}}_k(i,j,t)& |z_k(i,j,t)-z_{k-1}(i,j,t)|<T_{\mathrm{th}}\\ 0& |z_k(i,j,t)-z_{k-1}(i,j,t)|>T_{\mathrm{th}}\end{array}\right.---\left(7\right)$

$x_k(i,j,t)=y_k(i,j,t)-\hat{b}(i,j,t)---\left(8\right)$

其中，x_k(i,j,t)为温度t时第k帧探测元(i,j)自适应非均匀性校正后的输出，为非均匀性校正的偏置值；

所述的步骤三具体实现方法包括步骤3.1、3.2：

步骤3.1：根据非均匀性校正后的图像，从SRAM中读取双平台直方图映射表，压缩该图像；由FPGA内部控制单元在行有效期间，统计直方图信息，在场消隐期间，计算更新直方图映射表；采用乒乓结构按地址将一块SRAM划分为存储当前帧直方图、上一帧和当前帧的直方图映射表的三块区域，由FPGA内部控制单元生成SRAM控制器；

步骤3.2：由FPGA内部控制单元在行有效期间进行双边滤波分离出细节图像，自适应线性映射放大细节图像，与步骤3.1中的压缩图像形成合并形成输出图像，完成基于双平台直方图的双边滤波细节增强；由FPGA内FIFO首尾级联形成的双边滤波器的滤波窗口；

所述的细节图像分离方法为：

$B_k(i,j)=\frac{1}{w(i,j)}\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)x_k(i,j)---\left(9\right)$

D_k(i,j)＝x_k(i,j)-B_k(i,j)    (10)

$w(i,j)=\underset{(i,j)\&Element;S}{\mathrm{\&Sigma;}}\mathrm{\&phi;}(i,j)---\left(11\right)$

其中，φ为双边滤波器，w(i,j)为双边滤波权重，B_k(i,j)为第k帧探测元(i,j)的双边滤波器输出，D_k(i,j)为分离出的第k帧探测元(i,j)细节图像。

4.如权利要求3所述的一种红外自适应非均匀性校正及细节增强级联处理方法，其特征在于：所述的步骤2.2中由FPGA内2×S级FIFO首尾级联形成的(2×S+1)×(2×S+1)空域均值滤波器的滤波窗口，优选由FPGA内6级FIFO首尾级联形成的7×7空域均值滤波器的滤波窗口。

5.如权利要求3或4所述的一种红外自适应非均匀性校正及细节增强级联处理方法，其特征在于：具体实现装置包括基于FPGA实现的基于灰度相关的时域高通非均匀性校正模块、基于双平台直方图的双边滤波细节增强模块和外部存储器SRAM；基于FPGA实现的基于灰度相关的时域高通非均匀性校正模块和基于双平台直方图的双边滤波细节增强模块用于实现非均匀性校正和双边滤波细节增强级联处理；外部存储器SRAM用来存储直方图映射表和偏置校正矩阵。