问题标签 [netcdf]

在使用 Unidata IDV 通过 OPeNDAP 访问加尔维斯顿湾的 NOAA CO-OPS 海洋模型时：

http://opendap.co-ops.nos.noaa.gov/thredds/catalog/NOAA/GBOFS/MODELS/201302/catalog.html?dataset=NOAA/GBOFS/MODELS/201302/nos.gbofs_hsc.fields.nowcast。 20130201.t09z.nc

我收到一条错误消息“创建数据源时出错”。从他们的OPeNDAP 数据表中可以看出，许多NetCDF 气候和预报 (CF) 元数据约定都得到了满足，但显然不是全部。

有谁知道缺少什么或错误？

我正在检查 R 中的一批 netCDF 文件的元数据。由于 SO 似乎不允许文件上传，我在下面附上了一个简化示例 netCDF 文件的标题：

基本上，我试图将变量的几个属性，例如，“no2”的“standard_name”、“long_name”和“units”，一次放入一个向量中。但是，att.get.ncdf()在 Package 中ncdf似乎一次只允许一个。由于通过重复代码获得更多属性可能相当乏味和费力，我认为必须有一种更简洁的方法来做到这一点，但我只是没能在谷歌上找到它。有人可以和我分享他们的经验吗？

非常感谢！

摘要：我可以在 a 上显示经纬度图lattice::layerplot，并且可以在 a 上显示 Lambert 保角圆锥 (LCC) 图spam::image。如何在 上显示 LCC 地图lattice::layerplot？如何不这样做的示例如下 - 修复方面的帮助表示赞赏（甚至只是调试）。

细节：

我一直在使用格子图形（通过latticeExtra和rasterVis）成功地显示未投影的经纬度全球大气数据，效果很好（尽管我肯定对ggplot/很感兴趣ggmap）。特别是，我可以用地图覆盖这些地块，这对于我正在做的工作来说是必不可少的。但是，我目前无法lattice::layerplot用于某些区域数据投影 LCC：数据绘图，但我无法获得要叠加的地图。我可以通过更粗略的方式做到这一点，但更愿意知道如何在lattice/rasterVis或类似的情况下做到这一点（例如，ggplot/ggmap）。下面是两个几乎独立的示例……但是如果您已经知道如何执行此操作，请告诉我，我将跳过调试。（我在一个项目上落后了。）

netCDF 数据 ,ozone_lcc.nc随CRAN package=M3 ... 除了 M3 将其提供为

并且该文件扩展名 ( .ncf) 当前会导致CRAN package=raster出现问题（请参阅此帖子）。您可以重命名该文件（并将其放入某个当前工作目录），或者您可以只下载该文件(270 kB)，或者您可以从M3 压缩包中获取它（但请记住重命名它！）。

然后，您可以运行以下任何示例（前提是（IIRC）您没有运行 Windows，package=M3不会在其上构建（但 ICBW）），根据需要更改常量以适应您的系统。示例 1 生成了一个我知道（根据以前的经验）将与 a rasterin a一起使用的类型的映射levelplot；但是，在这种情况下，地图和数据/栅格的坐标值不匹配。示例 2 使用旧式基本图形，实际上绘制了数据和地图；不幸的是，我不知道如何制作它在levelplot. 因为我想让这段代码与一堆使用rasterand的其他代码一起工作，所以levelplot这是一个问题。

示例 1：产生类似的输出

示例 2：产生类似的输出

但是我看不到如何从返回的简单矩阵（等，但不多）中获得想要M3::get.map.lines.M3.proj 的SpatialLines东西sp::sp.lines，更不用说想要的东西了latticeExtra::layer。（我是一个新手，发现格子文档相当难以理解。）此外，我更愿意避免手动进行上面的 IOAPI 转换（尽管我当然更愿意这样做而不是跳过旧式图形的箍）。

我需要处理 netcdf 文件中的单个变量，该文件实际上包含许多属性和变量。我认为不可能更新 netcdf 文件（请参阅问题如何删除 Scientific.IO.NetCDF.NetCDFFile 中的变量？）

我的方法如下：

1. 从原始文件中获取要处理的变量
2. 处理变量
3. 将所有数据从原始 netcdf 但处理后的变量复制到最终文件
4. 将处理后的变量复制到最终文件

我的问题是编写第 3 步的代码。我从以下内容开始：

如何在少量代码中保存所有数据和属性，而无需重建整个数据结构？

我对 PyPy 很陌生，我使用 PyPy1.9 来解释我的 Python 代码。

如您所知，CPython 解释器可以导入打包在 .so 文件中的模块，该文件由 C/C++ 代码编译而成。

但是PyPy1.9解释器告诉我他在.so文件中找不到模块，我的代码有点慢，我想优化性能，求助！！

我有几个 netcdf 文件。每个 nc 文件都有几个变量。我只对两个变量“Soil_Moisture”和“Soil_Moisture_Dqx”感兴趣。

我想根据“Soil_Moisture_Dqx”过滤“Soil_Moisture”。每当相应的“Soil_Moisture_Dqx”像素的值大于 0.04 时，我想用 NA 替换“Soil_Moisture”中的值。

:这里是要下载的文件：

1- 我尝试了这个循环，但是当我输入 f[1] 或 f[2] 时，我得到了一些奇怪的东西，这意味着我的循环不正确。我很感激任何帮助我纠正我的循环。

2- 谁能告诉我再写一遍？

我们将海洋模型的每小时输出存储在一系列 netcdf 文件中，每个月一个。

我们将每月的第一个小时和最后一个小时存储在每个文件中。我们想使用 NcML 聚合这些文件，但我们不想在聚合中获得重复的时间值。

有没有办法做到这一点？

我知道有像wgrib2这样的软件可以将文件grib和grib2格式转换为NetCDF文件，但我需要另辟蹊径： from NetCDFto grib2，因为这里的当地气象局只能使用grib2格式的网格数据。

似乎一种解决方案可能是在 Python 中，使用NetCDF4-Python 库（或其他）来读取NetCDF文件并使用pygrib来编写grib2.

有没有更好的办法？

我在 GrADS 工作以生成模型运行图（netcdf 文件），我对这个程序非常陌生。我正在尝试找出如何将多个图像发送到元数据（所有来自动画的图像），而不必一次处理一个。

另一种选择是减慢屏幕上的动画速度。

尽我所能，我无法弄清楚如何让 plyr 工作。对于这个特定示例的任何帮助，以及任何解释您的示例为何有效的奖励积分，将不胜感激。

数据在这里，以防有人想让这个例子“可行”

设置：

此 for 循环有效，但速度慢且不并行：

我想做这样的事情：

当我运行它时，我收到以下错误消息：

看起来 plyr 正在将整个向量传递给函数，而不是单个值！非常感谢您对如何进行这项工作的建议以及对修复工作原因的解释！

提前致谢！