# 第一個問題(分行政區畫圖)
data <- read.csv ( "ntp.gis.csv" , fileEncoding = "UTF-8" , header = TRUE)
map <- get_map( location = "Taipei" , zoom = 11 , language = "zh-TW" , color = "bw" )
ggmap( map) +
geom_point( aes( x = zip_x, y = zip_y) , color= "red" , data = data) +
facet_wrap( ~ district)
# 第 2 & 3 個問題(自動取最適合的zoom-in尺度, 經緯度精確度問題)
# 板橋
banqiao <- subset( data, data$district== '板橋區' )
map <- get_map( location = c( lon = 121.453220 , lat = 25.006734 ) , zoom = 14 , language = "zh-TW" , color = "bw" )
ggmap( map) +
#geom_point(aes(x = zip_x, y = zip_y), color = "red", data = banqiao) +
geom_jitter( aes( x = zip_x, y = zip_y) , col= "red" , position = position_jitter( w= 0.0005 , h= 0.0005 ) , data= banqiao)
IyDnrKzkuIDlgIvllY/poYwo5YiG6KGM5pS/5Y2A55Wr5ZyWKQpkYXRhIDwtIHJlYWQuY3N2KCJudHAuZ2lzLmNzdiIsZmlsZUVuY29kaW5nID0gIlVURi04IiwgaGVhZGVyID0gVFJVRSkKCm1hcCA8LSBnZXRfbWFwKGxvY2F0aW9uID0gIlRhaXBlaSIsIHpvb20gPSAxMSwgbGFuZ3VhZ2UgPSAiemgtVFciLCBjb2xvciA9ICJidyIpCmdnbWFwKG1hcCkgKwpnZW9tX3BvaW50KGFlcyh4ID0gemlwX3gsIHkgPSB6aXBfeSksIGNvbG9yPSJyZWQiLCBkYXRhID0gZGF0YSkgKwpmYWNldF93cmFwKH4gZGlzdHJpY3QpCgojIOesrCAyICYgMyDlgIvllY/poYwo6Ieq5YuV5Y+W5pyA6YGp5ZCI55qEem9vbS1pbuWwuuW6piwg57aT57ev5bqm57K+56K65bqm5ZWP6aGMKQojIOadv+apiwpiYW5xaWFvIDwtIHN1YnNldChkYXRhLGRhdGEkZGlzdHJpY3Q9PSfmnb/mqYvljYAnKQptYXAgPC0gZ2V0X21hcChsb2NhdGlvbiA9IGMobG9uID0gMTIxLjQ1MzIyMCwgbGF0ID0gMjUuMDA2NzM0KSwgem9vbSA9IDE0LCBsYW5ndWFnZSA9ICJ6aC1UVyIsIGNvbG9yID0gImJ3IikKZ2dtYXAobWFwKSArCiNnZW9tX3BvaW50KGFlcyh4ID0gemlwX3gsIHkgPSB6aXBfeSksIGNvbG9yID0gInJlZCIsIGRhdGEgPSBiYW5xaWFvKSArCmdlb21faml0dGVyKGFlcyh4ID0gemlwX3gsIHkgPSB6aXBfeSksY29sPSJyZWQiLCBwb3NpdGlvbiA9IHBvc2l0aW9uX2ppdHRlcih3PTAuMDAwNSxoPTAuMDAwNSksZGF0YT1iYW5xaWFvKQo=