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)
