サンプルコード
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")