0

我正在尝试裁剪和屏蔽栅格数据并将栅格数据转换为矢量和带有 terra 的 sf 对象。但是，我总是收到以下错误消息： UseMethod 中的错误("st_as_sf")not applicable method for 'st_as_sf'applied to object of class "SpatVector"。对于克服这一障碍的任何帮助或建议，我将不胜感激。

我的代码如下所示

packs <- list("sf", "tidyverse", "terra")
lapply(packs, require, character.only=TRUE)

pop_count <- rast(raster.path)
adm_ner_sf <- st_read(shp.path)

pop_afg_sf <- pop_count %>%
  crop(vect(adm_afg_sf)) %>%
  mask(vect(adm_afg_sf)) %>%
  as.points() %>%
  st_as_sf() %>%
  rename(pop_count = gpw_v4_population_count_rev11_2010_2pt5_min)

我正在使用 SEDAC 人口密度栅格数据（正式类 SpatRaster），我试图通过使用阿富汗形状文件 (adm_afg_sf) 将其限制为阿富汗的形状。

任何反馈都会非常受欢迎。

4

1 回答 1

0

我认为您所做的与下面的最小自包含可重现示例等效，并且有效：

library(terra)
library(sf)

f <- system.file("ex/elev.tif", package="terra")
r <- rast(f)
p <- as.points(r)
s <- st_as_sf(p)

也许您有旧版本的sf和/或terra？尝试update.packages()

于 2021-08-01T14:53:10.730 回答