- 积分
- 57041
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-6-21
- 最后登录
- 1970-1-1
|
发表于 2015-7-17 23:40:50
|
显示全部楼层
本帖最后由 MeteoInfo 于 2018-12-17 16:33 编辑
 - f = addfile_micaps('D:/Temp/micaps/10101414.000')
- data = f['Temperature'][:]
- lon = f['Longitude'][:]
- lat = f['Latitude'][:]
- t = f.gettime(0)
- #To grid data
- x = arange(70, 140, 0.5)
- y = arange(15, 58, 0.5)
- gdata,gx,gy = griddata((lon, lat), data, xi=(x, y), method='idw')
- gdata.savegrid(x, y, 'D:/Temp/test/savegrid_test.txt')
- #Plot
- axesm()
- geoshow('country')
- geoshow('cn_province')
- levs = arange(0, 35, 2)
- layer = contourfm(x, y, gdata, levs)
- colorbar(layer)
- yticks([20,30,40,50])
- title('Temperature (' + t.strftime('%Y-%m-%d %H:00') + ')')
- xlim(72, 136)
- ylim(16, 55)
|
|