【技术实现步骤摘要】
本专利技术涉及遥感，具体涉及针对山地遥感植被指数的时序滤波算法、装置及存储介质。
技术介绍
1、山地地形高差大，天气变化快，容易受到云雨天气的干扰。同时，在非生长季时期，高海拔地区气温较低，常出现积雪覆盖。这就导致遥感影像的地表反射率受到影响，植被指数数据集在空间上的出现不连续，在云、雪覆盖区域出现异常值，对时间序列植被指数分析带来了较大的困难。
2、这类植被指数数据有一个共同点，即植被的季节效应，本专利技术意在研究一种算法，以消除云雪对于时间序列植被指数数据的干扰，从而获得高质量的空间连续数据。
技术实现思路
1、针对现有技术存在的不足，本专利技术提出针对山地遥感植被指数的时序滤波算法、装置及存储介质，用于消除云雪对于时间序列植被指数数据的干扰，从而获得高质量的空间连续数据。
2、第一方面，提供了针对山地遥感植被指数的时序滤波算法，包括：
3、获取植被指数数据集；
4、根据场景质量分类图栅格文件，获取掩膜栅格文件；
5、利用掩膜栅格文件标记无效像元，并对植被指数数据集进行计算，获取时序滤波重构结果。
6、进一步的，获取植被指数数据集，包括：
7、获取sentinel-2l2a级卫星数据；
8、对sentinel-2l2a级卫星数据采用植被指数计算公式进行计算，获取植被指数或其它衍生数据；
9、对植被指数或其它衍生数据进行波段叠加，获取植被指数数据集。
10、进一
11、根据sentinel-2l2a级卫星数据，获取场景质量分类波段；
12、根据场景质量分类波段，获取云、雪、云阴影与其他像元区分结果；
13、对云、雪、云阴影与其他像元区分结果采用重分类，获取掩膜栅格文件。
14、进一步的，利用掩膜栅格文件标记无效像元，并对植被指数数据集进行计算，获取时序滤波重构结果，包括：
15、根据掩膜栅格文件，分别标记数据集中每一波段的云雪及云阴影像元；
16、对标记的云雪及云阴影像元进行计算，获取无效像元的模拟植被指数值并进行替换；
17、对替换无效像元后的植被指数数据集采用sg滤波算法迭代，获取时序滤波重构结果。
18、进一步的，对标记的云雪及云阴影像元进行计算，获取无效像元的模拟植被指数值，包括：
19、通过计算
20、
21、获取无效像元的模拟植被指数值，其中，vi0(n,doy)为序号为doy的影像第n个像元计算后的模拟植被指数值，df为像元对应的掩膜编码值，doy2与doy1分别为向后和向前搜索到的相邻有效像元对应时相序号。
22、进一步的，对替换无效像元后的植被指数数据集采用sg滤波算法迭代，获取时序滤波重构结果，包括：
23、通过计算mse并与阈值a进行对比，当mse相对变化率大于a时，停止迭代获取时序滤波重构结果，
24、其中，mse相对变化率表示均方误差的相对变化率，msei表示第i次sg滤波后的均方误差，msei-1表示第i-1次sg滤波后的均方误差。
25、第二方面，提供了用于针对山地遥感植被指数的时序滤波算法的装置，包括：
26、植被指数数据集模块，被配置为获取植被指数数据集；
27、掩膜栅格文件模块，被配置为根据场景质量分类图栅格文件，获取掩膜栅格文件；
28、时序滤波重构结果模块，被配置为利用掩膜栅格文件标记无效像元，并对植被指数数据集进行计算，获取时序滤波重构结果。
29、第三方面，提供了用于针对山地遥感植被指数的时序滤波算法的装置，包括处理器和存储有程序指令的存储器，其特征在于，所述处理器被配置为在运行所述程序指令时，执行如前文任意一项所述的用于针对山地遥感植被指数的时序滤波算法。
30、第四方面，提供了一种存储介质，存储有程序指令，其特征在于，所述程序指令在运行时，执行如前文所述的用于针对山地遥感植被指数的时序滤波算法。
31、采用上述技术方案的专利技术，具有如下优点：
32、本专利技术基于sentinel-2卫星数据，结合线性内插与s-g滤波算法，建立一种植被指数数据集时序重构算法，该算法能够显著改善长时间序列山地植被指数数据集受云雪及云阴影影响情况，在保证数据真实性的基础上根据时序数据逐像元计算合适的模拟值，使数据集具有良好的空间连续性，显著提升数据集的质量与可用性，以精细化表征山地植被随时间变化趋势。
本文档来自技高网...【技术保护点】
1.针对山地遥感植被指数的时序滤波算法，其特征在于，包括：
2.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，获取植被指数数据集，包括：
3.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，根据场景质量分类图栅格文件，获取掩膜栅格文件，包括：
4.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，利用掩膜栅格文件标记无效像元，并对植被指数数据集进行计算，获取时序滤波重构结果，包括：
5.根据权利要求4所述的针对山地遥感植被指数的时序滤波算法，其特征在于，对标记的云雪及云阴影像元进行计算，获取无效像元的模拟植被指数值，包括：
6.根据权利要求4所述的针对山地遥感植被指数的时序滤波算法，其特征在于，对替换无效像元后的植被指数数据集采用SG滤波算法迭代，获取时序滤波重构结果，包括：
7.用于针对山地遥感植被指数的时序滤波算法的装置，其特征在于，包括：
8.用于针对山地遥感植被指数的时序滤波算法的装置，包括处理器和存储有程序指令的存储器，其特征在
9.一种存储介质，存储有程序指令，其特征在于，所述程序指令在运行时，执行如权利要求1至6任一项所述的用于针对山地遥感植被指数的时序滤波算法。
...【技术特征摘要】
1.针对山地遥感植被指数的时序滤波算法，其特征在于，包括：
2.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，获取植被指数数据集，包括：
3.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，根据场景质量分类图栅格文件，获取掩膜栅格文件，包括：
4.根据权利要求1所述的针对山地遥感植被指数的时序滤波算法，其特征在于，利用掩膜栅格文件标记无效像元，并对植被指数数据集进行计算，获取时序滤波重构结果，包括：
5.根据权利要求4所述的针对山地遥感植被指数的时序滤波算法，其特征在于，对标记的云雪及云阴影像元进行计算，获取无效像元的模拟植被指数值，包...
【专利技术属性】
技术研发人员：黄子洋，黄河，谭玲，刘中帅，何晓兵，张传霆，彭东明，龙春宇，
申请(专利权)人：招商局重庆交通科研设计院有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。