1

我在 R 中工作，尝试使用此处找到的数据（https://datadryad.org/resource/doi:10.5061/dryad.dk1j0；两个顶级文件）创建一个类似于此的表：[administrative_name, GDP2010, GDP2011 ....., GDP 2015]

据我所知，我需要从“admin_areas_GDP_HDI.nc”文件中提取行政单位的名称，并将它们与 GDP_per_capita_PPP_1990_2015.nc 文件中的年度数据结合起来。

使用 ncdf4 包，我设法打开档案，并获取所有属性和变量，但是我不知道如何访问数据并提取它。

我整天都在尝试访问数据，但我对 NetCDF 存档的经验有限，并且没有设法提取数据。任何指针都会帮助我！

4

1 回答 1

1

我喜欢使用raster包来处理 NetCDF 文件。它使用ncdf4包来读取文件，但提供了一些额外的工具来处理栅格。您没有提到要提取哪些数据，因此下面的示例显示了每个行政单位的平均 GDP。

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)
于 2019-04-05T15:19:59.317 回答