2016-03-22 17 views
0

我想加入從gDistance輸出到用於生成距離矩陣的形狀文件之一的距離,我無法弄清楚。我以爲我知道了,但意識到通過加入行名,距離被加入到錯誤的記錄中。這是我迄今爲止所做的。需要將gDistance的輸出連接到用於測量R中的距離的shape文件中

Plants_dist <- gDistance(Sites.rg, Plants.rg, byid=T, hausdorff=F, densifyFrac = NULL) 
Plants_dist.df <- as.data.frame(Plants_dist, row.names = TRUE) 

我感興趣的只是記錄16093.44米內...

Plants_dist.df[apply(Plants_dist.df, MARGIN = 1, function(x) any(x < 16093.44)),] 

這是它打破了...當我試圖合併上row.names記錄創建植物的距離到站點的輸出表我得到上面定義的適當距離值,但它將距離合併到錯誤記錄。

Dist_Site_Plants_list <- merge([email protected], Sites_Plants_list, by.x="row.names", by.y="row.names", all.x=F, sort=F) 

我想這是因爲該行的名字沒有適當的排序,但我無法弄清楚怎麼回事,要做到這一點

+0

歡迎來到Stack Overflow!如果您花費一些時間來使問題具有可重現性,您會發現您會得到更好的答案。請按照[這裏的指導](http://stackoverflow.com/a/5963610/1036500)。謝謝! – Ben

+0

感謝您的反饋本。將來我會對那些試圖回答它的人提出更加友好的問題。 – Collin

回答

0

我正在這樣複雜得多,它需要的是... cbind Plants_dist(gDistance輸出)與[email protected](shapefile屬性表)。然後根據連接列中的值(距離)進行子集合。