基于非抽样Contourlet变换的红外图像增强算法

随着红外技术的迅速发展，它已广泛应用于军事国防、遥感探测、无损检测等众多领域。由于红外图像的成像机理以及红外成像系统自身的原因，红外成像系统的成像效果不够理想，大多红外图像都有对比度低、图像模糊、灰度范围窄的缺点。实际应用中，为了提高红外图像的质量，需要对红外图像进行必要的增强处理。一般图像增强处理的方法是基于空间域和变换域的，前者主要包括直接灰度变换、空间滤波和直方图处理等；后者是将图像由时域变换到频域，再通过修正变换域内的系数达到增强图像的目的,它优于基于空间域的增强方法，代表性算法有小波变换的算法和基于Contourlet变换的算法等[1]。

2002年，DO和VETTERLI提出的Contourlet变换是多尺度几何分析方法中十分重要的一类。Contourlet能够实现一种“真正”二维的图像表示，能够提取在图像中非常重要的内在几何结构特征。类似于小波可以从滤波器角度考虑，Contourlet则利用不可分的滤波器建立了一个离散的多分辨率多方向率分析，实现灵活的多分辨率、局部的、具有方向性的图像表示。研究Contourlet变换的方向选择能力和非线性近似能力，体现Contourlet变换超越小波的优异表现。但是由于小波变换与Contourlet变换都缺乏平移不变性，图像增强结果会产生伪Gibbs失真，CUNHA A L等提出的非抽样Contourlet变换具有平移不变性，可以在一定程度上抑制这种失真。本文对其进行了改进，与基于小波变换、Contourlet变换的图像增强算法相比，该算取得了良好的增强效果[2]。

1 非抽样Contourlet变换

Contourlet变换也称金字塔形方向滤波器组PDFB(Pyramidal Direction Filter Bank)，其分解变换的实现可以分为两个步骤：拉普拉斯金字塔LP(Laplacian Pyramid)分解和方向滤波器组DFB(Directional Filter Bank)滤波。其变换思想是使用类似于线段的基函数去逼近原始图像，从而实现对图像信号的稀疏分离。因此，要实现Contourlet变换，首先需要对图像进行一个多尺度变换以检测不同尺度下的奇异点，然后再通过一个具有局部性的方向变换将同一尺度下相邻的奇异点连接成线段结构[3]。Contourlet变换的实现过程可以归纳为如下步骤：

(1)使图像通过类似于小波的多尺度变换以检测边缘上的奇异点；
(2)将步骤(1)所得到的图像通过局部化的方向变换完成轮廓线段的检测。

非抽样Contourlet变换NSCT(Nonsubsampled Contourlet Transform）由非抽样塔状滤波器NSP(Nonsubsampled Pyramid)将图像分解为低频部分和高频部分，然后由非抽样方向性滤波器组NSDFB(Nonsubasmpled Directional Filter Banks)将高频部分分解为若干个方向。

NSP是一种平移不变性的双通道滤波器结构，它使NSCT具有多尺度性质，并且下一层次的滤波器可以通过对上一层的滤波器抽样得到。其频域分解图如图1所示。

NSDFB由两通道的非采样滤波器组迭代构成，该滤波器组也没有进行采样，具有平移不变性。NSDFB可以将第一级变换所得到的高频部分分解为2的任意次幂个方向，每个方向上的高频部分与NSP得到的低频部分以及原始图像都有相同的大小。NSDFB分解是将信号在一组基函数上展开,它对应的基函数之间是冗余的[4]。频域分解图如图2所示。

NSP与NDFB可以保证信号完全重建的条件是滤波器必须满足等式：

其中，H0(z)、H1(z)表示分解滤波器，G0(z)、G1(z)表示重建滤波器。非抽样Contourlet变换在表达图像时具有Contourlet变换所具有的优点，还具有平移不变性[4]。

2 基于的非抽样的Contourlet变换图像增强算法

图像变换后，对变换系数分三种：强边缘、弱边缘和噪声。强边缘每个方向的系数值都较大；弱边缘在某一方向的系数大,但在其他方向上的系数小；噪声是指那些在所有方向上的系数都较小[5]。LAINE A F提出的增强函数：

x为输入原始图像的变换系数，0p1再次放大。此函数可以放大弱边缘的系数，保持强边缘的系数[6]。

基于Contourlet变换的图像增强算法的流程如图3所示，具体步骤如下:

(1)对图像进行非抽样Coutourlet变换，得到不同尺度不同方向的变换系数;
　 (2)按照上述原则对Contourlet变换系数进行处理;
(3)由修正后的变换系数重建增强图像。

3 实验结果与分析

本文采用客观评价标准信噪比SNR来衡量不同去噪方法去噪后图像的客观质量，SNR的定义如下：
　
实验结果如图4所示。原图是分辨率为640×480像素的红外热像仪拍摄的红外照片，建筑物为主要探测目标，建筑物后有少量植被。图4(b)为加噪图像，图4(c)为利用拉普拉斯变换处理的图像，图4(d)为利用Contourlet变换处理的图像，图4(e)为利用本文算法处理的图像。SNR值如表1所示。