作者giock18 (小武)
看板R_Language
标题[问题] 使用lmap::gdist找两组资料最近点的错误
时间Thu Oct 26 22:20:20 2017
[问题类型]:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
入门(写过其他程式，只是对语法不熟悉)
[问题叙述]:
我有两组经纬度座标，每组资料长得像这样
data$Lon data$Lat
100.7200 13.61500
100.6683 13.52000
100.5717 13.54500
100.5750 13.69333
100.5783 13.82333
100.3983 13.90833
另一组称为temp
temp$Lon & temp$Lat
其中一组是参考地标的经纬度，另外一组是实际的轨迹
我使用gdist来算出点和点之间的距离，并找最小值，
也就是最靠近这个地标的点。
因此我的程式码这样写
for (i in nrow(data)){
closest<-which.min(gdist(data$Lon[i],
OFPdata$Lat[i],temp$LONGITUDE,temp$LATITUDE,units = "nm",verbose = FALSE) )}
用回圈的方式去找每个地标对应到的最靠近轨迹
which.min虽然有用，但是应该是因为原本lmap::gdist的设计是输入单值，
但我却给他一栏资料，所以他会跳出警告如下，虽然程式还是可以work，
但想知道有没有更正确的写法？
Warning messages:
1: In while (abs(lamda - lamda.old) > 1e-11) { :
the condition has length > 1 and only the first element will be used
2: In while (abs(lamda - lamda.old) > 1e-11) { :
the condition has length > 1 and only the first element will be used
另外，为了提升程式的效率，想要将轨迹与地标的计算迭代缩小成
经度或纬度差5度以内的才计算，不知道有甚麽方法可以达成呢？
我有在stackoverflow发问，原文网址如下:
https://stackoverflow.com/questions/46843802/find-the-nearest-point-using-lmapgdist-function
先谢谢了
[环境叙述]:
R version 3.4.1 (2017-06-30)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)
[关键字]:
lmap,gdsit,nearest point,closer point,经纬度距离
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 118.165.215.45
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1509027623.A.1BF.html