请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
12
返回列表 发新帖
楼主: oppion123

grib2文件读取出来的数据是什么含义?

[复制链接]

新浪微博达人勋

发表于 2017-8-3 18:03:49 | 显示全部楼层
oppion123 发表于 2017-8-3 17:39
谢谢老师的解答。
文件头里面的信息如下:
DSAA

用MeteoInfoLab脚本读取数据会更加灵活方便,比如读取定点某个时间范围多个变量的值。可以参考下面的脚本程序:

  1. fn = 'C:/Temp/nww3.t00z.grib.grib2'
  2. f = addfile(fn)
  3. lon = 160.63
  4. lat = 50.62
  5. lonlim = str(lon)
  6. latlim = str(lat)
  7. tn = f.timenum()
  8. hline = 'Time,Lon,Lat,dwws,shcwwss'
  9. print hline
  10. for i in range(tn):
  11.     t = f.gettime(i)
  12.     dwws = f['Direction_of_wind_waves_surface'][i,latlim,lonlim]
  13.     shcwwss = f['Significant_height_of_combined_wind_waves_and_swell_surface'][i,latlim,lonlim]
  14.     data = t.strftime('%Y%m%d') + ',%.2f,%.2f,%.2f,%.2f' % (lon,lat,dwws,shcwwss)
  15.     print data


运行后输出:
Time,Lon,Lat,dwws,shcwwss
2017080100,160.63,50.62,150.12,1.64
2017080103,160.63,50.62,149.96,1.67
2017080106,160.63,50.62,149.77,1.70
2017080109,160.63,50.62,149.48,1.78
2017080112,160.63,50.62,148.91,1.88
2017080115,160.63,50.62,148.06,1.95
2017080118,160.63,50.62,147.03,1.97
2017080121,160.63,50.62,145.68,1.99
2017080200,160.63,50.62,143.92,2.09
2017080203,160.63,50.62,142.59,2.20
2017080206,160.63,50.62,142.03,2.22
2017080209,160.63,50.62,141.88,2.15
2017080212,160.63,50.62,141.85,2.05
2017080215,160.63,50.62,141.70,1.94
2017080218,160.63,50.62,141.43,1.83
2017080221,160.63,50.62,141.25,1.72
2017080300,160.63,50.62,141.64,1.64
2017080303,160.63,50.62,142.98,1.58
2017080306,160.63,50.62,146.00,1.54
2017080309,160.63,50.62,150.28,1.52
2017080312,160.63,50.62,154.97,1.50
2017080315,160.63,50.62,160.59,1.50
2017080318,160.63,50.62,164.16,1.48
2017080321,160.63,50.62,165.68,1.45
2017080400,160.63,50.62,166.99,1.42
2017080403,160.63,50.62,167.75,1.39
2017080406,160.63,50.62,168.23,1.35
2017080409,160.63,50.62,168.44,1.33
2017080412,160.63,50.62,168.51,1.31
2017080415,160.63,50.62,168.44,1.30
2017080418,160.63,50.62,168.50,1.29
2017080421,160.63,50.62,168.61,1.30
2017080500,160.63,50.62,169.09,1.33
2017080503,160.63,50.62,170.05,1.37
2017080506,160.63,50.62,171.23,1.44
2017080509,160.63,50.62,172.36,1.53
2017080512,160.63,50.62,173.39,1.65
2017080515,160.63,50.62,174.43,1.78
2017080518,160.63,50.62,175.25,1.88
2017080521,160.63,50.62,175.78,1.96
2017080600,160.63,50.62,176.04,2.01
2017080603,160.63,50.62,176.27,2.03
2017080606,160.63,50.62,176.57,2.05
2017080609,160.63,50.62,176.46,2.06
2017080612,160.63,50.62,176.21,2.15
2017080615,160.63,50.62,175.80,2.41
2017080618,160.63,50.62,175.20,2.76
2017080621,160.63,50.62,174.49,2.99
2017080700,160.63,50.62,173.66,3.04
2017080703,160.63,50.62,172.85,2.95
2017080706,160.63,50.62,171.86,2.82
2017080709,160.63,50.62,170.91,2.70
2017080712,160.63,50.62,169.97,2.59
2017080715,160.63,50.62,169.11,2.48
2017080718,160.63,50.62,168.61,2.39
2017080721,160.63,50.62,168.80,2.32
2017080800,160.63,50.62,169.70,2.26
2017080803,160.63,50.62,170.77,2.23
2017080806,160.63,50.62,171.13,2.20
2017080809,160.63,50.62,170.60,2.14
2017080812,160.63,50.62,169.60,2.07


密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-4 09:27:02 | 显示全部楼层
MeteoInfo 发表于 2017-8-3 18:03
用MeteoInfoLab脚本读取数据会更加灵活方便,比如读取定点某个时间范围多个变量的值。可以参考下面的脚本 ...

非常感谢老师的解答。我用您的程序读取了我自己下载下来的grib2文件,这个文件里只有2017年7月30~8月10号的浪和涌的有义波高,所以我把您例子中的浪向这个变量删除了。
表格文件保存后的数据如下:
DSAA
361 181
0.0 360.0
-90.0 90.0
0.0 12.699999809265137


经纬度160.63 50.62在(0~360,-90~90)范围内。我运行后的结果出错。我也试了其他几个GRIB2文件,均报错了。还得麻烦老师帮我看看错在哪里?

万分谢谢。

grib2文件和程序报错截图附上

程序报错截图

程序报错截图

grib2文件读取

grib2文件读取

HTSGW_02.t00z.grib2

2.24 MB, 下载次数: 0, 下载积分: 金钱 -5

grib2文件

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-4 09:56:08 | 显示全部楼层
oppion123 发表于 2017-8-4 09:27
非常感谢老师的解答。我用您的程序读取了我自己下载下来的grib2文件,这个文件里只有2017年7月30~8月10号 ...

在13行中删除最后的 %.2f 和前面的逗号
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-8-4 10:38:21 | 显示全部楼层
MeteoInfo 发表于 2017-8-4 09:56
在13行中删除最后的 %.2f 和前面的逗号

谢谢老师,我已经解决了。
这个程序是输出指定的经纬度的参数值,可不可以输出一段区间经纬度的参数值?
我试着将经纬度改成比如说(180-200  40-50)区间做个循环,但是运行报错,我不太了解这个脚本程序的一些函数,您能帮我看看吗?谢谢您。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-8-4 13:10:09 | 显示全部楼层
oppion123 发表于 2017-8-4 10:38
谢谢老师,我已经解决了。
这个程序是输出指定的经纬度的参数值,可不可以输出一段区间经纬度的参数值? ...

建议你加入MeteoInfo QQ群,这种技术问题交流起来方便一些。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-8-15 21:11:25 | 显示全部楼层
老师们好,请问用什么软件读MeteoInfo保存的grd文件,以及如何将数据输入到excel里,萌新万分感谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-4-19 09:45:43 | 显示全部楼层
你否把你读取Grib2是方法共享一下 给我学习
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表