如何从 R 中的 .tiff 文件中提取数据？

Question

我对地理空间分析完全陌生，不知道如何开始。我想知道 2000 年越南的性别比例是多少，按网格划分为https://hub.worldpop.org/geodata/summary?id=11943.

如何在 R 中下载和读取此数据，并提取 .tiff 文件的每个网格方块的男性和女性人口？

此外，我如何汇总数据以找到每个地区 20 年持有权的性别比例？

非常感谢。

Answer 1

使用

terra

，下载您想要比较的年龄。由于服务器似乎受到限制，这会很慢。使用2000，女/男

library(terra)
vnm_f_2k_20 <- rast('~/Downloads/vnm_f_20_2000.tif')
vnm_m_2k_20 <- rast('~/Downloads/vnm_m_20_2000.tif')
# plot(vnm_m_2k_20) you could but not much to look at 
vnm_20_ratio <- vnm_f_2k_20/vnm_m_2k_20 # female/male
writeRaster(vnm_20_ratio, '~/Downloads/vnm_20_2k_fm_ratio.tif')
plot(vnm_20_ratio)
viet_nam <- vect('~/Downloads/boundaryjson.geojson')
plot(viet_nam)

你有 shapefile，我有 geojson，我们都需要弄清楚哪个级别的详细信息，是在省级还是区级，以及所有涉及到的小南沙群岛。

为了考虑应用分析的级别，您确实应该回顾一下

worldpop

数据以及它是在哪个级别收集的。这就是空间分析中“支持”的概念。如果 worldpop 是省级的，则应使用省级，如果是地区，则应使用地区。为了演示的目的，如果我能弄清楚我的 geojson，让我们做省级。