裁剪和屏蔽 Rasta 数据并将其转换为矢量

问题描述

我正在尝试裁剪和屏蔽栅格数据，并将栅格数据转换为矢量和带有 terra 的 sf 对象。但是，我总是收到以下错误消息：UseMethod ("st_as_sf") 中的错误不适用于 'st_as_sf' 的方法应用于类 "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) 将其限制为阿富汗的形状。

任何反馈将不胜感激。

解决方法

在询问 R 问题时，请使用 R 附带的数据提供一个最小的独立可重现示例（请参阅帮助文件），并且仅包含重现错误所需的代码，显示错误消息发生的位置（不要使用管道），并且不要使用 library(tidyverse) 来保持示例最小。

我认为你所做的与下面的等效并且有效：

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()

geospatial r r terra