下载 WorldPop 人口统计网格

Question

我有兴趣从 WorldPop 下载 2020 年联合国调整后的刚果民主共和国人口统计，分辨率为 30 弧秒（参见 here ) as a GeoTIFF (.tif) file. I am trying to do this for multiple countries to I was exploring the wpgpDownloadR 包。

devtools::install_github("wpgp/wpgpDownloadR")
library(wpgpDownloadR)
wpgpListCountryDatasets(ISO3="COD")

我可以直接下载特定国家（例如刚果民主共和国、尼日利亚、比利时、法国等）的 WorldPop 人口统计 .tif 文件吗？

我有一个 R 代码，用于将此 GeoTIFF 栅格（.tif 文件）转换为 NetCDF（.nc 文件），如果我已经在本地预下载了 tif 文件，但如果我想要一个不同的国家/地区，我必须下载并重新执行此操作。

library(raster)

# Downloaded from above website
DRCWorldPop <- raster("cod_ppp_2020_1km_Aggregated_UNadj.tif") 

rnc <- writeRaster(DRCWorldPop, filename ='DRC_0000.nc', format = "CDF",  varname = "Susceptible", varunit = "Persons", longname = "Susceptible", overwrite = TRUE)

我正在尝试使用一个通用代码来使用 ftp 下载 任何国家 的光栅。你能给我指出正确的方向吗？谢谢你的时间。

Answer 1

你可以这样做

library(raster)
library(wpgpDownloadR)

iso3 <- "RWA"
cov <- "ppp_2000"
path <- dirname(tempdir())

f <- file.path(path, tolower(paste0(iso3, "_", cov, ".tif")))
if (!file.exists(f)) {
    f <- wpgpGetCountryDataset(ISO3 = "RWA", covariate = "ppp_2000", destDir=path) 
}
r <- raster(f)

这允许您遍历国家和协变量。它还可以避免下载您已有的文件。

Answer 2

我在下载 WorldPop 栅格时遇到了与 Ash 相同的问题，我认为指出 wpgpGetCountryDataset() 函数中存在错误可能会有所帮助。该函数将'covariate'参数转换为小写，我认为这导致它在后续调用其他内部函数时无法匹配必要的数据集，从而导致下载失败。

这仅适用于同时具有大写和小写文本大小写的协变量（例如，对于包含字符串 UNadj 的协变量，它会转换为 unadj）。这会导致错误消息，例如：

Entered Covariates: ppp_2020_1km_aggregated_unadj not present in WP

作为临时解决方法，我通过注释掉函数中的行 covariate <- tolower(covariate) 成功下载了数据。在 Github.

上作为 issue 提交

下载 WorldPop 人口统计网格

Downloading WorldPop Population Count rasters

tiff

aggregate

r

raster

netcdf