2つのポリゴンがあります。1つにはフィールド(X、Y、Z)が含まれ、もう1つには土壌タイプ(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は機能しません。あなたは土壌の種類について話しているのでそれが示唆されているけれども、あなたはあなたが地理座標を持っていることを指定しませんでした。
—
フーリエ