- 积分
- 4511
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-2-9
- 最后登录
- 1970-1-1
![[ksha_bonk] 粉丝数:70 微博数:639 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
发表于 2014-3-15 21:37:07
|
显示全部楼层
本帖最后由 沙颖凯 于 2014-3-15 21:43 编辑
添加经纬度和海岸线往往不是进一步的问题了,而是重头开始画,根据你的东西大致写了个例子。
%% An example of plotting geographical data in a given map projection (with coast line)
% latlim, lonlim is the range of lat/lon grid you would like to plot
% latitude, longitude is your lat/lon grid (LATnum*LONnum)
% sst, u, v, is the data imported from *.nc
% label is the interval of contour
label=[20 25 30 35];
latlim=[-90 90];
lonlim=[-180 180];
load coast
fig=figure('Color', 'w');
axesm('MapProjection', 'Eqdcylin',...
'MapLatLim', latlim, 'MapLonLim', lonlim)
framem('FLineWidth', 2.0)
gridm on
plabel on
mlabel('South')
tightmap
axis off
hold on
contourfm(latitude, longitude, sst, label, 'LineStyle', 'none')
Handle=quiverm(latitude, longitude, v, u, 1.5)
set(Handle, 'Color', 'k')
linem(lat, long, 'Color', 'k')
caxis([20 35]')
colorbar('EastOutside')
|
|