在應用年內的NDVI信息的時候，由于云的干擾焰薄，通常需要對NDVI產品如MOD13A2進行濾波處理進行重建拿诸，本文基于matlab平臺，結合SG濾波和質量控制文件對MOD13A2產品進行重新構建塞茅，基本思路是首先對23期數(shù)據(jù)進行濾波處理亩码，然后質量好的像元保持不變，而質量較差的像元為濾波后的值凡桥，構成新的NDVI序列蟀伸∈赐基本代碼如下

% @author yinlichang3064@163.com
[a,R]=geotiffread('I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_2017353.hdfout.1_km_16_days_NDVI.tif');
info=geotiffinfo('I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_2017353.hdfout.1_km_16_days_NDVI.tif');
[m,n]=size(a);
begin_year=2001;
end_year=2018;
qs=23;
%總共23期數(shù)據(jù)
for year=begin_year:end_year
    ndvisum=zeros(m*n,qs)+NaN;
    qcsum=zeros(m*n,qs)+NaN;
    k1=1;
    for k=1:16:366 
        if k<10
              filename=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),'00',int2str(k),'.hdfout.1_km_16_days_NDVI.tif'];
              filename_pr=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),'00',int2str(k),'.hdfout.1_km_16_days_pixel_reliability.tif'];
        else
           if k<100
              filename=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),'0',int2str(k),'.hdfout.1_km_16_days_NDVI.tif']; 
              filename_pr=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),'0',int2str(k),'.hdfout.1_km_16_days_pixel_reliability.tif']; 
            else
              filename=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),int2str(k),'.hdfout.1_km_16_days_NDVI.tif'];  
              filename_pr=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\NDVI_PR_',int2str(year),int2str(k),'.hdfout.1_km_16_days_pixel_reliability.tif'];  
            end 
        end
        ndvi=double(importdata(filename));%獲取ndvi
        ndvi=reshape(ndvi,m*n,1);
        ndvi(ndvi==-3000)=NaN;
        qc=importdata(filename_pr);%獲取ndvi的質量文件
        qc=reshape(qc,m*n,1);
        ndvisum(:,k1)=ndvi;
        qcsum(:,k1)=qc; 
        k1=k1+1;
    end
    for i=1:m*n
        data1=ndvisum(i,:);
        if min(data1)>=-2000
            qc=qcsum(i,:);
            data_sg=sgolayfilt(data1,3,5);
            data_sg(qc==0)=data1(qc==0);%進行還原
            ndvisum(i,:)=data_sg;
        end
    end
    for k=1:qs
        ndvi=ndvisum(:,k);
        ndvi=reshape(ndvi,m,n);
        ndvi(ndvi<-2000)=NaN;
        filename=['I:\MOD13A2\MOD13A2_yr_03-17_16DAY\基于質量文件和SG濾波插補\MOD13A2_',int2str(year),'_',int2str(k),'.tif'];
        geotiffwrite(filename,ndvi,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag)
    end
end

在用上述代碼的時候請引用以下文獻:
Yin, L.; Wang, X.; Feng, X.; Fu, B.; Chen, Y. A Comparison of SSEBop-Model-Based Evapotranspiration with Eight Evapotranspiration Products in the Yellow River Basin, China. Remote Sens. 2020, 12, 2528.
更多需求缅刽，請查看個人介紹