查找哪些地理编码点位于特定点的半径 25 英里的圆内

Question

我试图在某个地理编码点的 25 英里半径内绘制某些医疗设施。 ：

设施数据集如下所示：

Name        Lat          Long        Type        Color
A          42.09336    -76.79659      X          green
B          43.75840    -74.25250      X          green
C          43.16816    -77.60332      Y          blue

...

然而，设施列表跨越全国（美国），但我只想绘制圆圈内存在的设施。 缓冲圈的中心是坐标集(long =-73.857932, lat = 41.514096)和radius 25 miles 。

所以在我需要绘制的数据集中，我需要过滤设施列表、它们的纬度和经度、类型和颜色

我在这方面真的很新，而且期限很紧，所以如果有人能解释一下，那就太好了。

PS：我还想计算设施的类型（但我想这将是一个简单的 dplyr %>% n() 一旦创建过滤器，对吧？）

Answer 1

您可以使用该功能distHaversine从geosphere包（假设df是你的数据帧和43, -77是你的参考点的坐标）：

geosphere::distHaversine(c(43,  -77), df[, 2:3]) / 1609.334 <= 25
#[1]  TRUE FALSE FALSE

默认输出为米，因此除以 1609.334 将转换为英里。