[問題] 使用lmap::gdist找兩組資料最近點的錯誤

看板R_Language作者giock18 (小武)時間8年前 (2017/10/26 22:20)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串1/2 (看更多)

[問題類型]: 程式諮詢(我想用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://www.ptt.cc/bbs/R_Language/M.1509027623.A.1BF.html