MATLAB读取nc文件并转换为.tif格式

MATLAB读取nc文件并转换为.tif格式

1. 博文目的

（1）利用matlab读取nc文件，并将nc文件转换为.tif文件

2.时间

（1）开始时间：2017年04月08日

（2）结束时间：2017年04月09日

3.关键字

ncread, georasterref,geotiffwrite

4.MATLAB代码

%*************************************************************************%
 %程序目的：学习将nc文件转换为.tif文件并写出的方法
 %2017年04月09日
%*************************************************************************%

clc
clear all

InPath = 'F:\2科研\3书籍撰写\2_CSDN\3_GeoRasterref\1_TestData\';
InFile = strcat(InPath,'spei01.nc');
%获取nc文件的基本信息
ncdisp(InFile)

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')

5. 补充

5.1 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))])，增加程序的灵活性。

”