gadm_taiwan_map.R

ben — Jan 12, 2014, 5:31 AM

library(sp)
library(RColorBrewer)

# get spatial data for Taiwan on county level
con <- url("http://www.gadm.org/data/rda/TWN_adm2.RData")
print(load(con))

[1] "gadm"

close(con)

# plot Taiwan with random colors, length of the vector taiwan is 21
taiwan <- c("高雄市","台北市","彰化縣","嘉義縣","新竹縣","花蓮縣","宜蘭縣","高雄縣","基隆市","苗栗縣","南投縣","澎湖縣","屏東縣","台中縣","台中市","台南縣","台南市","新北市","台東縣","桃園縣","雲林縣")
gadm$NAME_2 <- as.factor(taiwan)
col = rainbow(length(levels(gadm$NAME_2)))
spplot(gadm, "NAME_2", col.regions=col, main="Taiwan Regions", colorkey = FALSE, lwd=.4, col=grey(.9))