%程序用于平滑处理数据,以及将间隔时间数据插值成每天
clc
clear
%光学卫星获取的数据存在受云影响的情况,通过滤波的方式将极低或极高等不合理值进行处理(以MODIS中的LAI为例)
[a,b]=xlsread('F:\LAI_4Day.xlsx','Sheet2');
LAI=a(:,7);%提取LAI列
LAI_smooth_gauss = smoothdata(LAI,'gaussian',[6 6]);%数据平滑,高斯加权移动平均值
LAI_smooth_movmean = smoothdata(LAI,'movmean',[2 8]);%数据平滑,移动平均值
%注:根据实际需要和实际数据调整平滑窗口大小和平滑方法
%%针对卫星的周期性观测,进行插值,得到每天的数据
x_1=1:4:897;%此为卫星每隔四天获取的周期性观测日期
x_2=1:1:897;%将四天插补到逐日
LAI_i=interp1(x_1,LAI_smooth_gauss,x_2, 'spline');%选择三次样条插值方法
转载本文请联系原作者获取授权,同时请注明本文来自李雷科学网博客。
链接地址:/blog-3436537-1241889.html
下一篇:arcgis模型构建器 批量处理矢量点位和矢量面