मैं ऑस्ट्रेलिया में दूरदराज के क्षेत्रों दिखा कोई शेपफ़ाइल, ऑस्ट्रेलियाई सांख्यिकी ब्यूरो से प्राप्त राशि शेपफ़ाइल बहुभुज रंग:आर - साजिश के साथ फंस() - एक स्लॉट मूल्य पर आधारित
http://www.abs.gov.au/AUSSTATS/[email protected]/DetailsPage/1270.0.55.005July%202011?OpenDocument
उसी यूआरएल पर एक पीडीएफ है "एएसजीएस रिमोटनेस स्ट्रक्चर एडिशन 2011 पीडीएफ मैप्स" - मैं इस पीडीएफ दस्तावेज़ से पहले मानचित्र को पुन: पेश करने की कोशिश कर रहा हूं।
मैं शेपफ़ाइल में पढ़ सकते हैं और data
स्लॉट के लिए रंग जानकारी जोड़ लिया है:
ra <- readShapeSpatial("RA_2011_AUST", delete_null_obj = TRUE)
[email protected]$COLOUR <- "#FFFFFF"
[email protected]$COLOUR[(as.numeric(as.character([email protected]$RA_CODE11)) %% 10) == 0] <- "#006837"
[email protected]$COLOUR[(as.numeric(as.character([email protected]$RA_CODE11)) %% 10) == 1] <- "#31A354"
[email protected]$COLOUR[(as.numeric(as.character([email protected]$RA_CODE11)) %% 10) == 2] <- "#78C679"
[email protected]$COLOUR[(as.numeric(as.character([email protected]$RA_CODE11)) %% 10) == 3] <- "#C2E699"
[email protected]$COLOUR[(as.numeric(as.character([email protected]$RA_CODE11)) %% 10) == 4] <- "#FFFFCC"
केवल एक चीज छोड़ दिया मुझे करने के लिए साजिश नक्शा है! यह वह जगह है जहाँ मैं अटक जाते हैं ...
[email protected]
35 बहुभुज, जिनमें से प्रत्येक एक स्लॉट ID
जो डेटा फ्रेम [email protected]
लिए एक सूचकांक है की एक सूची है। तो मुझे [email protected]$COLOUR[ID]
में रंग खोजने के लिए plot()
बताएं। खैर, काफी नहीं। 35 बहुभुजों में से प्रत्येक (वर्ग "बहुभुज") में बहुभुज की अपनी सूची है (वर्ग "बहुभुज"); कुल मिलाकर 6902 बहुभुज हैं !!!
plot()
की मेरी समझ यह है कि मुझे इसे उसी क्रम में रंगों का एक वेक्टर पास करना होगा क्योंकि बहुभुज प्लॉट किया जाएगा। इसलिए मेरा मानना है कि मुझे संबंधित पॉलीगॉन के लिए रंग मान रखने वाले प्रत्येक तत्व के साथ लंबाई 6902 का वेक्टर बनाना होगा। मैं अब तक कैसे कर रहा हूँ?
बहुभुज क्रम में प्लॉट किए जाने पर यह आसान होगा, लेकिन वे नहीं हैं। 35 बहुभुजों में से प्रत्येक में स्लॉट plotOrder
है जो एक पूर्णांक वेक्टर है, इसलिए रंग वेक्टर, संभवतः, इन वैक्टरों में से प्रत्येक के मानों द्वारा आदेश दिया जाना चाहिए।
इस बिंदु पर यह सब कुछ जटिल लगता है। क्या मैं यहां पूरी तरह से ट्रैक बंद कर रहा हूँ?
आपकी सलाह के लिए धन्यवाद!
आपको rgdal पैकेज से readOGR का उपयोग करके आकारफाइल पढ़ना चाहिए, और यदि आकारफाइल के साथ .prj फ़ाइल है तो प्रक्षेपण डेटा में लोड होगा। – Spacedman