爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13157|回复: 4

[求助] 关于nc数据绘制hgt图像

[复制链接]

新浪微博达人勋

发表于 2019-3-29 13:51:19 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
根据前辈们的程序,我改成了hgt绘图,但是现在遇到一个问题,每个时次的绘图,他的值是变化的,以500hpa为例,等高线一会儿是100的间距,一会儿是40的间距,就差了一个时次,这是为什么呢,有没有办法固定间距。谢谢。 1.png 2.png
另附程序:
from netCDF4 import Dataset
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
hgt=Dataset('./hgt.2004.nc')
hgts=hgt.variables['hgt'][530][5][:]
m=Basemap(projection='cyl',llcrnrlat=0,urcrnrlat=70,llcrnrlon=50,urcrnrlon=150)
m.drawcountries()
m.drawcoastlines()
lat=hgt.variables['lat'][:]
lon=hgt.variables['lon'][:]
index1=np.logical_and(lon>=50,lon<=150);index2=np.logical_and(lat>=0,lat<=70)
lons=lon[index1];lats=lat[index2]

hgt1=hgts[index2,:]
hgt2=hgt1[:,index1]
nx,ny=np.meshgrid(lons,lats)
x,y=m(nx,ny)
m.contour(x,y,hgt2,fontsize=10,colors='k')


hgt3=m.contour(x,y,hgt2,colors='k')
plt.clabel(hgt3, fontsize=8,fmt='%.0f')

plt.title('Height and Temperature ',size=20)

plt.show()

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-3-29 16:14:49 | 显示全部楼层
自己定义每个level的数值就可以了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-3-29 21:50:59 | 显示全部楼层
Wetter 发表于 2019-3-29 16:14
自己定义每个level的数值就可以了

找了半天都没找到怎么定义level这个数值。。。请问怎么定义
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-3-29 21:53:15 | 显示全部楼层
简单的世界 发表于 2019-3-29 21:50
找了半天都没找到怎么定义level这个数值。。。请问怎么定义

直接搜索啊
https://matplotlib.org/api/_as_gen/matplotlib.pyplot.contour.html
参见level一栏
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-3-29 22:04:42 | 显示全部楼层
Wetter 发表于 2019-3-29 21:53
直接搜索啊
https://matplotlib.org/api/_as_gen/matplotlib.pyplot.contour.html
参见level一栏

谢谢!看懂了。加一个levels命令在后面就能让数据更加丰富了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表