두 개의 다각형이 있습니다. 하나는 들판 (X, Y, Z)을 포함하고 다른 하나는 토양 유형 (A, B, C, D)을 포함합니다. 모든 분야에서 어떤 종류의 토양이 포함되어 있는지 알고 싶습니다. 나는 다음을 시도했다.
library(rgdal)
library(rgeos)
Field<-readOGR("./","Field")
Soil<-readOGR("./","Soil")
Results<-gIntersects(Soil,Field,byid=TRUE)
rownames(Results)<-Field@data$FieldName
colnames(Results)<-Soil@data$SoilType
> Results
A B C D
Z TRUE FALSE FALSE FALSE
Y FALSE TRUE TRUE FALSE
X TRUE TRUE TRUE TRUE
어떤 필드에 어떤 토양 유형이 들어 있는지 알려주는 좋은 결과를 얻었습니다. 그러나 대신 영역을 어떻게 얻습니까?
1
참고로, 포인트가 위도 및 경도이면 st_intersection이 작동하지 않습니다. 토양 유형에 대해 이야기하고 있기 때문에 힌트를 받았지만 지리적 좌표가 있음을 지정하지 않았습니다.
—
푸리에