0

我想在 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)

我该如何解决?

感谢您的任何意见和建议。

4

1 回答 1

0

这是否给了您您的期望?

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

于 2018-11-29T21:00:03.607 回答