Matlab处理SPEI全球数据集

SPEI数据集下载网址：https://digital.csic.es/handle/10261/202305

SPEI定义及计算方法参见我的《R语言spei包计算标准化降水蒸散指数》

目前该数据时间尺度为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个月数据
        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