SPEI數(shù)據集下載網址：https://digital.csic.es/handle/10261/202305
SPEI定義及計算方法參見我的簡書《R語言spei包計算標準化降水蒸散指數(shù)》
目前該數(shù)據時間尺度為1901-2018，空間分辨率0.5°羔挡，其中計算蒸散的方法采用彭曼公式
代碼
clc;
clear;
ncdisp('I:\data\Climate--nc文件\SPEI\spei12.nc')
data=ncread('I:\data\Climate--nc文件\SPEI\spei12.nc','spei');
for year=1901:2018
data1=data(:,:,1+12*(year-1901):12*(year-1900)); %得到每年的12個月數(shù)據
data3=sum(data1,3)/12;
data4=rot90(data3);
data_yearly=flipud(data4);
filename=strcat('I:\data\SPEI\yearly\spei12_',num2str(year),'.tif');
georef=georasterref('RasterSize',size(data_yearly),'Latlim',[-90 90],'Lonlim',[0 360]);
geotiffwrite(filename,data_yearly,georef);
for month=1:12
data2=data1(:,:,month);
data4=rot90(data2);
data_monthly=flipud(data4);
filename_mon = strcat('I:\data\SPEI\monthly\',num2str(year),'\',num2str(year),num2str(month,'%02d'),'.tif');
georef=georasterref('RasterSize',size(data_monthly),'Latlim',[-90 90],'Lonlim',[0 360]);
geotiffwrite(filename_mon,data_monthly,georef);
end
end
