[問題] 反向地理編碼(用座標查行政區)

[BUQ (固執又暴燥)]

[問題類型]:

程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)


[軟體熟悉度]:
請把以下不需要的部份刪除
一個月新手(沒寫過程式，R 是我的第一次)

[問題敘述]:

小弟有41*86=3526筆資料
每個資料都是 WGS84座標系統
要如何知道每個座標所在的行政區...例如：新北市 土城區

各方高手有什麼解法比較快的....
小弟現在是一個解法都沒有
或者   相關的套件   關鍵字  也行


座標格式如下  一行經度  一行緯度
[1] 121.998
[1] 21.52036



[程式範例]:



[環境敘述]:

請提供 sessionInfo() 的輸出結果，

R version 3.3.3 (2017-03-06)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)



[關鍵字]: 反向地理編碼


--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.135.58.147
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1492143808.A.CC7.html
→ cywhale: 不知道有無較快方法? download政府開放圖資取得shp file
→ cywhale: http://data.gov.tw/node/7441 then importShapefile
→ cywhale: 利用sp package point.in.polygon()即可取得所在polygon
→ cywhale: 有了所在polygon index可以回頭查你import shape file中
→ cywhale: countycode, countyname, towncode, townname等you want
→ BUQ: 解壓縮有好多個檔案哪一個？要怎麼匯進R？
→ BUQ: 謝謝C大..已解決...晚上貼解法
推 Edster: 用 point in polygon 等應該已經是最快的了.
→ BUQ: 我是用別的方法解的 但是很想要知道高手們說的辦法..
→ BUQ: 可是我看不懂...沒有沒詳細一點的講解...<(_ _)>先謝謝了
→ BUQ: 請問匯入是用哪個套件的哪個函數？
→ locka: 期待原po大大分享不用比對台灣的區界地圖的解法～XD
→ BUQ: 好喔....
今早才在本板找到一個文章教我
library(sos)#超強找套件 只要下關鍵字就好
findFn("Reverse geocoding")
結果 跑出一個網址
http://localhost:24955/session/file102864b750ad.html

我就選了 第一個來試用  (分數很高)
library(googleway)
API="這請打上自己向google map API key"    #申請是免費的 但有訪問上限
kk=google_reverse_geocode(location =c(23.48742,120.00024),language="zh-TW",key=API)
kk$results$address_components[[1]]$long_name

[1] "海埔北路" "鰲鼓村"   "東石鄉"   "嘉義縣"   "台灣"     "614"

最後補上API申請網址

https://developers.google.com/maps/documentation/javascript/get-api-key?hl=zh-tw

※ 編輯: BUQ (1.162.241.154), 04/14/2017 23:45:49
※ 編輯: BUQ (1.162.241.154), 04/14/2017 23:46:58
推 locka: 感謝分享，原來也可以這樣！cool~
推 cywhale: 這解法第一次看到，感謝分享~
→ BUQ: 我也是運氣好找到的 板上大部分文章基本上我是都看不懂太高端