[发明专利]陆面蒸散发遥感估算的时空尺度匹配方法及装置

摘要

本发明提供了一种陆面蒸散发遥感估算的时空尺度匹配方法及装置，方法包括：S1获取多空间、多时间分辨率的遥感数据集；S2根据获取的遥感数据集进行时空数据融合，得到同时具有高时间和高空间分辨率的遥感数据集；S3利用预设数据模型将S1和S2生成的数据集进行ET参数反演，得到不同时间和空间分辨率的遥感参数集；根据S3获取的不同时间和空间分辨率的遥感参数，分别进行空间及时间尺度匹配判断；若空间尺度匹配判断结果和时间尺度匹配判断结果满足ET估算的精度要求，则将步骤S4中满足ET估算精度要求的遥感数据作为匹配的数据集。本发明能够解决现有ET估算中时间、空间尺度匹配不明确而导致的ET反演不确定性的问题。

主权项

1.一种陆面蒸散发遥感估算的时空尺度匹配方法，其特征在于，包括如下步骤：S1、利用传感器获取多空间、多时间分辨率的遥感数据集；S2、根据步骤S1获取的遥感数据集进行时空数据融合，得到同时具有高时间和高空间分辨率的遥感数据集；S3、利用预设数据模型将步骤S1和S2生成的数据集进行ET参数反演，得到ET参数集，得到不同时间和空间分辨率的遥感参数集；S4、根据步骤S3获取的不同时间和空间分辨率的遥感参数集，分别进行空间尺度匹配判断以及时间尺度匹配判断；S5、若空间尺度匹配判断结果和时间尺度匹配判断结果满足ET估算的精度要求，则将步骤S4中满足ET估算精度要求的遥感数据作为匹配的数据集；其中，根据步骤S3获取的不同时间和空间分辨率的遥感参数集，进行空间尺度匹配判断，具体包括：S41、获取不同分辨率遥感影像生成的土地利用图，将土地利用图矢量化，将相同分辨率遥感图像得到的矢量化的土地利用图与栅格的遥感土地利用图叠加，计算空间松散度指数：SP＝N_i/S_i；式中，SP为空间松散度指数；N_i为预设空间分辨率下生成的土地利用斑块数；S_i为上述土地利用的面积；S42、计算下垫面异质性指数：式中，A为下垫面异质性指数，A_a为被选择的遥感数据像元总数；M为混合像元判定值，当为混合像元时，值为1，否则为0；S_m为对应像元相似度系数，为混合像元中所包含的地物类型的种类数；S43、当SP≥a，同时A≤b时，判断为空间尺度高度匹配；当只满足SP≥a，或只满足A≤b时，判断为空间尺度中度匹配；当二者均不满足时，判断为空间尺度不匹配；式中，a和b分别为对应的预设阈值；其中，根据步骤S3获取的不同时间和空间分辨率的遥感参数集，进行时间尺度匹配判断，具体包括：S41’、获取在t‑t₀区间不同蒸散发参量的生育期曲线数据，并判断所述曲线数据与卫星实测的曲线之间的相关性是否满足：R²≥c，RMSE≤d，P≤0.01；式中，c，d为预设的阈值；t‑t₀为时间区间；R为相关系数；RMSE为均方根误差；P值为显著性水平判别值；S42’、按照植被生长期或作物生育期进行分段，在此分段的基础上，计算不同的数据时间段之间可获取的参数稳定性指标VP_z：VP_z＝VAR_tb‑tc(P_ztb,P_ztb+1,P_ztb+2,......,P_ztc)，并判断VP_z取值是否趋近于0；其中，P_z为参数值，tb，…tc为不同的时间点数值，VAR_tb‑tc为时间区间tb到tc的方差计算模型，VP_z为tb到tc的方差值；S43’、若S41’和S42’中有一项满足条件，则判断为时间尺度匹配，满足ET估算的精度要求；其中，若空间尺度匹配判断结果和时间尺度匹配判断结果满足ET估算的精度要求，则将步骤S4中满足ET估算精度要求的遥感数据作为匹配的数据集，具体包括：若空间尺度匹配结果为高度匹配，时间尺度匹配结果为匹配，则认为满足ET估算精度要求，此时将步骤S4中满足ET估算精度要求的遥感数据作为匹配的数据集；若空间尺度匹配结果为中度匹配，时间尺度匹配结果为匹配，则将空间匹配以及时间匹配的参数集输入到ET模型中，若得到模拟的ET估算结果与实际观测的ET值关系满足条件：R²≥e，RMSE≤f，P≤g时，则认为满足ET估算精度要求，此时将步骤S4中满足ET估算精度要求的遥感数据作为匹配的数据集；式中，e，f，g为对应的预设阈值。