Chuyển đổi một đối tượng đa giác không gian sang khung dữ liệu bằng R

Mục tiêu của tôi là sửa đổi một shapefile hiện có bằng cách hợp nhất một số đa giác nhất định.

Sau khi nhập shapefile và sử dụng lệnh UnionSpatialPolygons, tôi nhận được phác thảo đa giác mà tôi muốn.

Tuy nhiên, đây hiện là đối tượng SpatialPolygons chứ không phải SpatialPolygonsDataFrame, vì vậy tôi không thể xuất nó sang shapefile bằng writeOGR.

Làm thế nào tôi có thể khắc phục vấn đề này?

Như tên đã nói, SpatialPolygonsDataFrame về cơ bản chỉ là một đối tượng SpatialPolygons với dữ liệu được đính kèm (bảng thuộc tính). Dữ liệu phải có ít nhất nhiều hàng như có các tính năng

library(rgdal)
ob <- SpatialPolygons(..)# Your SpatialPolygons Object
spp <-     SpatialPolygonsDataFrame(ob,data=as.data.frame("yourData"),proj4string=CRS("+proj=    aea > +ellps=GRS80 +datum=WGS84"))
writeOGR(spp,"shapes","testShape",driver="ESRI Shapefile",)

----BIÊN TẬP----

Nếu bạn muốn chuyển đổi SpatialPolygonsDataFramelưng của bạn thành một SpatialPolygonsđối tượng, bạn chỉ cần giải quyết cấu trúc đối tượng trong R

ob <- SpatialPolygons(spp@polygons,proj4string=spp@proj4string)

Các vấn đề:

1: kết quả của UnionSpatialPolygons là một đa giác không gian

2: chuyển đổi kết quả trở lại thành khung dữ liệu đa giác không gian là một nỗi đau thực sự

-a. bạn cần một khung dữ liệu rất chính xác để gắn vào đa giác không gian

-b. dữ liệu bạn đã sử dụng cho UnionSpatialPolygons có nhiều hàng hơn đầu ra và không được định dạng theo cách cần thiết.

Giải pháp (xấu xí) của tôi:

### Coerce into spatial polygon data frame with id and row name of spatial polygon

# Make a data frame that meets the requirements above:

df<- data.frame(id = getSpPPolygonsIDSlots(your.spatialpolygon))
row.names(df) <- getSpPPolygonsIDSlots(your.spatialpolygon)

# Make spatial polygon data frame
spdf <- SpatialPolygonsDataFrame(your.spatialpolygon, data =df)

# Then don't forget to make sure the projection is correct
# XXXX is your SRID

proj4string(spdf) <- CRS("+init=epsg:XXXX");
spdf <- spTransform(spdf , CRS("+init=epsg:XXXX"));