基于小波包的ETM遥感图像融合方法研究.pdfVIP

基于小波包的ETM+遥感图像融合方法研究 张登荣俞乐刘辅兵 (浙江大学理学院空间信息技术研究所，杭卅『310027) 摘要：本文在总结了像索级图像融台方法的基础上，介绍了小波包变换的方法及其在遥感 图像融合中的应用，讨论了基于小波包变换的遥感图像融合中的小波基选择。采用熵、联合 熵、偏差指数和边缘指数等作为定量指标，对融合结果进行了质量评价。通过ETM+实验表 明基于小波包的融合方法在光谱信息保持和纹理信息增强上相较于其他小波基方法有一定 优势。 关键词：图像融合z小波包；小波基； 融合质量评价 0引言 随着遥感技术的快速发展，遥感数据资源日益丰富，充分利用不同传感器获得的不同波 段，不同空间分辨率，不同时间分辨率，不同光谱分辨率的多源遥感信息关系到遥感应用的 深入发展。与单源遥感影像相比，多源遥感图像融合可以减少系统的多义性和不确定性，增 加特征提取的精确性，能形成对被感知对象的完整一致的信息描述“J。当前，常用的像素级 J， 遥感影像融合的方法有：IHS变换融合方法”“，HPF融合方法…，PCA变换融合方法”．7 小波变换融合方法[8-11]以及各种方法问的相互配合。H-IS方法是最常用的一种图像融合方法， 能大大提高图像的可判读性，但需要同时对3个波段进行融合，由于不同波段数据的光谱曲 线不同，融合后的图像扭曲了原始光谱特性，会产生光谱退化现象9】。胂方法是一种光谱 信息丢失较少的方法，但融合效果对滤波器的依赖很大，很难较好地保持高分辨率遥感影像 的纹理…。PCA方法适用于多波段融合，但单纯的用I：)121来代替多波段会导致原图像物理 特征的丢失，产生光谱退化16J。基于小波变换的融合方法能最大限度地保留光谱特性，能满 足遥感应用的要求，同时可以对单个和多个波段进行融合操作【9】。但是小波变换方法本身只 对低频信号进行剖分，没有对高频信号进一步分解。很多信息从而得不到充分利用。小波包 变换[121则能同时对低频和高频部分进行分解，具有更好的时频特性，在图像融合领域已经 有研究证明了基于小波包变换的融合方法的有效性[13-151．本文以ETM+多光谱影像和全色影 像为例，通过基于小波包的融合方法与传统方法之问的比较实验表明基于小波包融合方法既 很好的保留了光谱信息又提高了空问分辨率。 1小波包变换 小波包概念最早足由R．R．Coifman和M．VWickerhause等人在小波变换的基础上提出 来的。小波包的定义为： 若W。(f)满足下式： I％。(f)=√2∑h(k)W．(2t-k) 1 w2。(f)=压∑g(I)％(2卜t) L I 确定的小波包。若≯(f)为正交尺度函数，妒(f)为正交小波函数，则称为正交小波包。 ·530· 为了便于说明，用，0)表示二维空问，C表示低频子空间，D表示高频子空问。下标 代表分解层数，上标h，v，d分别代表高频的水平，垂直，对角分量。相对于图1的小波 分析，小波包可以表示为图2。 阳塞由 图1小波变换示意图 广]Tr] 『mnl 鲁皇E蚓 。一一孑j耀≮、、～。 Et才：—吲 累、 景 只J ≯J j竺竺】竺【竺匡蕉至歪童删，碰g…m一～咖叫d 图2小波包变换示意图 小渡包变抉算法，即小波包的分解和重构算法，可以表示为： 设f(t)在子空问U雌中的系数(c‰}。：为 c肚n=￡邝)2叫2瓦(21一k)dt cn，j进一步分解为 fc鬻=∑h(1-2k)c；： k2n址=∑g(t-2k)c7．，