matlab变异系数CV代码

clear
[a,R]=readgeoraster('F:\LIN\npp\MOD17A3NPP\2000_npp.tif');  
info=geotiffinfo('F:\LIN\npp\MOD17A3NPP\2000_npp.tif');
[m,n]=size(a);
datasum=zeros(m*n,21)+NaN; 
for year=2000:2020 
    filename=['F:\LIN\npp\MOD17A3NPP\',int2str(year),'_npp.tif'];
    data=importdata(filename);
    data=reshape(data,m*n,1);
    datasum(:,year-1999)=data;
end

CV=zeros(1,m*n)+nan;
for i=1:length(datasum)
    a=datasum(i,:);
    if min(a)>=0
        mean0=mean(a);
        b=0;
        for j=1:21
            c=(a(j)-mean0)^2;
            b=b+c;
            
        end
        b=sqrt(b/20);
        CV(i)=b/mean0;  
    end
end
CV=reshape(CV,m,n);
geotiffwrite('F:\LIN\AA_figure\figure代码\npp变异系数.tif',CV,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
disp('OK!')