- 积分
- 352
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-9-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 happymgj 于 2017-8-21 16:18 编辑
guojie=shaperead('bou2_4l.shp');%读取国界shp文件的内容bou1_4lx=[guojie(:).X];%提取经度信息
bou1_4ly=[guojie(:).Y];%提取纬度信息
plot(bou1_4lx,bou1_4ly,'-k','LineWidth',0.8);%绘国界
hold on;
latmax=leftup_lat/100;
lonmin=leftup_lon/100;
latmin=rightdown_lat/100;
lonmax=rightdown_lon/100;
%tbb=contourf(lonmin:0.1:lonmax,latmax:-0.1:latmin,data1);
contourf(lonmin:0.1:lonmax,latmax:-0.1:latmin,data1,[-80:2:-52],'LineStyle','none');
tbbcolor=[1 0 1;0.6902 0.18824 0.37647;1 0 0;1 0.75294 0.79608;0.72156 0.52549 0.04314;0.8549 0.64706 0.12549;1 1 0;0.13333 0.5451 0.13333;0.60392 0.80392 0.19608;0.19608 0.80392 0.19608;0.39216 0.58431 0.92941;0 1 1;0.52941 0.80784 0.92157;0.68627 0.93333 0.93333];
colormap(tbbcolor);
tbbyanse=colorbar;
set(tbbyanse,'ytick',[-80:2:-52]);
axis([117 123 36 41]);
set(gca,'Box','off');
|
|