1.本发明涉及时空谱遥感数据融合技术领域,尤其是一种面向高光谱图像和多光谱图像的时空谱融合方法。
背景技术：2.近年来随着高光谱成像技术的快速发展,高光谱遥感数据在各领域应用研究中取得了良好的发展与突破。高光谱遥感数据具有波段多，光谱分辨率高的优点，但高光谱数据也存在着空间分辨率低，卫星种类少而导致高光谱数据覆盖频次低的缺点，单一传感器无法同时满足同类数据的高空间、高光谱、高时间特性，大大制约了高光谱数据的利用程度。多传感器的数据融合能够有效解决这一问题，有效提高光谱数据的利用效率。
3.时空谱遥感数据融合是将两个或者多个影像数据的时间、空间、光谱信息进行协同组合。数据融合的概念和方法发展于20世纪70年代，有文献记录的遥感数据融合起源于20世纪90年代，21世纪以来，众多学者开展了融合技术和应用研究，国内外研究文献近万篇，也有众多学者结合经典方法和新兴高光谱数据进行了测试验证，对融合方法选择起到了很好的指导意义。目前的融合方法按指标提升可分为三类：空间维提升(成分替换法、多分辨率分析法)、光谱维提升(稀疏解混法、人工智能法)和时间维提升(权重函数法、线性优化分解法、人工智能法)，从而衍生了空谱、时空、时空谱融合的研究。空谱融合多采用模型优化的方法解决低分辨率像元的混合像元问题，包括基于锐化的方法(gram-smidt变换、brovey算法、颜色空间变换法、主成分分析法、空谱采样建模法如稀疏表达和神经网络)，基于成像模型的方法(基于混合像元分解、张量分解、耦合非负矩阵分解(cnmf)、低秩分解(lrf))，基于张量的融合方法，基于深度网络的方法如引入残差注意力网络提升空间和光谱的融合信息等，然而目前研究高光谱与其他数据融合的研究并未考虑高光谱反演特征的时序变化。时空融合主要考虑将时间序列影像的变化信息注入基准影像，发展至今已涌现100多种方法，主要包括基于空间加权starfm、基于空间解混stdaf、多种组合的方法)、改进glp时空融合方法，以及一系列衍生研究的方法，目前的时空研究主要集中在多光谱与多光谱图像之间(如modis与landsat、gf1 pms与gf1 wfv)的时空融合模拟，高光谱图像空间时间尺度提升的研究应用文献记录还比较少。
4.高光谱影像通常具有较低的空间分辨率，混合像元是不可忽视的问题，因此研究高光谱的融合，需考虑混合像元分解问题。结合加权、解混、插值思想的fsdaf(flexible spatial-temporal data fusion)模型利用权重函数将解混结果与空间插值相结合，能获得较好的鲁棒性，适用于异质性高、变化大的场景中。本发明基于fsdaf方法进行高光谱图像的时空融合预测，基于t0时刻的gf5 ashi和gf1 wfv图像，以及t1、t2时刻的gf1 wfv图像，融合得到t1、t2时刻的gf5 ashi模拟图像和分类结果，并与传统方法的融合和分类结果进行比较，解决了高光谱卫星数据获取困难、时间频次低等问题，为能充分发挥高光谱特色的时空动态预测研究提供了方法。
技术实现要素：5.针对上述现有技术中存在的不足，本发明的目的在于提供了一种面向高光谱图像和多光谱图像的时空谱融合方法，为解决高光谱卫星数据获取困难、时间频次低等问题，为充分发挥高光谱特色的时空动态预测研究提供了方法。
6.为解决上述问题，本发明公开了一种面向高光谱图像和多光谱图像的时空谱融合方法，在时空谱融合中，gf5 ashi提供光谱信息，gf1wfv提供时间和空间信息，将参考时刻和预测时刻分别定义为t1、t2，gf5 ahsi和gf1 wfv影像分别定义为hsi、msi；包括以下步骤：
7.s1.对t1时刻hsi图像的混合像元分解，获得端元类别和丰度图；
8.s2.估计每一类影像的从t1到t2的时间变化；
9.s3.假设t1到t2端元类别不发生变化，将时间变化残差加入到t1时刻的hsi图像上，得到t2时刻的hsi图像预测值；
10.s4.利用薄板样条插值法tps对t2时刻msi图像的每个波段进行插值，预测得到t2时刻的hsi图像，并将残差分配给预测的hsi图像；
11.s5.使用邻域信息得到t2时刻的hsi影像。
12.作为优选，s1的具体方法如下：
13.预定义(xi,yi)为每个像元的索引；i为msi的像元索引；j为每个hsi像元中对应msi图像的像元索引，j＝1...n；m1(xi,y
i,b
)与m2(xi,y
i,b
)分别表示t1和t2时刻，多光谱图像在b波段在(xi,yi)处的像元值；h1(x
ij
,y
ij,b
)与h2(x
ij
,y
ij,b
)分别表示t1和t2时刻，多光谱图像在b波段在(xi,yi)处的像元值；hc(xi,yi)表示高光谱图像上，端元类别c在(xi,yi)处的像元值；δm(xi,y
i,b
)表示t1和t2时间，多光谱图像在b波段在(xi,yi)像元处的变化值；δh(c
,b
)表示t1和t2时间，地物端元类别c在高光谱图像上的变化值；
14.高光谱图像的低空间分辨率特性导致混合像元的产生，首先对t1时刻的hsi图像混合像元分解，得到地物端元类别m和以及像元内每类地物的丰度hm(xi,yi)：
15.定义h(xi,yi)表示为gf5 ashi的第i个像元，n为h(xi,yi)中包含的端元个数，则第i个像元的每类丰度表示为：
[0016][0017]
式中，aj为地物cj在f(xi,yi)中的丰度，ε为模型误差，该误差可能是由带宽、太阳几何形状和观察角度的差异引起的；
[0018]
通过hsi中每个像元中所包含的类别数量，可计算出一个msi图像中的端元类别：
[0019]hm
(xi,yi)＝nc(xi，yi)/m
ꢀꢀꢀ
(2)
[0020]
nc(xi,yi)为(xi,yi)处的msi内属于c类的hsi图像像素数。
[0021]
作为优选，s2的具体方法如下：预测时间的变化过程中，假设地物类型没有变化,那么hsi中地物类别的时间变化可表示为δh(c，b),msi图像的时间变化为：
[0022]
δm(xi,yi,b)＝mn(xi,yi,b)-mm(xi,yi,b)
ꢀꢀ
(3)
[0023]
结合混合像元分解模型，又可表示为：
[0024]
[0025]
式中，选择n(n》l)个高光谱像元组成混合像元方程组对δh(class，b)求解，尽量选择每个类别中的高纯度像元。
[0026]
作为优选，s3的具体方法如下：
[0027]
预测地物类型发生变化的事件：假设t0-tp的土地类型不发生变化，将时间变化残差添加至t0多光谱图像上，得到tp时刻多光谱图像的时间预测值h
tp
；
[0028]hn-tp
(x
ij
,y
ij
,b)＝h1(x
ij
,y
ij
,b)+δh(class,b)
ꢀꢀ
(5)
[0029]
式中，
[0030]
每个msi的像素等于hsi像素及系统误差φ的和,t0、tp时刻的msi图像表示为：
[0031][0032][0033]
式中，每个msi图像像素包含m个hsi图像像素。
[0034]
作为优选，s4的具体方法如下：残差计算和分配，前面预测了时间变化，但未考虑土地类型变化较大的情况，因此需要引入残差r来表征类别差异；
[0035][0036][0037]
空间变化信息从tp时刻的msi图像中获得，对每个波段图像b使用薄板样条法进行插值：
[0038][0039]
其中，r
i2
＝(x-xi)2+(y-yi)2，n是已知像元个数，
[0040]
通过最小化来优化上式中的系数，
[0041][0042]
tp时刻的hsi图像空间预测(表示为下标sp)表示为：
[0043]hn-sp
(x
ij
,y
ij
,b)＝f
tps-b
(x
ij
,y
ij
)
ꢀꢀ
(12)
[0044]
如果图像的变化是同质的，那么预测的hsi图像空间要素与真是图像之间可认为没有误差，基于时间变化预测hsi图像的误差表示为：
[0045]eho
(x
ij
,y
ij
,b)＝h
n-sp
(x
ij
,y
ij
,b)-h
n-tp
(x
ij
,y
ij
,b)
ꢀꢀ
(13)
[0046]
如果图像地类是异质的，假设每个msi像元对应的hsi像元有相同的误差，则
[0047]ehe
(x
ij
,y
ij
,b)＝r(xi,yi,b)
ꢀꢀ
(14)
[0048]
引入hi作为表征同质程度的银子，若窗口内的第k个hsi像元与中心像元(x
ij
,y
ij
)的土地类型一致，则令ik＝1，反之为0；
[0049]ehe
(x
ij
,y
ij
,b)＝r(xi,yi,b)
ꢀꢀ
(15)
[0050]
将同质、异质误差结合，则综合权重weight可表示为：
[0051]
weight(x
ij
,y
ij
,b)＝e
ho
(x
ij
,y
ij
,b)
×
hi(x
ij
,y
ij
)+e
he
(x
ij
,y
ij
,b)
×
(1-hi(x
ij
,y
ij
))
ꢀꢀ
(16)
[0052]
将weight归一化得到weight，则分配到hsi像元上的残差为
[0053]
r(x
ij
,y
ij
,b)＝m
×
r(xi,yi,b)
×
weight(x
ij
,y
ij
,b)
ꢀꢀ
(17)
[0054]
hsi像元的全部变化表示为：
[0055][0056]
作为优选，s5的具体方法如下：最终预测，综合考虑邻近像元k与目标中心像元之间的空间距离、光谱相似度，得出权重wk，则最终hsi像元预测公式表达为：
[0057][0058]
作为优选，还包括以下步骤：
[0059]
s6.模型验证：通过分析建立了高光谱与多光谱影像的时空谱融合模型，对高分一号多光谱和高分五号高光谱数据进行了处理，并利用高分一号多光谱影像作为验证图像，对预测图像的信噪比、图像质量、信息熵进行验证。
[0060]
由于采用上述技术方案，本发明的有益效果如下：
[0061]
一种面向高光谱图像和多光谱图像的时空谱融合方法，基于多光谱高空间分辨率和高光谱图像的空间和光谱综合优势，结合混合像元分解模型进行数据融合，面向目前全球高光谱卫星少，数据难获取和全面覆盖等问题，计算出需求时刻的高光谱图像，解决各行业高光谱遥感数据应用需求，也是高光谱发展和时空融合的技术趋势，可解决环境变化监测、农业发展、地表类型变化等问题。
附图说明
[0062]
图1本发明的技术流程图；
[0063]
图2本发明的数据预处理流程图；
[0064]
图3高光谱图像和多光谱图像的像元图。
具体实施方式
[0065]
以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
[0066]
如图1所示，本发明提供一种面向高光谱图像和多光谱图像的时空谱融合方法，在时空谱融合中，gf5 ashi提供光谱信息，gf1wfv提供时间和空间信息，将参考时刻和预测时刻分别定义为t1、t2，gf5 ahsi和gf1 wfv影像分别定义为hsi、msi；基于fsdaf方法的高光谱时空预测模型修正，主要目标是获取hsi图像的高光谱信息，以及msi图像的高空间和时间尺度信息，包括以下步骤：
[0067]
s1.对t1时刻hsi图像的混合像元分解，获得端元类别和丰度图；
[0068]
混合像元：是构成遥感影像的基本单元，若其中包含多种覆盖类型，称为混合像元；
[0069]
端元：每个混合像元的基本组成单元；
[0070]
丰度：每个端元在其所在像元中所占比例，hsi图像的空间分辨率低，一个像元包含的地物类别较多，每种类别即为单元，每种类别在一个像元中所占的比例即为丰度；
[0071]
s2.估计每一类影像的从t1到t2的时间变化；
[0072]
s3.假设t1到t2端元类别不发生变化，将时间变化残差(各种端元类别的丰度比例变化实际值与观测值的残差)加入到t1时刻的hsi图像上，得到t2时刻的hsi图像预测值；
[0073]
s4.利用薄板样条插值法(tps：是一种常用的2d空间插值方法，经常用于图像变形等，通过少量的控制点就可以驱动图像进行变化。)对t2时刻msi图像的每个波段进行插值，预测得到t2时刻的hsi图像，并将残差分配给预测的hsi图像；
[0074]
s5.使用邻域信息得到t2时刻的hsi影像；
[0075]
s6.模型验证：通过分析建立了高光谱与多光谱影像的时空谱融合模型，对高分一号多光谱和高分五号高光谱数据进行了处理，并利用高分一号多光谱影像作为验证图像，对预测图像的信噪比、图像质量、信息熵等进行验证。
[0076]
进一步地，s1的具体方法如下：
[0077]
预定义(xi,yi)为每个像元的索引；i为msi的像元索引；j为每个hsi像元中对应msi图像的像元索引，j＝1...n；m1(xi,yi,b)与m2(xi,yi,b)分别表示t1和t2时刻，多光谱图像在b波段在(xi,yi)处的像元值；h1(x
ij
,y
ij
,b)与h2(x
ij
,y
ij
,b)分别表示t1和t2时刻，多光谱图像在b波段在(xi,yi)处的像元值；hc(xi,yi)表示高光谱图像上，端元类别c在(xi,yi)处的像元值；δm(xi,yi,b)表示t1和t2时间，多光谱图像在b波段在(xi,yi)像元处的变化值；δh(c
,b
)表示t1和t2时间，地物端元类别c在高光谱图像上的变化值；
[0078]
高光谱图像的低空间分辨率特性导致混合像元的产生，首先对t1时刻的hsi图像混合像元分解，得到地物端元类别m和以及像元内每类地物的丰度hm(xi,yi)：
[0079]
如图3所示，左-gf5 ashi一个像元，右-gf1 wfv对应混合像元，定义h(xi,yi)表示为gf5 ashi的第i个像元，n为h(xi,yi)中包含的端元个数，则第i个像元的每类丰度表示为：
[0080][0081]
式中，aj为地物cj在f(xi,yi)中的丰度，ε为模型误差，该误差可能是由带宽、太阳几何形状和观察角度的差异引起的。
[0082]
通过hsi中每个像元中所包含的类别数量，可计算出一个msi图像中的端元类别：
[0083]hm
(xi,yi)＝nc(xi，yi)/m
ꢀꢀꢀ
(2)
[0084]
nc(xi,yi)为(xi,yi)处的msi内属于c类的hsi图像像素数，m为地物端元类别。
[0085]
s2的具体方法如下：预测时间的变化过程中，假设地物类型没有变化。
[0086]
那么hsi中地物类别的时间变化可表示为δh(c，b)。msi图像的时间变化为：
[0087]
δm(xi,yi,b)＝mn(xi,yi,b)-mm(xi,yi,b)
ꢀꢀ
(3)
[0088]
结合混合像元分解模型，又可表示为：
[0089][0090]
式中，选择n(n》l)个高光谱像元组成混合像元方程组对δh(class，b)求解，尽量选择每个类别中的高纯度像元。
[0091]
s3的具体方法如下：空间预测(预测地物类型发生变化的事件)：假设t0-tp的土地类型不发生变化，将时间变化残差添加至t0多光谱图像上，得到tp时刻多光谱图像的时间预测值h
tp
。
[0092]hn-tp
(x
ij
,y
ij
,b)＝h1(x
ij
,y
ij
,b)+δh(class,b)
ꢀꢀ
(5)
[0093]
式中，
[0094]
每个msi的像素等于hsi像素及系统误差φ的和。t0、tp时刻的msi图像表示为：
[0095][0096][0097]
式中，每个msi图像像素包含m个hsi图像像素；
[0098]
s4的具体方法如下：残差计算和分配
[0099]
前面预测了时间变化，但未考虑土地类型变化较大的情况，因此需要引入残差r来表征类别差异。
[0100][0101][0102]
空间变化信息从tp时刻的msi图像中获得，对每个波段图像b使用薄板样条法进行插值：
[0103][0104]
其中，r
i2
＝(x-xi)2+(y-yi)2，n是已知像元个数。
[0105]
通过最小化来优化上式中的系数。
[0106][0107]
tp时刻的hsi图像空间预测(表示为下标sp)表示为：
[0108]hn-sp
(x
ij
,y
ij
,b)＝f
tps-b
(x
ij
,y
ij
)
ꢀꢀ
(12)
[0109]
如果图像的变化是同质的，那么预测的hsi图像空间要素与真是图像之间可认为没有误差。基于时间变化预测hsi图像的误差表示为：
[0110]eho
(x
ij
,y
ij
,b)＝h
n-sp
(x
ij
,y
ij
,b)-h
n-tp
(x
ij
,y
ij
,b)
ꢀꢀ
(13)
[0111]
如果图像地类是异质的，假设每个msi像元对应的hsi像元有相同的误差，则
[0112]ehe
(x
ij
,y
ij
,b)＝r(xi,yi,b)
ꢀꢀ
(14)
[0113]
引入hi作为表征同质程度的银子，若窗口内的第k个hsi像元与中心像元(x
ij
,y
ij
)的土地类型一致，则令ik＝1，反之为0，
[0114]ehe
(x
ij
,y
ij
,b)＝r(xi,yi,b)
ꢀꢀ
(15)
[0115]
将同质、异质误差结合，则综合权重weight可表示为：
[0116]
weight(x
ij
,y
ij
,b)＝e
ho
(x
ij
,y
ij
,b)
×
hi(x
ij
,y
ij
)+e
he
(x
ij
,y
ij
,b)
×
(1-hi(x
ij
,y
ij
))
ꢀꢀ
(16)
[0117]
将weight归一化得到weight。则分配到hsi像元上的残差为
[0118]
r(x
ij
,y
ij
,b)＝m
×
r(xi,yi,b)
×
weight(x
ij
,y
ij
,b)
ꢀꢀ
(17)
[0119]
hsi像元的全部变化表示为：
[0120][0121]
s5的具体方法如下：最终预测。综合考虑邻近像元k与目标中心像元之间的空间距离、光谱相似度，得出权重wk，则最终hsi像元预测公式表达为：
[0122][0123]
以上所述仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
技术特征：1.一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，在时空谱融合中，gf5 ashi提供光谱信息，gf1wfv提供时间和空间信息，将参考时刻和预测时刻分别定义为t1、t2，gf5 ahsi和gf1 wfv影像分别定义为hsi、msi；包括以下步骤：s1.对t1时刻hsi图像的混合像元分解，获得端元类别和丰度图；s2.估计每一类影像的从t1到t2的时间变化；s3.假设t1到t2端元类别不发生变化，将时间变化残差加入到t1时刻的hsi图像上，得到t2时刻的hsi图像预测值；s4.利用薄板样条插值法tps对t2时刻msi图像的每个波段进行插值，预测得到t2时刻的hsi图像，并将残差分配给预测的hsi图像；s5.使用邻域信息得到t2时刻的hsi影像。2.根据权利要求1所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，s1的具体方法如下：预定义(x
i
,y
i
)为每个像元的索引；i为msi的像元索引；j为每个hsi像元中对应msi图像的像元索引，j＝1...n；m1(x
i
,y
i,b
)与m2(x
i
,y
i,b
)分别表示t1和t2时刻，多光谱图像在b波段在(x
i
,y
i
)处的像元值；h1(x
ij
,y
ij,b
)与h2(x
ij
,y
ij,b
)分别表示t1和t2时刻，多光谱图像在b波段在(x
i
,y
i
)处的像元值；h
c
(x
i
,y
i
)表示高光谱图像上，端元类别c在(x
i
,y
i
)处的像元值；δm(x
i
,y
i,b
)表示t1和t2时间，多光谱图像在b波段在(x
i
,y
i
)像元处的变化值；δh(c
,b
)表示t1和t2时间，地物端元类别c在高光谱图像上的变化值；高光谱图像的低空间分辨率特性导致混合像元的产生，首先对t1时刻的hsi图像混合像元分解，得到地物端元类别m和以及像元内每类地物的丰度h
m
(x
i
,y
i
)：定义h(x
i
,y
i
)表示为gf5 ashi的第i个像元，n为h(x
i
,y
i
)中包含的端元个数，则第i个像元的每类丰度表示为：式中，a
j
为地物c
j
在f(x
i
,y
i
)中的丰度，ε为模型误差，该误差可能是由带宽、太阳几何形状和观察角度的差异引起的；通过hsi中每个像元中所包含的类别数量，可计算出一个msi图像中的端元类别：h
m
(x
i
,y
i
)＝n
c
(x
i
，y
i
)/m
ꢀꢀꢀꢀ
(2)n
c
(x
i
,y
i
)为(x
i
,y
i
)处的msi内属于c类的hsi图像像素数。3.根据权利要求2所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，s2的具体方法如下：预测时间的变化过程中，假设地物类型没有变化,那么hsi中地物类别的时间变化可表示为δh(c，b),msi图像的时间变化为：δm(x
i
,y
i
,b)＝m
n
(x
i
,y
i
,b)-m
m
(x
i
,y
i
,b)
ꢀꢀꢀꢀ
(3)结合混合像元分解模型，又可表示为：式中，选择n(n＞l)个高光谱像元组成混合像元方程组对δh(class，b)求解，尽量选择每个类别中的高纯度像元。
4.根据权利要求3所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，s3的具体方法如下：预测地物类型发生变化的事件：假设t0-tp的土地类型不发生变化，将时间变化残差添加至t0多光谱图像上，得到tp时刻多光谱图像的时间预测值h
tp
；h
n-tp
(x
ij
,y
ij
,b)＝h1(x
ij
,y
ij
,b)+δh(class,b)
ꢀꢀꢀꢀ
(5)式中，每个msi的像素等于hsi像素及系统误差φ的和,t0、tp时刻的msi图像表示为：每个msi的像素等于hsi像素及系统误差φ的和,t0、tp时刻的msi图像表示为：式中，每个msi图像像素包含m个hsi图像像素。5.根据权利要求4所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，s4的具体方法如下：残差计算和分配，前面预测了时间变化，但未考虑土地类型变化较大的情况，因此需要引入残差r来表征类别差异；空间变化信息从tp时刻的msi图像中获得，对每个波段图像b使用薄板样条法进行插值：其中，r
i2
＝(x-x
i
)2+(y-y
i
)2，n是已知像元个数，通过最小化来优化上式中的系数，tp时刻的hsi图像空间预测(表示为下标sp)表示为：h
n-sp
(x
ij
,y
ij
,b)＝f
tps-b
(x
ij
,y
ij
)
ꢀꢀꢀꢀ
(12)如果图像的变化是同质的，那么预测的hsi图像空间要素与真是图像之间可认为没有误差，基于时间变化预测hsi图像的误差表示为：e
ho
(x
ij
,y
ij
,b)＝h
n-sp
(x
ij
,y
ij
,b)-h
n-tp
(x
ij
,y
ij
,b)
ꢀꢀꢀꢀ
(13)如果图像地类是异质的，假设每个msi像元对应的hsi像元有相同的误差，则e
he
(x
ij
,y
ij
,b)＝r(x
i
,y
i
,b)
ꢀꢀꢀꢀ
(14)
引入hi作为表征同质程度的银子，若窗口内的第k个hsi像元与中心像元(x
ij
,y
ij
)的土地类型一致，则令i
k
＝1，反之为0；e
he
(x
ij
,y
ij
,b)＝r(x
i
,y
i
,b)
ꢀꢀꢀꢀ
(15)将同质、异质误差结合，则综合权重weight可表示为：weight(x
ij
,y
ij
,b)＝e
ho
(x
ij
,y
ij
,b)
×
hi(x
ij
,y
ij
)+e
he
(x
ij
,y
ij
,b)
×
(1-hi(x
ij
,y
ij
))
ꢀꢀꢀꢀ
(16)将weight归一化得到weight，则分配到hsi像元上的残差为r(x
ij
,y
ij
,b)＝m
×
r(x
i
,y
i
,b)
×
weight(x
ij
,y
ij
,b)
ꢀꢀꢀꢀ
(17)hsi像元的全部变化表示为：6.根据权利要求5所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，s5的具体方法如下：最终预测，综合考虑邻近像元k与目标中心像元之间的空间距离、光谱相似度，得出权重wk，则最终hsi像元预测公式表达为：7.根据权利要求1所述的一种面向高光谱图像和多光谱图像的时空谱融合方法，其特征在于，还包括以下步骤：s6.模型验证：通过分析建立了高光谱与多光谱影像的时空谱融合模型，对高分一号多光谱和高分五号高光谱数据进行了处理，并利用高分一号多光谱影像作为验证图像，对预测图像的信噪比、图像质量、信息熵进行验证。
技术总结本发明涉及时空谱遥感数据融合技术领域,尤其是一种面向高光谱图像和多光谱图像的时空谱融合方法，通过对高分一号多光谱(GF1WFV)图像和高分五号高光谱(GF5AHSI)图像的预处理，并对光谱图像空间分辨率和高光谱图像光谱分辨率的分析，将参考时刻的高光谱分辨率和预测时刻的多光谱高空间分辨率进行保留，构建时空谱综合预测模型，得到集合新时相的高空间高光谱图像，并利用预测时刻的多光谱图像进行验证。本发明的有益效果是在某时刻高光谱图像难以获取时，可用过去时刻的多光谱图像和高光谱图像进行预测，得到一幅能够满足需求时刻高空间分辨率的高光谱图像，解决数据源的问题，从而服务于众多高光谱应用领域。而服务于众多高光谱应用领域。而服务于众多高光谱应用领域。
技术研发人员：李雁 李卫海 甘淑 袁希平 张泉 王艳英 王琳 胡琳
受保护的技术使用者：云南省军民融合发展研究院（云南省国防科技工业局综合研究所、云南省卫星遥感技术应用工程中心）
技术研发日：2022.06.22
技术公布日：2022/11/1
