- 积分
- 226
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2020-2-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我根据meteoinfo脚本汇总中的四川盆地三维地形图脚本,修改了经度参数后画了一张图,脚本是这样:
f = addfile('E:/dixing-china.nc')
elev = f['z']['26:34','98:108']/1000.0
lchina = shaperead('G:/BaiduNetdiskDownload/MeteoInfo\map\cn_province.shp')
#Plot
figure(figsize=[1300,950],newfig=False)
ax = axes3d(tickfontsize=20,position=[0.1, 0.13, 0.75, 0.85],bbox=False)
levs = arange(0,6.1,0.1)
cols = makecolors(len(levs) + 1, cmap='temp_diff_18lev',reverse=False)
ls = ax.plot_surface(elev, levs, colors=cols, edge=False)
colorbar(ls,ticks=arange(0,6.1,1),orientation='horizontal',label='Altitude(km)',labelloc='right',fontsize=24,shrink=0.8,aspect=40,yshift=35)
ylabel('Latitude (degrees)',fontsize=20)
xlabel('Longitude (degrees)',fontsize=20)
zlabel('Altitude (km)',fontsize=20)
zlim(0, 6)
xlim(98, 108)
ylim(26,34)
但是我画出来的图看不到中间部分,想请教一下该怎么改。
|
|