refvec=[1 60 70];
worldmap([20 60],[70 140])%描绘出中国地图的轮廓,其中的China是描绘中国地图,要是换做其他地图,就可以换成‘World’,或者'North Pole', 'South Pole', 'Pacific'
%陆地水格网数据显示geoshow(gca,data/24,refvec,'DisplayType','texturemap');%data是陆地水格网数据
%海岸线边界显示
S = shaperead('landareas','UseGeoCoords',true);
geoshow([S.Lat], [S.Lon],'Color','black');
ylabel(colorbar,'cm/yr','fontsize',10);title('Trend map of Terrestrial water change from GLDAS')
%中国边界显示
% L = shaperead('D:\MATLAB\toolbox\fengweiigg-GRACE_Matlab_Toolbox-6e91b05\GRACE_data\boundary_files\China\中华人民共和国.shp','UseGeoCoords',true);
% geoshow([L.Lat], [L.Lon],'Color','black');
B = shaperead('D:\BaiduNetdiskDownload\区划\区划\省.shp','UseGeoCoords',true);
geoshow([B.Lat], [B.Lon],'Color','black');