专利名称：基于压缩感知的多光谱与全色图像超分辨融合方法
技术领域：
本发明属于图像处理技术领域，更进一步涉及遥感图像处理技术领域中的基于压缩感知的多光谱与全色图像超分辨融合方法。本发明可以应用在土地利用、城区识别和森林资源调查等遥感领域。本发明是将低空间分辨率的多光谱图像与高空间分辨率的全色图像进行融合，同时利用压缩感知技术进行超分辨重建，以得到一幅同时具有高空间分辨率和光谱信息的融合图像。
背景技术：
多光谱图像和全色图像融合是指将从多光谱探测器获得的同一场景的多谱图像的信息综合利用，利用它们在时空上的相关性及信息上的互补性，得到对景物更全面、清晰的描述，更有利于人眼的识别和机器的自动探测。多光谱图像具有丰富的光谱信息，但空间分辨率较低；全色图像具有丰富的空间细节信息，空间分辨率较高，但光谱信息较少。要获得同时具有高的光谱和空间分辨率的图像，可通过两个途径实现一是采用更新型的传感器，这样能获取更高质量的图像，但对新型传感器的研制将耗费大量时间，且成本高；二是进行多光谱图像和全色图像的融合，这样能在现有传感器技术的基础上获取更高质量的图像，且技术成本低。目前的多光谱图像融合技术主要有两类一类是IHS色彩变换法、主成分分析法(PCA)等；另一类是小波变换法、Laplacian金字塔分解等多尺度变换方法，但这些方法在多光谱图像和全色图像融合中存在着一个共同的问题均采用“注入”的方法或运用分辨率变换提取全色图像的细节信息进行融合，利用全色图像简单的替代多光谱图像的亮度分量或当抽取全色图像的高频分量与多光谱图像的亮度分量的高频分量不一致时，会产生严重的颜色失真和光谱扭曲；另一方面，融合后图像的空间分辨率受限于融合前全色图像的分辨率。湖南大学申请的专利“一种基于压缩传感理论的卫星遥感图像融合方法”(申请号 201010283310.0,申请公布号CN 101996396A)中，该专利申请的方法是把压缩感知技术运用到多光谱图像的融合中，先向量化低空间分辨率多光谱图像和高空间分辨率全色图像； 然后构造高空间分辨率多光谱图像块的稀疏表示字典；再用基追踪算法求解高空间分辨率多光谱图像在该字典下的稀疏表示；再将稀疏表示与预先设定的过完备字典相乘得到高空间分辨率多光谱图像块向量；最后将向量表示转化为图像块得到融合图像。该方法能够较好的克服颜色失真和光谱扭曲，但是仍然存在的不足是，融合后的多光谱图像的分辨率取决于融合前全色图像的分辨率，只能得到与融合前全色图像相同分辨率的图像。
发明内容
本发明的目的在于克服上述已有技术的缺点，提出一种基于压缩感知的多光谱图像与全色图像的超分辨融合方法，以提高空间分辨率，同时降低颜色失真和光谱扭曲。为实现上述目的，本发明在压缩感知的框架下实现全色图像与多光谱图像的融合，其技术方案是首先将全色图像与多光谱图像作为压缩采样后的观测值，将高分辨率多
4光谱图像与低分辨率多光谱图像构造的观测矩阵，与高分辨率多光谱图像和高分辨率全色图像构造的观测矩阵进行联合；然后利用KSVD (K-Singular Value Decomposition)方法对高分辨率多光谱图像随机抽块进行字典学习，得到训练字典；最后利用压缩感知基追踪算法得到高分辨率多光谱图像。本发明的具体步骤如下(1)分别输入低分辨率多光谱图像和高分辨率全色图像；(2)获得高分辨率联合采样值矩阵；2a)对低分辨率多光谱图像和高分辨率全色图像进行分块列向量化；2b)将分块列向量化的高分辨率全色图像矩阵连接到分块列向量化的低分辨率多光谱图像矩阵后，得到高分辨率联合采样值矩阵；(3)对高分辨率的多光谱图像和低分辨率的多光谱图像构造观测矩阵，获得高-低观测矩阵；(4)对高分辨率的多光谱图像和高分辨率的全色图像构造观测矩阵，获得高-高观测矩阵；(5)将高-低观测矩阵连接到高-高观测矩阵后，得到高分辨率联合观测矩阵；(6)利用KSVD方法对高分辨率多光谱图像进行字典学习，得到高分辨率过完备字
Ffti. ·
ZN (7)利用基追踪算法求解高分辨率过完备字典下的高分辨率稀疏系数；(8)重构高分辨率多光谱图像；8a)将高分辨率过完备字典与高分辨率稀疏系数相乘，得到列向量化的多光谱图像；8b)将列向量化的多光谱图像转化为图像块，恢复到原图像的位置，得到超分辨率的多光谱图像。本发明与现有的技术相比具有以下优点第一，本发明将多光谱图像融合引入到压缩感知领域来，克服了现有技术中颜色失真和光谱的扭曲的缺陷，使得本发明可以保持空间分辨率和光谱信息，有利于后期对图像的处理和识别。第二，本发明引入了矩阵联合和字典学习的方法，克服了现有技术中融合图像分辨率受限于融合前全色图像分辨率的缺陷，使得本发明得到的融合图像分辨率比现有技术的融合前全色图像分辨率明显提高。第三，本发明引入了压缩感知基追踪算法，克服了现有技术中传感器和融合过程中噪声的影响，使得本发明比现有技术的峰值信噪比大大提高。本发明的技术过程和效果可结合以下附图详细说明。
图1为本发明的流程图；图2为本发明的仿真图。具体实施方法参照附图1，本发明具体实施方式
如下
步骤1，分别输入低分辨率多光谱图像和高分辨率全色图像；分别读入待处理的低分辨多光谱图像和高分辨率的全色图像。本发明实施例中输入的低分辨多光谱图像大小为125X 125X4，分辨率为細；高分辨率的全色图像大小为500 X 500，分辨率为lm。步骤2，获得高分辨率联合采样值矩阵；将低分辨率多光谱图像的红、绿、蓝、近红外四个颜色通道图像分为5X5的小块然后列向量化，将第二块构成的列向量连接到第一块构成的列向量后，将第三块构成的列向量连接到第二块构成的列向量后，以此类推，低分辨率的多光谱图像变换为一个行数为 100的矩阵，本发明实施例中低分辨率的多光谱图像变换为一个大小为100X625的矩阵。将高分辨率全色图像分为20X20的小块然后列向量化，将第二块构成的列向量连接到第一块构成的列向量后，将第三块构成的列向量连接到第二块构成的列向量后，以此类推，高分辨率的全色图像变成一个行数为400的矩阵，本发明实施例中高分辨率的多光谱图像变换为一个大小为400X625的矩阵。将高分辨率全色图像矩阵连接到低分辨率多光谱图像矩阵后，得到一个行数为500的高分辨率联合采样值矩阵，本发明实施例中高分辨率联合采样值矩阵大小为 500X625o步骤3，对高分辨率的多光谱图像和低分辨率的多光谱图像构造观测矩阵，获得高-低观测矩阵；将低分辨率多光谱图像的红、绿、蓝、近红外四个颜色通道图像分为5X5的小块然后列向量化，将高分辨率多光谱图像的红、绿、蓝、近红外四个颜色通道图像分为40X40 的小块然后列向量化；将高分辨率多光谱图像每8个像素取平均得到低分辨率多光谱图像的1个像素， 将高分辨率多光谱图像的8X8的一个小块取平均得到低分辨率多光谱图像的1 X 1的一个小块，高分辨率多光谱图像与低分辨率多光谱图像的下采样比为8 1；本发明实施例中按下式求解高-低观测矩阵Y^ =MX^gh+VpIf*'为低分辨率的多光谱图像，X^a为高分辨率的多光谱图像，M为待求解的高低观测矩阵，Vp为噪声向量，本发明实施例中高-低观测矩阵大小为100X6400。步骤4，对高分辨率的多光谱图像和高分辨率的全色图像构造观测矩阵，获得高-高观测矩阵；将高分辨率多光谱图像的红、绿、蓝、近红外四个颜色通道图像分40X40的小块然后列向量化，将高分辨率全色图像分为40X40的小块然后列向量化，将低分辨率全色图像分为20X20的小块然后列向量化；将高分辨率多光谱图像的红、绿、蓝、近红外四个颜色通道的每一个像素加权求和，得到高分辨率全色图像对应位置的一个像素，本发明实施例中按下式求解权值矩阵yPAN = M2x+v2Ypan为高分辨率的全色图像，Y^a为高分辨率的多光谱图像，M2为待求解的权值矩阵，Vp为噪声向量。本发明实施例中红、绿、蓝、近红外每一谱段所取的权值系数大小分别为0. 1139,0. 2315,0. 2308,0. 4239，权值矩阵大小为 1600X6400。将高分辨率全色图像每2个像素取平均得到低分辨率全色图像的1个像素，将高分辨率全色图像的2X2的一个小块取平均得到低分辨率全色图像的1X1的一个小块，高分辨率多光谱图像与低分辨率多光谱图像的下采样比为2 1，本发明实施例中按下式求解高分辨率全色图像下采样矩阵
权利要求
1.一种基于压缩感知的多光谱与全色图像超分辨融合方法，包括以下步骤(1)分别输入低分辨率多光谱图像和高分辨率全色图像；(2)获得高分辨率联合采样值矩阵；2a)对低分辨率多光谱图像和高分辨率全色图像进行分块列向量化； 2b)将分块列向量化的高分辨率全色图像矩阵连接到分块列向量化的低分辨率多光谱图像矩阵后，得到高分辨率联合采样值矩阵；(3)对高分辨率的多光谱图像和低分辨率的多光谱图像构造观测矩阵，获得高-低观测矩阵；(4)对高分辨率的多光谱图像和高分辨率的全色图像构造观测矩阵，获得高-高观测矩阵；(5)将高-低观测矩阵连接到高-高观测矩阵后，得到高分辨率联合观测矩阵；(6)利用KSVD方法对高分辨率多光谱图像进行字典学习，得到高分辨率过完备字典；(7)利用基追踪算法求解高分辨率过完备字典下的高分辨率稀疏系数；(8)重构高分辨率多光谱图像；8a)将高分辨率过完备字典与高分辨率稀疏系数相乘，得到列向量化的多光谱图像； 8b)将列向量化的多光谱图像转化为图像块，恢复到原图像的位置，得到超分辨率的多光谱图像。
2.根据权利要求1所述的基于压缩感知的多光谱与全色图像超分辨融合方法，其特征在于步骤(3)所述构造高-低观测矩阵的方法具体步骤如下第一步，对低分辨率的多光谱图像和高分辨率的多光谱图像进行分块列向量化； 第二步，对高分辨率多光谱图像8 1下采样，得到低分辨率的多光谱图像； 第三步，由高分辨率的多光谱图像和低分辨率的多光谱图像构造高-低观测矩阵。
3.根据权利要求1所述的基于压缩感知的多光谱与全色图像超分辨融合方法，其特征在于步骤(4)所述构造高-高观测矩阵的方法具体步骤如下第一步，对高分辨率的多光谱图像、高分辨率的全色图像、低分辨率的全色图像进行分块列向量化；第二步，对高分辨率多光谱图像的每一谱段求权值之和，得到高分辨率全色图像； 第三步，对高分辨率全色图像2 1下采样，得到低分辨率的全色图像； 第四步，由高分辨率的多光谱图像和高分辨率的全色图像构造高-高观测矩阵。
4.根据权利要求1所述的基于压缩感知的多光谱与全色图像超分辨融合方法，其特征在于步骤(6)所述KSVD方法训练字典的具体步骤如下第一步，对输入的与待融合图像统计特性相近且比全色图像高一倍分辨率的多幅高分辨率多光谱图像进行分块列向量化，随机抽取高分辨率多光谱图像的25000个图像块进行字典训练；第二步，按照下式求解训练字典D
5.根据权利要求1所述的基于压缩感知的多光谱与全色图像超分辨融合方法，其特征在于步骤(7)所述的基追踪算法如下 min Il α H1 Subject to Ι^-Φα^ <ε其中，minll α ||工为最小化α的一范数，α为待求解的高分辨率稀疏系数，Subject 切表示||少-0>0^为求解时11|| α II1W限制条件，||广0)<为求解y-Φα的二范数，y为高分辨率联合采样值矩阵，Φ为高分辨率联合观测矩阵与高分辨率训练字典的乘积，ε为重构误差。
全文摘要
本发明公开了一种基于压缩感知的多光谱与全色图像超分辨融合方法。其步骤为(1)分别输入低分辨率多光谱图像和高分辨率全色图像；(2)获得高分辨率联合采样值矩阵；(3)获得高-低观测矩阵；(4)获得高-高观测矩阵；(5)获得高分辨率联合观测矩阵；(6)学习高分辨率过完备字典；(7)获得高分辨率稀疏系数；(8)重构高分辨率多光谱图像。本发明将压缩感知技术引入到多光谱图像融合领域来，克服了现有技术中颜色失真和光谱扭曲的缺陷，融合得到的多光谱图像比高分辨率全色图像高一倍分辨率。
文档编号G06T7/00GK102542549SQ20121000172
公开日2012年7月4日 申请日期2012年1月4日 优先权日2012年1月4日
发明者侯彪, 刘芳, 杨丽霞, 杨淑媛, 焦李成, 王爽, 陈璞花, 靳焕庭, 马文萍 申请人:西安电子科技大学