一种基于结构-纹理分解的遥感图像变分融合方法技术

本发明专利技术公开了一种基于结构‑纹理分解的多光谱与全色图像变分融合方法，基于遥感领域中PAN图像和HRMS图像之间结构分量和纹理分量的相似性，分别对全局信息和局部空间细节(高阶信息)进行保留。通过梯度稀疏度来确保PAN图像和HRMS图像的结构相似度，根据低轶属性来描述其纹理相似度，提出我们的变分模型。然后通过交替方向乘子法求解该模型，从而得到具有高分辨率的多光谱图像。实验表明：本发明专利技术提出的方法考虑到PAN图像和HRMS图像之间结构分量和纹理分量的相似性，使用不同的模拟数据集和真实数据集进行融合实验，与对比方法相比，我们的融合图像的各项指标均最接近于理想值，视觉上最接近于参考图像。

A variational fusion method of remote sensing image based on structure texture decomposition

【技术实现步骤摘要】
一种基于结构-纹理分解的遥感图像变分融合方法
本专利技术属于图像融合领域，涉及一种基于结构-纹理(cartoon-texture)分解的多光谱与全色图像变分融合方法，适用于由低分辨率多光谱图像与高分辨率全色图像融合得到高分辨率多光谱图像的应用场景。
技术介绍
在遥感领域，通过卫星获取的地球图像的数据具有重要的应用价值，例如应用于土地覆盖分类，灾害检测和植被绘图等。由于卫星的存储与传输负担等限制，通常使用不同的传感器分别获得同一画面的光谱信息丰富但空间分辨率低的多光谱图像(LRMS)和光谱信息单一但空间分辨率高的全色图像(PAN)，并通过对其进行融合来获得同时具有丰富的光谱信息和高分辨率的多光谱图像(HRMS)，以便有利于对信息的分析和理解，该图像融合技术也被称为全锐化技术。基于不同的融合原理，全锐化方法通常被分为4类，包括成分替换方法(componentsubstitution，CS)，多分辨率分析法(multi-resolutionanalysis，MRA)，基于深度卷积神经网络(deepconvolutionalneuronnetwork-based，DCNN)的全锐化方法，以及变分法。CS方法将MS图像映射到一个新的空间，用PAN图像替换其中的空间分量后逆变换到原空间。MRA方法是通过多分辨率分解提取出PAN图像的细节然后将其注入MS波段中，这两类全锐化方法均可以得到分辨率较高的多光谱图像，但通常遭受严重的光谱失真或局部伪像。随着深度学习理论的不断发展，研究人员们开始探索基于深度学习的全锐化方法(deeplearning-basedpansharpeningmethod，DLP)，并展现出优异的性能。但这类方法需要大量的计算资源和训练数据，因此在资源受限的情况下缺乏灵活性。近年来，变分方法依靠其优异的融合效果成为了一个广受欢迎的融合方法，这类方法的核心为基于某些假设和先验信息建立目标能量函数。例如，C.Chen等提出引入L2,1约束来鼓励动态梯度稀疏性和群体稀疏性的SIRF方法，较好的提升了融合质量。XFu等提出一种基于局部梯度约束(LocalGradientConstraint，LGC)的全锐化变分模型，该方法考虑PAN和HRMS图像在不同局部块和光谱波段中的梯度差异来进行局部梯度约束，实现了高质量的锐化。在已提出的众多变分方法中，考虑到PAN图像和MS图像取自同一场景，它们具有共同的几何特征，因此通常利用梯度相似性来对取自同一场景的PAN图像和MS图像的几何结构一致性进行约束。然而梯度算子只能获得图像的一阶边缘信息，因此融合过程中全色图像的高阶空间信息并未得到有效的利用。
技术实现思路
本专利技术针对现有变分方法的不足，提供一种新的将卡通-纹理分解与二阶导变分模型结合的变分全锐化方法。本专利技术所采用的技术方案是：一种结构-纹理分解与二阶导变分模型结合的变分全锐化方法。图1展示了一张全色图像的梯度图像以及结构-纹理分量，其中结构分量包含了全局结构信息，而纹理分量中包含着一些梯度图像中不存在的图案细节(即高阶信息)，这对于全锐化过程中的细节信息提取更有利。基于PAN图像和HRMS图像之间结构分量和纹理分量的相似性，我们分别使用梯度稀疏性和低轶性进行先验约束来保留空间信息，通过数据保真项保留MS图像的光谱信息，通过交替方向乘子法(alternativedirectionmultipliermethod，ADMM)制定并求解变分模型，从而获得具有高分辨率的多光谱图像。步骤1，选取同一时间同一画面的经过配准处理的多光谱图像和全色图像，低分辨率多光谱图像的记为M，全色图像记为P，融合得到的高分辨率多光谱图像记为F，全色图像与多光谱图像之间的分辨率之比为c；步骤2，构建数据保真项J(F)来保留低分辨率多光谱图像M的光谱信息，构建先验约束项K(F)以确保融合图像F与全色图像P之间的空间结构一致性，从而获得目标能量函数：argminFJ(F)+K(F)(1)步骤3：基于步骤2的目标能量函数，通过交替方向乘子法来求解优化问题，从而得到融合图像。进一步的，步骤2的具体实现包括以下子步骤，步骤2.1，构建数据保真项J(F)，来保留低分辨率多光谱图像M的光谱信息，低分辨率多光谱图像M被视为高分辨率多光谱图像F的下采样模糊版本，即M＝ΨF+n，其中Ψ表示下采样运算符，n表示均值为0的高斯噪声；因此，为了保留LRMS图像的光谱信息，构造保真项为：其中表示Frobenius范数；步骤2.2，基于全色图像和相应的高分辨率多光谱图像之间结构分量和纹理分量的相似性，在结构-纹理空间构建先验约束项K(F)；具体原理如下：融合图像F表示为纹理分量R和结构分量S的组合：F＝R+S(3)其中结构分量S代表的是图像中的全局结构信息，即大面积平滑区域和锐利的边缘，纹理分量R代表的是图像中的小尺寸细节内容；选择使用核范数来对纹理分量进行约束，因此，先验约束项K(F)：其中，和表示将U和V分别复制到b个波段，U和V为全色图像P的结构分量和纹理分量，由图像分解算法获得；R和S分别为F的结构和纹理分量，为梯度算子，||·||2,1表示群体稀疏；||·||*表示核范数，ρ和τ为平衡各项的参数；为块映射函数，用于将矩阵无重叠的分割为r×r的块并将每个块重新排列成长度为r2的列向量；步骤2.3：基于子步骤2.1和2.2，提出的目标能量函数模型为：进一步的，步骤3的具体实现包括以下子步骤，步骤3.1：对于步骤2得到的目标能量函数模型，通过交替方向乘子法来求解该优化问题，通过引入辅助变量Y，(5)的增广拉格朗日函数表示为：其中λ为拉格朗日参数，(6)式的优化问题被分解成4个子问题，分别进行求解；步骤3.2：求解F子问题：其中下标t表示迭代次数，上述问题相当于：其中I是单位矩阵，等式(8)通过高斯消去法求解；步骤3.3：求解R子问题：等式(9)通过快速迭代收缩阈值算法来解决；步骤3.4：求解S子问题：等式(10)是低轶优化问题，通过奇异值阈值算法解决；步骤3.5：求解Y子问题，该子问题由下式直接更新：Yt+1＝Yt+λ(Ft-Rt-St)(11)。本专利技术主要针对由低分辨率多光谱图像和高分辨率全色图像融合得到高分辨率多光谱图像的应用需求。考虑到在图像全锐化过程中，梯度算子只能获得图像的一阶边缘信息，导致融合过程中全色图像的高阶空间信息并未得到有效的利用。而图像的结构-纹理分解的理论中，结构分量包含了全局结构信息，而纹理分量中包含着一些梯度图像中不存在的图案细节(即高阶信息)，这对于全锐化过程中的细节信息提取更有利。因此，基于PAN图像和HRMS图像之间结构分量和纹理分量的相似性，我们提出一种卡通-纹理分解与二阶导变分模型结合的变分全锐化方法。我们的模型分别使用梯度稀疏性和低轶性进行先验约束来保留空间信本文档来自技高网...

【技术保护点】
1.一种基于结构-纹理分解的遥感图像变分融合方法，其特征在于，包括以下步骤：/n步骤1，选取同一时间同一画面的经过配准处理的多光谱图像和全色图像，低分辨率多光谱图像的记为M，全色图像记为P，融合得到的高分辨率多光谱图像记为F，全色图像与多光谱图像之间的分辨率之比为c；/n步骤2，构建数据保真项J(F)来保留低分辨率多光谱图像M的光谱信息，构建先验约束项K(F)以确保融合图像F与全色图像P之间的空间结构一致性，从而获得目标能量函数：/narg min

【技术特征摘要】
1.一种基于结构-纹理分解的遥感图像变分融合方法，其特征在于，包括以下步骤：
步骤1，选取同一时间同一画面的经过配准处理的多光谱图像和全色图像，低分辨率多光谱图像的记为M，全色图像记为P，融合得到的高分辨率多光谱图像记为F，全色图像与多光谱图像之间的分辨率之比为c；
步骤2，构建数据保真项J(F)来保留低分辨率多光谱图像M的光谱信息，构建先验约束项K(F)以确保融合图像F与全色图像P之间的空间结构一致性，从而获得目标能量函数：
argminFJ(F)+K(F)(1)
步骤3：基于步骤2的目标能量函数，通过交替方向乘子法来求解优化问题，从而得到融合图像。


2.如权利要求1所述的一种基于结构-纹理分解的遥感图像变分融合方法，其特征在于：步骤2的具体实现包括以下子步骤，
步骤2.1，构建数据保真项J(F)，来保留低分辨率多光谱图像M的光谱信息，低分辨率多光谱图像M被视为高分辨率多光谱图像F的下采样模糊版本，即M＝ΨF+n，其中Ψ表示下采样运算符，n表示均值为0的高斯噪声；因此，为了保留LRMS图像的光谱信息，构造保真项为：



其中表示Frobenius范数；
步骤2.2，基于全色图像和相应的高分辨率多光谱图像之间结构分量和纹理分量的相似性，在结构-纹理空间构建先验约束项K(F)；具体原理如下：
融合图像F表示为纹理分量R和结构分量S的组合：
F＝R+S(3)
其中结构分量S代表的是图像中的全局结构信息，即大面积平滑区域和锐利的边缘，纹理分量R代表的是图像中的小尺寸细节内容；
选择...

【专利技术属性】
技术研发人员：田昕，陈月荣，杨晋陵，李松，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42