MeteoInfoLab脚本示例：读取BUFR数据

GIS-Yuan

1. fn = r'D:\soft\MeteoInfo\sample\2170_Z_SURF_C_BCCD_20200608080112_O_AWS_FTM_PQC.BIN'
   
2. f = addfile(fn, keepopen=True)
   
3. obs = f['obs']
   
4. lon = obs.member_array('Longitude_high_accuracy')
   
5. lat = obs.member_array('Latitude_high_accuracy')
   
6. lon = (lon - 1.8E7) * 1.E-5
   
7. lat = (lat - 9.E6) * 1.E-5
   
8. alt = obs.member_array('Height_of_station_ground_above_mean_sea_level')
   
9. alt = (alt - 4000) * 1.E-1
   
10. st_name = obs.member_array('Station_or_site_name')
    
11. st_id = obs.member_array('WMO_station_number')
    
12. year = obs.member_array('Year')
    
13. month = obs.member_array('Month')
    
14. day = obs.member_array('Day')
    
15. hour = obs.member_array('Hour')
    
16. seq1 = obs.member_array('seq1')
    
17. pres = seq1.member_array('Pressure') * 1e1
    
18. seq2 = obs.member_array('seq2')
    
19. temp = seq2.member_array('Temperature-air_temperature') * 1e-1
    
20. seq5 = obs.member_array('seq5')
    
21. windspeed = seq5.member_array('Wind_speed')
    
22. #print('lon: {}; lat: {}'.format(lon, lat))
    
23. print('Pressure: {}; Temperature: {}; windspeed:{}'.format(pres, temp,windspeed))
    
24. 
    
25. f.close()

复制代码

王老师，解析出来的数值有误。

1. Pressure: 6130.0; Temperature: 368.2; windspeed:1762
   

复制代码

MeteoInfo

GIS-Yuan 发表于 2020-6-12 09:08
王老师，解析出来的数值有误。

需要看bufr编码格式文档，相关变量的单位、比例因子、基准值等信息来协助获取真实数据。

Dahoe

谢谢，王老师，是我读取的文件有误

棉花糖heyuan

请问大佬，怎么写入bufr数据？

shuying0218

请问如果在python中直接应用这段代码的话，需要导入什么库包呀？

来自安卓APP客户端

王宇涛

能否发下完整的代码

MeteoInfo

王宇涛 发表于 2022-4-19 16:03
能否发下完整的代码

为什么会认为代码不完整？

不吃鱼的小虾米

老师您好，我的一个BUFR数据可以成功读出来，但是我有很多个数据，请问怎么能批量读取呢

MeteoInfo

不吃鱼的小虾米 发表于 2022-10-28 16:14
老师您好，我的一个BUFR数据可以成功读出来，但是我有很多个数据，请问怎么能批量读取呢

可以用循环

不吃鱼的小虾米

MeteoInfo 发表于 2022-10-29 09:40
可以用循环

老师您好，我还是这个软件的小白，自己写循环的时候总是语法错误，请问有没有可以学习基础语法的网址呢？我是想写一个数组，数组里是所有文件的命名，因为文件很多且命名无规则，NCL和matlab里可以用systemfunc和dir来实现，请问MeteoInfoLab是否也有类似的函数呢？谢谢您！