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

悠悠 2024-04-18 15:46 15阅读 0赞

# 一. 环境配置 #

## 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=88918  
**GrADS基础绘图**：https://wenku.baidu.com/view/180c01e6793e0912a21614791711cc7931b778f6.html  
**GrADS绘图软件使用手册**：https://max.book118.com/html/2018/0523/168074252.shtm  
**Grads 里地图消失的问题**：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**

**![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3NoeWpoeXAxMQ_size_16_color_FFFFFF_t_70][]**

**![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3NoeWpoeXAxMQ_size_16_color_FFFFFF_t_70 1][]**

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3NoeWpoeXAxMQ_size_16_color_FFFFFF_t_70]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/18/b379d79c4a504215b5760295dcc461c8.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3NoeWpoeXAxMQ_size_16_color_FFFFFF_t_70 1]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/18/f89dff715010475192b9a00f472cebd2.png