Неправильные размеры в геопространственной NetCDF - PullRequest
0 голосов
/ 29 ноября 2018

Я хотел бы загрузить следующий геопространственный файл в 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)

Как мне это исправить?

Спасибо за любые комментарии и предложения.

1 Ответ

0 голосов
/ 30 ноября 2018

Дает ли это то, что вы ожидаете?

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

...