问题标签 [netcdf4]

我在一个文件夹中有许多 NCDF 文件。我尝试使用raster和ncdf4包将它们提取到光栅砖中。如果我分别提取每个 NCDF 文件，它就可以工作。但是，我尝试使用 for 循环提取所有文件，然后它给了我错误。

as.integer(nl) 中的错误：无法将类型“S4”强制转换为“整数”类型的向量

如果我分别看，r它们r1似乎具有相同的程度，并且都是光栅砖类型：

请帮忙。

我不得不在我的 Linux（Ubuntu 15.10）上更新一些 Anaconda/Python 包，结果我现在收到以下我以前从未见过的错误消息：

在 Python 脚本中调用基于 Fortran 的程序时会发生该错误。但是，在没有 Python 的情况下运行相同的程序时，一切正常...

更新后我也重新编译了程序，但没有任何改变。

有谁知道这里发生了什么？

我有一个带有以下内容的 netCDF 文件。ncdump：

如何提取data_array组中的最后一个变量？如果我使用ncks，我可以像这样提取所有内容data_array：

ncks -v data_array test_nc.nc output_nc.nc

但是，我只想提取对应于 N=3 in 的变量data_array。任何使用 python netCDF4 或 nco 工具的解决方案都对我有用（但不使用 cdo），也需要在 Windows 上工作。

我正在尝试从数据集模型编写 Netcdf-4 文件，遵循此http://www.unidata.ucar.edu/software/thredds/current/netcdf-java/tutorial/NetcdfWriting.html导致以下代码：

但是它产生的只是Permission denied 错误。有了这个块，我希望从填充的数据集中编写 .nc4 文件。别的不知道怎么写。我正在使用 Netcdf nc_inq_libvers=4.4.1.1

谢谢！

使用python netcdf4库，我想为我的 netcdf4 读取脚本编写一个测试数据集。但是，我无法生成所需的输出。

这是目前我的写脚本：

但这会返回一个TypeError：

TypeError：需要一个整数

我应该如何更改输入或编写脚本以获得所需的结果？

我有一个 netCDF 文件 (.nc)，其中包含 16 年（1998 - 2014）的每日降水量（5844 层）。3 个维度是时间（大小 5844）、纬度（大小 19）和经度（大小 20）在 R 中是否有一种直接的方法来计算每个光栅单元：

• 月平均和年平均
• 累积比较（例如 1 月至 3 月与所有 1 月至 3 月的平均值相比）

到目前为止，我有：

我的第一个挑战是计算每个栅格单元的月平均值。我不确定在牢记最终目标（累积比较）的同时如何最好地进行。如何轻松访问某个月份的几天？

希望我把我的问题说清楚了，我们将不胜感激朝正确方向迈出的第一步。此处的示例数据

我在组合多个 for 循环时遇到问题。我将举两个例子，我想结合起来。如果我知道如何用两个来做，我也可以用多个循环来做。

如果有人知道如何将其编写为 lapply 函数，那也很好。

我如何尝试结合但没有奏效：

我已经解决了我的问题。见下文。但我仍然很好奇如何解决上述两个 for 循环，所以我将把这个问题悬而未决。

这是我的解决方案：

我认识到这是我的安装失败，很抱歉将这个无趣且无关紧要的问题放在你的脚下，但对于我的生活，我无法弄清楚出了什么问题，我已经没有想法。我希望有人能够迅速指出显而易见的事情。

我正在尝试分析一个 python 脚本（使用 Kern 的 line_profiler），并且该脚本需要加载 netCDF4 模块。我已经用 pip 安装了 line_profiler 和 netCDF4。当我为已安装软件包列表排队 pip 时，两者都报告为存在并更新。

在不使用分析器的情况下，我的脚本运行没有问题，这意味着 netCDF4 模块已正确加载。但是，如果我从“myscript”目录运行“kernprof -l -v myscript.py”，我会收到以下错误：

我从 /opt/local/bin/python 的安装中运行 Python，它在我的 PATH 中首先列出。

所以，无论如何，如果我设置的默认 Python 版本与我的 PATH 中最先出现的版本相同，并且该默认版本能够访问 netCDF4 模块，为什么 line_profiler 不是呢？

对于我的脑筋急转弯，我已经在文档和邮件列表档案中搜索了一段时间，并且很难将处理此聚合所需的步骤放在一起。

CFSR 1 小时数据文件数据来自这里：http ://rda.ucar.edu/datasets/ds094.0/

每小时文件声明了 2 个时间维度，一个设置了边界，另一个没有设置。

问题是，当我逐步创建每个数据集时，不同的每小时文件会将名称交换为 2 个时间维度名称。因此AggregationExisting无法找到某些文件的维度名称“时间”，例如在 0300 文件中的 u-component_of_wind_isobaric 变量上，因为它被声明为 time1。

我正在调用的代码：

时间维度名称问题在 2 个文件中说明：

这是我第一次使用 NetCDF 库，所以我正在购买一些预处理工具来合并这些具有这种怪癖的数据集。例如，我可以将所有变量移动到同一时间维度并重命名吗？即使是指向我错过的示例的链接也会有所帮助。否则我猜我会考虑手动删除尺寸并使用 readDataSlice() 手动将数据复制到新的合并文件中。

我试图通过将其变量的时间分辨率除以 2 来减小 netcdf 文件的大小，方法如下：

代码运行没有错误，并且新文件确实包含时间维度为原始时间维度一半的变量。

原始文件的大小是 1.1 Go，但新文件的大小是 1.4 Go。我希望生成一个大小约为原始文件一半的文件，因为我将时间分辨率减半。我不确定这是怎么发生的。

你能解释一下吗？