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

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

terra r-raster raster r

katefull06 · 技术社区 · 8 月前

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

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 回复 | 直到 8 月前

Andrea Manica 8 月前

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

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

Robert Hijmans 8 月前

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

这是一个变通办法 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')

推荐文章

Amp · 使用R ggplot2删除geom_radial中axis.line和panel.border之间的空格

8 月前

Hard_Course · 用另一列中的值替换行的最后一个非NA条目

8 月前

Mark R · 使用geom_sf()删除地球仪上不需要的网格线

8 月前

Joe · 根据对工作日和本周早些时候的日期的了解,找到一个日期

8 月前

Ben · 统计向量中的单词在字符串中出现的频率

8 月前

TheCodeNovice · R中符号格式的尾随零和其他问题[重复]

8 月前

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

8 月前

dez93_2000 · 在R管道子功能中引用管道对象的当前状态

8 月前

accibio · 在ggplot2中为同一变量创建两个连续的颜色渐变比例

8 月前

Mankka · 如何在Ggplot2中绘制均匀的径向图

8 月前