R中两个栅格的匹配范围？

Question

我正在尝试使用R中的“grainchanger”包将分辨率更高的光栅聚合为分辨率更高的光栅。

我有10公里的英国网格轮廓，我已经从shapefile转换为光栅。

我还有一张英国的栅格土地覆盖图。

我需要查看土地覆盖图中每10公里见方的土地覆盖%s。

当我尝试使用grainchanger包聚合此内容时，会出现以下错误：

Error in .local(x, y, ...) : extents do not overlap

我想这是因为我的10km栅格光栅是正方形，而土地覆盖地图没有填充。

如何更改土地覆盖地图的范围，使其与10km栅格光栅匹配？

这是我的代码：

library(raster)

# continuous landscape
show_landscape(cont_ls)

# load the coarse resolution raster
g_raster <- raster("raster.tif", package = "grainchanger")

show_landscape(g_raster)

# direct aggregation using range
dda <- nomove_agg(coarse_dat = g_raster,
              fine_dat = cont_ls, 
              agg_fun = var_range)
#> aggregation assumes all cells are rectangular
#> ● set `is_grid = FALSE` if coarse_dat is not a grid

Answer 1

此错误消息表明您可能存在投影问题，但如果没有更多信息，很难知道。

您可以检查范围是否重叠

library(raster)
ext.ras <- extent(6e+05, 709800, 5690220, 5800020)
ext.pol <- extent(499386.6, 517068.2, 6840730, 6857271)


plot(ext.ras, xlim = c( 499386.6,709800), ylim= c(5690220,6857271), col="red")
plot(ext.pol, add=T, col="blue")

摘自此答案。

如果它们没有重叠，并且您希望它们重叠，则可以使用spTransform（）更改投影。