Jak przekonwertować wielokąt przestrzenny na SpatialPolygonsDataFrame i dodać kolumnę do tabeli atrybutów

coast<-readShapeSpatial("coastline.shp")
landc<-readShapeSpatial("landcover.shp")
ro<-readShapeSpatial("roads.shp")
bc<-gBuffer(ro,width=100)
landc$ratings=1
landc$ratings[landc$LANDUSE_ID==4]=0 

Powyżej biorę dowolną kategorię, która ma 4, i w nowej kolumnie wpisuję ją jako 0.

W tym momencie chcę również nazwę kolumny ratingsdla bc, gdzie zajmie 0, jeśli jest w buforze i 1, jeśli jest na zewnątrz. Problem polega na tym, że bcjest SpatialPolygonsi nie zawiera tabeli atrybutów.

Oczywiście, aby dodać kolumnę do SpatialPolygonobiektu, musisz przekonwertować ją na obiekt SpatialPolygonsDataFrame, ale nie wiem jak.

Próbowałem tego:

buf_df<-as.data.frame(bc)
s_po<-SpatialPolygonsDataFrame(bc,buf_df)
s_po$ratings=0

ale pojawia się ten błąd:

row.names of data and Polygons IDs do not match 

Co mają wspólnego z twoim problemem obiekty „coast”, „ro” i „bc”? Problem może polegać na tym, że używasz „readShapeSpatial”. Czy próbowałeś readOGR w rgdal? Jeśli czytasz plik kształtu wielokąta, readOGR spowoduje utworzenie obiektu SpatialPolygonsDataFrame.

Jeśli faktycznie masz obiekt SpatialPolygons i chcesz przekonać się do SpatialPolygonsDataFrame, określona ramka danych będzie musiała dopasować nazwy runa do identyfikatorów wielokątów w gnieździe wielokątów. Oto szybki przykład.

library(sp)

# create some SpatialPolygons with ID's "2" and "3"
( p <- SpatialPolygons(list(Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "2"),
     Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "3"))) )
class(p)    

# Create a dataframe and display default rownames
( p.df <- data.frame( ID=1:length(p)) ) 
rownames(p.df)

# Try to coerce to SpatialPolygonsDataFrame (will throw error)
p <- SpatialPolygonsDataFrame(p, p.df) 

# Extract polygon ID's
( pid <- sapply(slot(p, "polygons"), function(x) slot(x, "ID")) )

# Create dataframe with correct rownames
( p.df <- data.frame( ID=1:length(p), row.names = pid) )    

# Try coersion again and check class
 p <- SpatialPolygonsDataFrame(p, p.df)
 class(p) 

# Now we can add a column
p@data$ratings <- 1:2 

# Or modify an existing one
p[p$ID < 2 ,] <- 5

Próbować:

#Code taken from the question:
s_po <- SpatialPolygonsDataFrame(bc, buf_df, match.ID = F) 

match.ID unika wymogu nazw rown dla identyfikatora wielokątów dopasowania

To jest dość proste:

library("rgdal")
polygons <- readOGR('path_to/file.shp',
                      layer = 'file')
class(polygons)
>[1] "SpatialPolygonsDataFrame"
>attr(,"package")
>[1] "sp"

poly_df <- as.data.frame(polygons)
# do some staff with "poly_df" that doesn't support SpatialPolygonsDataFrame
# then convert it to SPDF back again
s_poly <- SpatialPolygonsDataFrame(polygons, poly_df)
# add new column to SPDF:
s_poly$new_column <- "some data" 

Gdy pojawia się błąd: „Wiersz.nazwy danych i identyfikatory wielokątów nie pasują”, to rozwiązanie wydaje się pomocne: zmień nazwy identyfikatorów ramki danych, aby dopasować identyfikatory wielokątów:

newdata <- data.frame(whatever you want in here)
row.names(newdata) <- (however the new polygons are labeled)
polygons <- SpatialPolygonsDataFrame(polygons, newdata)

Uważam, że następujące rozwiązanie ogólnie działa.

Najpierw utwórz pustą ramkę danych o polu ID jako:

df <- data.frame(ID=character(), stringsAsFactors=FALSE )

Następnie uzyskaj identyfikatory wielokąta przestrzennego bc:

for (i in bc@polygons ) { df <- rbind(df, data.frame(ID=i@ID, stringsAsFactors=FALSE))  }
# and set rowname=ID
row.names(df) <- df$ID

Następnie użyj df jako drugiego argumentu do funkcji konwersji przestrzennej ramki danych:

spatial_df <- SpatialPolygonsDataFrame(bc, df)

Ponieważ dfi spatial_dfsą obiektami ramki danych, kolumny można łatwo dodawać