我正在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存档的经验有限,并且还没有设法提取数据。任何指针都可以帮助我!
答案 0 :(得分:0)
我喜欢使用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)