0

我们正在尝试读取包含 13 个变量的 grib2 天气预报,其中一个压力水平为 6 lat/long 一天。grib2 文件为 62mb。要读取一天的每小时数据变量,大约需要一个小时,尤其是 ReadGrib。

我们正在使用 wgrib.exe,版本 v2.0.3。我们在 Windows 10 中使用 R 版本 3.4。

如果阅读需要这么长时间,还是我们做错了什么。请指教。

我们使用的代码是:

file <- "./xxx.grib2"

### define the geographical coordinated to be extracted
lat_list <- c(36.8200, 43.2800, 52.7200)
lon_list <- c(85.8000, 87.3500, 95.8000)

var_list <- c("APCP", "DSWRF", "PRES")
lev_list = c("surface")

### read data
data <- ReadGrib(file, variable = var_list,  levels = lev_list)

for (i in 1: length(lat_list)) {
   cat("\n", "         ... extracting point",  i,  "/",  length(lat_list))
   data_extract <- BuildProfile(data, spatial.average = F,
                    lat = lat_list[i], lon = lon_list[i])
}
4

1 回答 1

0

您可以尝试在地理兴趣范围内定义域参数:

domain=c(85.8000,95.8000,52.7200,36.8200)

### read data
data <- ReadGrib(file, variable = var_list,  levels = lev_list,domain=domain)
于 2018-06-27T12:19:09.040 回答