从坐标获取大致的大陆信息

Question

我有一些坐标，想知道它们属于哪个大陆。

coords <- data.frame(lat=29.51974, long=-94.88704)

然而，有时由于某些测量误差，坐标有点偏离。例如，上面给出的坐标在Galveston湾，因此并不总是被标记为“北美洲”。

一些类似的答案使用rworldmap，例如这里，但大陆边界非常精确，因此会显示缺失的大陆。其他版本（如使用google API），例如这里，对于给定的示例确实有效，但对于一些更偏离的示例则不太适用（例如几英里进入海洋）。

我找到了matlab的答案这里，看起来很完美。我考虑导出由matlab代码生成的shapefile，并将其导入R，但我没有拥有Matlab。

我的问题是：是否有一种方法在R中具有广泛的边界来标记坐标属于大陆，就像在这里（与matlab解决方案相同的链接）给出的那样。

- Apriconer

---

请查看 rgeos::gBuffer。也就是说，如果一个点不符合所有给定的大陆，请扩展大陆直到找到一个点适合其中的一个。 - r2evans

如果你想要快速解决这个问题，最简单的方法似乎是从你提供的 MATLAB 解决方案中抓取数据并编写自己的函数。 - s_baldur

Answer 1

你可以将matlab example转换为R列表，然后使用maps::map.where()函数（注意我手动关闭了多边形）

NAm=data.frame(lat=c(90,       90,  78.13,   57.5,  15,  15,  1.25,  1.25,  51,  60,    60, 90),
           lon=c(-168.75 ,-10 ,-10     ,-37.5 ,-30 ,-75 ,-82.5 ,-105  ,-180 ,-180 ,-168.75, -168.75))

NAm2 = data.frame(lat=c(51,    51,  60, 51),
          lon=c(166.6, 180, 180, 166.6))

SAm = data.frame(lat=c(1.25,   1.25,  15,  15, -60, -60, 1.25),
         lon=c(-105, -82.5,  -75, -30, -30, -105, -105))

europe=data.frame(lat=c(90,   90,  42.5, 42.5, 40.79, 41, 40.55, 40.40, 40.05, 39.17, 35.46, 
            33,   38,  35.42, 28.25, 15,  57.5,  78.13, 90),
          lon=c(-10, 77.5, 48.8, 30,   28.81, 29, 27.31, 26.75, 26.36, 25.19, 27.91,
                27.5, 10, -10,  -13,   -30, -37.5, -10, -10))

africa=data.frame(lat=c(15,  28.25 ,35.42 ,38 ,33   ,31.74 ,29.54 ,27.78 ,11.3 ,12.5 ,-60 ,-60, 15),
          lon=c(-30 ,-13   ,-10 ,10 ,27.5 ,34.58 ,34.92 ,34.46 ,44.3 ,52    ,75 ,-30, -30))

australia=data.frame(lat=c(-11.88, -10.27, -10 ,-30    ,-52.5 ,-31.88, -11.88),
            lon=c(110,      140  ,145 ,161.25 ,142.5  ,110, 110))

asia=data.frame(lat=c(90   ,42.5 ,42.5 ,40.79 ,41 ,40.55 ,40.4  ,40.05 ,39.17 ,35.46 ,33   ,
              31.74 ,29.54 ,27.78 ,11.3 ,12.5 ,-60 ,-60 ,-31.88 ,-11.88 ,-10.27 ,33.13 ,51    ,60  ,90, 90),
        lon=c(77.5 ,48.8 ,30   ,28.81 ,29 ,27.31 ,26.75 ,26.36 ,25.19 ,27.91 ,27.5 ,
              34.58 ,34.92 ,34.46 ,44.3 ,52   ,75  ,110  ,110   ,110    ,140    ,140   ,166.6 ,180 ,180, 77.5))

asia2=data.frame(lat=c(90    ,90      ,60      ,60, 90),
         lon=c(-180 ,-168.75 ,-168.75 ,-180, -180))

antarctica=data.frame(lat=c(-60, -60, -90, -90, -60),
              lon=c(-180, 180, 180, -180, -180))

continents=list(
        y=c(NAm$lat, NA, NAm2$lat, NA, SAm$lat, NA, europe$lat,NA,africa$lat,NA,
            australia$lat,NA,asia$lat,NA,asia2$lat,NA,antarctica$lat),
        x=c(NAm$lon, NA, NAm2$lon, NA, SAm$lon, NA,europe$lon,NA,africa$lon,NA,
            australia$lon,NA,asia$lon,NA,asia2$lon,NA,antarctica$lon),
        names=c("North America", "North America:2", "South America", "Europe",
            "Africa","Australia","Asia","Asia:2","Antarctica"))
class(continents) <- "map"

通过这个，你可以使用地图包：

maps::map.where(continents, x=-94.88704, y=29.51974)

- Alex Deckmyn