如何解决裁剪和屏蔽 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()
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。