Saya punya dua poligon. Satu berisi bidang (X, Y, Z) dan yang lainnya berisi jenis tanah (A, B, C, D). Saya ingin tahu area mana dari setiap bidang yang berisi jenis tanah apa. Saya mencoba yang berikut ini:
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
dan mencapai hasil yang baik dengan itu memberi tahu saya bidang mana yang berisi jenis tanah mana. Namun, bagaimana cara mendapatkan area itu?
1
Sebagai catatan, st_intersection tidak akan berfungsi jika poin Anda adalah garis lintang dan bujur. Anda tidak menentukan bahwa Anda memiliki koordinat geografis, meskipun itu mengisyaratkan karena Anda berbicara tentang jenis tanah.
—
Fourier