我有两个多边形。一个包含字段(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将不起作用。尽管有人在谈论土壤类型,但暗示您没有指定地理坐标。
—
傅里叶