遥感影像像素级融合方法概述

  • 格式:docx
  • 大小:28.26 KB
  • 文档页数:3

遥感影像像素级融合方法概述

摘要:本文介绍了6种不同的像素级遥感影像融合方法,并较为深入的探讨了几种融合方法的原理、实现过程和优缺点,并对几种融合方法进行了综合评价,得到了针对多光谱影像融合较好的方法,表明了不同传感器影像融合的潜力,同时论文也对融合影像的应用推广进行了初步的对比分析。

关键词:融合 原理 评价 方法

0 引言

影像融合按照其水平和特点可以分为像素级、特征级和决策级融合[1],本文讨论的是像素级的影像融合方法,像素级融合是对传感器的原始信息及预处理的各个阶段上产生的信息分别进行融合处理。对不同卫星和不同传感器的遥感影像进行融合,可以取长补短,产生新的融合影像。比较常见的影像融合是高分辨率全色影像和低分辨率多光谱影像的融合,生成的新的融合影像既具有全色影像的高分辨率,也一定程度上保持了多光谱影像的光谱信息。

目前比较成功的高分辨率商业卫星如QuickBird、IKONOS和SPOT5,都搭载有高分辨率全色和低分辨率的多光谱相机,能同时提供全色和多光谱影像,并提供两者的融合产品。QuickBird、IKONOS和SPOT5的全色和多光谱影像分辨率之比都为1:4或1:2,而CBERS-02B星的全色和多光谱影像分辨率之比大约为1:8,虽然CBERS-02B星HR与CCD的融合影像能一定程度上保持HR的高空间分辨率和CCD的高光谱分辨率,但1:8的大比率给CBERS-02B星的全色和多光谱影像融合也带来了一些困难和问题,如全色和多光谱影像之间的配准精度难已保证、融合影像出现斑块效应等等。

1 像素级影像融合方法

1.1 HSV融合方法 HSV融合方法属于一种颜色变换的融合方法,首先介绍一下HSV颜色变换,HSV颜色变换是把标准的RGB图像变换到为色度H (Hue)、饱和度S (Saturation)和亮度V (Value)图像。HSV融合方法流程是对多光谱影像3个波段使用HSV颜色正变换为H、S和V三幅图像,然后用高分辨率影像替代H图像,最后对H、S和V图像实施HSV颜色变换的逆变换得到融合影像。

1.2 GS融合方法 GS (Gram-Schmidt) 融合方法是商用遥感图像处理软件ENVI中独有的影像融合方法,具体算法不清楚。

1.3 HPF融合方法 HPF(High Pass Filter) 融合方法基本思想是用一个高通滤波器对高分辨影像滤波得到高频信息,然后按照一定的权重注入到多光谱波段中生成融合影像。

HPF融合方法的基本流程如下:①计算多光谱影像和高分辨率影像分辨率之比,以此来确定高通滤波器的卷积滤波核的大小;② 对高分辨率影像进行高通滤波得到高频影像,同时对多光谱影像重采样到高分辨影像像元大小;③把高频信息加入到每个多光谱影像中得到融合影像。

1.4 Brovey融合方法 Brovey变换融合方法是一种通过归一化后的三个波段多光谱影像与高分辨率影像乘积的融合方法,它可以应用于不同传感器的影像数据融合,如TM多光谱和SPOT全色影像的融合。Brovey变换既一定程度上保持多光谱数据的完整性又增强了图像的细节。

Brovey变换融合CBERS-02B CCD多光谱数据和HR全色数据的计算

1.5 SFIM融合方法 Liu[2](2000) 提出了基于光滑滤波的亮度调节 (Smoothing

Filter-based Intensity Modulation, SFIM) 影像融合方法,SFIM融合方法的主要流程是,首先对高分辨率影像进行平滑滤波 (低通滤波) 生成均值平滑影像,然后再进行代数运算

滤波核的大小一般定为低分影像对高分影像空间分辨率比的大小。

1.6 小波变换融合方法 Mallat和Atrous算法是二维离散小波变换影像融合中最常用的两种算法。Mallat算法是由S.Mallat在1988年提出的,它是正交小波的构造方法和正交小波变换的快速算法。基于小波变换的融合方法有两种实现途径:叠加方法(Additive Method)或取代方法(Substitution Method)。除了简单的叠加方法和取代方法之外还有结合IHS和PCA变换的叠加和取代方法,考虑到Nunez[3]表明用Mallat算法进行小波分解,通过叠加和取代方法得到的融合影像的质量基本一致,本文中小波变换融合采用基于Mallat算法的简单叠加方法而不考虑取代方法。

2 融合方法对比

2.1 定性评价 对于融合方法的定性评价一般是目视判读融合结果图像,主要对比空间细节和光谱保持。目视判读的定性评价缺点在于有较大的主观性,不同的人会有不同的评价结果,因此一般还要进行融合结果的定量评价。

2.2 定量评价 定量评价一般采用3个指标:对比度、高频信息的相关系数已经融合结果的相关系数。

对比度和图像熵这两个指标都是基于灰度共生矩阵计算的,图像的灰度共生矩阵表达了图像灰度关于方向、相邻间隔、变化幅度的综合信息。在灰度共生矩阵的基础上抽取出图像的纹理特征参数,称为二次统计量,对于这二次统计量的分析,可以很好地理解图像的纹理特性。

经过正规化处理的灰度共生矩阵是抽取二次统计量纹理特征系数的基础。

图像的对比度是用于评价图像纹理的参数,可理解为图像的清晰程度,即纹理的清晰程度。图像中,纹理的沟纹越深,其对比度越大,图像就越清晰,视觉效果越好。相邻像素间的灰度值差异在同一幅图像中相同的灰度对的多少最能影响结果。熵是图像所具有的信息量的度量,纹理的复杂度越高就意味着图像信息量越大,其熵也越大。

高频信息的相关系数由Zhou[4](1998) 提出,考虑到影像的空间细节信息主要保留在影像的高频部分,对融合影像和高分辨率全色影像高通滤波,然后计算滤波后两影像之间的相关系数,相关系数理想值为1,相关系数值越高表示融合影像的空间细节越接近于全色影像。

遥感影像融合后,融合结果的光谱特性应与原始多光谱数据尽可能保持一致,这样才能保证融合前基于光谱特性可分的目标在融合后仍然可分。

其中fi,j和gi,j分别是像元(i,j)在影像f和g中的灰度值,ef和eg分别是影像f和g影像的灰度值均值。融合影像和原始影像相关系数的理想值为1,值越大表明光谱质量越好。

3 结论

综合分析,就空间细节保持方面,较好的方法依次是小波、SFIM和HPF,而同时考虑空间细节保持和光谱保持,较好的方法依次是HPF、小波和SFIM。同时考虑空间细节和光谱保持,较好的方法依次是HPF、GS和Brovey。

文中提到的6种方法得到的融合影像效果均有其优缺点,融合的结果既有高分辨率影像的空间信息,也有多光谱影像的光谱信息,影像相对全色和多光谱有较好的目视效果,但是相对与HR与XS的融合影像还是有些不足,融合影像的颜色不如HR与XS融合影像自然丰富,某些地方斑块效应比较明显,这主要是由于1:8的分辨率之差所导致的。

参考文献:

[1]POHL C, VAN GENDEREN J L. Multisensor Image Fusion in Remote Sensing: Concept,

Methods and Applications [J]. International Journal of Remote Sensing, 1998, 19(5):

823-854.

[2]Liu, J.G., 2000, Smoothing Filter Based Intensity Modulation: A Spectral Preserving

Imaging Fusion Technique for Improving Spatial Details. International Journal of

Remote Sensing, 21(18): 3461-3472.

[3]NUNEZ J, OTAZU X, FORS O, et al. Multiresolution-based image fusion with additive

wavelet decomposition. IEEE Transactions on Geoscience and Remote Sensing, 1999,

37(3): 1204-1211.

[4]ZHOU J, CIVCO D L, SILANDER J A. A Wavelet Transform Method to Merge Landsat

TM and SPOT Panchromatic Data [J]. International Journal of Remote Sensing, 1998,

19(4): 743-757.