CRU气象数据处理

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