Rで地理空間データ分析

事前準備

  • spdep,maptoolsをinstall.packages()
  • libraryで読み込み

基礎操作

  • .shpファイルの読み込み
building<-readShapePoly("building.shp")
  • ポリゴンデータの代表点を抽出
building.cor<-coordinate(building)
  • ポリゴンデータをプロット
plot(building,col="gray",border="black")

日本語データを読み込む時は注意 name<-read.csv(file("buildingname.csv",encoding="Shift_JIS"),stringsAsFactors=F,header=T)

  • テキストを図に追加する
par(family = "HiraKakuProN-W3")
text(building.cor[,1],building.cor[,2],name[,2],cex=0.5,col="black")

shpファイルの入手先

qiita.com

参考資料