这确实是我无法解释的奇怪行为。
library(terra)
r1 = rast(nrow = 2, ncol = 2)
values(r1) = c(0, 1, 2, 3)
r2 = r1
ext(r2) = ext(r2)/2
ext(r2)
ext(r1)
您会看到“r1”的范围也减少了一半。为什么?顺便说一句,您可以更改“r2”的值而不影响“r1”的值。据我所知,它仅适用于范围。
使用“栅格”包执行相同过程时也不会发生这种情况。
library(raster)
r1 = raster(nrow = 2, ncol = 2)
values(r1) = c(0, 1, 2, 3)
r2 = r1
extent(r2) = extent(r2)/2
extent(r2)
extent(r1)