基于光谱重构的三步高光谱图像融合方法
【专利摘要】基于光谱重构的三步高光谱图像融合方法，即通过波段选择、光谱重构的处理方法对主成分分析图像融合进行完善，得到更好的融合效果。该方法如下：1）通过人机交互接口模块获得待融合的高光谱图像数据和高空间分辨率图像数据，并初始化相关参数；2）通过高光谱图像波段选择模块，设计波段选择矩阵，对初始高光谱图像进行波段选择；3）通过高光谱图像-高空间分辨率图像融合模块，对输入的图像进行主成分分析融合；4）通过光谱重构模块，对初始融合图像进行高精度光谱重构；5）通过融合结果输出模块，输出最终的融合结果。本发明的优点在于减少融合数据量，降低光谱畸变，提升光谱保持性，增强图像质量，该方法融合质量好，使用范围广。
【专利说明】基于光谱重构的三步高光谱图像融合方法
【技术领域】
[0001]本发明涉及一种基于光谱重构的三步高光谱图像融合方法，适用于高光谱图像处理系统，属于高光谱遥感图像处理领域。
【背景技术】
[0002]随着20世纪80年代初期成像光谱概念的出现，光学遥感进入了 一个崭新的阶
段-高光谱遥感(Hyperspectral Remote Sensing)阶段。高光谱遥感指集探测器技术、
精密光学机械、微弱信号检测、计算机技术、信息处理技术于一体的，能够获取具有高的光谱分辨率图像的遥感科学和技术，是当前模式识别及遥感信息处理研究领域中的一个研究执占。
[0003]与常规遥感图像相比，高光谱遥感图像具有“图谱合一”的突出特点，既能表达空间信息，又能表达光谱信息，以三维数据立方体的形式呈现。在遥感应用中，有时会要求图像同时具有高的空间分辨率和高的光谱分辨率。然而，由于当前技术条件的限制，现有高光谱成像仪器很难提供这样的数据。高光谱图像拥有较高的光谱分辨率，可达10nm，甚至更高，但其空间分辨率较低，即空间的细节表现能力较差；而普通遥感图像(如全色图像等)具有相对较高的空间分辨率，但不能表达光谱信息。因此，将高光谱遥感图像和对应的具有高的空间分辨率的遥感图像进行融合，可优化图像质量，得到新的高光谱图像，既具有较高的空间细节表现能力，同时又保留图像的光谱特性。
[0004]目前很多传统多源遥感数据融合方法只适用于普通遥感图像间的融合以及多光谱图像与全色数据的融合，由于高光谱图像的巨大数据量，使得传统融合方法效率低、难适用。在图像融合领域，使用多源遥感数据进行融合的方法主要包括直接加权法、主成分分析(Principal Component Analysis, PCA)法、高通滤波(High-Pass Filter, HPF)算法和小波变换(Wavelet Translation, WT)法等。尽管前三种方法提高了高光谱数据各波段图像的空间分辨率，但融合后的图像光谱失真严重，容易增加噪声，甚至影响视觉效果；而小波变换法保持了原始光谱特性，但空间分辨率效果较差。因此有必要采用新的方式解决高光谱图像融合问题。
[0005]PCA是一种统计相关分析技术，在多数实际问题中，由于指标较多及指标间有一定的相关性，势必增加分析问题的复杂性，PCA就是设法将原有指标重新组合成一组新的相互无关的几个综合指标来代替原来指标，是研究如何通过少数几个主成分来解释多变量的协方差结构的一种方法。PCA变换后的新子空间中，图像间的相关性基本消除，其中的第一主分量包含了绝大部分的信息。
[0006]压缩感知是当前信号重构研究的热门理论，该理论指出，在一定条件下，可用少量压缩信号高精度重构原始信号，由此衍生了一系列信号重构方法，如正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法等,利用最小二乘求解的贪婪方法快速求解原始信号。
[0007]本发明将压缩感知波段重构技术引入融合过程进行图像优化，首先对高光谱图像进行波段选择，再进行高光谱图像与高空间分辨率图像的PCA融合，最后再对其进行波段维的重构过程，得到同时包含高空间分辨率空间信息和光谱信息的图像，形成一套新的高光谱图像融合方法体系。
【发明内容】
[0008]本发明的目的在于提供一种基于光谱重构的三步高光谱图像融合方法，发明了一种运用压缩感知光谱重构技术进行波段选择，从而以少量数据进行快速、精确融合，最终再高质量重构光谱维信息的多源图像融合方法。本方法既提高融合过程的效率，又增强了高光谱图像融合的光谱保持性，融合图像性能好，使用范围广。
[0009]本发明所涉及的方法流程包括以下步骤:(1)获得待融合的高光谱图像和高空间分辨率图像数据，初始化相关参数；(2)利用波段选择矩阵对高光谱图像进行波段选择，以减少融合数据量；(3)计算高光谱图像的主成分，并用高空间分辨率图像与其进行PCA融合；(4)通过压缩感知光谱重构方法进行光谱的高精度重构；(5)输出最终的高光谱图像融合结果。
[0010]下面对该方法流程各步骤进行详细说明:
[0011](1)获得待融合的高光谱图像和高空间分辨率图像数据，初始化相关参数:待融合数据为波段数为n的高光谱图像X= (X1, X2,, Xn)T和高空间分辨率图像X，Xi (i=l, 2，…，n)表示第i波段的图像，设置合适的波段采样比r。
[0012](2)利用波段选择矩阵对高光谱图像进行波段选择:按照合理的方式和采样比对高光谱图像进行波段选择，通常选择等间隔波段选择，采样比为大于I的正整数，即r=n/m=2, 3，…，其中m为选择后的高光谱图像波段数，通过波段选择矩阵O对高光谱图像进行波段选择，即得到s=Ox。
[0013](3)计算高光谱图像的主成分，并用高空间分辨率图像与其进行PCA融合:对采样后的高光谱图像s进行PCA变换，得到第一主成分PC;，将X与PC;作直方图匹配得到PCnew,用PCnew代替/5C1 ):与其余主成分一起做PCA逆变换，得到初始融合图像S。
[0014](4)通过压缩感知光谱重构方法进行光谱的高精度重构:
[0015]
【权利要求】
1.基于光谱重构的三步高光谱图像融合方法，其特征在于:通过波段选择和重构提高融合过程的效率，并增强光谱保持性、稳定性，具有人机交互接口模块、高光谱图像波段选择模块、高光谱图像-高空间分辨率图像融合模块、光谱重构模块、融合结果输出模块这五个功能模块。包括如下步骤: (1)获得待融合的高光谱图像数据和高空间分辨率图像数据，初始化相关参数:假设X为待融合的高空间分辨率图像，X= (X1, X2, , xn)T为待融合的低空间分辨率的高光谱图像，其中n是高光谱数据的波段总数，Xi, i=l,…，n代表第i波段的图像响应数值，完成初始参数设置步骤； (2)设置合适的采样比r，根据设置的采样比和起始波段设计得到波段选择矩阵O的表达，通过计算式S=Ox获得降维后的高光谱图像S，完成波段选择步骤； (3)利用降维的高光谱图像s与空间分辨率较高的遥感图像X进行PCA融合:对采样后的高光谱图像s进行PCA变换，得到第一主成分，将高空间分辨率图像X与PC: ?作直方图匹配得到PCnew,用PCnew代替尸C1后与其余主成分一起做PCA逆变换，得到初始融合图像S ； (4)基于压缩感知光谱重构方法对图像S进行重构，根据初始参数设计正交的稀疏变换基屯，通过S=OX’=O Wy求解图像的稀疏表示y，并经过正交反变换得到重构图像
【文档编号】G06T5/50GK103714572SQ201310656056
【公开日】2014年4月9日 申请日期:2013年12月6日 优先权日:2013年12月6日
【发明者】尹继豪, 孙建颖, 姜志国, 余万科, 朱红梅 申请人:北京航空航天大学