问题标签 [ncdf4]

我正在尝试在 r 中读取多个 nc4 文件。下面是我用来执行这个任务的代码；

但我一直收到这个错误：

R_nc4_open 中的错误：没有这样的文件或目录。

但是文件列表被正确识别为 chr [1:1440]，如全局环境值中所示。

有人可以帮我解决我做错了什么吗？

我需要在 R 中创建一个三维数组，其中包含分辨率为 538x907 像素的栅格数据。我一个月内每个小时都有这个光栅，所以一月份有 744 个光栅文件。我必须通过 R 更改一些值，然后想将它们汇总回一个可以由包 ncdf4 处理的数组。因此，我需要创建一个看起来像的三维数组prectdata[1:538, 1:907, 1:744]（第一个和第二个用于光栅的 x 和 y 维度，第三个用于时间维度）。我如何将 744 个光栅矩阵连接到包 ncdf4 的三维数组？

我想使用函数“filled.contour”，但我的数据文件是 ncdf4 文件。

如何将ncdf4文件转换为R中的data.frame？

一些忠告？

谢谢

我需要在 R 3.4.0 上安装 ncdf4 包，但它失败并告诉我错误“找不到 -lnetcdf”。如何成功安装 ncdf4。

我是 NetCDF 文件的新手。我有来自 NOOA 的AVHRR Pathfinder 数据版本 5.3的 1981-2014 年全球海面温度的每日数据。该数据集由每年的 365x2 .nc 文件组成：给定年份中每一天的夜间和日间温度各一个。

我希望将每年的各种每日 .nc 文件合并为一个每年的年度 .nc 文件。这样做最有效的方法是什么？我已经nco使用 Anaconda 探索了这个包，但我遇到了问题。任何帮助将不胜感激！

所以我一直在 R 中使用以下教程来习惯 NetCDF（我使用 netcdf4 包）：http ://geog.uoregon.edu/bartlein/courses/geog607/Rmd/netCDF_01.htm

但是，我无法在我的数据集中为任何特定年份创建漂亮的图像，可以在这里找到：https ://www.ncdc.noaa.gov/paleo-search/study/19419

这是我的代码，我在其中尝试获取公元 100 年的地图：

特别是，当我尝试运行最后一行时，我收到以下错误消息：

image.default（lon，lat，drried.slice，col = rev（brewer.pal（10，：z 的尺寸不是长度（x）（-1）乘以长度（y）（-1））中的错误

有人可以解释我做错了什么吗？

顺便说一句，为了让您对我的数据结构有所了解，请看：

我正在尝试从 .nc 文件中提取数据。由于我的文件中有 7 个变量，我想使用 foreach 循环 ncvar_get 函数通过所有 7 个变量。

这是我的代码：

错误是：

你能解释一下这段代码有什么问题吗？我认为这与集群找不到从文件中获取哪个变量这一事实有关，因为“从 C 调用返回的错误”通常来自 ncvar_get 变量参数。

当我尝试使用 writeOGR 和循环来保存我的 shapefile 时，它​​什么也没做，只是给我一条错误消息：

writeOGR(plot.locationsSP_DROUGHT, dsn, layer1, driver = "ESRI Shapefile") 出错：层存在，使用新层名

本质上，我将每个对象转换为 CSV 文件，然后转换为 shapefile，并希望同时保存 CSV 文件和 shapefile。这是我的代码片段：

这是我正在使用的完整代码：

帮助将不胜感激。我可能正在做一些非常愚蠢的事情。

我有一个R 网格文件，其中包含 1981 年的月度温度数据，我读入并尝试使用以下代码写入 NetCDF：

这会写入NetCDF 文件，但它似乎只有一个月（我不确定是哪个月），而不是我使用 panoply 检查它时一年中的 12 个月。

是否可以编写 NetCDF 文件并尽可能多地保留 R-grid 文件中的数据？尤其是每个月的数据！

编辑：

新的工作代码：

我有一个以下格式的栅格堆栈，我使用此方法将其转换为 NetCDF。这可行，但纬度和经度变量以“米”为单位，这是光栅文件的范围，我需要它们以十进制度为单位。R Grid文件在这里，格式如下：

理想情况下，我想将其转换为具有以下格式的 NetCDF 文件以供模型读取：

我尝试使用此方法将光栅砖范围转换为纬度/经度：

但这给出了这个错误：

编辑：也试过

但是当我在 panoply 中检查这个时，纬度和经度坐标都是 0.0。

我不知道从这里去哪里。