Tôi có hai đa giác. Một cái chứa các trường (X, Y, Z) và cái kia chứa các loại đất (A, B, C, D). Tôi muốn biết khu vực nào của mọi lĩnh vực đều chứa loại đất nào. Tôi đã thử như sau:
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
và đạt được kết quả tốt với nó cho tôi biết lĩnh vực nào chứa loại đất nào. Tuy nhiên, làm thế nào để tôi có được khu vực thay thế?
1
Lưu ý, st_intersection sẽ không hoạt động nếu điểm của bạn là vĩ độ và kinh độ. Bạn đã không xác định rằng bạn có tọa độ địa lý, mặc dù điều đó được gợi ý vì bạn đang nói về các loại đất.
—
Fourier