一、数据的下载
数据详细解释文档:GMD - The global water resources and use model WaterGAP v2.2d: model description and evaluation
二、数据读取
下载对应水文分量的”nc4”文件,以TWS为例,利用如下matlab代码读取
clearvars
addpath('E:\RSE\WGHM\data');
file='watergap_22d_WFDEI-GPCC_histsoc_tws_monthly_1901_2016.nc4';
lon=ncread(file,'lon');lat=ncread(file,'lat');
tws=rot90(ncread(file,'tws'));
index_str=1;index_end=12;num_file=index_end-index_str+1;
TWS=tws(:,:,index_str:index_end);%截取1901年的数据
for kk=1:num_file
tep(1,1)=TWS(1,1,1);temp(:,:)=TWS(:,:,kk);
for jj=1:size(lon,1)
for ii=1:size(lat,1)
if temp(ii,jj)==tep
temp(ii,jj)=0;%将数据中的缺省值替换为0
end
end
end
TWS_new(kk,:,:)=temp(:,:);
end
TWS_mean=mean(TWS_new);%去平均
for i=1:num_file
TWS_new(i,:,:)=TWS_new(i,:,:)-TWS_mean;
end
三、绘制全球格网图
利用冯伟老师的“GRACE_Matlab_Toolbox-master”工具包绘制全球格网图
for j=1:num_file
gws_1(:,:)=TWS_new(j,:,:);
gws_1Q=gws_1(:,1:360);gws_1H=gws_1(:,361:720);
gws_2=[gws_1H,gws_1Q]; %gldas_2(find(isnan(gldas_2)==1))=0;%找到缺失值NaN赋值为0
TWS(:,:,j)=gws_2(:,:);
end
for kk=1:num_file
temp(:,:)=TWS(:,:,kk);
WGHM(:,:,kk)=flipud(temp);
end
figure(1);
temp(:,:)=WGHM(:,:,1);
gmt_grid2map((temp)*0.1,-30,30,0,0,'cm',['WGHM TWS 1901.01']);