一种融合生成高分辨率多光谱图像的方法技术

本发明专利技术提供一种融合生成高分辨率多光谱图像的方法，其特征在于包括以下步骤：对来自光学遥感卫星系统的原始低分辨率多光谱图像和与其空间配准的高分辨率全色图像，利用低分辨率多光谱波段和高分辨率全色波段之间的统计关系，合成出仿真的低分辨率全色图像，再利用高分辨率全色图像和仿真的低分辨率全色图像确定不同分辨率图像之间的细节图像，然后按比例将细节图像添加于原始低分辨率多光谱图像中生成高分辨率多光谱图像。采用该方法能够有效消除或减少融合的高分辨率多光谱图像所存在的光谱畸变，获得真正的高分辨率多光谱图像。

【技术实现步骤摘要】

本专利技术涉及遥感图像处理数据融合
，特别是光学遥感卫星同一传感器低分辨率多光谱波段和高分辨率全色波段融合生成高分辨率多光谱图像的技术，尤其是，适用于对艾科诺斯卫星等遥感图像处理。
技术介绍
许多光学遥感卫星系统同时提供了空间分辨率低的多光谱图像和空间分辨率高的全色图像，如艾科诺斯卫星同时提供了空间分辨率是1米的全色图像和空间分辨率是4米的包含蓝、绿、红、近红外四个波段的多光谱图像。为了信息的优化利用，尤其是利用全色波段的高分辨率信息，许多图像融合技巧努力将全色波段的高空间分辨率信息与多光谱信息结合起来，以生成空间分辨率是1米的多光谱图像，同时尽可能地保持原始多光谱波段的光谱信息。图像融合的前提是待融合的图像是空间配准的，配准误差在一个像元以内。目前比较常用的图像融合方法是高分辨率全色波段的分量替换方法，如PCA(主成分分析)、HSI方法(色度-饱和度-亮度，彩色空间变换)等。分量替换的融合方法如图1所示，以PCA(主成分分析)方法为例，包括以下步骤①图1中的步骤2对原始低分辨率多光谱图像各波段进行重采样操作，获得超象素多光谱图像。超象素是指多光谱图像一个像元的地面覆盖面积与对应的全色波段一个像元的地面覆盖面积相同。重采样的比例是原始低分辨率多光谱图像的像元空间分辨率与原始的全色高分辨图像像元空间分辨率的比。②图1中的步骤4将低分辨率超象素多光谱图像的所有波段进行主成分变换，变换后生成第一主成分(PC1)、第二主成分(PC2)，第三主成分(PC3)，第四主成分(PC4)。③图1中的步骤6将高分辨率全色波段进行直方图变换，使高分辨率全色波段与第一主成分分量的灰度特征趋向一致。④图1中的步骤8用经过直方图变换的高分辨率全色波段替换第一主成分。⑤图1中的步骤10主成分反变换，将经过直方图变换的高分辨率全色波段与第二主成分(PC2)，第三主成分(PC3)，第四主成分(PC4)一起进行主成分反变换，变换得到的各分量就是融合的高分辨率的多光谱图像。图1中各步骤的功能如下 对原始低分辨率多光谱图像各波段重采样是为了获得超象素多光谱图像。超象素多光谱图像一个像元的地面覆盖面积与对应的全色波段一个像元的地面覆盖面积相同，因而超象素多光谱图像与高分辨率全色图像有完全相同的地面覆盖区域。主成分变换是可逆变换，是为了获得与高分辨率全色波段相近似的新的分量。高分辨率全色波段的直方图匹配是为了使其灰度特征与第一主成分尽可能相似，以便将它替换第一主成分进行主成分反变换后，生成的高分辨率的多光谱图像尽可能保持原始低分辨率多光谱图像的光谱信息。主成分反变换将图像变换为原来的空间，获得融合了高分辨率全色波段细节的高分辨率多光谱图像。上述融合方法有不足之处高分辨率全色波段与第一主成分的数据特性有很大差别。直方图匹配就是为了缩小二者之间的差别，但不能消除二者的差别，从而导致融合的高分辨率多光谱图像存在光谱畸变。
技术实现思路
本专利技术针对现有技术中存在的缺陷或不足，提供，采用该方法能够有效消除或减少融合的高分辨率多光谱图像所存在的光谱畸变，获得真正的高分辨率多光谱图像。本专利技术的技术方案如下，其特征在于包括以下步骤对来自光学遥感卫星系统的原始低分辨率多光谱图像和与其空间配准的高分辨率全色图像，利用低分辨率多光谱波段和高分辨率全色波段之间的统计关系，合成出仿真的低分辨率全色图像，再利用高分辨率全色图像和仿真的低分辨率全色图像确定不同分辨率图像之间的细节图像，然后按比例将细节图像添加于原始低分辨率多光谱图像中生成高分辨率多光谱图像。所述细节图像，通过计算高分辨率全色图像与仿真的低分辨率全色图像的差而获得。所述细节图像添加比例是不同的，不同的比例分别对应于原始低分辨率多光谱图像的各不同波段。所述合成出仿真的低分辨率全色图像包括以下步骤对原始低分辨率多光谱图像各波段进行重采样操作，获得超像素多光谱图像；将超像素多光谱图像的各个波段作为自变量，将全色的高分辨率图像作为因变量进行多元线性回归，合成出仿真的超像素低分辨率全色图像。所述多元线性回归的线性项系数为多光谱波段细节与全色波段细节的比，并作为细节图像的添加比例。所述空间配准的配准误差在一个像元以内。本专利技术的技术效果如下本专利技术利用低分辨率多光谱图像和高分辨率全色波段的统计关系获得仿真的低分辨率全色波段，然后利用高分辨率全色波段和仿真的低分辨率全色波段，获得图像的细节(高频)信息，最后只将细节信息按比例合理添加到多光谱图像各波段中，从而生成高分辨率的多光谱图像。本专利技术包括重采样低分辨率多光谱图像、合成仿真低分辨率全色图像、获得高频细节图像、高频细节按比例添加到低分辨率多光谱图像生成高分辨率多光谱图像四个部分。重采样低分辨率多光谱图像，获得超象素多光谱图像；将高分辨率全色图像和超象素多光谱图像实施多元线性回归，利用超象素多光谱图像合成仿真低分辨率超象素全色图像；然后利用高分辨率全色波段和仿真的低分辨率全色波段，获得图像的细节(高频)信息，最后只将细节信息接比例合理添加到多光谱图像各波段中，从而生成高分辨率的多光谱图像。这样做的理由是在同样的空间分辨率条件下，构成全色图像的多光谱图像各波段的信息比例是不同的，从而导致构成全色图像多光谱各波段高频细节的信息比例也不同。本专利技术与现有技术相比有以下特点第一，将高分辨率全色图像和超象素多光谱图像实施多元线性回归，利用超象素多光谱图像合成仿真低分辨率超象素全色图像；如此得到的低分辨率超象素全色图像比通过主成分变化获得的第一主成分与全色图像的特性更相近。第二，利用高分辨率全色波段和仿真的低分辨率全色波段，获得图像的细节(高频)信息。然后，将细节信息按比例合理添加到多光谱图像各波段中，从而生成高分辨率的多光谱图像。这样获得的高分辨率的多光谱图像与低分辨率多光谱图像相比只是增加了细节。附图说明图1为分量替换数据融合的方法框图，以PCA(主成分分析)方法为例。图2为本专利技术多光谱波段和全色波段数据融合的流程图。图3为本专利技术图2中图像融合单元16的融合过程图。图4为本专利技术图3中细节融合单元20的多光谱波段与高频细节图像叠加的实现流程图。具体实施例方式下面结合附图对本专利技术作进一步的详细说明。本专利技术方法实现艾科诺斯数据融合的流程如图2所示。图2包括重采样单元12，合成仿真单元14和图像融合单元16。单元12对原始低分辨率多光谱图像按波段进行重采样，得到超象素多光谱图像。重采样方法为双三次卷积，按行和列依次进行。双三次卷积方法产生的图像在保持图像信息和图像的光滑性方面有优越性，能够避免块状和锯齿效应。重采样。对原始低分辨率多光谱图像按波段进行重采样，得到超象素多光谱图像。然后进入合成仿真单元14。单元14分别将超象素多光谱四个波段的图像像素按行排列，形成四个向量x1＝{blueBand(i)|i＝0，1，…，n-1}；x2＝{greenBand(i)|i＝0，1，…，n-1}；x3＝{redBand(i)|i＝0，1，…，n-1}；x4＝{nirBand(i)|i＝0，1，…，n-1}，同样将高分辨率全色波段像素按行排列，形成向量。超象素多光谱四个波段向量分量个数与全色波段像素个数相同，是图像像素的总数；将超象素多光谱图像的各个波段作为自变量，将全色的高分辨率图像作为应本文档来自技高网...

【技术保护点】
一种融合生成高分辨率多光谱图像的方法，其特征在于包括以下步骤：对来自光学遥感卫星系统的原始低分辨率多光谱图像和与其空间配准的高分辨率全色图像，利用低分辨率多光谱波段和高分辨率全色波段之间的统计关系，合成出仿真的低分辨率全色图像，再利用高分辨率全色图像和仿真的低分辨率全色图像确定不同分辨率图像之间的细节图像，然后按比例将细节图像添加于原始低分辨率多光谱图像中生成高分辨率多光谱图像。

【技术特征摘要】

【专利技术属性】
技术研发人员：唐娉，
申请(专利权)人：中国科学院遥感应用研究所，
类型：发明
国别省市：11[中国|北京]