मेरे पास दो बहुभुज हैं। एक में फ़ील्ड (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 काम नहीं करेगा। आपने निर्दिष्ट नहीं किया कि आपके पास भौगोलिक निर्देशांक थे, हालांकि यह संकेत दिया गया है क्योंकि आप मिट्टी के प्रकारों के बारे में बात कर रहे हैं।
—
फूरियर