10.2 地図データの入手とRでの表示

サンプルコード


install.packages("maptools")
library(maptools)

map <- readShapeSpatial("C:/rvis/japan_ver71/japan_ver71.shp",
                         proj4string=CRS("+init=epsg:4612"))

install.packages("rgdal")
library(rgdal)
EPSG <- make_EPSG()
EPSG[grep("JGD2000", EPSG$note), 1:2]
EPSG[262,]

plot(map,axes=T)
grid()

names(map)

map_fuk <- map[map$KEN=="福岡県",]
plot(map_fuk, axes=T)
install.packgaes("geosphere")
library(geosphere)
text(centroid(map_fuk), labels=map_fuk$SHIKUCHOSON, col="#666666",cex=0.6)
grid()

map_fuk <- map[floor(map$JCODE/1000)==40,]

install.packages("gpclib")
library(gpclib)

gpclibPermit()
map_union <- unionSpatialPolygons(map, map$KEN)
plot(map_union, axes=T)
grid()
map_fuk_union <- unionSpatialPolygons(map_fuk, map_fuk$KEN)
plot(map_fuk_union, axes=T)
grid()

river_fuk <- readShapeSpatial("C:/rvis/W05-07_40/W05-07_40_GB02.shp",
                                proj4string=CRS("+init=epsg:4612"))
plot(map_fuk_union, axes=T)
grid()
names(river_fuk)

plot(river_fuk[river_fuk$W05_003 %in% 1:3,], add=T, col="skyblue")




  • 最終更新:2013-05-23 18:45:47

このWIKIを編集するにはパスワード入力が必要です

認証パスワード