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