使用gContains测试lat点是否位于地图中

Question

我在dismo包中有一个从gmap函数得到的区域的地图。然后我试着测试某些纬度和经度点是否在这张地图上。

以下是我的可复制的例子：

library(dismo)
library(sp)
library(rgdal)
library(XML)
library(rgeos)


spatP <- SpatialPoints(gmap("Palo Alto, CA", lonlat = TRUE), 
    proj4string = CRS("+proj=longlat"))

point <- data.frame( x =  -122.129610, y = 37.399761)
pointSp <- SpatialPoints(point, proj4string = CRS("+proj=longlat"))

gContains(spatP, pointSp, byid = FALSE)

这一点就在地图上，但我显然遗漏了一些东西。

Answer 1

想出了解决办法。如果我将SpatialPoints数据帧spatP转换为凸多边形，它就能工作。

polyHull <- gConvexHull(spatP)
gContains(polyHull, pointSp)
[1] TRUE