r - R - terra::distance() 等效于 raster::gridDistance(..., origin = x, omit = y)

Question

更新- 这个问题与 相关terra 1.4-1，现在已经过时terra 1.5-12（或更早？）。

---

我正在寻找terra相当于raster::gridDistance(..., origin = my_origin, omit = my_omit).

terra::gridDistance我从terra v0.2-8 这里发现了一个旧网页，但我可以收集terra::distance到的是当前的替换raster::gridDistance（替换功能列表terra's）。

但是，我不知道omit = my_omit在terra::distance. 从文档页面看来，任何非 NA 都被视为origin，但没有提及 ，也没有omit将 更改origin为特定值的选项（与raster::gridDistance下面的示例不同）。

这是来自的示例raster::gridDistance：

library(raster)

# world lon/lat raster
r <- raster(ncol=10,nrow=10, vals=1)
r[48] <- 2     # this will be the origin
r[66:68] <- 3  # this will be the area that can't be traversed
plot(r)

d <- gridDistance(r,origin=2,omit=3) 
plot(d)

任何人都可以使用复制这个例子terra::distance吗？

我可以使用 做我需要做的事情raster，但我仍在学习如何正确处理空间栅格数据，并且正在努力学习新terra包。

Answer 1

感谢@lovalery 在评论中指出这一点。

在我使用的原始问题中terra 1.4-1。截至 2022 年 1 月，terra 1.5-12已发布（2022 年 1 月 13 日），它现在包含一个与该功能terra::gridDistance()非常相似的raster::gridDistance()功能。

为了完整起见，这是来自的示例?terra::gridDistance：

#world lon/lat raster
r <- rast(ncol=10,nrow=10, vals=1)
r[48] <- 2
r[66:68] <- 3
d <- gridDistance(r,origin=2,omit=3) 
plot(d)

#UTM small area
crs(r) <- "+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
d <- gridDistance(r,origin=2,omit=3) 
plot(d)