如何读取GRIB数据?快看Python大神整理的干货! 橙子心法 百家号17-11-0116:30 GRIB是WMO开发的一种用于交换和存储规则分布数据的二进制文件格式,主要用来表示数值天气预报的产品资料.现行的GRIB 码版本有GRIB1 和GRIB2 两种格式.由于小编只接触过GRIB1数据,所以本文仅讲解如何读取GRIB1数据. 使用Python读取GRIB数据要根据在何种平台分别讨论. 一.Linux及Mac OS系统下 Linux.Mac OS系统下使用Python语言实现Grib数据

之前写了一篇<基于Python的GRIB数据可视化>的文章,好多博友在评论里问我Windows系统下如何读取GRIB数据,在这里我做一下说明. 一.在Windows下Python为什么无法读取GRIB 大家在windows系统不能读取GRIB数据的主要原因是,GRIB_API在Windows下无法编译安装,从而导致pygrib安装失败.我曾经也为这个问题苦恼了很久,也到ECMWF论坛里找了很久,也给ECMWF发了邮件,回应我没有做Windows版本的打算,所以在Windows下直接用pygri

一.环境准备(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为例,下载好了之后在我的文件夹下是这样的: 可以看到里面有一个应用程序,那么从命

from osgeo import gdal import numpy as np path = './data/201912/anl_surf125.2019120100' dataset = gdal.Open(path) number_bands = dataset.RasterCount for i in range(1, number_bands): band = dataset.GetRasterBand(i) data_array = band.ReadAsArray() prin

基于Python的Grib数据可视化           利用Python语言实现Grib数据可视化主要依靠三个库——pygrib.numpy和matplotlib.pygrib是欧洲中期天气预报中心(ECMWF)的GRIG API C库的Python接口,通过这个库可以将Grib数据读取出来:numpy是Python的一种开源的数值计算扩展,这种工具可用来存储和处理大型矩阵:matplotlib是python著名的绘图库,它提供了一整套和matlab相似的命令API,十分适合交互式地进行制图:

参考 https://www.ecmwf.int/assets/elearning/eccodes/eccodes2/story_html5.htmlhttps://confluence.ecmwf.int/display/ECC/GRIB+exampleshttps://confluence.ecmwf.int/download/attachments/97363968/eccodes_grib_python_2018.pdf 关于Python读取GRIB格式数据,Kallan写过一篇“基于P

1.wgrib的使用 在cmd命令行下键入wgrib后即可察看wgrib相关命令参数,简要介绍如下: l        Inventory/diagnostic–output selections 详细目录或诊断的输出选择 -s                 short inventory 简短目录 -v                 verbose inventory 详细目录 -V                 diagnostic output <not inventory> 诊断

该文章来自博客:http://windforestwing.blog.163.c ... 412007103084743804/如有错误 ,大家及时指出啊!ps:meteoinfo可以直接处理grib数据=========================== grib格式转换心得 1.wgrib的使用 在cmd命令行下键入wgrib后即可察看wgrib相关命令参数,简要介绍如下: l        Inventory/diagnostic–output selections 详细目录或诊断的输出

参考:XML DOM 参考手册(w3school) 参考:python专题-读取xml文件 参考:请问用python怎么修改xml的节点值? 1. 读取标签内的文本(Python) 如下的 xml 文件一部分,读取里面的经纬度 <?xml version="1.0" encoding="UTF-8" ?> <ProductMetaData> <SatelliteID>GF1</SatelliteID> <Sens

LFD,非官方的Windows二进制文件的Python扩展包 LFD,非官方版本.32和64位.Windows.二进制文件.科学开源.Python扩展包 克里斯托夫·戈尔克(by Christoph Gohlke),LFD(荧光动力学实验室),加利福尼亚大学,Irvine, 网址:http://www.lfd.uci.edu/~gohlke/pythonlibs/ ------------------------------ A APSW,另一个Python Sqlite封装. Aspell-p

对于低分数据来说,常用的校正方式就是给定数据的经纬度查找表来进行校正.在GDAL中,这种校正方式叫Geolocation array.常用的数据有国外的MODIS数据,国内的如风云系列(FY)和海洋系列卫星(HY)都是使用这种方式来进行处理. GDAL库读取MODIS数据(HDF格式)的时候,会自动读取里面的经纬度查找表,并构造一个叫做GeoLocation的元数据段.如下所示: 只要有这个信息,使用GDAL就可以来对数据进行校正,但是对于国内的风云系列以及海洋系列卫星的数据来说,在数据里面会有

一.首先在百度api注册获得ak密钥 二.新建js文件,我命名为loadBMap.js,里面创建script,代码如下: /** * 加载地图 * @param {Function} callback 回调函数名称,回调函数将会挂载到window上,例如:window.initBaiduMapScript,在方法中就能拿到BMap对象 */ export function loadBMap(callback) { var script = document.createElement('scri

全球长时间序列中尺度气旋数据(http://cioss.coas.oregonstate.edu/eddies/)有netCDF格式,散点数据类型,只有一个很大的维Nobs = 2590938.尝试读取了部分经纬度和track数据,并绘制散度图.脚本程序: fn = 'D:/Temp/nc/tracks.20130125.nc' f = addfile(fn) n = 1000 lon = f['lon'][0:n] lat = f['lat'][0:n] track = f['track'][

项目结构 源码:ImageInfo.zip 第一步:添加需要的架包metadate-extractor.jar 架包下载地址:https://code.google.com/p/metadata-extractor/downloads/list 或者去Maven仓库下载 http://search.maven.org/#search%7Cga%7C1%7Cmetadata-extractor 第二步:编写解析代码 package com.drew.metadata; import java.io

如果读取图片经纬度,使用 String latValue = exifInterface.getAttribute(ExifInterface.TAG_GPS_LATITUDE); String lngValue=exifInterface.getAttribute(ExifInterface.TAG_GPS_LONGITUDE); 直接读取经纬度的话会返回null 要 String latValue = exifInterface.getAttribute(ExifInterface.TAG_

通过第三方库exifread读取照片信息.exifread官网:https://pypi.org/project/ExifRead/ 一.安装exifreadpip install exifread 二.读取照片信息,以及根据经纬度通过百度地图API获取位置 import exifread import json import urllib.request # Open image file for reading (binary mode) f = open('001.jpg', 'rb')

1.获取秘钥 调用百度地图API实现得申请百度账号或者登陆百度账号,然后申请自己的ak秘钥.链接如下:http://lbsyun.baidu.com/apiconsole/key?application=key 2.调用API将经纬度信息解析成json信息 3.可以自行通过json格式选择自己想要的数据,比如国家.省份.市区等. 4.批量转换信息 先要批量读取经纬度信息,我是将经纬度信息存在excel表中,然后再通过调用openpyxl中的函数读取经纬度信息. 读取完经纬度信息之后,再for循环

示例数据: #格式为txt文本 ID,L,B 001,116.5154,45.1154 001,116.5160,45.1153 ... 002,xxx,xxx ... 目标:建立轨迹数据结构,即dict{key=ID1:value=[[L,B],[L,B].....],key=ID2:value=[[L,B],[L,B].....]},即一个轨迹由一个ID和一串轨迹点组成,其中ID即为Key,轨迹点为Value,格式为list 其中L为大地经度,B为大地纬度: vec_dict = {} ke

参考 https://www.ecmwf.int/assets/elearning/eccodes/eccodes2/story_html5.htmlhttps://confluence.ecmwf.int/display/OPTR/ecCodes%3A+GRIB+data+decoding+and+encoding+software+2018 基本解码流程 1. 指定打开方式(“读”或“写”),打开一个或多个GRIB文件; 2. 根据不同加载方式,加载一个或多个GRIB messages到内存

pygrib-2.0.3/docs/index.html 导入pygrib模块 >>> import pygrib 打开grib文件,获取grib消息迭代器 >>> grbs = pygrib.open('sampledata/flux.grb') 和打开正常的python文件对象一样,pygrib使用seek, tell, read, readline, 和close方法打开grib文件.唯一的不同的是,偏移量用grib消息衡量,而不是用字节. >>>