matlab圆形数据集,Matlab处理SPEI全球数据集

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