spatialpolygonをSpatialPolygonsDataFrameに変換し、属性テーブルに列を追加する方法

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 

上記では、4のカテゴリを取り上げ、新しい列に0を入力しています。

この時点で、私が指定したカラムたいratingsためだけでなくbc、それは外であれば、それはそれは、バッファ内にある場合は0を取り、1になります。問題があるということbcでSpatialPolygons、それは属性テーブルが含まれていません。

明らかにSpatialPolygonオブジェクトに列を追加するにはSpatialPolygonsDataFrame、それをに変換する必要がありますが、その方法はわかりません。

私はこれを試しました：

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

しかし、このエラーはポップアップします：

row.names of data and Polygons IDs do not match 

「海岸」、「ro」、および「bc」オブジェクトは、問題とどのような関係がありますか？問題は、「readShapeSpatial」を使用していることにある可能性があります。rgdalでreadOGRを試しましたか？ポリゴンシェープファイルを読み込んでいる場合、readOGRはSpatialPolygonsDataFrameオブジェクトになります。

実際、SpatialPolygonsオブジェクトがあり、SpatialPolygonsDataFrameに強制的に変換する場合は、指定されたデータフレームの行名がポリゴンスロットのポリゴンIDに一致する必要があります。以下に簡単な例を示します。

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

試してください：

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

match.ID 一致するポリゴンIDの行名の要件を回避します

とても簡単です：

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" 

エラー：「データの行名とポリゴンIDが一致しない」が発生した場合、この解決策が役立つと思われます。データフレームのIDをポリゴンのIDに一致するように変更します。

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

次の解決策が一般的に機能することがわかります。

最初に、IDをフィールドとして空のデータフレームを作成します。

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

次に、空間ポリゴンのIDを取得します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

次に、dfを空間データフレーム変換関数の2番目の引数として使用します。

spatial_df <- SpatialPolygonsDataFrame(bc, df)

ようdfとspatial_dfデータフレームのオブジェクトであり、列は簡単に追加することができます