我正在尝试获取两个 shapefile(位于某些大都市区边界内的人口普查区)的交集。我能够成功获得相交特征,但是当我尝试将 sf_intersect 的输出转换为 SpatialPolygonsDataframe 时,出现错误:
"Error in as_Spatial(from) : conversion from feature type sfc_GEOMETRY to sp is not supported"
这是我的代码:
library(sf)
library(dplyr)
library(tigris)
library(sp)
#download shapefiles corresponding to metro areas
metro_shapefiles<-core_based_statistical_areas(cb = FALSE, year = 2016)
#convert to sf and filter
metro_shapefiles<-st_as_sf(metro_shapefiles)%>%filter(GEOID==31080 )
#Data for California
census_tracts_california<-tracts(state="CA",year=2016)
census_tracts_california<-st_as_sf(census_tracts_california)
#INTERSECT AND CONVERT BACK TO SP
census_tracts_intersected1<-st_intersection(census_tracts_california,
metro_shapefiles)
#back to spatial
census_tracts_intersected1<-as(census_tracts_intersected1,"Spatial")
最佳答案
错误消息告诉您不能将 sfc_GEOMETRY
转换为 Spatial
对象。没有 sp
等效对象。
在您的交集结果中,您混合了多种几何图形(因此,您将返回一个 sfc_GEOMETRY
作为您的“几何图形”)。您可以在此处查看所有几何图形:
types <- vapply(sf::st_geometry(census_tracts_intersected1), function(x) {
class(x)[2]
}, "")
unique(types)
# [1] "POLYGON" "MULTILINESTRING" "MULTIPOLYGON"
如果需要,您可以提取每种类型的几何图形,并将它们分别转换为 SP:
lines <- census_tracts_intersected1[ grepl("*LINE", types), ]
polys <- census_tracts_intersected1[ grepl("*POLYGON", types), ]
spLines <- as(lines, "Spatial")
spPolys <- as(polys, "Spatial")
附加信息
我在评论中提到您可以使用 st_join
。但是,这可能不会给您想要的结果。 sf
库中有几何二元谓词,例如 ?st_intersects
,以及几何运算,例如 ?st_intersection
谓词返回一个稀疏(默认)或密集矩阵,告诉您 x 的每个几何图形与 y 的哪个几何图形相交。如果您在 st_join
中使用它,它将返回相交的(原始)几何图形,而不是稀疏矩阵。
而操作(例如 st_intersection
)将计算交点,并返回新的几何图形。
使用示例
谓词 (st_intersects
) 可以在 st_join
中使用,它们将返回“相交”的原始几何图形
sf_join <- sf::st_join(census_tracts_california, metro_shapefiles, join = st_intersects)
在这种情况下,这给出了一个单一的type
对象
types <- vapply(sf::st_geometry(sf_join), function(x) {
class(x)[2]
}, "")
unique(types)
# [1] "MULTIPOLYGON"
## so you can convert to a Spatial object
spPoly <- as(sf_join, "Spatial")
但是您需要确定 st_intersect
的结果是否是您所追求的结果,或者您是否需要 st_intersection
给出的新几何图形。
进一步阅读
感谢用户@lbussett 对 st_intersect
和 st_intersection
之间差异的描述
关于r - st_intersection 到 spatialpolygon 数据框不起作用,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/48796828/