我没有使用你正在使用的方法,但我在 sp 中有一个在你的点和多边形上完美运行 .
我挑选了你的代码并将 lat 和 lon 作为向量,将点坐标作为值来满足函数要求 .
但是你可以很容易地创建一个数据框并明确地将列用作lat / lon值 .
以下是它的要点:
require(sp)
## Your polygon
lat
lon
## Your Point
lng=-71.05609
lt=48.40909
# sp function which tests for points in polygons
point.in.polygon(lt, lng, lat, lon, mode.checked=FALSE)
这是输出:
[1] 1
从文档中解释这个:
整数数组值是:
0点严格地在多边形外部
1点严格地是多边形的内部
2点位于多边形边缘的相对内部
3点是多边形的顶点
因为你的点是基于此的1,它应该完全在你的 Map 显示的多边形内!使用这些类型的数据获得良好输出的关键是以正确的格式提供变量 .
您可以轻松地将数据框 df 与 df$lat 和 df$lon 作为两个多边形变量以及测试框架 test 与 test$lat 和 test$lon 作为一系列点 . 您只需将等式中的每一个替换为等式:
point.in.polygon(df$lat, df$lon, test$lat, test$lon, mode.checked=FALSE)
并且它将返回0,1的2和3的向量