मैं आर पैकेज के लिए नया हूँ और आरस्टर का उपयोग कर रहा हूँ। मुझे मौजूदा रैस्टर फ़ाइल से पॉलीगॉन निकालने में समस्या है। अगर मैं उपयोग करता हूं
extract(raster, poly_shape)
रेखापुंज पर कार्य यह हमेशा डेटा के साथ एक सूची बनाता है। मैं वास्तव में चाहता हूं कि एक और रेखापुंज फ़ाइल निकाली जाए जिसे मैं आर्कगिस के साथ फिर से लोड कर सकूं। थोड़ा और पढ़ने के बाद मुझे लगता है कि फसल का कार्य वही है जिसकी मुझे वास्तव में आवश्यकता है। लेकिन जब मैं इस फ़ंक्शन का उपयोग करने का प्रयास करता हूं
crop(raster, poly_shape)
मुझे यह त्रुटि मिलती है:
Error in .local(x, y, ...) : extents do not overlap
In addition: Warning message:
In intersect(extent(x), extent(y)) : Objects do not overlap
फ़ाइलें raster और poly_shape दोनों फ़ंक्शन के लिए समान हैं। क्या आप बता सकते हैं कि यहां क्या गलत हो सकता है? क्या यह भी सही है कि फसल कार्य एक और रेखापुंज बनाता है और सूची नहीं?
संपादित करें : हद () फ़ंक्शन मेरे लिए काम नहीं करता है। मुझे अब भी वही त्रुटि मिलती है। लेकिन मुझे यकीन है कि 2 डेटासेट ओवरलैप हैं! उसके साथ
extract(raster, poly_shape)
मुझे इससे सही डेटा मिलता है। बस एक सूची के रूप में और एक रेखापुंज के रूप में नहीं, जैसे कि मैं यह करना चाहता हूं। मैंने पहले आर्कगिस में डेटासेट लोड किया था और वे बहुत अच्छी तरह से फिट थे इसलिए मैंने प्रोजेक्शन की जांच नहीं की। अब मैंने कोशिश की
projection(raster) # "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs"
projection(poly_shape) # "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs"
और आप देख सकते हैं कि अनुमान फिट नहीं हैं। निकालने का कार्य सही तरीके से फ़ाइलों को स्वचालित रूप से बदलने में सक्षम प्रतीत होता है। मुझे पता है कि क्योंकि मैंने निम्नलिखित किया है:
- मैं बहुभुज का सटीक हिस्सा काटता हूं जिसे मैंने आरजीसीआईएस में भी आर में निकाला था
- मैंने निकाले गए आर बहुभुज (सूची) के सभी मूल्यों की राशि की गणना की
- मैंने उन सभी रैस्टर कोशिकाओं के योग की गणना की, जिन्हें मैंने आर्कगिस में काटा
2 का एक ही परिणाम है, इसलिए मुझे लगता है कि निष्कर्ष यह होना चाहिए कि निकालने का कार्य सही था। अब मेरे पास 2 विकल्प हैं जो मुझे लगता है:
- मुझे फिर से निकाली गई सूची से एक रेखापुंज प्राप्त करने का एक तरीका चाहिए या
- 2 डेटासेट (रेखापुंज + पॉली_शैप) को एक ही आपत्ति का उपयोग करने की आवश्यकता होती है और फसल कार्य करना चाहिए
आप यहां क्या करने का सुझाव देंगे?