[发明专利]基于子字典稀疏重构的全色与多光谱遥感图像融合方法

权利要求书

1.一种基于子字典稀疏重构的全色与多光谱遥感图像融合方法，其特征在于，包括以下步骤：

记YPAN和YMS分别表示已知的全色图像和多光谱图像，XMS表示未知高空间分辨率的多光谱图像，其中YMS和XMS的第b个波段的图像分别记为YMS_b和XMS_b，b＝1,2,...,B，B表示光谱波段总数目，记全色图像与多光谱图像的空间分辨率之比为γ:1；

步骤一，高空间分辨率多光谱图像字典的构造

步骤1.1，将B个波段的YMS_b分别上采样成与YPAN大小相同的图像，记为MS_b；

步骤1.2，从YPAN以及其对应位置处的MS_b中随机选M个大小为图像块，并将图像块拉直成向量依次排列，组成第a个波段的子字典，a∈b；即：

D_a＝[y_MS_b,1,y_MS_b,2,...,y_MS_b,M,yPAN₁,yPAN₂,...,yPAN_M] 式1

式1中y_MS_b,j与yPAN_j分别表示MS_b与YPAN图像中的第j个图像块拉直成的向量，j＝1,2,...,M；

步骤1.3，采用步骤1.2的方法依次构造B个波段的子字典D_b；

步骤1.4，子字典的级联构成了高空间分辨率多光谱图像字典D，即D＝[D₁,D₂,...,D_b,...D_B]^T；

步骤二，基于稀疏重构的图像融合

步骤2.1，对YMS_b、YPAN分别以大小为的滑动窗口依次按照从左上到右下的顺序提取图像块，提取图像块时的步长均为1，再将图像块拉直并依次排列组成矩阵yMS_b与yPAN；

步骤2.2，令其中和为单位矩阵，1为长度为γ的全1向量；令M₂＝(w₁I,w₂I,...,w_bI,...,w_BI)，其中为单位矩阵，w_b表示每个波段对应线性组合的权重，并且满足

步骤2.3，建立模型，令：

y＝M×xMS+v 式2

其中，xMS表示XMS_b以的滑动窗口依次按照从左上到右下的顺序提取图像块，并将图像块拉直依次排列组成的矩阵；v1与v2分别表示未知高空间分辨率多光谱图像退化成多光谱图像与全色图像的噪声；

步骤2.4，采用BP算法求解下面方程

其中，Φ＝MD，D表示步骤1.4中的字典，α为稀疏表示系数；ε表示容许误差，||α||₀表示α的l₀范数，||y-Φα||₂表示y-Φα的l₂范数；

步骤2.5，将xMS进行重构：

xMS＝D·α 式4

步骤2.6，将xMS_b矩阵中的每一列排列成大小的块，再将这些块按照从左上到右下的顺序放到XMS_b的对应位置处，并取平均，从而得到融合图像XMS_b。