作者BUQ (固执又暴燥)
看板R_Language
标题[问题] 反向地理编码(用座标查行政区)
时间Fri Apr 14 12:23:25 2017
[问题类型]:
程式谘询(我想用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://webptt.com/cn.aspx?n=bbs/R_Language/M.1492143808.A.CC7.html
1F：→ cywhale: 不知道有无较快方法? download政府开放图资取得shp file 04/14 14:11
3F：→ cywhale: 利用sp package point.in.polygon()即可取得所在polygon 04/14 14:13
4F：→ cywhale: 有了所在polygon index可以回头查你import shape file中 04/14 14:15
5F：→ cywhale: countycode, countyname, towncode, townname等you want 04/14 14:16
6F：→ BUQ: 解压缩有好多个档案哪一个？要怎麽汇进R？ 04/14 14:38
7F：→ BUQ: 谢谢C大..已解决...晚上贴解法 04/14 15:19
8F：推 Edster: 用 point in polygon 等应该已经是最快的了. 04/14 15:40
9F：→ BUQ: 我是用别的方法解的 但是很想要知道高手们说的办法.. 04/14 15:54
10F：→ BUQ: 可是我看不懂...没有没详细一点的讲解...<(_ _)>先谢谢了 04/14 15:55
11F：→ BUQ: 请问汇入是用哪个套件的哪个函数？ 04/14 22:53
12F：→ locka: 期待原po大大分享不用比对台湾的区界地图的解法～XD 04/14 23:17
13F：→ BUQ: 好喔.... 04/14 23:24
今早才在本板找到一个文章教我
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
14F：推 locka: 感谢分享，原来也可以这样！cool~ 04/14 23:50
15F：推 cywhale: 这解法第一次看到，感谢分享~ 04/15 10:02
16F：→ BUQ: 我也是运气好找到的 板上大部分文章基本上我是都看不懂太高端 04/15 10:35