{{tag>R GIS windows}}
====== RでGIS ======
GISというか地図表示
[[http://www.okada.jp.org/RWiki/?ShapeFile%A5%E9%A5%A4%A5%D6%A5%E9%A5%EA|RjpWikiの記事]]を参考（というかほとんどそのまま）．下記データを使用した場合に書き換えただけ．

<del>よりmaptoolsに特化した記事は以下．
http://www.okada.jp.org/RWiki/?maptools%20%A4%C7%C3%CF%BF%DE%C9%BD%BC%A8</del>無くなってた．

http://www.okada.jp.org/RWiki/?Rmap%A4%F2%BB%C8%A4%C3%A4%BF%C3%CF%BF%DE%C9%BD%BC%A8にRを使った地図表示についてまとめられている．

===== 使用パッケージ =====
[[http://cran.r-project.org/web/packages/maptools/index.html|maptools]]
と[[http://cran.r-project.org/web/packages/sp/index.html|sp]]
===== データの入手 =====
[[http://www.globalmap.org/download/kanni01.html|みんなの地球地図ホーム ＞ ダウンロードするには（簡易版）＞ Shape形式データ]]
http://www.globalmap.org/download/data/shp/1_1/shape_jp_1_1.zip
を「C:/GIS_mapdata/」に解凍した場合を記す．
すなわち，「C:/GIS_mapdata/shape_jp_1_1」に色々データが入っている状態．

===== 表示 =====
行政域を表示
  library(maptools)
  setwd("C:/GIS_mapdata/shape_jp_1_1")
  getwd()
  
  jpn <- read.shape("bnda_1_1.shp")
  jpn_poly <- Map2poly(jpn)
  plot(jpn_poly)

行政域以外の地図を表示させたい場合，適宜
  jpn <- read.shape("bnda_1_1.shp")
の.shpファイルを変えればよい．

===== maptoolsでShapeFileを切り出す方法 =====
日本全体は不要なので，北海道だけの.shpファイルが欲しい．
  library(maptools)
  setwd("C:/GIS_mapdata/shape_jp_1_1")
  getwd()
  
  jpn <- read.shape("bnda_1_1.shp")
  jpn_poly <- Map2poly(jpn)
  #plot(jpn_poly)
  jpn_df <- jpn$att.data
  
  hokkaido_poly <- subset(jpn_poly,jpn_df$nam=="HOKKAIDO")
  hokkaido_df <- subset(jpn_df,jpn_df$nam=="HOKKAIDO")
  
  write.polylistShape(hokkaido_poly,hokkaido_df,"hokkaido")

北海道だけのシェープファイルを作成後，
  library(maptools)
  setwd("C:/GIS_mapdata/shape_jp_1_1")
  hokkaido <- read.shape("hokkaido.shp")
  hokkaido_poly <- Map2poly(hokkaido)
  
  plot(hokkaido_poly)
で北海道だけが表示される．

札幌市ならsubsetの条件をlaa=="SAPPORO-SHI"にすればよい．