刚好有做过同样的事情 可以使用sp package 的 SpatialPointsDataFrame()跟over() (诚如cywhale大大提到的，先汇入台湾行政区shape file档案) df为所求经纬度的data frame twn_shp为汇入的shape file (class为SpatialPolygonsDataFrame) spdf <- SpatialPointsDataFrame(coords=(df$lng,df$lat), data=df, proj4string = CRS("+proj=longlat +ellps=GRS80 +no_defs")) result <- over(spdf, twn_shp) over会回传一个新的SpatialPointsDataFrame 其中 result$Substitute为村里 result$T_Name为行政区 result$C_Name为都市... 之前研究这个的时候才发现sp的over很强大 (help over会看到一整份pdf档...有很多有趣的用法) 这应该是更快更直接的方法，分享给大家~ ps. 是说有人知道可以帮忙解释proj4string那段吗 哈哈 那里是我上网看别人的范例这样下参数， 应该跟指定座标系统有关 不过不清楚其含义 ※ 引述《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), 来自: 61.230.94.38 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1492160213.A.C7C.html ※ 编辑: locka (61.230.94.38), 04/14/2017 17:01:53 rgdal::readOGR ※ 编辑: locka (61.230.94.38), 04/14/2017 23:18:52