बहुभुज डेटा को रेखापुंज डेटा परिवर्तित करना

मैं एक रेखापुंज डेटा को बहुभुज डेटा में बदलने की कोशिश कर रहा हूं। पहले मैं रैपर पर अलग-अलग लैट / लोन पर पॉइंट डेटा बनाने के लिए शेपफाइलहोम का उपयोग कर रहा था। लेकिन यह एक बड़ा डेटा तैयार करेगा और इसे क्वेरी करना एक बड़ा मुद्दा होगा।

इसलिए, यह बेहतर होगा कि मैं रैस्टर डेटा को बहुभुज डेटा में परिवर्तित करने में सक्षम हूं। मैं सागा को खोजने में सक्षम था , लेकिन ऐसा लगता है कि यह केवल .grd फ़ाइलों को परिवर्तित कर सकता है । क्या सामान्य रूप से बहुभुज में GeoTIFF प्रारूप या रेखापुंज डेटासेट को बदलने का एक तरीका है?

अद्यतन: gdal_polygonize.py का उपयोग करना

मैं निम्नलिखित लाइन ऑफ कोड निष्पादित कर रहा हूं,

gdal_polygonize.py ~/development/Biodiversity/biodiversity/GeoTIFF/bio11.tif -f "ESRI Shapefile" bio11.shp

लेकिन मुझे निम्नलिखित त्रुटि मिल रही है। मेरा गदल संस्करण 1.9.0 है।

Traceback (most recent call last):
  File "/usr/bin/gdal_polygonize.py", line 34, in <module>
    import gdal, ogr, osr
  File "/usr/lib/python2.7/dist-packages/gdal.py", line 2, in <module>
    from osgeo.gdal import deprecation_warn
  File "/usr/lib/python2.7/dist-packages/osgeo/__init__.py", line 21, in <module>
    _gdal = swig_import_helper()
  File "/usr/lib/python2.7/dist-packages/osgeo/__init__.py", line 17, in swig_import_helper
    _mod = imp.load_module('_gdal', fp, pathname, description)
ImportError: libjvm.so: cannot open shared object file: No such file or directory

raster gdal polygon-creation

— Sam007
स्रोत

यदि आप टिफ़्स के साथ भूमि वर्गीकरण कर रहे हैं, तो लिंक देखें ( gis.stackexchange.com/q/29570/8104 ) जहां अनपेक्षित अधिकतम संभावना वर्गीकरण का उपयोग करके टिफ़ को पॉलीगॉन में बदलने की विधि है।

— आरोन

धन्यवाद @Aaron, लेकिन मैं इस परियोजना पर खुले स्रोत पर काम करने की जरूरत है

— Sam007

वहाँ एक बार में 200 tifs कन्वर्ट करने के लिए एक रास्ता है?

— बेन्जो

Gdal_polygonize.py का प्रयास करें । उदाहरण के लिए, mylayer.shpइनपुट से शेपफाइल बनाने के लिए myraster.tif:

gdal_polygonize.py myraster.tif -f "ESRI Shapefile" mylayer.shp

या सीधे पोस्टगिस डेटाबेस में आउटपुट करने के लिए ( GDAL के लिए PostgreSQL / PostGIS ड्राइवर विवरण देखें ):

gdal_polygonize.py myraster.tif -f PostgreSQL PG:"dbname='postgis' user='postgres'" mylayer

("ERROR 1: कोई फ़ील्ड परिभाषा जो 'मायलेयर' के लिए नहीं मिली है, उसे अनदेखा करें, क्या यह एक तालिका है?" संदेश)

GDALPolygonize के नीचे पायथन बाइंडिंग हैं , यदि आपको इस उपकरण को अनुकूलित करने की आवश्यकता है।

— माइक टी
स्रोत

धन्यवाद माइक। मुझे इसकी कोशिश करने दीजिए। इसके अलावा, इसे आकार देने के रूप में आउटपुट करने के बजाय, क्या मैं इसे सीधे PostGIS में असाइन कर सकता हूं और इसे जियोसेवर से एक्सेस कर सकता हूं।

— सैम ००।

@ Sam007 yup, PostGIS आउटपुट भी काम करता है, अद्यतन देखें

— माइक टी।

आकृति बनाने में यह त्रुटि क्यों हुई। पोस्टगिस डेटाबेस में नहीं जोड़ा गया। dpaste.de/L8tSK क्या यह समस्या है यदि GeoTIFF दूसरे डायर में है?

— सैम ००

यह एक विचलित करने वाली त्रुटि की तरह दिखता है, जो पायथन द्वारा उपयोग किए गए SWIG बाइंडिंग से संबंधित है। GDAL का कौन सा संस्करण आपके पास है, और आपने इसे कैसे स्थापित किया?

— माइक टी।

@ Sam007 कोशिश करें sudo apt-get install python-gdal, लेकिन अगर वह काम नहीं करता है, तो हाँ एक नया प्रश्न पोस्ट करें। मैं अभी भी परीक्षण के आसपास नहीं मिला है कि यह OpenGeo के साथ कैसे काम करता है या नहीं।

— माइक टी