参考链接：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 

来源：CSDN 

原文：https://blog.csdn.net/yangjh1991/article/details/69788778 

版权声明：本文为博主原创文章，转载请附上博文链接！