问题标签 [ncdf4]

我有一个将 nc 文件的一部分转换为 csv 文件的脚本。脚本本身可以工作，但问题是我需要指定确切的目录，包括文件名和输出 csv。我有兴趣为文件夹 test1 中的所有 nc 文件运行脚本，并将其转换为文件夹 test2 中同名的 csv。我尝试修改脚本，但没有奏效。这是我的脚本。

我修改脚本的失败尝试是

以前，它本来是

其中 example1 是 nc 文件名，exampleresult 是 csv 名称。

我正在尝试从简单海洋数据同化 (SODA) 2.2.4 访问 NetCDF 文件。我想将它们下载为栅格图层，包含温度、年份、深度、纬度和经度数据。

这些文件位于 DOD（分布式海洋数据系统）/OPeNDAP 中，因为它们太大而无法作为 netcdf 文件下载。我以前没有使用过这种格式，因为我一直都能下载 .nc 文件。在该站点上，它说我应该能够像使用 R 脚本中的 .nc 文件一样使用此链接，以便转换为栅格图层并下载。但是，当我尝试这样做时，它失败了。但是，我可以在 R 中查看有关文件的属性，如下所示。

该文件在我的计算机上不存在，但应该存在于此链接中？帮助？

关于将 .nc 文件转换为 .tiff 文件，我遇到了丢失像素地理信息的问题。我知道其他用户也遇到了同样的问题，并试图通过 kotlin 解决但失败了。我更喜欢使用 R 的解决方案。有关 kotlin 方法 URL，请参见此处：https ://gis.stackexchange.com/questions/259700/converting-sentinel-3-data-netcdf-to-geotiff

我下载了 ESA 的免费 Sentinel-3 数据（网址：https ://scihub.copernicus.eu/dhus/#/home ）。不幸的是，此数据以 .nc 格式提供，因此我想将其转换为 .tiff 格式。我已经尝试了各种方法，但都失败了。到目前为止我已经尝试过：

我非常感谢任何保留像素地理信息的解决方案！

非常感谢，Explorer

我正在尝试ncdf使用rCurl. 任何人都可以就为什么这不起作用提供任何建议吗？

R_nc4_open 中的错误：NetCDF：未知文件格式

nc_open(files[1]) 中的错误：

nc_open 尝试打开文件 A20021722014263.L3m_SCSU_NSST_sst_4km.nc 时出错

我正在尝试将文件转换为.nc文件以.csv在 R 中进行进一步分析，因为我习惯于使用.csv.

基本上我认为要解决我的问题（下面有更多详细信息），我需要_FillValue在.nc文件中添加一个，但我尝试过的一切都不起作用。

按照http://geog.uoregon.edu/bartlein/courses/geog490/week04-netCDF.html#replace-netcdf-fillvalues-with-r-nas.nc中采取的步骤，我已经成功地为许多文件做到了这一点，直到第 3.4.3 节。

但是，我最近获得了对另一个.nc文件的访问权限，并且相同的过程无法正常工作。我想我已经把它缩小到_FillValue新.nc文件中没有的事实。

看起来_FillValue应该是“9.97e+36”。我尝试使用将此数字添加为缺失值

这似乎添加missing_value:9.97e+36到.nc文件中。但是，当我运行时：tmp_array <- ncvar_get(ncin,dname)tmp_array 仍然有 9.97e+36。

我希望 tmp_array 已经将 9.97e+36 替换NA为它对它工作的文件所做的那样。

有没有办法可以将 _FillValue 添加到我的文件中，以便用 _FillValue 替换这些值NA？

如果需要，这是不工作的文件的信息：

我想计算 R 中每月降雨量栅格的标准化降水指数 (SPI)。我在 SO（R 标准化降水指数 .nc 文件）中看到了一个与此相关的问题，但它没有任何答案。所以，我提出一个新问题。我从以下来源下载了每月降雨量数据 https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.03/cruts.1905011326.v4.03/pre/

我已经使用raster包在 R 中阅读了它。

我知道它可以使用包的calc功能来完成，raster并且对于 SPI 计算，R 中有许多包可用（例如spi，、、SPEI）precintcon。我能够计算点数据的 SPI。但我无法使用该calc函数进行栅格数据的 SPI 计算。现在我的问题是如何在不同的时间尺度（1、2、3、6、12、24）上计算这个栅格上的 SPI？

我有一个从 1979 年到 2019 年具有 3 个维度（经度、纬度、时间）的月度温度 NetCDF 文件。我从这个网站下载了该文件：https ://cds.climate.copernicus.eu/cdsapp#!/dataset /reanalysis-era5-single-levels-monthly-means?tab=form。我想创建一个函数，允许我选择纬度高于 23 度的所有网格单元，对于这些网格单元，只需使用例如 12 月、1 月、2 月的月份来计算年平均值。任何想法将不胜感激。

到目前为止，我只是设法单独执行任务，但我正在寻找更多类似的东西：

我正在寻找类似的东西：

到目前为止，我只是设法对每年的这些月份进行了子集化，但现在我需要合并上述条件。太感谢了！

我正在尝试安装一个名为“ncdf4”的 R 包。我尝试将其安装在 R-studio 界面的“包”部分，但也尝试进入install.packages("ncdf4")控制台。

这是我在两种情况下得到的输出：

有人能够从我提供的代码中破译出什么问题吗？谢谢。

我正在尝试.nc从文件中写入.csv文件。

这是我.csv桌子的一部分 ( df)

如您所见，我有两个维度和一个变量。

这是我正在使用的代码：

当我运行最后一行时，我收到以下错误消息：

有关我的尺寸的更多信息：

yvals独特的价值观：

xvals独特的价值观：

并且elev列有 959 个唯一值。

另外，我不确定重新排序尺寸是否会给我long,lat和之间的正确对应关系elev。无论如何，我完全不知道如何解决这个问题。

感谢大家！

我有这个 shapefile，我从中获取前两列 (left和top) 来访问另一个数据库并获取elevation.

我访问了数据库并编写了一个 netCDF 文件（StackOverflow Question）。我使用了相同的 959 个地理点。

但是，当我与我一起绘图时，plot(RasterLayer)我得到了这个：

为什么栅格图有这种差距，为什么没有与 shapefile 图相同的形状？