- 积分
- 453
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
根据前辈们的程序,我改成了hgt绘图,但是现在遇到一个问题,每个时次的绘图,他的值是变化的,以500hpa为例,等高线一会儿是100的间距,一会儿是40的间距,就差了一个时次,这是为什么呢,有没有办法固定间距。谢谢。
、
另附程序:
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()
|
|