代码之家 › 专栏 › 技术社区 › katefull06

在R中使用terra修改范围时,会为单独的SpatRaster重写范围

terra r-raster raster r

katefull06 · 技术社区 · 6 月前

这真是奇怪的行为,我无法解释。

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”的范围也减少了一半。为什么??? 顺便说一句,您可以在不影响“r1”值的情况下更改“r2”的值。据我所知,它只适用于某种程度。

使用“光栅”包执行相同过程时也不会发生这种情况。

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)

1 回复 | 直到 6 月前

Andrea Manica 6 月前

一个简单的任务 terra 创建浅拷贝,这意味着如果编辑派生对象,则可以更改原始对象。如果你想避免这种行为,解决方案是使用 deepcopy .文件 深度拷贝 显示了类似的行为 time 正如你所发现的 extent : https://rspatial.github.io/terra/reference/deepcopy.html

编辑 :注意@Robert Hijmans的回答,他说 dev 版本 terra 现在表现不同,不需要 深度拷贝

Robert Hijmans 6 月前

这是一个现已在开发版本中修复的错误。

这是一个变通办法 as.vector() :

library(terra)
r1 = rast(nrow = 2, ncol = 2)
values(r1) = c(0, 1, 2, 3)
r2 = r1
ext(r2) = as.vector(ext(r2)/2)
ext(r2)
#SpatExtent : -90, 90, -45, 45 (xmin, xmax, ymin, ymax)
ext(r1)
#SpatExtent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)

在开发版本中,它按预期工作

library(terra)
#terra 1.8.24
r1 = rast(nrow = 2, ncol = 2)
values(r1) = c(0, 1, 2, 3)
r2 = r1
ext(r2) = ext(r2)/2
ext(r2)
#SpatExtent : -90, 90, -45, 45 (xmin, xmax, ymin, ymax)
ext(r1)
#SpatExtent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)

您可以使用安装开发版本

install.packages('terra', repos='https://rspatial.r-universe.dev')

推荐文章

dwiz · R中具有焦函数和parlappy的光栅并行处理列表

7 年前

krishnab · 谷歌地球引擎:寻找正确的图像尺寸以获得卫星图像的最大分辨率

7 年前

Hugo · 具有随机森林的图像分类(光栅堆栈)(package ranger)

7 年前

B. Davis · 创建包含特定光栅值的多边形

8 年前

MapUser · 如何从光栅中提取多边形?

8 年前

corvus · R: 从density()获取数字输出

9 年前

Joke O. · 使用R中的savitzky golay/sgolayflt(信号)平滑rasterstack

9 年前

N.Varela · R: 如何通过在基本R图中用光栅覆盖地图来获得正确的分辨率(光栅)?

9 年前

Geo-sp · 在水平图面板区域外添加文本

9 年前

user32882 · 将光栅文件作为阵列处理?

11 年前