แปลงวัตถุรูปหลายเหลี่ยมเชิงพื้นที่เป็นกรอบข้อมูลโดยใช้ R

เป้าหมายของฉันคือการแก้ไข shapefile ที่มีอยู่โดยการผสานรูปหลายเหลี่ยมบางอย่าง

หลังจากนำเข้ารูปร่างไฟล์และใช้คำสั่ง UnionSpatialPolygons ฉันได้รับโครงร่างรูปหลายเหลี่ยมที่ฉันต้องการ

อย่างไรก็ตามตอนนี้เป็นวัตถุ SpatialPolygons ไม่ใช่ SpatialPolygonsDataFrame ดังนั้นฉันจึงไม่สามารถส่งออกไปยังรูปร่างไฟล์โดยใช้ writeOGR

ฉันจะแก้ไขปัญหานี้ได้อย่างไร

ดังที่ชื่อกล่าวไว้ SpatialPolygonsDataFrame นั้นเป็นเพียงวัตถุ SpatialPolygons ที่มีข้อมูลแนบอยู่ (ตารางแอตทริบิวต์) ข้อมูลจะต้องมีแถวอย่างน้อยที่สุดเท่าที่มีคุณสมบัติ

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

---- ---- แก้ไข

หากคุณต้องการแปลงSpatialPolygonsDataFrameกลับไปเป็นSpatialPolygonsวัตถุคุณเพียงแค่ต้องระบุโครงสร้างวัตถุภายใน R

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

ปัญหา:

1: ผลลัพธ์ของ UnionSpatialPolygons เป็นรูปหลายเหลี่ยมเชิงพื้นที่

2: การแปลงผลลัพธ์กลับเป็นกรอบข้อมูลรูปหลายเหลี่ยมเชิงพื้นที่เป็นความเจ็บปวดที่แท้จริง

-a คุณต้องการ data frame ที่แน่นอนมากเพื่อแนบไปกับรูปหลายเหลี่ยมเชิงพื้นที่

-b ข้อมูลที่คุณใช้สำหรับ UnionSpatialPolygons มีแถวมากกว่าเอาท์พุทและไม่ได้จัดรูปแบบตามที่ต้องการ

โซลูชันของฉัน (น่าเกลียด):

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