R : 하나의 SpatialLines (DataFrame) 내에서 모든 선의 교차점에서 모든 xy 좌표 가져 오기

Question

개체 또는 SpatialLinesDataFrame 내의 모든 교차점의 xy 좌표를 얻는 방법을 찾고 있습니다. 나는 rgeos의 함수 gIntersect을 찾았지만 두 개의 데이터 세트 사이의 교차점 만 보았습니다. 50 만 개가 넘는 라인의 데이터 세트로 작업하고 있기 때문에 모든 라인을 별도의 파일로 만들고 어떤 라인이 다른 라인과 교차하는지 확인하는 데 너무 많은 시간이 걸립니다. ArcMap에는 몇 초 만에 할 수있는 Intersect 함수가 있으며 R에 그러한 함수가 있는지 궁금합니다. 고마워요! 당신이 spatstat에서 psp 객체로 개체를 변환하는 경우

Answer 1

당신은 spatstat 기능 selfcrossing.psp을 사용할 수 있습니다. 그러나 교차 수가 잠재적으로 엄청날 수 있기 때문에 50 만 라인에 어떻게 대처할 지 확신하지 못합니다. 아래 코드는 무작위 세그먼트 패턴을 생성하고 자체 교차점을 찾습니다.

이 코드는 많은 메모리를 차지하고 R을 죽일 수 있으므로주의하십시오. 따라서 50 만 줄을 처리하기 전에 점진적으로 증가하는 예제를 사용해보십시오. 아래의 코드는 내 5 년 된 랩톱에서 꽤 많은 메모리를 사용하고 5 초가 걸렸습니다.

set.seed(42) 
N <- 1e4 
x <- psp(runif(N), runif(N), runif(N), runif(N), owin(), check=FALSE) 
y <- selfcrossing.psp(x)