R中的接近分析

时间:2017-06-20 14:56:36

标签: r gps gis

我在R中有一组GPS坐标,我想将其视为对应于感兴趣的患者的另一组GPS坐标的“曝光”。如何在R中进行接近度分析,将患者分为两组:在暴露坐标的x米范围内和远处的患者?

1 个答案:

答案 0 :(得分:0)

我认为这样的事情应该有效:

library(sp)
exp <- data.frame(lat= 40.741895,long = -73.989308)
patients <- data.frame(lat = rnorm(10,exp$lat,0.1),long = rnorm(10,exp$long,0.1))
coordinates(patients) <- ~ long + lat
coordinates(exp) <- ~ long + lat
d <- spDistsN1(coordinates(patients),coordinates(exp),longlat = TRUE)

强烈启发:http://handlebarsjs.com/block_helpers.html