作家
登录

基于Python的Grib数据可视化

作者: 来源: 2017-10-31 11:12:56 阅读 我要评论

3、安装basemap

  1. python3 setup.py install 

二、grib数据攫取

(三)pygrib安装

  • pygrib依附
    • Jasper
    • GRIB API
    • numpy
    • pyproj
  • 安装过程

因为之前已经安装了numpy和pyproj,这里只需安装Jasper和GRIB API即可安装pygrib

1、安装Jasper

  1. ./configure 
  2.  
  3. make && date 
  4.  
  5. sudo make install && date  

2、安装GRIB API

  1. ./configure --with-jasper='/usr/local/' 
  2.  
  3. make && date 
  4.  
  5. sudo make install && date  

3、安装pygrib

安装pygrib之前起重要根据本身的实际情况修改文件目次下的setup.cfg文件,最重要的就是修改grib_api_dir和jasper_dir,这两个是方才安装的Jasper和GRIB API的路径,如不雅这两个地址不精确安装会报错

修改好就可以正常安装了

  1. python3 setup.py install 

固然我做的器械和蔼象沾边,然则我本身并不是气候专业出身,所有这些器械都是我慢慢研究揣摩出来的,所以有些方面可能讲的比较外行,有纰谬的处所迎接大年夜家留言斧正。

(一)导入pygrib模块

  1. >>> import pygrib 

(二)打开Grib文件

2、安装pyproj

  1. >>> grbs = pygrib.open('/Users/Kallan/Documents/data/echhae50.082'

(七)图像展示

  1. >>> grbs.seek(0) 
  2.  
  3. >>> for grb in grbs: 
  4.  
  5. grb 
  6.  
  7. 1:Geopotential Height:gpm (instant):regular_ll:isobaricInhPa:level 500:fcst time 24 :from 201507081200  

信息解读

1 :数据列表的行号,有的文件可能包含多个数据

Geopotential Height:数据的名称

regular_ll:惯例数据,其拭魅这个字段我也不清跋扈

isobaricInhPa:这个字段表示的是数据属性,此处表示是以hPa为单位的等压面

fcst time 24 :预告时效

from 201507081200 :起报时光

综合膳绫擎的信息可以得出,这个文件是大年夜2015年7月8日12时开端的24小时后500hPa等压面高度场数据

(四)导出文件数据

  1. >>> grb = grbs.select(name='Geopotential Height')[0] 
  2.  
  3. >>> data = grb.values 
  4.  
  5. >>> print(data.shape,data.min(),data.max()) 
  6.  
  7. (37, 37) 5368.6796875 5941.0390625 
  8.  
  9. >>> lat,lon=grb.latlons() 

      推荐阅读

      华硕也出妖板:双路Xeon竟能超频 史上第一

    2017年架构师最重要的48个小时 | 8折倒计时 都说自古要搬出华擎,比如比来搞了全球第一款Mini-ITX迷你规格的发烧级X299,但老大年夜哥华硕今天也发了块妖板,为双路工作站/办事器打造,却支>>>详细阅读


    本文标题:基于Python的Grib数据可视化

    地址:http://www.17bianji.com/lsqh/38292.html

关键词: 探索发现

乐购科技部分新闻及文章转载自互联网,供读者交流和学习,若有涉及作者版权等问题请及时与我们联系,以便更正、删除或按规定办理。感谢所有提供资讯的网站,欢迎各类媒体与乐购科技进行文章共享合作。

网友点评
自媒体专栏

评论

热度

精彩导读
栏目ID=71的表不存在(操作类型=0)