使用经纬度计算距离时出错

Question

我有一个包含 4 列的数据集：

start_lat 41.90687、41.94367、41.93259、41.89076

start_lng -87.62622, -87.64895, -87.63643, -87.63170

end_lat 41.90672、41.98404、41.93650、41.91831

end_lng -87.63483, -87.66027, -87.64754

我想添加名为 distance 的第 5 列，我使用以下公式，但出现以下错误：

trips_202007 <- trips_202007 %>%
  rowwise()
  mutate(distance = distm(c(trips_202007["start_lng"], trips_202007["start_lat"]), c(trips_202007["end_lng"], trips_202007["end_lat"]), fun=distHaversine))

.pointsToMatrix(x) 中的错误：“list”对象不能被强制输入“double”

有没有办法避免错误？

非常感谢您的建议。

这是我在 dput(trips_202007, file="trips_202007.txt"): https://drive.google.com/file/d/1lnXeNqIRCDad0WotgoUQhrfPg0AXApLr/view?usp=sharing之后得到的文件

Answer 1

您遇到以下错误：

Error in .pointsToMatrix(x) : 'list' object cannot be coerced to type 'double'

这是因为您的函数正在尝试计算矩阵。您问 R 是否会将结果放入一个单元格中，而 R 不希望这样。对于这种情况，我发现以下功能正在工作：

    dt.haversine <- function(lat_from, lon_from, lat_to, lon_to, r = 6378137){
          radians <- pi/180
          lat_to <- lat_to * radians
          lat_from <- lat_from * radians
          lon_to <- lon_to * radians
          lon_from <- lon_from * radians
          dLat <- (lat_to - lat_from)
          dLon <- (lon_to - lon_from)
          a <- (sin(dLat/2)^2) + (cos(lat_from) * cos(lat_to)) * (sin(dLon/2)^2)
          return(2 * atan2(sqrt(a), sqrt(1 - a)) * r)
}

先看那个函数我找到了here，然后用这段代码重新开始我们的计算：

trips_202007 <- trips_202007 %>%
          mutate(distance = dt.haversine(start_lat, start_lng, end_lat, end_lng))
trips_202007