参考链接:https://blog.csdn.net/yangjh1991/article/details/69788778
Lon = ncread(InFile,'lon'); %读取经度数据
Lat = ncread(InFile,'lat');
SPEI = ncread(InFile,'spei',[1,1,7],[720,360,1]); %读取1900年7月份的SPEI数据
%地理参考系的建立
GeoRef = georasterref('Rastersize',[360,720],'Latlim',[-90,90],'Lonlim',[-180,180]);
%数据的写出
SPEI_Tif = strcat(InPath,'190007SPEI01.tif');
geotiffwrite(SPEI_Tif,flip(rot90(SPEI)),GeoRef)
disp('finished')
% georasterref函数
% georasterref函数用于生成tif文件的坐标系统,包括三个参数'rastersize','Latlim'和'Lonlim',在形成地理坐标系时,可以在前面先用ncread读取相应变量,然后通过size,min和max等函数实现坐标系统的产生,例如:4中代码的“GeoRef = georasterref('Rastersize',[360,720],'Latlim',[-90,90],'Lonlim',[-180,180])”可以改写为“
% GeoRef = georasterref('Rastersize',size(SPEI),'Latlim',[double(min(Lat)),double(max(Lat))],'Lonlim',[double(min(Lon)),double(min(Lon))]),增加程序的灵活性。
---------------------
作者:木易1991