gleam數(shù)據(jù)3.6a

鏈接：https://pan.baidu.com/s/1xrOXYby6sYJMAEi6bltHXw?pwd=gnwd?

提取碼：gnwd?

nc to tif mtlab代碼

%%

clc

clear

%% 讀取文件信息

ncFile = 'SMsurf_1980-2021_GLEAM_v3.6a_MO.nc'; %設(shè)定文件

ncdisp(ncFile)

outpath='E:\DATA\MWSNDVI\gleam\MONTH\E\'

%% 讀取數(shù)據(jù)

time=1980;

y=1;

lat=ncread(ncFile,'lat');

lon=ncread(ncFile,'lon');

E=ncread(ncFile,'E'); % Actual evaporation from GLEAM v3.6a

for i=1:1:504

? ?if y==13

? ? ? ?time=time+1;

? ? ? ?y=1;

? ?end

outname=[int2str(time),int2str(y),'gleam'];

out=E(:,:,i);

out=rot90(out);

%out = fliplr(out);

R=georasterref('RasterSize', size(out),'Latlim', [-89.875 89.875], 'Lonlim', [-179.875 179.875]);

geotiffwrite([outpath,outname], out, R);

? ?y=y+1;

? ?disp(outname)

end

?代碼可以根據(jù)實(shí)際情況修改，這是以E作為例子；

如果要循環(huán)文件夾的話，可以把變量名稱從 文件名稱中取 ；在將文件夾寫入循環(huán) 嵌套三個(gè)循環(huán)；命名循環(huán)比較簡(jiǎn)單；