[问题] 反向地理编码(用座标查行政区)

楼主: BUQ (固执又暴燥)   2017-04-14 12:23:25
[问题类型]:
程式咨询(我想用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)
[关键字]: 反向地理编码
作者: cywhale (cywhale)   2017-04-14 14:11:00
不知道有无较快方法? download政府开放图资取得shp filehttp://data.gov.tw/node/7441 then importShapefile利用sp package point.in.polygon()即可取得所在polygon有了所在polygon index可以回头查你import shape file中countycode, countyname, towncode, townname等you want
楼主: BUQ (固执又暴燥)   2017-04-14 14:38:00
解压缩有好多个档案哪一个?要怎么汇进R?谢谢C大..已解决...晚上贴解法
作者: Edster (Edster)   2017-04-14 15:40:00
用 point in polygon 等应该已经是最快的了.
楼主: BUQ (固执又暴燥)   2017-04-14 15:54:00
我是用别的方法解的 但是很想要知道高手们说的办法..可是我看不懂...没有没详细一点的讲解...<(_ _)>先谢谢了请问汇入是用哪个套件的哪个函数?
作者: locka (locka)   2017-04-14 23:17:00
期待原po大大分享不用比对台湾的区界地图的解法~XD
楼主: BUQ (固执又暴燥)   2017-04-14 23:24:00
好喔....
作者: locka (locka)   2017-04-14 23:50:00
感谢分享,原来也可以这样!cool~
作者: cywhale (cywhale)   2017-04-15 10:02:00
这解法第一次看到,感谢分享~
楼主: BUQ (固执又暴燥)   2017-04-15 10:35:00
我也是运气好找到的 板上大部分文章基本上我是都看不懂太高端

Links booklink

Contact Us: admin [ a t ] ucptt.com