请问MICAPS第三类资料怎么根据站点来抽取资料

James_z

参考了王老师的http://bbs.06climate.com/forum.p ... 1290&extra=page%3D1  这篇贴子的抽取MICAPS第一类站点代码，
但是我在自己使用MICAPS第三类降水站点的资料有一步报错
就是这一步，df = f1.read_dataframe()，
报错代码是：  File "<iostream>", line 8, in <module>
  File "D:\meteoinfo\MeteoInfo\pylib\mipylib\dataset\dimdatafile.py", line 158, in read_dataframe
    df = self.dataset.getDataInfo().readDataFrame()
AttributeError: 'org.meteoinfo.data.meteodata.micaps.MICAPS3DataInf' object has no attribute 'readDataFrame'
但是我看了下第三类跟第一类资料都是站点资料呀，为啥用不了
这是print第三类资料的结果
Description: diamond 3 2018年07月05日00时00分00秒一小时降水
Time: 2018-07-05 00:00
File Name: F:/micapsdata/20180705/SURFACE_N/r1/18070500.000
Dimensions: 1
        station = 2401;
Global Attributes:
        : data_format = "MICAPS 3"
Variations: 5
        String Stid(station);
        float Longitude(station);
        float Latitude(station);
        float Altitude(station);
        float Var1(station);

请问有什么办法可以根据站点抽取MICAPS第三类资料的降水数据吗？谢谢大家！

MeteoInfo

micaps 3目前不支持readDataFrame，可以直接用变量来提取数据，参考下面的脚本：

1. f = addfile('D:/Temp/micaps/18070500.000')
   
2. stid = f['Stid'][:]
   
3. data = f['Var1'][:]
   
4. lon = f['Longitude'][:]
   
5. lat = f['Latitude'][:]
   
6. t = f.gettime(0)
   
7. #Plot
   
8. axesm(bgcolor=(204,255,255))
   
9. geoshow('country', edgecolor=None, facecolor=(255,251,195))
   
10. geoshow('cn_province', edgecolor='k')
    
11. levs = arange(0., 15., 1.)
    
12. levs[0] = 0.01
    
13. layer = scatter(lon, lat, data, levs, cmap='WhBlGrYeRe')
    
14. layer.legend().getLegendBreaks()[0].setSize(1)
    
15. colorbar(layer)
    
16. yticks([20,30,40,50])
    
17. title('Precipitation (' + t.strftime('%Y-%m-%d %H:00') + ')')
    
18. xlim(72, 136)
    
19. ylim(16, 55)

f = addfile('D:/Temp/micaps/18070500.000')<br /> stid = f['Stid'][:]<br /> data = f['Var1'][:]<br /> lon = f['Longitude'][:]<br /> lat = f['Latitude'][:]<br /> t = f.gettime(0)<br /> #Plot<br /> axesm(bgcolor=(204,255,255))<br /> geoshow('country', edgecolor=None, facecolor=(255,251,195))<br /> geoshow('cn_province', edgecolor='k')<br /> levs = arange(0., 15., 1.)<br /> levs[0] = 0.01<br /> layer = scatter(lon, lat, data, levs, cmap='WhBlGrYeRe')<br /> layer.legend().getLegendBreaks()[0].setSize(1)<br /> colorbar(layer)<br /> yticks([20,30,40,50])<br /> title('Precipitation (' + t.strftime('%Y-%m-%d %H:00') + ')')<br /> xlim(72, 136)<br /> ylim(16, 55)

James_z

MeteoInfo 发表于 2022-11-28 16:33
micaps 3目前不支持readDataFrame，可以直接用变量来提取数据，参考下面的脚本：

感谢王老师耐心回答！我还有个小问题，请问MICAPS第三类资料能根据站点抽取数据吗，比如我想要提取青海省清水河站（编号：56034）的降水数据，但是MICAPS第三类资料用不了.readDataFrame()这条命令，目前有什么函数可以实现嘛

MeteoInfo

James_z 发表于 2022-11-28 18:02
感谢王老师耐心回答！我还有个小问题，请问MICAPS第三类资料能根据站点抽取数据吗，比如我想要 ...

读取数据后可以利用stid和data数组构造一个DataFrame，例如：
df = DataFrame(data=data, index=stid, columns=['pre'])
下一个MeteoInfo版本会支持micaps 3数据文件的 read_dataframe 方法。

James_z

MeteoInfo 发表于 2022-11-28 20:32
读取数据后可以利用stid和data数组构造一个DataFrame，例如：
df = DataFrame(data=data, index=stid, c ...

谢谢王老师！刚刚试了这句代码，能成功输出站点降水数据了！感谢您！！