- 积分
- 2
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-5
- 最后登录
- 1970-1-1

|
发表于 2019-7-16 18:38:01
|
显示全部楼层
老师,您好,我是新手小白,看不太懂代码,看您发帖,模仿着想试试画图,我把shp文件替换成自己的,但是还是有问题,不知道什么原因,你可以帮忙看一下吗?
fn = 'G:/lianxi/dixing-china.nc'
f = addfile(fn)
data = f['z'][:,:]
#Plot
axesm(tickfontsize=16)
world = shaperead('C:/MI/MeteoInfo/map/country.shp')
lchina = shaperead('C:/MI/MeteoInfo/map/china.shp')
geoshow(world)
geoshow(lchina,edgecolor='k',size=0.3)
levels = arange(0, 5200, 400)
layer = imshowm(data,levels,cmap='temp_diff_18lev')
masklayer(lchina, [layer])
xlim(70,136)
ylim(15,55)
xlabel('Longitude',fontsize=17)
ylabel('Latitude',fontsize=17)
colorbar(layer,shrink=1,label='Height(m)') |
|