[发明专利]一种普适的高分辨率遥感图像融合方法

权利要求书

1.一种普适的高分辨率遥感图像融合方法，其步骤为：

第一步、由全色影像波谱范围与多光谱影像波谱范围之间的关系和公式1-3计算出基于辐射亮度radiance的、由多光谱波段影像模拟的全色影像Pan_HSyn；

PanHSyn=XSH1+XSH2+...+XSHn+(XSH1λ1,2-λ1,1+XSH2λ2,2-λ2,1)×λ2,1-λ1,22---(1)]]>

+(XSH2λ2,2-λ2,1+XSH3λ3,2-λ3,1)×λ3,1-λ2,22+...+(XSH(n-1)λ(n-1),2-λ(n-1),1+XSHnλn,2-λn,1)×λn,1-λ(n-1),22]]>

也可以写成：

其中参数由下式获得：

其中，XS_Hi表示多光谱影像第i个波段的辐射亮度；λ_i，1和λ_i，2是多光谱影像的第i个波段的两端波长数值，n是多光谱总的波段数目；

第二步、用波谱范围比值来修正Pan_HSyn，以弥补全色与多光谱数据波谱范围的差异，即：

设：

p=dλpandλmul---(4)]]>

其中，p为影像波谱范围比，为全色波段的波谱范围，d_λmul为多光谱波段的波谱范围；

当全色影像的波谱范围和多光谱影像的波谱范围相一致的时候，p＝1；

令

PanHSyn*=PanHSyn*p---(5)]]>

第三步、在原始全色波段影像和之间寻找一一对应的像素点，这些点分别代表图像上的几种不同类型地物，例如：植被、水体、建筑物、裸地等；每一种地物类别分别取5-10个点作为样本点；

第四步、计算原始全色波段影像和样本点之间的线性回归系数k、b；

第五步、由公式

PanHSyn**=kPanHSyn*+b---(6)]]>

计算得到新的模拟全色的值

第六步、将代入到式7中

XSPi=PanHXSHiPanHSyn**---(7)]]>

其中XSP_i为融合后的高分辨率图像第i个波段的灰度值；Pan_H是全色图像，XS_Hi是将多光谱图像重采样为全色图像相同像元大小后，第i个波段的灰度值。