Cách chuyển đổi một spatialpolygon thành SpatialPolygonsDataFrame và thêm một cột vào bảng thuộc tính

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 

Ở trên, tôi lấy bất kỳ danh mục nào có 4 và trong cột mới đặt nó là 0.

Tại thời điểm này, tôi muốn cột cũng được đặt tên ratingscho bc, trong đó nó sẽ lấy 0 nếu nó nằm trong bộ đệm và 1 nếu nó ở bên ngoài. Vấn đề là các bclà SpatialPolygonsvà nó không chứa các bảng thuộc tính.

Rõ ràng để thêm một cột vào một SpatialPolygonđối tượng bạn phải chuyển đổi nó thành một SpatialPolygonsDataFrame, nhưng tôi không biết làm thế nào.

Tôi đã thử điều này:

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

nhưng lỗi này bật lên:

row.names of data and Polygons IDs do not match 

Các đối tượng "bờ biển", "ro" và "bc" có liên quan gì đến vấn đề của bạn? Vấn đề có thể nằm ở chỗ bạn đang sử dụng "readShapeSpatial". Bạn đã thử readOGR trong rgdal chưa? Nếu bạn đang đọc một shapefile đa giác, readOGR sẽ dẫn đến một đối tượng SpatialPolygonsDataFrame.

Nếu bạn trong thực tế, hãy có một đối tượng SpatialPolygons và muốn ép buộc vào SpatialPolygonsDataFrame, khung dữ liệu được chỉ định sẽ cần phải có các tên gọi phù hợp với ID đa giác trong khe đa giác. Dưới đây là một ví dụ nhanh.

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

Thử:

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

match.ID tránh yêu cầu của các tên gọi cho ID đa giác khớp

Nó khá đơn giản:

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" 

Khi Lỗi: "hàng. Tên dữ liệu và ID đa giác không khớp" phát sinh giải pháp này có vẻ hữu ích: đổi tên ID của khung dữ liệu để khớp với ID của đa giác:

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

Tôi tìm thấy giải pháp sau đây thường hoạt động.

Đầu tiên, tạo một khung dữ liệu trống với ID là trường:

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

Sau đó lấy ID của đa giác không gian 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

Sau đó sử dụng df làm đối số thứ hai cho hàm chuyển đổi khung dữ liệu không gian:

spatial_df <- SpatialPolygonsDataFrame(bc, df)

Vì dfvà spatial_dflà các đối tượng khung dữ liệu, các cột có thể dễ dàng được thêm vào