PTT數位生活區 / R_Language

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

看板R_Language作者 (固執又暴燥)時間8年前 (2017/04/14 12:23), 8年前編輯推噓3(3推 0噓 13→)
留言16則, 4人參與, 最新討論串1/2 (看更多)
[問題類型]: 程式諮詢(我想用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
→
04/14 14:11, , 1F
不知道有無較快方法? download政府開放圖資取得shp file
04/14 14:11, 1F
→
04/14 14:12, , 2F
http://data.gov.tw/node/7441 then importShapefile
04/14 14:12, 2F
→
04/14 14:13, , 3F
利用sp package point.in.polygon()即可取得所在polygon
04/14 14:13, 3F
→
04/14 14:15, , 4F
有了所在polygon index可以回頭查你import shape file中
04/14 14:15, 4F
→
04/14 14:16, , 5F
countycode, countyname, towncode, townname等you want
04/14 14:16, 5F
→
04/14 14:38, , 6F
解壓縮有好多個檔案哪一個？要怎麼匯進R？
04/14 14:38, 6F
→
04/14 15:19, , 7F
謝謝C大..已解決...晚上貼解法
04/14 15:19, 7F
推
04/14 15:40, , 8F
用 point in polygon 等應該已經是最快的了.
04/14 15:40, 8F
→
04/14 15:54, , 9F
我是用別的方法解的 但是很想要知道高手們說的辦法..
04/14 15:54, 9F
→
04/14 15:55, , 10F
可是我看不懂...沒有沒詳細一點的講解...<(_ _)>先謝謝了
04/14 15:55, 10F
→
04/14 22:53, , 11F
請問匯入是用哪個套件的哪個函數？
04/14 22:53, 11F
→
04/14 23:17, , 12F
期待原po大大分享不用比對台灣的區界地圖的解法～XD
04/14 23:17, 12F
→
04/14 23:24, , 13F
好喔....
04/14 23:24, 13F
今早才在本板找到一個文章教我 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
推
04/14 23:50, , 14F
感謝分享，原來也可以這樣！cool~
04/14 23:50, 14F
推
04/15 10:02, , 15F
這解法第一次看到，感謝分享~
04/15 10:02, 15F
→
04/15 10:35, , 16F
我也是運氣好找到的 板上大部分文章基本上我是都看不懂太高端
04/15 10:35, 16F
‣ 更多 BUQ 的文章
文章代碼(AID): #1Oy4x0p7 (R_Language)
R_Language 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 BUQ 的文章
文章代碼(AID): #1Oy4x0p7 (R_Language)