请教meteoInfoLAB数据处理的问题！

道狼52203

王老师您好！又要麻烦你请教两个小问题。
1、读取micaps3类数据以后，怎么取某一个区域的数据，我看了脚本用下面方法出错了
f = addfile_micaps(os.path.join(datadir, '000.000'))
pr = f.stationdata('Var1')
polygon1 = [(90.71, 40.82), (90.71, 45.4), (96.5, 45.4), (96.5, 40.82), (90.71, 40.82)]
pr1=pr.maskout(polygon1)

Traceback (most recent call last):
  File "<iostream>", line 27, in <module>
  File "D:\MeteoInfo\pylib\mipylib\dimarray.py", line 574, in maskin
    return PyStationData(self.data.maskin(polygon))
at org.meteoinfo.global.MIMath.getExtent(MIMath.java:454)
at org.meteoinfo.geoprocess.GeoComputation.pointInPolygons(GeoComputation.java:233)
at org.meteoinfo.data.StationData.maskin(StationData.java:966)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
at java.lang.reflect.Method.invoke(Unknown Source)
java.lang.ClassCastException: java.lang.ClassCastException: org.python.core.PyTuple cannot be cast to org.meteoinfo.shape.Shape

2、怎么批量取站点数据

stid=（51463，51076,52203，...）
st=pr.filter（stid）这样好像不行

道狼52203

用shp 文件取区域文件不提示错误，结果还是跟以前一样的？？
hmlk=shaperead('D:/MeteoInfo/map/hm.shp')
mdata = pr.maskout(hmlk)
print len(mdata)
tmin = mdata.min()
tmax = mdata.max()
print tmin,tmax

MeteoInfo

1、下载最新的MeteoInfo版本：http://yun.baidu.com/share/link?shareid=669776748&uk=51062435

2、参考下面的脚本：

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

复制代码