在R中处理grib2文件的RNomad花费的时间太长
在R中处理grib2文件的RNomad花费的时间太长
提问于 2018-05-03 12:53:51
我们正在尝试读取一个包含13个变量的grib2，天气预报在一个气压水平上持续6个纬度/经度，持续一天。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])
}回答 1
Stack Overflow用户
发布于 2018-06-27 20:19:09
您可以尝试在地理兴趣范围内定义domain参数：
domain=c(85.8000,95.8000,52.7200,36.8200)
### read data
data <- ReadGrib(file, variable = var_list, levels = lev_list,domain=domain)页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/50147178
复制相关文章
相似问题
腾讯云开发者
