怎样MeteoInfoLAB取数据问题？

道狼52203

王老师您好！
用Micaps第三类格式做最高气温分布图是，可以用f.max和f.min取出最大和最小值，怎样取它的站名、区站号、经纬度等信息？
谢谢！！

MeteoInfo

你是指最大值、最小值所在的站点吗？

道狼52203

对，最大值、最小值所在的站点。

MeteoInfo

道狼52203 发表于 2015-7-17 11:14
对，最大值、最小值所在的站点。

更新了软件（还没有打包发布），给PyStationData类增加了minloc和maxloc函数，并重载了__getitem__。下面是一个例子，不知道能否满足你的需求？

1. f = addfile_micaps('D:/Temp/micaps/13042408.000')
   
2. pr = f.stationdata('Var1')
   
3. axesm()
   
4. mlayer = shaperead('D:/Temp/map/country1.shp')
   
5. geoshow(mlayer)
   
6. layer = scatterm(pr, s=4, edge=False)
   
7. title('MICAPS 3')
   
8. xlim(70, 140)
   
9. ylim(15, 55)
   
10. colorbar(layer, orientation='horizontal')
    
11. minv, minidx = pr.minloc()
    
12. maxv, maxidx = pr.maxloc()
    
13. min_st, min_lon, min_lat = pr[minidx]
    
14. max_st, max_lon, max_lat = pr[maxidx]
    
15. print 'Min: %.2f, Stid: %s, Lon: %.2f, Lat: %.2f' %(minv, min_st, min_lon, min_lat)
    
16. print 'Max: %.2f, Stid: %s, Lon: %.2f, Lat: %.2f' %(maxv, max_st, max_lon, max_lat)

f = addfile_micaps('D:/Temp/micaps/13042408.000') <br /> pr = f.stationdata('Var1') <br /> axesm() <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> layer = scatterm(pr, s=4, edge=False) <br /> title('MICAPS 3') <br /> xlim(70, 140) <br /> ylim(15, 55) <br /> colorbar(layer, orientation='horizontal') <br /> minv, minidx = pr.minloc() <br /> maxv, maxidx = pr.maxloc() <br /> min_st, min_lon, min_lat = pr[minidx] <br /> max_st, max_lon, max_lat = pr[maxidx] <br /> print 'Min: %.2f, Stid: %s, Lon: %.2f, Lat: %.2f' %(minv, min_st, min_lon, min_lat) <br /> print 'Max: %.2f, Stid: %s, Lon: %.2f, Lat: %.2f' %(maxv, max_st, max_lon, max_lat)

道狼52203

谢谢！！ 可以下载了吗？

MeteoInfo

道狼52203 发表于 2015-7-17 15:52
谢谢！！ 可以下载了吗？

已经打包上传：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435