[发明专利]一种基于离散余弦变换的全色与多光谱图像融合方法

说明书

本发明公开了一种基于离散余弦变换的全色与多光谱图像融合方法，可解决现有全色与多光谱图像融合方法的缺点与不足。该基于离散余弦变换的全色与多光谱图像融合方法通过引入图像压缩标准中广泛采用的DCT变换来提取全色图像中的高频信息，并通过DCT反变换将其成功注入低分辨率多光谱图像中。与现有的技术相比，本发明在保证空间细节失真和谱失真均较小的前提下，融合速度明显提高，从而为全色与多光谱图像融合提供了更有力的技术支持。

技术领域

本发明涉及遥感技术领域，具体的说是涉及一种基于离散余弦变换的全色与多光谱图像融合方法。

背景技术

高分辨率遥感卫星如QuickBird、IKONOS，GeoEye等能够同时提供全色(Panchromatic, Pan) 波段图像和多光谱(Multispectral, MS)波段图像.全色图像具有高空间分辨率，反映了目标地物的空间结构信息，能够详尽地表达地物的细节特征，但光谱分辨率不足；多光谱图像光谱信息丰富，利于对地物的识别与解译，但空间分辨率低． 然而, 在遥感应用如目视解译、变化检测等任务中常常需要同时具有高空间和高光谱分辨率的MS图像. 全色与多光谱图像融合技术（Pan-sharpening）就是用来将全色图像的空间细节信息注入多光谱波段图像以提高其空间分辨率, 同时尽可能地保持多光谱波段图像的光谱特性的一种遥感图像处理技术。

传统的全色与多光谱图像融合技术可分为两大类：分量替换法（Componentsubstitution, CS），多分辨率分析法(Multiresolution-analysis, MRA). 分量替换法的主要思想是将MS图像按点做谱段变换，获得一个类似与全色图像的主分量，然后用全色图像替换该主分量并进行按点反变换，所得结果即为既具有高空间分辨率又具有高光谱分辨率的MS图像。根据所采用的谱段变换方式不同，典型方法有基于HIS、基于PCA、基于GA等等的方法。CS 方法具有计算效率高, 融合图像空间失真小的优点，但缺点是融合图像的光谱失真较大。多分辨率分析法主要利用小波变换的多分辨特性，通过多级小波变换提取全色图像的空间细节信息，然后将其注入多光谱图像，从而获得既具有高空间分辨率又具有高光谱分辨率的MS图像。多分辨率分析法具有融合图像的光谱失真较小的优点，但缺点是计算复杂度较高，且融合图像空间失真较大。

发明内容

本发明的目的在于解决现有技术所存在的问题，提供了一种基于离散余弦变换的快速有效的全色与多光谱图像融合方法。

为了实现上述目的，本发明采用以下技术方案：

一种基于离散余弦变换的全色与多光谱图像融合方法，包括以下步骤：

（1） 输入高分辨全色图像P和低分辨多光谱图像MS_k

（2）将高分辨全色图像P分为N×N大小的不重叠块

（3）对每个N×N块做离散余弦变换（DCT， Discrete Cosine Transform）

（4）用低分辨多光谱图像的像点值代替对应块的DC分量

（5）对步骤（4）的结果做离散余弦反变换（ICDT，）所得结果即为高分辨多光谱图像。

进一步的，所述步骤（1）中对高分辨全色图像和低分辨多光谱图像的尺寸比例无2的幂次倍要求限制，可以是任意倍。

进一步的，所述步骤（2）中N的取值为全色图像分辨率和高光谱图像分辨之比。即全色图像的总块数与多光谱图像的总像素点数相同，每个N×N大小块对应低分辨率多光谱图像的一个像素点。

进一步的，所述步骤（4）中低分辨多光谱图像的像点值代替对应块的DC分量的实现公式为：