求助距离问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有以下两组数据a1和a2，其两列数据都分别表示经度和纬度，这里只列举一小部分
a1
x1      y1      z1
43.46    65.18    2.29
43.54    65.51    1.45
43.55    66.05    1.52
44.01    65.56    2.24
44.04    65.55    1.7
44.08    65.56    1.54
44.13    77.08    2.03
44.14    77.15    2.36
44.15    79.55    1.54
a2
x2      y2      z2
43.46    83.42    2.04
43.46    66.18    1.68
43.47    79.79    1.68
43.48    81.6    2.02
43.48    77.26    2.02
43.49    73.32    2.2
43.5    66.89    1.66
43.5    62.95    1.66
43.5    60.43    1.74
43.51    72.61    1.74
43.52    84.55    2.61
43.52    67.61    2.77
43.53    71.89    2.21
43.54    76.55    2.25
43.54    68.32    1.94
43.55    79.08    2.15
43.55    71.18    2.15
43.55    69.04    2.07
43.56    70.46    2.65
在a2中找出与a1中每个（x1，y1）距离最短的点（x2，y2），把匹配出来的每一行（x2，y2，z2）添加到a1相应的行后边形成（x1，y1，z1，x2，y2，z2）

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：小部分 距离

没看懂，你打算要做什么呢

给你一个小函数，big radios dintance距离的

1. getDistance <- function(lat1,lon1,lat2,lon2){
   
2.     R <- 6378.137
   
3.     #lat weidu lon jingdu
   
4.     rad <- function(d){d*pi/180}
   
5.     d <- acos(sin(rad(lat1))*sin(rad(lat2)) + cos(rad(lat1))*cos(rad(lat2)) * cos(rad(lon2)-rad(lon1))) * R
   
6.     return(d)
   
7. }

复制代码

就是输入两个点之间的纬度经度，就会给出两点之间的最短距离。计量单位为公里。