R аварийно завершает работу при открытии файла netcdf

Я скачал файл netcdf и пытаюсь открыть его в R. Вот мой код

download.file("https://data.giss.nasa.gov/impacts/agmipcf/agmerra/AgMERRA_1980_prate.nc4",destfile = "AgMERRA_1980_prate.nc4", method="libcurl")

Я хочу открыть файл netcdf с помощью R

library(ncdf4)
my.file <- nc_open("AgMERRA_1980_prate.nc4")

Однако каждый раз, когда я это делаю, R падает.

введите здесь описание изображения

Что-то не так в моем коде или что-то не так с R studio?

sessionInfo() R версия 3.5.0 (23 апреля 2018 г.) Платформа: x86_64-w64-mingw32/x64 (64-разрядная версия) Работает под: Windows 7 x64 (сборка 7601) Пакет обновления 1

Матричные продукты: по умолчанию

ИЗМЕНИТЬ

Если я вручную скачаю файл, я смогу его открыть. Так что должно быть что-то не так в том, как я загружаю его. Какие-либо предложения?


r raster netcdf4
person 89_Simple    schedule 18.08.2018    source источник

Ответы (2)


arrow_upward
0
arrow_downward

Я не уверен, что здесь происходит, возможно, это связано с Windows. Пробовал скачивать без аргумента method="libcurl" и вроде работает.

download.file("https://data.giss.nasa.gov/impacts/agmipcf/agmerra/AgMERRA_1980_prate.nc4",
destfile = "AgMERRA_1980_prate.nc4")



library(ncdf4)
my.file <- nc_open("AgMERRA_1980_prate.nc4")


File AgMERRA_1980_prate.nc4 (NC_FORMAT_NETCDF4):

     1 variables (excluding dimension variables):
        short prate[longitude,latitude,time]   (Chunking: [1440,720,1])  (Compression: level 9)
            _FillValue: 32767
            description: Precipitation Rate
            units: mm/day
            add_offset: 0
            scale_factor: 0.100000001490116
            vMin_original_data: 0
            vMax_original_data: 457.399993896484
            vRange: 457.399993896484

     3 dimensions:
        time  Size:366   *** is unlimited ***
            units: days since 1980 01-01-01 12:00:00
        latitude  Size:720
            units: degrees_north
        longitude  Size:1440
            units: degrees_east

    4 global attributes:
        history: Tue Aug 12 16:42:13 EDT 2014
        source: AgMIP / Alex Ruane
        title: AgMERRA v1.1 Precipitation Rate
        center: NASA GISS

Информация о моем сеансе -

R version 3.5.0 (2018-04-23)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS High Sierra 10.13.6
person Suhas Hegde    schedule 18.08.2018

arrow_upward
3
arrow_downward

Я подозреваю, что этот вопрос дублирует Загрузка файлов NetCDF с R: Вручную работает, download.file выдает ошибку. Предложение @Luis использовать mode = "wb" вместо mode = "w" по умолчанию помогло мне избежать сбоев nc_open() с R 4.0.2, RStudio 1.3.959 и ncdf 1.17. wb сообщает download.file(), чтобы обрабатывать файл как двоичный, соответствующий формату netCDF.

Для данных, представляющих интерес здесь, это будет

download.file("https://data.giss.nasa.gov/impacts/agmipcf/agmerra/AgMERRA_1980_prate.nc4", destfile = "AgMERRA_1980_prate.nc4", method = "libcurl", mode = "wb")
person Todd West    schedule 06.08.2020
comment
Очень полезный! С model = "wb" у меня работает любой из способов на винде. - person Danny; 21.07.2021