在應(yīng)用年內(nèi)的NDVI信息的時(shí)候，由于云的干擾，通常需要對(duì)NDVI產(chǎn)品如MOD13A2進(jìn)行濾波處理進(jìn)行重建，本文基于matlab平臺(tái)，結(jié)合SG濾波和質(zhì)量控制文件對(duì)MOD13A2產(chǎn)品進(jìn)行重新構(gòu)建，基本思路是首先對(duì)23期數(shù)據(jù)進(jìn)行濾波處理，然后質(zhì)量好的像元保持不變，而質(zhì)量較差的像元為濾波后的值，構(gòu)成新的NDVI序列?；敬a如下

% @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的質(zhì)量文件
        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);%進(jìn)行還原
            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\基于質(zhì)量文件和SG濾波插補(bǔ)\MOD13A2_',int2str(year),'_',int2str(k),'.tif'];
        geotiffwrite(filename,ndvi,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag)
    end
end

在用上述代碼的時(shí)候請(qǐng)引用以下文獻(xiàn):
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.
更多需求，請(qǐng)查看個(gè)人介紹