2015年9月19日 星期六

用R繪製shapefile (SHP) 資料練習

library(ggmap)
library(rgdal)
library(ggplot2)

mydb <- dbConnect(MySQL(), username="root", host="localhost", dbname="CSD", password="p9996t90k", client.flag = CLIENT_MULTI_STATEMENTS)
dbSendQuery(mydb,"SET NAMES utf8")
dbClearResult(dbListResults(mydb)[[1]])

tpn = readOGR("Taipei","G97_63000_U0200_2015") %>% spTransform(CRS("+proj=longlat +datum=WGS84"))

L302 <- subset(cust_geo_data, cust_geo_data$cod_loc == "L302")
include_town <- unique(L302$village)
area_town <- c("北投區", "士林區", "內湖區", "中山區")
L302.mp <- tpn[(tpn$TOWN %in% area_town), ]
town.f = L302.mp %>% fortify(region = 'VILLAGE')

NHNeighbourhoods  = merge(town.f, L302.mp@data, by.x = 'id', by.y = 'VILLAGE')

nhpostcodes <- unique(NHNeighbourhoods$id)
nhvalues = data.frame(id = c(nhpostcodes),
                      value = c(runif(171 ,5.0, 25.0)))

nh = merge(NHNeighbourhoods, nhvalues, by.x='id')
tpMap = map = get_map(location = "北投區", zoom = 12)
ggmap(tpMap) +
  geom_polygon(aes(fill = value, x = long, y = lat, group = group),
               data = nh,
               alpha = 0.8,
               color = "lightyellow",

               size = 0.2)