Извлечение данных из файла NetCDF - PullRequest
0 голосов
/ 05 апреля 2019

Я работаю в R, пытаясь использовать данные, найденные здесь (https://datadryad.org/resource/doi:10.5061/dryad.dk1j0; два верхних файла), чтобы создать таблицу, подобную этой: [административное_имя, GDP2010, GDP2011 ....., GDP 2015]

Насколько я вижу, мне нужно извлечь названия административных единиц из файла "admin_areas_GDP_HDI.nc" и объединить их с годовыми данными в файле GDP_per_capita_PPP_1990_2015.nc.

С пакетом ncdf4 мне удалось открыть архивы и получить все атрибуты и переменные, однако я не знаю, как получить доступ к данным и извлечь их.

Я пытался получить доступ к данным весь день, но у меня ограниченный опыт работы с архивами NetCDF, и мне не удалось извлечь данные.Любые указатели помогут мне!

1 Ответ

0 голосов
/ 05 апреля 2019

Мне нравится использовать пакет растр для работы с файлами NetCDF.Он использует пакет ncdf4 для чтения в файлах, но предлагает некоторые дополнительные инструменты для обработки растров.Вы не упомянули, какие данные вы хотите извлечь, поэтому в приведенном ниже примере показан средний ВВП для каждой административной единицы.

library(raster)

#Read in NetCDF files
ad -> brick('admin_areas_GDP_HDI.nc')
gdp -> brick('GDP_per_capita_PPP_1990_2015_v2.nc')

#Calculate mean GDP using admin zones
zoneMean -> zonal(gdp, ad[[1]], fun='mean', na.rm=T)
...