2012-06-17 11 views
5

मुझे लगता है कि मुझे अपने भयानक शब्द पिछले प्रश्न (इसे हटाया गया) को फिर से लिखना होगा। यहाँ एक और कोशिश है। मैं आसन्न क्षेत्रों में शामिल होना चाहता हूं, इस तरह कि उनकी आम सीमा गायब हो जाती है और केवल उनकी बाहरी रेखा देखी जा सकती है।आर (समेकित स्थानिक डेटा) में पड़ोसी क्षेत्रों को मर्ज करें?

यहाँ एक प्रतिलिपि प्रस्तुत करने योग्य उदाहरण है:

require(shapefiles) 
require(sp) 

xx <- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1], 
        IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) 

# show all the subregions 
plot(xx) 

enter image description here

अब केवल उन्हीं क्षेत्रों क्षेत्रों 3 और 5

plot(xx[c(3,5),]) 

कैसे मैं सिर्फ इन क्षेत्रों समेकित कर सकते हैं पर विचार करें। व्यावहारिक रूप से मैं जो करना चाहता हूं वह पूरे महाद्वीप का नक्शा रखने जैसा है जो सभी देशों को दिखाता है और नक्शा तैयार करता है जो उत्तरी अमेरिका और दक्षिण अमेरिका को दिखाता है।

मेरे लिए यह एक बहुत ही आम काम जैसा दिखता है लेकिन मुझे अब तक ऐसा करने के लिए सही कार्य नहीं मिल रहा है। क्या मुझे बस एक फंक्शन याद आती है या क्या मैं इसे मैन्युअल रूप से कर सकता हूं?

enter image description here

+2

'Maptools' पैकेज से 'UnionSpatialPolygons()' का उपयोग करें: http://gis.stackexchange.com/questions/21360/merging-2-polygon-adjacent-to-each-other-using-r –

+0

मुझे पता है कि आर में आपका प्रश्न कहता है, लेकिन स्थानिक मैनिपुलेशन लाइब्रेरी उत्कृष्टता पोस्टग्रेस/पोस्टगिस है। आप पोस्टग्रेज़ में आर संग्रहीत प्रक्रियाओं को भी लिख सकते हैं, http://www.bostongis.com/PrinterFriendly.aspx?content_name=postgresql_plr_tut01 देखें, जिसका अर्थ है कि आपको स्थानिक इंडेक्सिंग और पोस्टगिस से कार्य और आर के आंकड़े मिलते हैं। आयात कार्य को पोस्ट करने के लिए एक shp साथ ही पोस्टग्रेस के साथ जहाजों। –

उत्तर

9

rgeos पैकेज Spatial* डेटा, कि इस मामले में इस्तेमाल किया जा सकता से निपटने के लिए उत्कृष्ट उपकरणों का एक नंबर प्रदान करता है।

उदाहरण के लिए:

library(rgeos) 
regionOfInterest <- gUnion(xx[3,], xx[5,]) 

यह भी एक ही परिणाम है, और एकाधिक बहुभुज के लिए और अधिक उपयोगी हो सकता है:

regionOfInterest <- gUnionCascaded(xx[c(3,5), ]) 

plot(regionOfInterest) से परिणाम:

enter image description here

संबंधित मुद्दे