[英]How do you convert from a multipolygon to a spatial polygons data frame
你如何在R
从MULTIPOLYGON
转换为SpatialPolygonsDataFrame
?
我在网上找不到任何其他资源,我试图离散化作为GEOMETRYCOLLECTION
、 XY
和sfg
类的同一个文件,但是,这导致了在不访问单个点的情况下引用同一个多边形的连续循环。 我很高兴提供任何额外的说明,并感谢任何见解。
显示格式的多面Ref_circles
输出之一的示例如下:
Ref_circles[1]
[[1]]
[[1]][[1]]
[,1] [,2]
[1,] 51.62730 4.340600
[2,] 51.62549 4.343550
[3,] 51.61800 4.357353
[4,] 51.61124 4.371529
[5,] 51.60523 4.386039
[6,] 51.59998 4.400845
...
[311,] 51.63570 4.322473
[312,] 51.62894 4.336649
[313,] 51.62730 4.340600
我对您的第一个问题是:为什么您要从sf
对象(MULTIPOLYGON)转换为sp
对象(SpatialPolygonDataFrame),因为sf
取代了sp
?
可能有一些方法可以在sf
库中实现您的最终目标,而无需进行此转换。
如果你还想这样做,它会沿着
library(sf)
library(sp)
## using a MULTIPOLYGON data set supplied with library(sf)
nc <- sf::st_read(system.file("shape/nc.shp", package="sf"))
## convert the geometry of the `sf` object to SpatialPolygons
spd <- sf::as_Spatial(st_geometry(nc), IDs = as.character(1:nrow(nc)))
class(spd)
# [1] "SpatialPolygons"
# attr(,"package")
# [1] "sp"
## grab the data from the sf object
df <- nc
df$geometry <- NULL
df <- as.data.frame(df)
## create the SpatialPolygonsDataFrame
spd <- sp::SpatialPolygonsDataFrame(spd, data = df)
class(spd)
# [1] "SpatialPolygonsDataFrame"
# attr(,"package")
# [1] "sp"
head(spd@data)
# AREA PERIMETER CNTY_ CNTY_ID NAME FIPS FIPSNO CRESS_ID BIR74 SID74 NWBIR74 BIR79 SID79 NWBIR79
# 1 0.114 1.442 1825 1825 Ashe 37009 37009 5 1091 1 10 1364 0 19
# 2 0.061 1.231 1827 1827 Alleghany 37005 37005 3 487 0 10 542 3 12
# 3 0.143 1.630 1828 1828 Surry 37171 37171 86 3188 5 208 3616 6 260
# 4 0.070 2.968 1831 1831 Currituck 37053 37053 27 508 1 123 830 2 145
# 5 0.153 2.206 1832 1832 Northampton 37131 37131 66 1421 9 1066 1606 3 1197
# 6 0.097 1.670 1833 1833 Hertford 37091 37091 46 1452 7 954 1838 5 1237
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.