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")