[問題類型]:
程式諮詢(我想用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,經緯度距離