ओपन सोर्स टूल्स का उपयोग करके सैटेलाइट इमेज की इमेज बाउंड्री / फुटप्रिंट की गणना?

Question

मुझे कई एकल परत रेखापुंज छवियों की बहुभुज रूपरेखाओं को बनाने की आवश्यकता है, न कि हद / बाउंडिंग बॉक्स, बल्कि नोदाटा मानों के बिना क्षेत्र, जैसा कि यहाँ दिखाया गया है: रेखापुंज के पदचिह्न दिखाते हुए आकार आकृति बनाना? ।

ऊपर दिए गए प्रश्न के उत्तर में, छवि सीमा प्लगइन का उल्लेख किया गया था, हालांकि, मैं इसे QGIS 1.8.0 लिस्बोआ में नहीं मिला, जो उबंटू में चल रहा है।

उपकरण अभी भी उपलब्ध है?

यदि नहीं, तो क्या ऐसा करने का एक तरीका है ओपनसोर्स टूल, या तो गाल्ड, आर, क्यूजीआईएस, ग्रैस, या इसी तरह का उपयोग करना, आर्कपैक के साथ नहीं?

Comments

संबंधित: gis.stackexchange.com/questions/9074/…

— कॉमरेड चे

Answer 1

आप आधिकारिक रिपॉजिटरी में पुराने रिपॉजिटरी को प्लग इन कर सकते हैं।

http://pyqgis.org/repo/contributed

छवि सीमा प्लगइन में वैध पिक्सेल की गणना के लिए विकल्प है, लेकिन, छवि को पूर्ण दृश्य होने की आवश्यकता है, कैसे सीबर्स या लैंडसैट, जहां यह प्रक्रिया पहले 4 कोनों की गणना करती है। इमेज को बिना इमेजिंग (सैटेलाइट सेंसर द्वारा स्कैन किए गए) क्षेत्र के लिए शून्य मान के साथ नोडा की आवश्यकता है।

छवि सीमा के लेखक

Comments

हाय @ लमोट्टा, रिपॉजिटरी लिंक के लिए धन्यवाद। इस समय मैं जिन छवियों पर काम कर रहा हूं वे सभी एसएआर चित्र हैं, एनए मान (या क्यूजी में मूल्य उपकरण के अनुसार -3.4e + 38)। यदि कोई एकमात्र विकल्प के रूप में शून्य होने के बजाय टूल में मैन्युअल रूप से नोडाटा मूल्य जोड़ सकता है, तो यह बहुत अच्छा होगा, अन्यथा, मैं अपने एनए मूल्यों को शून्य का नया मान देने के लिए पहले आर में एक स्क्रिप्ट चला सकता हूं। यह एक महान उपकरण है, हालांकि, धन्यवाद!

— spib

क्या छवि सीमा के लिए एक नया भंडार है? लिंक उपलब्ध कराने 404. है हो सकता है कि निम्नलिखित ?: github.com/lmotta/imagefootprint_plugin

— हारून

क्या कहीं भी हम छवि सीमा के लिए प्रलेखन पा सकते हैं?

— लूनुह

Answer 2

मुझे आखिरकार ऐसा करने का एक तरीका मिल गया:

चरण 1: gdalwarp -dstnodata 0 -dstalpha -of GTiff foo1 foo2

यह दो महत्वपूर्ण चीजें करता है: यह गंतव्य डेटा (सीमा के बाहर) को 0 पर सेट करता है , और यह एक अल्फा बैंड बनाता है।

चरण 2: gdal_polygonize.py foo2 -b 2 -f "ESRI Shapefile" foo3

दूसरा चरण 1 में बनाए गए अल्फा बैंड (बैंड 2) का उपयोग करता है, और उस बैंड से एक आकृति बनाता है।

इसके बाद आसानी से बैश स्क्रिप्ट में स्क्रिप्ट की जा सकती है, यदि आपके पास कई चित्र हैं, तो सटीक रूपरेखा बनाने के लिए।  

Comments

शानदार, यह समाधान किसी भी प्लगइन्स या सॉफ़्टवेयर पर निर्भर नहीं करता है। NoData के प्रति संवेदनशील होने वाले रेखापुंज पदचिह्न बनाना उन जीआईएस कार्यों में से एक है जो उतना सरल नहीं है जितना लगता है कि यह होना चाहिए।

— चार्ली Parr

जब मैं इस पद्धति का उपयोग करता हूं तो मैं बहुभुजों से भरी आकृति के साथ समाप्त होता हूं जो मुझे लगता है कि एक साथ समूहीकृत विभिन्न रंगों का प्रतिनिधित्व करते हैं। मैं इससे कैसे बच सकता हूं और केवल कुछ बहुभुज हां-डेटा का प्रतिनिधित्व कर सकता हूं ??

— लूनुह

Answer 3

मैं प्रयोग किया जाता है gdal_translateके रूप में सुझाव दिया GDAL परियोजना द्वारा।

gdal_translate -b mask -of vrt -a_nodata 0 test.tif test.vrt
# Note the  -a_nodata 0 doesn't seem to work when the mask is input, so do another pass
gdal_translate -b 1 -of vrt -a_nodata 0 test.vrt test2.vrt
gdal_polygonize.py -q  -8 test2.vrt -b 1 -f "ESRI Shapefile" testdata.shp

Comments

GDAL ट्रंक संस्करण के साथ अब स्रोत छवि osgeo-org.1560.x6.nabble.com/ ... के खिलाफ सीधे gdal_polygonize का उपयोग करना संभव है । हालांकि, बिना डेटा के बहुभुज को हटाने के लिए कुछ सफाई की आवश्यकता होती है।

— user30184

Answer 4

छवि सीमा प्लगइन मेरे लिए भी काम नहीं करता था, इसलिए मैंने GDAL के साथ एक ही दृष्टिकोण का उपयोग किया। फिर भी यह मेरे लिए पहला कदम बदलने के बाद ही काम आया:

चरण 1 : gdalwarp -srcnodata 0 -dstalpha -of GTiff foo1 foo2

मैं लैंडसेट 8 बैंड (जहां कोई डेटा नहीं है 0) के साथ काम कर रहा हूं और जब -dstnodataमुझे मिलता है तो फंक्शन का उपयोग करते हैं:

band1 साथ कोई डेटा = 'कोई डाटा नहीं'
band2 पूरे दृश्य / हद तक के लिए (अल्फा बैंड) = '255'

-srcnodataफ़ंक्शन के साथ मुझे मिलता है:

band1 साथ कोई डेटा = 'कोई डाटा नहीं'
band2 कोई डेटा = 'कोई डाटा नहीं' और मान्य डेटा = '255' जो तब मान्य डेटा के लिए ही निकालने बहुभुज की अनुमति देता है के साथ (अल्फा बैंड)।

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

Answer 5

छवि सीमा अप्रचलित रह गई। आप फ़ुटप्रिंट प्लगइन का उपयोग कर सकते हैं। http://pt.slideshare.net/LuizMotta3/qgis-ibama-imagefootprint