绘制grib2格式大气数据(windows下 GrADS wgrib2 g2ctl perl)

一. 环境配置 1.下载安装ActivePerl.exe编译环境

(或者StrawberryPerl编译环境也可以)         (包含命令：perl, 若cmd中不能使用这些命令，请配置对应目录到path环境变量中)         下载路径：https://www.activestate.com/products/activeperl/downloads/

2.下载g2ctl.pl文件

(是处理grib2的数据的  grib2ctl.pl是处理grib1的数据的)         下载路径：ftp://ftp.cpc.ncep.noaa.gov/wd51we/g2ctl/g2ctl         帮助说明：https://www.cpc.ncep.noaa.gov/products/wesley/g2ctl.html          (前两步为了用.grib2文件生成对应的.ctl文件，供GrADS绘图使用)    

3.安装OpenGrADS2.1.a2

(含有wgrib2.exe，可用)，为了大气数据的绘图。(包含命令：wgrib2,gribmap,grads, XLaunch.exe; 若cmd中不能使用这些命令，请配置对应目录到path环境变量中)         下载路径：https://sourceforge.net/projects/opengrads/         3.0 配置环境变量(如D:\ProgramFiles\OpenGrADS\Contents\Resources)到path,为了找到绘制的笔刷等文件         3.1  配置启动 安装目录下的 XMing https://jingyan.baidu.com/article/6079ad0e9a348528ff86dba8.html         3.1' (或者 安装XMing运行环境,并确保XMing运行)

4.附加信息

（使用可单独运行的工具wgrib2.exe，读取查看grib2数据）         下载路径: ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/Windows_64/         帮助文档：ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/intro_grib2.pdf         (下载安装 wgrib2.exe 以及相应的dll放在同一个文件夹下，配置该路径到环境变量path)

 

二. 绘制步骤

在grib数据文件目录下执行：(g2ctl.pl文件 数据文件、ctl文件、index文件最好都放入同一个目录中

1. 生成ctl管理文件

cmd窗口下执行 "perl g2ctl.pl mydata.grib2 >mydata.ctl" 生成GrADS的管理文件ctl 特别注意: .ctl文件的编码方式不能是 UCS-2 LE BOM ,可以是UTF-8,以便GrADS识别，可能需要手动修改文件的编码格式，

2. 生成索引文件idx

cmd窗口下执行 "gribmap –v –i mydata.ctl",生成索引文件 mydata.idx,以便GrADS软件使用

3. 启动XLaunch.exe

(路径如:D:\ProgramFiles\OpenGrADS\Contents\Resources\Xming\XLaunch.exe)

4. 绘制大气数据图形 4.1 命令窗口交互方式绘图

            cmd下运行"grads"命令 进入"ga-> " 交互窗口，输入命令绘图             常用命令：                 "open mydata.ctl"      打开数据管理文件                 "q file"                           显示文件中的变量（或者称“要素”）信息                 "d 变量名称"                 绘制指定变量对应的图形

4.2 文件脚本形式绘图

            将上述命令 放入文件中如 (draw_my_data.gs文件)             cmd或者使用代码调用 "grads -cl draw_my_data.gs"             grads命令参数说明及帮助：可使用"grads -help"可查看

 

三.参考链接

Grads 处理 【GRIB】 格式文件的准备：http://blog.sina.com.cn/s/blog_86cf6ab60101qujf.html利用grads处理Grib2数据的过程体会：http://bbs.06climate.com/forum.php?mod=viewthread&tid=88918GrADS基础绘图：https://wenku.baidu.com/view/180c01e6793e0912a21614791711cc7931b778f6.htmlGrADS绘图软件使用手册：https://max.book118.com/html/2018/0523/168074252.shtmGrads 里地图消失的问题：http://blog.sina.com.cn/s/blog_4e1fe80801013cmk.html 及 http://blog.sina.com.cn/s/blog_4d5c3ba901000au9.html#comment1

四. 其他 wgrib2命令

打印grib数据时，精确到分钟和秒： 

.\wgrib2.exe mygrb2file   -S 

获取第一行的详细数据，并打印输出到csv文件：

 .\wgrib2.exe mygrb2file -d 1 -csv line1.csv   仅仅获取温度数据：.\wgrib2.exe mygrb2file  -match ":(TMP:200 mb):"

以目录中的"行"为单位，逐行输出具体数据 格式为.asc栅格格式：.\wgrib2.exe mygrb2file  -AAIG     输出指定行的具体数据 格式为.asc栅格格式：     如: .\wgrib2.exe mygrb2file -d 1 -AAIG     如：.\wgrib2.exe mygrb2file  -match ":(TMP:200 mb):" -AAIG