r语言 - st_intersection至spatialpolygon dataFrame无法正常工作

Question

我试图获得两个Shapefiles的交叉点(属于某些大都市地区边界内的人口普查区域(。我能够成功获得相交功能，但是当我尝试将sf_intersect的输出转换为pacialpolygonsdataframe时，我会得到错误：

" as_spatial中的错误(来自(：从功能类型转换 不支持SP_的SFC_CEOMETRY"

">

这是我的代码：

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相交的每个几何几何的几何形状。如果您在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博客上。

• 在Wikipedia上所做的事情的空间谓词和示例(带有一些很好的插图(

---

信用 @lbussett的描述st_intersect和st_intersection

Answer 2

转换为Spatial对象无法处理混合线和多边形。交叉路口后，您可以使用以下方式仅提取多边形(并丢弃任何行(：
st_collection_extract("POLYGON")

您的示例似乎不再失败了，因此我创建了一个新示例，将两个多边形与共享方面相交。这导致多边形和线的交点输出。在第二次尝试中，我在成功转换为Spatial对象之前，通过ST_COLLECTION_EXTRACT((函数将交叉点输送。

library(sf)
library(dplyr)
library(sp)
#Create some boxes
BoxA <- st_polygon(list(cbind(c(0,0,2,2,0),c(0,2,2,0,0))))
BoxB <- st_polygon(list(cbind(c(1,1,3,3,1),c(1,3,3,1,1))))
BoxC <- st_polygon(list(cbind(c(2,2,4,4,2),c(0,2,2,0,0))))
#Create a funny shaped union to help demonstrate the intersection issue
BoxAB <- st_union(BoxA,BoxB)
plot(BoxAB)
plot(BoxC,add=TRUE,border="blue")

与Intersect的示例多边形

#Intersect of BoxAB with BoxC results in a line and a polygon
BoxIntersects<-st_intersection(BoxAB,BoxC)
plot(BoxIntersects)

由多边形和线组成的交点

#back to spatial fails
SpatialVersionOfIntersects<-as(BoxIntersects,"Spatial")

.as_spatial中的错误(来自，cast，ids(：不支持从功能类型SFC_GEOMETRY转换为SP

#Intersect again, but this time extract only the polygons
BoxIntersects<-st_intersection(BoxAB,BoxC) %>% st_collection_extract("POLYGON")
#back to spatial suceeds!
SpatialVersionOfIntersects<-as(BoxIntersects,"Spatial")