У меня есть два полигона. Одно содержит поля (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 не будет работать, если ваши точки широты и долготы. Вы не указали, что у вас есть географические координаты, хотя на это намекают, поскольку вы говорите о типах почв.
—
Фурье