Question

我想在R中加载以下地理空间文件：ftp://ftp.nodc.noaa.gov/pub/data.nodc/icoads/1930s/1930s/ICOADS_R3.0.0_1930-10.nc。问题在于，即使我应该获得三个维度，使用后续代码也只能获得一个维度：

require("raster")
require("ncdf4")

nc_data <- nc_open("ICOADS_R3.0.0_1930-10.nc")

id.array <- ncvar_get(nc_data, "ID")
dim(id.array)

我该如何解决？

感谢您的任何评论和建议。

Answer 1

这能给您您所期望的吗？

library(tidync)
library(magrittr)

tfile <- tempfile(fileext = ".nc")
download.file("ftp://ftp.nodc.noaa.gov/pub/data.nodc/icoads/1930s/1930s/ICOADS_R3.0.0_1930-10.nc", tfile) 

id <- tidync(tfile) %>% activate("ID") %>% hyper_tibble()
dim(id)

[1] 69779 3

tidync仅在Github上：https://github.com/hypertidy/tidync

地理空间NetCDF中的尺寸错误

1 个答案: