5_批量读取nc并计算月均LAI_批量下载NC号_

clc; %清屏 clear; %清空 datadir='G:\zjj\2007\5\'; %指定批量数据所在的文件夹 filelist=dir([datadir,'*.nc']); %指定批量数据的类型 k=length(filelist); b=0; for s=1:k filename=[datadir,filelist(s).name]; M=ncread(filename,'LAI'); M=M'; A(:,:,s)=sum(M,3); A(find(isnan(A)==1))=0; b=b+A(:,:,s); C = GetLccMap(0.05); %0.05是该项目数据的度数，GetLccMap就是用来计算mask的 D =(C==13|C==15|C==16|C==17); C(D)=0; mask=C==0; b(mask)=0; end b=b./5; q=imagesc(b); set(q,'alphadata',b~=0)%将背景值变成白色