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

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

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

我还有一张英国的土地覆盖图，是光栅形式的

我需要查看土地覆盖图中每个10km正方形的土地覆盖%s

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

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

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

如何更改土地覆盖图的范围，使其与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

库（光栅）
#连续景观
展示景观（续）
#加载粗分辨率光栅
g_光栅● 如果粗略数据不是网格，则设置“is_grid=FALSE”

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

您可以检查范围是否与重叠

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")

库（光栅）
ext.ras此错误消息表明您可能存在投影问题，但如果没有更多信息，很难知道

您可以检查范围是否与重叠

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")

库（光栅）
分机ras