python 读取grib \grib2-程序员宅基地

技术标签: python  操作系统  

一、环境准备
(1)、python3环境

(2)、wgirb工具(用于读取grib1文件),下载地址:  ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib

(3)、wgirb2工具(用于读取grib2文件),下载地址:ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2

二、环境测试
工具下好了,windows下无需安装,怎么看是否能用呢?很简答,以wgrib为例,下载好了之后在我的文件夹下是这样的:

 

可以看到里面有一个应用程序,那么从命令行进入该文件夹下,运行命令:wgrib.exe就可以查看一些使用帮助了,如下图所示

 

可以看到很多提示信息,也就是教我们如何使用这个工具。wgrib2也是同样的查看。

接下来就可以使用python来读取grib1和grib2数据了(其实就是在python里面调用cmd命令而已)

三、读取数据
grib1和grib2数据我了解不多,说错的地方还请大家指正,它们是两种气象数据,在我看来,它相当于一个大仓库,里面是一个个“小房间”,每个房间都有编号,每个“小房间”存放的是一个气象变量,例如温度、湿度、压力等等,那么如何读取我们需要的气候变量呢?下面分别以grib1和grib2数据为例,来读取温度变量的值。二者读取方式相似,只是使用的工具不同。

刚刚说了,每一个“小房间”存放的是一个气象变量且“小房间”有编号,那么我们有2中方法进行读取,一种是根据“小房间”编号进行读取,另一种是根据变量名称进行匹配读取,一般来说后者使用的较多,毕竟使用名称是肯定错不了的。下面进行分别介绍。

1、按编号读取
既然要根据编号读取,那么首先要知道编号才行,还好wgrib工具中也提供了这和功能。查看“小房间”编号的代码如下:

import os
# 将当前目录转到wgrib工具目录
os.chdir("G:\mygrib\wgrib1")
# 执行cmd命令,查看各个“小房间”具体信息
os.system(r"wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -v")
运行得到的结果如下图所示:


如上图所示,我们可以得到很多有用的信息,例如:2m高的温度变量所在的“小房间”编号是82,接下来就可以根据这个编号82来进行读取2m温度变量,读取代码如下:

import os
# 将当前目录转到wgrib工具目录
os.chdir("G:\mygrib\wgrib1")
# 执行cmd命令,查看各个“小房间”具体信息
# os.system(r"wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -v")
# 知道了编号82,那么就可以根据这个编号把数据读到一个txt里面,注意:wgirb工具不支持读到csv文件,而wgrib2工具支持
os.system("wgrib G:\mete_data\W_NAFP_C_ECMF_20170713190330_P_C3E07131200071312001 -d 82 -h -text -o 2m_tmp.txt")
读取结果如下图所示:

 

除此之外,wgrib工具还支持“变量名称匹配读取”,在这里就不介绍了,可以看下一部分的利用wgrib2工具进行“变量名称匹配读取”二者是类似的。

2、按名称匹配读取
一般来说,用名称匹配读取这个方法的占多数,因为虽然grib文件可以看做“仓库”,一个个气象变量可以看做“小房间”,但是不同的grib文件内部的“小房间”分布可能是不一样的,这样的话同一气候变量的编号也就不确定,无法单靠一个编号批量读取很多文件。而利用名称匹配读取就不存在这个问题,因为气象变量的名称是不会变的。以读取grib2文件里的2m处的温度为例,这里使用wgrib2工具,wgrib2工具功能更强大,可以直接读出为csv文件,代码如下:

import os
# 将当前目录转到wgrib2工具目录
os.chdir("D:\mydownload\wgrib2")
# 进行气象变量名称匹配,读取2m处温度这个气象变量到2m_tmp.csv文件中
os.system("wgrib2 G:\mete_data\Z_NAFP_C_BABJ_20170709000000_P_CNPC-T639-GMFS-HNEHE-03900.grib2 -match ':TMP:2 m' -csv 2m_tmp.csv")
读取结果如下图所示:

 

 

csv格式的数据更便于我们处理,倒数第一列就是温度值。

同理,其他的气候变量同样可以使用匹配的方式读取出来,只是改变下代码-match后面的气象变量即可,大家可能会问这些气象变量的简写如何得知,这些都是固定的,目前为止我只用了温度和降水,其他的我也不知道,有需要的话可以去查一查。

结束
这两个工具的功能还有很多,目前我用到的只是一丢丢,有其他需要的童鞋可以使用wgrib.exe和wgrib2.exe命令来进行查看帮助,发掘更多功能。在linux下同样可以使用这两个工具进行读取,只是在安装和命令上有一点小小的差异,大家感兴趣的可以一起讨论下,水平有限,肯定有写的不对的地方,欢迎指正!
---------------------
作者:全村的希望呀
来源:CSDN
原文:https://blog.csdn.net/m0_37600626/article/details/79311129?utm_source=copy
版权声明:本文为博主原创文章,转载请附上博文链接!

posted on 2018-10-17 14:16 秦瑞It行程实录 阅读( ...) 评论( ...) 编辑 收藏

转载于:https://www.cnblogs.com/ruiy/p/9803849.html

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_33859665/article/details/93875352

智能推荐

java中的volatile关键字的功能详解_volatile 关键字,-程序员宅基地

文章浏览阅读4.6k次。Cookie的应用场景:1,记录用户的登陆状态,提示用户是否记住密码;2,购物车购物功能;我们知道,在web开发过程中,我们都要和cookie打交道,有时候离开了cookie还真玩不转。cookie最典型的应用莫过于登陆提示,最近在做一个小项目,正好要用到cookie的知识,在这里顺便做一下总结。_volatile 关键字,

通过 ICMP 协议实现 Ping Tunnel 建立可穿透网络隧道-程序员宅基地

文章浏览阅读7.1k次。Twitter via Ping Tunnel周四 Cola 没去幼儿园,中午带着他去 KFC 吃东西。回来的时候小林指着西总布胡同说走这条路回去还是原路返回,他说还是..._ping tunnel

基于springboot+vue.js的名城小区物业管理系统(附带文章和源代码设计说明文档ppt)-程序员宅基地

文章浏览阅读817次,点赞18次,收藏20次。博主介绍:CSDN深耕的技术专家、博客专家、有着常年的工作经验、全栈领域优质创作者,博客之星、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java、小程序技术领域和毕业项目实战文末获取源码+数据库。

Arthas使用教程 阿里巴巴开源项目、史上最强java线上诊断工具-程序员宅基地

文章浏览阅读1.4w次,点赞31次,收藏263次。什么是 Arthas摘录一段官方 Github 上的简介Arthas 是Alibaba开源的Java诊断工具,深受开发者喜爱。当你遇到以下类似问题而束手无策时,Arthas 可以帮助你解决:这个类从哪个 jar 包加载的?为什么会报各种类相关的 Exception?我改的代码为什么没有执行到?难道是我没 commit?分支搞错了?遇到问题无法在线上 debug,难道只能通过加日志再重新发布吗?线上遇到某个用户的数据处理有问题,但线上同样无法 debug,线下无法重现!是否有一个全局视角来._arthas

java 方法注释格式_JAVA注释方法及格式-程序员宅基地

文章浏览阅读2.7k次。2019独角兽企业重金招聘Python工程师标准>>>JAVA注释方法及格式1、单行(single-line)--短注释://……单独行注释:在代码中单起一行注释, 注释前最好有一行空行,并与其后的代码具有一样的缩进层级。如果单行无法完成,则应采用块注释。注释格式:/* 注释内容 */行头注释:在代码行的开头进行注释。主要为了使该行代码失去意义。注释格式:// 注释内容行尾注释:..._方法注释

egret4.X版本项目无法与egret 5.X项目共存解决_egret 4.x老项目升级-程序员宅基地

文章浏览阅读367次。在编译egret 5.X 项目项目中执行egret clean_egret 4.x老项目升级

随便推点

机器学习实战 | 自动化特征工程工具Featuretools应用-程序员宅基地

文章浏览阅读1.3w次,点赞7次,收藏70次。本篇讲解使用自动化特征工程工具Featuretools,对数据进行自动化特征工程的方法,并借助于BigMart Sales数据集来演示自动化特征工程的相关应用。_featuretools

eNSP 常用指令_stp region configuration-程序员宅基地

文章浏览阅读4.2k次,点赞21次,收藏80次。作为菜鸟记录下学习eNSP的笔记按照学校讲课进度记录,以网络学习的指令进行补充分类指令进入退出保存system-view 进入配置模式ctrl+z /quit 退出端口配置模式save 保存当前配置记得选Y命名和改名sysname SWC 给设备命名sys name j..._stp region configuration

【光伏功率预测】基于EMD-PCA-LSTM的光伏功率预测模型(Matlab代码实现)_计及相似日的 lstm 光伏出力预测模型研究-程序员宅基地

文章浏览阅读917次。随着电厂规模的不断扩增,电厂的数据量也呈爆炸式的增长,传统的神经网络光伏功率预测模型[7-10]一方面受电厂来源数据的制约,忽略了部分环境因素对光伏功率的影响[11] ,缺乏对多元环境序列信息的有效利用;另一方面,由于光伏功率与多元环境序列信息呈非线性变化,随着网络输入变量的增多,会导致模型收敛速度减慢[12-14] ,并出现过拟合问题;因此,要提高光伏功率预测模型的准确性,不仅要充分利用影响光伏功率的关键环境因素,也要进一步挖掘光伏功率预测与关键环境因素随时间变化的本质特征。长短期记忆神经网络;_计及相似日的 lstm 光伏出力预测模型研究

Flutter开发学习实战之商城项目(后端+Flutter端)_flutter 商城-程序员宅基地

文章浏览阅读4.7k次,点赞3次,收藏39次。Flutter_mall 商城项目引言:此Flutter工程项目是在学习 youxinLu 大佬写的一个商城项目:作者项目简介:Flutter_Mall是一款Flutter开源在线商城应用程序,是基于litemall基础上进行开发,litemall包含了Spring Boot后端 + Vue管理员前端 + 微信小程序用户前端 + Vue用户移动端感兴趣的同学可以自行研究部署,Flutter_Mall基本上包含了litemall中小程序的功能。学习的最好方法就是动手实践:于是将 youxinL_flutter 商城

Leetcode 104. 二叉树的最大深度(DAY 2)-程序员宅基地

文章浏览阅读380次。原题题目代码实现/** * Definition for a binary tree node. * struct TreeNode { * int val; * struct TreeNode *left; * struct TreeNode *right; * }; */int maxDepth(struct TreeNode* root){ if(!root) return 0; if(!root->left

记一次RATEL脱壳配合Il2CppDumper解密完成的样本分析-程序员宅基地

文章浏览阅读2.4k次,点赞2次,收藏8次。今天分析一个网赚类的样本,业务方要我搞清楚为什么每次完成任务出现红包时,总会随之弹出一个插屏广告。大概就是下图这样,说实话,弹几层广告确实挺恶心,关起来都麻烦先在测试机装上应用,发现有root检测,试过几个通用的过root检测的frida脚本,发现不行看了下代码,是网易的壳正好,最近在看Ratel,刚好来练练手,试试在非root环境下Ratel的脱壳怎么样Rate..._il2cppdumper

推荐文章

热门文章

相关标签