这个问题已经以不同的标题多次问了，老实说，我没有遇到任何想要的解决方案。
我已下载ecmwf上的控制预测数据，其中包含6个小时（时间步长）的NCEP模型（来源和时间）的2019年7月（月）的总降水量（参数）。数据采用Grib2格式。我尝试了几种不同的方法，但都没有成功以数字形式提取每天的预测数据。我可以看到有31个波段，我认为这意味着总共31天。以下是一个起始代码（我没有包括使用控制台测试其他功能时使用的所有其他命令）。我尝试使用rNOMADS和raster功能，但没有成功。 rNOMADS软件包需要级别，我不确定在哪里以及如何从我的数据中找到这些级别。这是我正在处理的link到Grib数据。
rm(list = ls())
graphics.off()
install.packages("rNOMADS")
library(weathercan)
library(rNOMADS)
library(rgdal)
library(raster)
library(ncdf4)
library(fields)
AreaOfInterest = extent(c(-104.8, -99.54, 47.65, 50.52))
model.data <- readGDAL("PCP6hours.grib")
Info = GDALinfo("PCP6hours.grib", silent = TRUE)
Row = Info[["rows"]]
z=geometry(model.data)
plot(z)
我想要什么？
我只想将AreaofInterest的每日天气预报保存到data.frame中。