一种星载全色与多光谱图像融合方法技术

本发明专利技术提供了一种星载全色与多光谱图像融合方法，首先利用高斯导向滤波联合全色图像完成对原始多光谱图像的重采样放大，通过该联合上采样过程将全色图像的空间细节信息注入至放大的多光谱图像中，再对原始多光谱图像使用超球彩色变换提取其光谱信息，在超球彩色空间中将光谱信息和全色图像的空间信息融合，经HCT反变换得到具有高空间分辨率的多光谱图像。本发明专利技术提高了重采样多光谱图像的空间分辨率，可以直接拓展至高光谱图像的融合处理，普遍性和实用性更强。

A fusion method for spaceborne panchromatic and multispectral images

The invention provides a method for fusion of satellite-borne panchromatic and multispectral images. Firstly, the original multispectral image is re-sampled and enlarged by Gaussian-guided filtering combined with panchromatic image, and the spatial detail information of the panchromatic image is injected into the enlarged multispectral image through the joint up-sampling process, and then the original multispectral image is enlarged. The hypersphere color transform is used to extract the spectral information of the image. The spectral information and the spatial information of the panchromatic image are fused in the hypersphere color space. The multispectral image with high spatial resolution is obtained by inverse HCT transform. The method improves the spatial resolution of the resampled multi-spectral image, and can be directly extended to the fusion processing of the hyperspectral image, thus having stronger universality and practicability.

【技术实现步骤摘要】
一种星载全色与多光谱图像融合方法
本专利技术涉及星载遥感图像的可视化增强处理技术，尤其是针对星载全色图像与多光谱图像的图像融合方法。
技术介绍
由于成像传感器自身的物理局限性以及星地间数据传输技术的壁垒，目前的星载成像传感器仅能提供光谱分辨率高但空间分辨率低的多波段多光谱图像或者空间分辨率高但光谱信息贫乏的单波段全色图像。在轨的星载传感器，例如GeoEye-1、WorldView-2/3/4、高景1号等能够同时提供多波段多光谱遥感图像和单波段高空间分辨率的全色图像。通过融合多光谱图像和全色图像，可以在保持多光谱图像光谱信息的同时提高其空间分辨率，融合结果可服务于高精度地物分类、城市规划、地理信息系统等多个遥感应用。全色图像空间分辨率往往高于多光谱图像，因此融合时首先要将多光谱图像重采样放大至与全色图像同尺寸，然后再将放大后的多光谱图像和全色图像输入至融合模型，输出得到融合的多光谱图像。在对多光谱图像的重采样处理上，目前国内外通用的做法是选用经典的插值方法包括最近邻插值法、双线性插值法或者双三次插值法对多光谱图像的各波段独立地进行插值放大(参见文献IEEETransactionsonGeoscienceandRemoteSensing,45(10):3012-3021,2007)。事实上，由上述插值法得到的多光谱图像会出现不同程度的光谱失真与空间信息失真，从而进一步影响融合图像的信息保持(参见文献IEEEGeoscienceandRemoteSensingLetters,4(1):27-31,2007；IEEEJournalofSelectedTopicsinSignalProcessing,5(3):446-453,2011)。此外，经典的插值方法不能够提高重采样放大图像的空间分辨率，因为在插值过程中没有利用外部有用的高空间分辨率细节信息，所以重采样放大图像的清晰度并没有得到提高。
技术实现思路
为了克服现有技术的不足，本专利技术针对星载多光谱图像和全色图像提出一种图像融合方法，首先利用高斯导向滤波(Gaussianguidedfiltering)联合全色图像完成对原始多光谱图像的重采样放大，通过该联合上采样过程将全色图像的空间细节信息注入至放大的多光谱图像中，再对原始多光谱图像使用超球彩色变换(hypersphericalcolortransform,HCT)提取其光谱信息，在超球彩色空间中将光谱信息和全色图像的空间信息融合，经HCT反变换得到具有高空间分辨率的多光谱图像。本专利技术解决其技术问题所采用的技术方案包括以下步骤：第一步、对于包含N个波段M1,M2,…,MN的原始多光谱图像，对应的全色图像记为PH，通过卷积运算对全色图像PH进行低通滤波，低通滤波器记为LPF，输出图像PL＝PH*LPF；对PL的像素点进行隔行取样和隔列取样，输出结果PL＝(PL)↓2；第二步、设置图像PL为高斯导向滤波的指导图像，记为I，任意第j个多光谱波段Mj为高斯导向滤波的输入图像，记为P，高斯导向滤波输出图像记为Q；假设滤波输出图像Q与指导图像I之间存在局部线性关系，输出图像Q在像素点i的像素值其中像素点i位于以像素点k为中心的窗口wk中，|w|为窗口wk内所有像素的个数，该窗口半径为r，窗口大小为(2r+1)×(2r+1)像素，窗口wk内的高斯权系数定义为Xi为窗口内像素点i的像素值，Xk为窗口内像素点k的像素值，ak和bk为窗口系数，其中Pi为输入图像P在i点的像素值，ε是规则化参数且ε∈(0,1)；将Q表示成矩阵形式，得到Q＝A×I+B，其中A和B是窗口系数矩阵；第三步、利用最近邻插值法对A和B分别做插值放大，放大至与全色图像同尺寸，得到对应的新窗口系数矩阵和第四步、利用全色图像PH和和计算出大尺寸的输出图像重复第二步至第四步对原始多光谱图像的各波段进行滤波放大，得到对应的滤波输出图像第五步、利用超球彩色变换到超球彩色空间，得到一个强度分量IM和N-1个角度分量θ1，θ2，…，θN-1，利用双线性插值法对N-1个角度分量θ1,θ2,…,θN-1分别做插值放大，放大至与全色图像同尺寸，分别记为构建新强度分量第六步、对第五步得到的N-1个角度分量和新强度分量INEW做HCT反变换，得到的N波段输出图像记为F1,F2,…,FN；第七步、计算高斯导向滤波输出图像与全色图像PH的均值与方差，得到匹配后的全色图像j＝1,2,…,N，其中μP与μj分别为全色图像PH和图像的均值，σP与σj分别为全色图像PH和图像的标准差；计算各波段的待注入空间细节计算第j个波段的融合结果Sj＝Fj+Dj，得到所有波段的融合结果(S1,S2,…,SN)。本专利技术的有益效果是：先后包含了两个过程去提高多光谱图像的空间分辨率，首先使用高斯导向滤波联合全色图像，将全色图像的空间细节信息部分注入至重采样放大的多光谱图像中，提高了重采样多光谱图像的空间分辨率；然后利用空间细节的再注入(第八步)过程进一步提高HCT反变换输出图像的空间分辨率。此外，本专利技术利用HCT变换有效地分离出多光谱图像的原始光谱信息，经HCT反变换将光谱信息融入至重采样的多光谱图像中，有利于提高融合结果的光谱保真能力。本专利技术与专利：一种基于导向滤波的遥感图像融合方法(申请号：2017102421174，公开号：CN107146212A)主要的区别之处在于：1)本专利技术使用的高斯导向滤波不同于传统的导向滤波，能够克服传统导向滤波在滤波窗口输出沿水平与垂直方向产生的偏执、不均匀现象；2)IHS(Intensity-Hue-Saturation)变换每次只能输入三个波段的多光谱图像，而本专利技术使用HCT变换可以突破对输入多光谱图像波段数量的限制，可以直接拓展至高光谱图像的融合处理；3)本专利技术不受全色图像与多光谱图像光谱响应匹配关系的影响，在对各波段注入空间细节时，不需要对待注入的空间细节进行专门的调制，普遍性和实用性更强。附图说明图1是本专利技术的原理示意图；图2是本专利技术的流程图。具体实施方式下面结合附图和实施例对本专利技术进一步说明，本专利技术包括但不仅限于下述实施例。假设原始多光谱图像包含N个波段，记为M1,M2,…,MN，各波段图像尺寸均为R1行×C1列，全色图像(记为PH)的尺寸为R2行×C2列。第一步、下采样全色图像：通过卷积运算对全色图像PH进行低通滤波，低通滤波器记为LPF，输出图像记为PL，*为卷积运算符PL＝PH*LPF(1)对PL作下2抽样处理，即对像素点隔行取样和隔列取样，输出结果记为PL，其尺寸减小至原全色图像尺寸的一半，即R2/2行×C2/2列PL＝(PL)↓2(2)假定全色图像的空间分辨率是原始多光谱图像空间分辨率的2倍，即R2/2＝R1，C2/2＝C1，公式(2)输出的下采样全色图像PL与原始多光谱图像同尺寸。第二步、高斯导向滤波设置图像PL为高斯导向滤波的指导图像，记为I，任意第j个多光谱波段(Mj)为高斯导向滤波的输入图像，记为P，高斯导向滤波输出图像记为Q。假设滤波输出图像Q与指导图像I之间存在局部线性关系，即其中像素点i位于以像素点k为中心的邻域(或窗口)wk中，该窗口半径为r，窗口大小为(2r+1)×(2r+1)像素(通常取半径r≥1)。窗口wk内的高斯权系数定义为其本文档来自技高网...

【技术保护点】
1.一种星载全色与多光谱图像融合方法，其特征在于包括下述步骤：第一步、对于包含N个波段M1,M2,…,MN的原始多光谱图像，对应的全色图像记为PH，通过卷积运算对全色图像PH进行低通滤波，低通滤波器记为LPF，输出图像PL＝PH*LPF；对PL的像素点进行隔行取样和隔列取样，输出结果PL＝(PL)↓2；第二步、设置图像PL为高斯导向滤波的指导图像，记为I，任意第j个多光谱波段Mj为高斯导向滤波的输入图像，记为P，高斯导向滤波输出图像记为Q；假设滤波输出图像Q与指导图像I之间存在局部线性关系，输出图像Q在像素点i的像素值

【技术特征摘要】
1.一种星载全色与多光谱图像融合方法，其特征在于包括下述步骤：第一步、对于包含N个波段M1,M2,…,MN的原始多光谱图像，对应的全色图像记为PH，通过卷积运算对全色图像PH进行低通滤波，低通滤波器记为LPF，输出图像PL＝PH*LPF；对PL的像素点进行隔行取样和隔列取样，输出结果PL＝(PL)↓2；第二步、设置图像PL为高斯导向滤波的指导图像，记为I，任意第j个多光谱波段Mj为高斯导向滤波的输入图像，记为P，高斯导向滤波输出图像记为Q；假设滤波输出图像Q与指导图像I之间存在局部线性关系，输出图像Q在像素点i的像素值其中像素点i位于以像素点k为中心的窗口wk中，|w|为窗口wk内所有像素的个数，该窗口半径为r，窗口大小为(2r+1)×(2r+1)像素，窗口wk内的高斯权系数定义为Xi为窗口内像素点i的像素值，Xk为窗口内像素点k的像素值，ak和bk为窗口系数，其中Pi为输入图像P在i点的像素值，ε是规则化参数且ε∈(0,1)；将Q表示成矩阵形式，得到Q＝A×I+B，...

【专利技术属性】
技术研发人员：李旭，潘煜，高昂，李立欣，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61