- 积分
- 42414
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-5
- 最后登录
- 1970-1-1
|
发表于 2013-7-2 11:41:12
|
显示全部楼层
本帖最后由 royalosyin 于 2013-7-2 11:56 编辑
试一试下面的脚本,直接用matlab自带的map工具箱。我用下面的脚本直接绘制NCEP或者CRU数据没有问题。
clear; clc
%% 读入你自己的数据,包括图数据和经纬度。
indata = ;
lat = ;
lon = ;
[x,y] = meshgrid(lon, lat);
maxiu = max(max(indata));
miniu = min(min(indata));
figure
%%开始绘图
project2; %设置头像
pcolorm(y, x, indata);
caxis([miniu maxiu]);
cmap = colormap( jet(20) );
colormap(cmap);
% colormap(flipud (colormap) ) %颜色反转了
shading interp
colorbar( 'horiz' );
%%画国家边界
load coast
plotm(lat, long, 'k')
%% 设置 title
tightmap %if you do not use this command, the pcolorm will fail.
title(['SPGI Date: ' num2str(time(itime,1)) '-' num2str(time(itime,2))...
'-' num2str(time(itime,3)) ], 'FontSize', 12); %'FontWeight','bold'
------- 下面是 project2.m 主要用于设置投影--------------------------------------------
------- 极地投影可以参考 http://www.mathworks.com.au/help/map/axes-for-drawing-maps.html
axesm(...
'MapProjection','wagner4',...
'MapLatLimit',[-48 -34],...
'Maplonlimit',[166 179],...
'labelformat','compass',...
'grid','off',...
'fontsize',8,...
'origin',[0 180 0]);
%'parallellabel','on',...
%'meridianlabel','on',...
%'mlabelparallel','south',...
|
|