RによるGISのテーブル結合

ここを参考に。
元shpデータ内テーブルと結合したいテーブル内にマッチングさせるIDが必要です。

まず、もとshpをインポート
A_shp=readShapePoly("A.shp", IDvar="OBJECTID")

テーブルの確認
A_shp@data
summary(A_shp)

結合したいテーブルデータをインポート
B_table=read.table(“B_table.txt”,header=T)

元shpを上書きしたくないので新しいデータをコピー
new_shp=A_shp

A_IDとB_IDを紐付けて結合
A_shp@data=data.frame(A_shp@data,table[match(A_shp@data$A_ID,B_table@data$B_ID),])

こんだけ！

ちなみにmaptoolsのspCbindではエラーが出て失敗しました。↑のやり方のほうがラクです。