%程序用于平滑处理数据,以及将间隔时间数据插值成每天

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');%选择三次样条插值方法