最近有个项目需要开发个气象信息API,可以通过经纬度查找未来几天的气象信息。
经过几天的研究,现在简单总结一下。
1、数据来源
数据来源采自美国国家环境预报中心的GFS(全球预报系统),该系统每天发布4次全球范围的气象数据,分辨率最高可达到0.25° x 0.25°,精度还是比较可观的。
2、数据下载
GFS数据提供FTP下载方式:http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/
每次发布的数据保存在命名为gfs.YYYYMMDDHH的文件夹中。本次需要的数据精度为0.25°(0p25),所以数据的文件名为:gfs.t{HH}z.pgrb2.0p25.f{XXX}
其中HH表示发布的时间,XXX表示未来几小时的预报数据。例如gfs.t00z.pgrb2.0p25.f001 表示0时发布的未来1小时气象数据信息。
那么问题来了,这个文件太大了,因为分辨率比较高,每个文件高达200多MB,不仅下载速度慢,处理慢,对存储空间也是一个考验。
于是找到了一个可以在线过滤要下载数据的地址:
http://nomads.ncep.noaa.gov/cgi-bin/filter_gfs_0p25.pl
这个网站允许用户仅下载过滤后的文件,并提供下载网址,这样文件就会小很多。
例如,在本项目中,仅需要下载降雨数据(PWAT),并且经纬度范围为特定区域,每次发布的数据全部下载下来也才8M多点。
3、wgrib2
GFS的数据格式为GRIB,那么什么是GRIB呢? 官方的解释如下&#