1.CRU气象数据
http://data.ceda.ac.uk/badc/cru/data/cru_ts/cru_ts_4.02/data/
2.matlab提取代码
clear;clc;close all
infile='Z:\cru_ts4.02.1901.2017.tmp.dat.nc';
outpath = 'Z:\1_720_month_TEMave\';
latlim = [-90 90];
lonlim = [-180 180];
rasterSize = [360 720];
%define the geo reference
R = georefcells(latlim,lonlim,rasterSize,'ColumnsStartFrom','south','RowsStartFrom','west');
start = [1,1,1];
count = [Inf,Inf,Inf];
stride = [1,1,1];
IData = ncread(infile,'tmp',start,count,stride);
n=size(IData,3);
for i=1:n
Data = IData(:,:,i);
Data = Data';
Data =flipud(Data);
Data(find(isnan(Data)==1)) = -999;
Month = rem(i,12);
if Month == 0
Month = 12;
else
Month;
end
Year = 1900 + ceil(i/12);
outname=strcat(outpath,'TEMave_',num2str(Year*100+Month))
geotiffwrite([outname, '_720.tif'],Data,R);
end