[发明专利]一种全色多光谱影像融合方法

权利要求书

1.一种全色多光谱影像融合方法，其特征在于包括以下步骤：

(1)通过以下公式对全色影像Pan降质，得到模拟全色影像Pan^*：

其中，m和n分别表示全色影像Pan的长度和宽度，表示经过均值滤波后位于模拟全色影像Pan^*的(x,y)处的灰度值，f(i,j)表示全色影像Pan的(i,j)处的灰度值；

(2)将原始多光谱影像MS重采样，得到多光谱影像MS′，使得多光谱影像MS′在长和宽具有和全色影像Pan相同数量的采样点；

(3)确定模拟全色影像Pan^*中各像素的基础融合系数；

(4)根据多光谱影像判断模拟全色影像Pan^*中各像素的地物类型；

(5)根据地物类型通过以下公式对基础融合系数进行区别性拉伸，得到最终融合系数；

λ(i,j)＝ωλ_basic(i,j)+1-ω

其中，λ(i,j)表示融合影像上(i,j)处的最终融合系数，ω为根据地物类型设置的权重值，若地物类型为建筑用地则ω设置为ω_building，若地物类型为植被则ω设置为ω_vegetation，若地物类型为水体则ω设置为ω_water，ω_building＞ω_vegetation＞ω_water；

(6)根据最终融合系数通过以下公式进行融合：

其中，表示融合影像在(i,j)处的灰度值，MS'_band(i,j)表示重采样得到的多光谱影像MS′在(i,j)处的灰度值。

2.根据权利要求1所述的全色多光谱影像融合方法，其特征在于：步骤(2)中，使用双线性内插法对原始多光谱影像MS进行重采样，以得到多光谱影像MS′：

首先，通过以下公式计算重采样后位于多光谱影像MS′中(dstY,dstX)处的像素在原始多光谱影像MS中的位置(srcY,srcX)：

srcX＝dstX*(srcWidth/dstWidth)

srcY＝dstY*(srcHeight/dstHeight)

其中srcHeight和srcWidth分别表示原始多光谱影像MS的长与宽，dstHeight和dstWidth分别表示重采样后的多光谱影像MS′的长与宽；

然后，将srcY分为整数部分i₂与小数部分u，将srcX分为整数部分j₂与小数部分v，通过以下公式对重采样进行描述：

其中，表示重采样得到的多光谱影像MS′在(dstY,dstX)处的像素灰度值，f′(srcY,srcX)表示原始多光谱影像MS中(dstY,dstX)的灰度值。

3.根据权利要求1所述的全色多光谱影像融合方法，其特征在于：步骤(3)中，通过以下公式确定基础融合系数：

其中，μ₁和μ₂为设置的系数，0<μ₁<3，1<μ₂<3，λ_basic(i,j)表示融合影像在(i,j)处的基础融合系数，pan(i,j)表示全色影像Pan在(i,j)处的灰度值，pan^*(i,j)表示模拟全色影像Pan^*在(i,j)处的灰度值。

4.根据权利要求3所述的全色多光谱影像融合方法，其特征在于：μ₁和μ₂均设置为2。