使两个栅格具有相同的维度

问题描述

我有两个栅格（第一个是从 www.worldpop.org 下载的网格化人口计数 GeoTIFF 栅格，另一个是 GADM 1 级数据（这是国家/省/地区的索引），它是转换为 R 中的栅格）。

假设我在从 here 下载的 30 弧秒分辨率 tif 文件中选择了以色列和该国 2020 年联合国调整后的人口数。

Router

library(raster)
library(rgdal)
library(terra)

rm(list = ls())

#----------------------------------------------------------------#
# Source 1: WorldPop UN-Adjusted Population Count GeoTIFF raster #
#----------------------------------------------------------------#

# Import the GeoTIFF file into R workspace

url <- "https://data.worldpop.org/GIS/Population/Global_2000_2020_1km_UNadj/2020/ISR/isr_ppp_2020_1km_Aggregated_UNadj.tif"

fname <- basename(url)

if (!file.exists(fname)) download.file(url,fname,mode="wb")

WorldPop <- raster("fname")

WorldPop <- replace(WorldPop,is.na(WorldPop),0)

WorldPop

dim(WorldPop)

plot(rast("fname"),main = "2020 UN-Adjusted Population Count for Israel",col=rev(rainbow(10,end=0.7)),breaks=c(0,10,25,50,100,250,1000,100000))

inhabitableRaster <- WorldPop
values(inhabitableRaster) <- ifelse(values(WorldPop) > 0,1,0)
names(inhabitableRaster) <- "Inhabitable"
inhabitableRaster

dim(inhabitableRaster)
table(as.matrix(inhabitableRaster))

plot(inhabitableRaster,main = "Inhabitable Status (0 = uninhabitable / 1 = inhabitable)")

#---------------------------------------#
# Source 2: From GADM: Level1Identifier #
#---------------------------------------#

Level1Identifier <- getData("GADM",level=1,country="ISR")
Level1Identifier$NAME_1 # List of all states/provinces/regions

Level1Raster <- raster(Level1Identifier,resolution = res(WorldPop)[1])
Level1Raster <- rasterize(Level1Identifier,Level1Raster) 

Level1Raster <- crop(Level1Raster,WorldPop)
names(Level1Raster) <- "Level1Identifier"
table(values(Level1Raster))

plot(Level1Raster)

我希望 Level1Raster 与 WorldPop 栅格具有相同的维度。这怎么办？

我可以将 Level1Raster 的原点设置为 WorldPop 的原点，但这不会将 GADM 栅格的维度（nrow x ncols）更改为 匹配 WorldPop 人口计数栅格。

crs(Level1Raster)
crs(WorldPop)

origin(Level1Raster)
origin(WorldPop)

extent(Level1Raster)
extent(WorldPop)

res(Level1Raster)
res(WorldPop)

dim(Level1Raster)
dim(WorldPop)

我也尝试了以下

origin(Level1Raster) <- origin(WorldPop)

dim(Level1Raster)
dim(WorldPop)

以上代码行将 WorldPop 的列/行从 (474,200) 更改为 Level1Raster 的 (464,196)。但是，我希望反过来，我希望 Level1Raster 维度与 WorldPop 相同。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）