参考链接: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
版权声明:本文为博主原创文章,转载请附上博文链接!
转载本文请联系原作者获取授权,同时请注明本文来自张乐乐科学网博客。
链接地址:http://blog.sciencenet.cn/blog-1103122-1155424.html
上一篇:MATLAB创建GIMMS NDVI批量下载链接
下一篇:MATLAB将日期显示为数值