我有3个栅格,在我再次重新分类之前,我必须乘以它们进行重新分类,但在我进行乘法之前,我必须使3个栅格的范围相同。距离光栅
的xmax是不同的ext(dem)
ext(slope)
ext(distance)
ext(dem)
SpatExtent : 1089831.75, 1721631.75, 4748091, 5516991 (xmin, xmax, ymin, ymax)
ext(slope)
SpatExtent : 1089831.75, 1721631.75, 4748091, 5516991 (xmin, xmax, ymin, ymax)
ext(distance)
SpatExtent : 1089831.75, 1721031.75, 4748091, 5516991 (xmin, xmax, ymin, ymax)
我已经尝试了其他问题的答案,但仍然失败,我已经尝试了以下的坡度和dem的范围内种植距离
distance <- crop(slope, distance)
distance <- crop(dem, distance)
但是在运行它并再次检查ext之后,它仍然没有修复它。我检查了投影,它们也都是一样的。
您不知道show(dem)
和show(distance)
,因此很难给出具体的答案。我们所能看到的是distance
的x范围比其他光栅的x范围小600 (m)。在将来,请通过在代码中创建spatraster(参见as.character(dem)
开始)提供一个可重复的示例,如
library(terra)
dem <- rast(ncol=, nrow=, xmin=, ....)
crop
只能删除整个行和列(您不能拥有栅格单元的一小部分)。我怀疑你所期望的程度不符合那个要求。但它可能会让你更接近。如果非常接近(也就是说,差异远远小于空间分辨率,并且行数和列数匹配),则可以执行
ext(distance) <- ext(slope)
否则,考虑resample
distance2 <- resample(distance, slope)