Question

我正试图获得两个shapefile的交集（人口普查区属于某些大都市区的边界）。我能够成功获得交叉功能，但是当我尝试将sf_intersect的输出转换为SpatialPolygonsDataframe时，我得到错误：

“as_Spatial（from）中的错误：从要素类型转换不支持sfc_GEOMETRY到sp“

这是我的代码：

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

Answer 1

错误消息告诉您无法将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给出的新几何图形。

进一步阅读

每个联接的信息都在sf blog上。
空间谓词以及不同操作在wikipedia上的作用示例（附有一些好的插图）

感谢用户@lbussett对st_intersect和st_intersection

之间差异的描述

st_intersection到spatialpolygon数据帧不起作用

1 个答案:

其他信息

使用示例

进一步阅读