각 점의 지정된 거리 내에있는 점 찾기

Question

UTM (Universal Transverse Mercator) 좌표로 장소 목록과 지리적 위치가있는 데이터 프레임이 있습니다.

Place X_UTM Y_UTM 
    1 574262.0 6140492 
    2 571251.2 6141669 
    3 570841.9 6142535 
    4 570233.8 6141213 
    5 578269.2 6140304 
    6 575067.1 6137444 

내가 확인하고 싶은 다른 장소의 주어진 유클리드 거리 내에있는 각 장소합니다 (dataframe의 각 행)에 대한 :

는

그것은 다음과 같이 보입니다. 이 경우, 나는 어느 장소가 1 킬로미터보다 더 가까운지를 찾고 싶다. 나는 이런 식으로 뭔가를 시도했습니다

: 1,000m보다 더 가까운 지점을 반환하는 함수이어야한다

foo <- function(x, y) dist(c(x, y), method = "euclidian") < 1000 

. 그런 다음 :

x <- lapply(df(,c(i, x, y)), FUN = foo) 

i이 "Place"입니다

, x는 "X_UTM"하고 y는 "Y_UTM"입니다. 이것은 전혀 작동하지 않습니다.

내가 (위의 숫자에서 얻을 수없는) 같은 것을 보일 것입니다 후 오전 출력 :

# Place Closest 
#  1 2, 5 
#  2  1 
#  3  NA 
#  4  5 
#  5 1, 4 
#  6  NA 

Answer 1

sp::spDists을 사용하여 거리 매트릭스를 반환하고 조건에 맞는 각 열/행의 요소를 찾을 수도 있습니다. 예를 들어

: 내가 잘못 될 경우

d <- read.table(text='Place X_UTM  Y_UTM 
1  574261.98 6140492.13 
2  571251.23 6141669.26 
3  570841.92 6142534.86 
4  570233.75 6141212.5 
5  578269.25 6140303.78 
6  575067.07 6137444.36', header=TRUE) 

library(sp) 
i <- apply(spDists(as.matrix(d[, c('X_UTM', 'Y_UTM')])), 2, 
      function(x) paste(which(x < 1000 & x != 0), collapse=', ')) 

data.frame(Place=d$Place, Closest=i) 

## Place Closest 
## 1  1   
## 2  2  3 
## 3  3  2 
## 4  4   
## 5  5   
## 6  6   

Answer 2

난 당신이 십자가에 필요할 것 생각은 장소 좌표 가입 할 수 있습니다. 그 이유는 모든 쌍의 장소가 가장 가까운 이웃이 될 수 있고 특정 쌍을 배제 할 수있는 선험적 정보가 없다고 가정 할 때 모든 쌍을 검사해야합니다.

십자가가 데이터 프레임 df의 조인 얻는 한 가지 방법은 merge에 매개 변수로 by = NULL 설정 자체와 병합하는 것입니다

: 이제 당신이해야 할 모든이의 최소 거리를 찾을 수있다

df.cross <- merge(x = df, y = df, by = NULL) 
df.cross$distance <- apply(df.cross[, c('X_UTM.x', 'X_UTM.y', 'Y_UTM.x', 'Y_UTM.y')], 
     1, 
     function(x) dist(x[1], x[2], x[3], x[4])) 

을 각 쌍의 장소.