基于HIS变换和对比度塔形分解图像融合方法.docVIP

基于HIS变换和对比度塔形分解图像融合方法 　　［摘 要］本文采用了HIS变换和对比度塔形分解相结合的方法对彩色图像和全色波段图像进行融合。该方法首先对彩色图像进行HIS变换，然后对I分量图像和全色波段图像分别进行对比度塔形分解，接着对于各分解层采用基于区域能量的融合规则及其相应的融合算子，最后通过对比度金字塔重建和HIS逆变换得到最终的融合图像。经实验对比表明，该方法具有较好的融合效果。?? 　　［关键词］图像融合 HIS变换 对比度塔形?? 　　［中图分类号］TP ［文献标识码］A ［文章编号］1009－5489(2008)05－0193－03?お? 　　 　　Edward Waltz和James Llinas对数据融合给出了如下的定义：将不同类型传感器获取的图像数据经预处理后，采用一定的算法将各幅图像中所包含的信息优势或互补性信息有机地结合起来，以产生新的数据，来获得对同一事物或目标的更客观、更本质的认识。从而大大提高融合图像的信息含量并使其在特征提取、分类、目标识别以及目视效果等方面更为有效。?? 　　目前图像融合的方法有很多，主要分为三类：像素级、特征级和决策级(或称符号级)的融合。而讨论得最多的是像素级图像融合。本文选用了具有高光谱分辨率而空间分辨率较低的多光谱影像与由全色波段的传感器获得的具有高空间分辨率但光谱分辨率较低的全色影像，进行了基于HIS变换和对比度塔形分解的图像融合。与其他常用的方法如HIS变换，小波变换相比，它不仅保留了图像的光谱信息，又很好地保留了空间细节信息，而且计算量也不是很大。最后，将融合后的与其他融合方法所得图像进行了对比。?? 　　 　　一、本文融合方法的流程图?? 　　 　　二、HIS变换?? 　　HIS变换也称彩色变换或蒙塞尔变换，在多源图像融合中，常常需要把RGB空间转换到色度H、亮度I及饱和度S所在的空间，是其中在IHS空间融合不同分辨率的数据，即基于IHS变换的多源图像融合，以便准确定量地描述颜色特征。先将多光谱影像的三个波段映射到RGB三个通道，再进行HIS变换，将全色影像与多光谱影像的I分量进行直方图匹配后替换I分量，最后用HIS逆变换获得融合图像。?? 　　假定RGB值已归一化到范围?В?0，1］?В?角度?Е泉?是点与HIS空间的红色轴之间的夹角。其正变换如下：?? 　　若??B≤G?В???H=θ?В蝗绻???B＞G?В???H=360°－θ?В华? 　　其中?Е?=??arccos??{1/2〔(R－G)+(R－B)〕〔(R－G)??2+(R－B)(G－B)〕????1/2??}?И? 　　??S=1－3(R+G－B)〔??min??(R，G，B)〕；I=13(R+G+B)；?И? 　　逆变换公式如下：?? 　　当??0≤H＜120°?В?在这个领域内：?? 　　??B=I(1－S)；R=I〔1+〕S??cos??H??cos??(60°－H)；G=3I－(R+B)；?И? 　　?У?120°≤H＜240°，H=H－120°；?И? 　　??R=I(1－S)；G=I〔1+S??cos??(H－120°)??cos??(180°－H)〕；B=3I－(R+G)；?И? 　　?У?240°≤H≤360°，H=H－240°；?И? 　　??G=I(1－S)；B=I〔1+S??cos??(H－240°)??cos??(300°－H)〕；R=3I－(G+B)。?И? 　　该方法的特点是高频信息丰富，提高了融合图像的空间细节表现能力。但光谱信息损失较大。本文方法不是简单的用全色波段影像图代替I分量图，而是分别对它们进行了塔形分解，并融合后产生新的I′分量。对原来的HIS变换是一种改进。?? 　　三、遥感图像的对比度塔形分解?? 　　图像的对比度塔形分解是一种图像的多尺度、多分辨率的金字塔形，建立步骤如下：?? 　　1.建立图像的高斯塔形分解?? 　　设原图像??G??0??为，以??G??0?ё魑?高斯金字塔的零层，高斯金字塔的第l层图像??G??l?д庋?构造：?? 　　??G??l=∑2m=－2∑2n=－2w(m,n)G????l－1??(2i+m,2j+n)，0＜l≤N,0≤i≤C??l,0≤j＜R??l?И? 　　式中：??N?А?―高斯金字塔顶层的层号；?お?C??l?А?―高斯金字塔第l层图像的列数；?? 　　??R??l?А?―高斯金字塔第l层图像的行数；?お?w(m,n)?А?―为5×5的窗口函数。?? 　　这样，由??G??0?В???G??l?В?…，??G??n?Ь凸钩闪送枷竦母咚菇鹱炙?，其中??G??n??为金字塔的顶层。?? 　　2.由高斯金字塔建立图像的对比度金字塔?? 　　先将??