การแปลงข้อมูล Raster เป็นข้อมูลรูปหลายเหลี่ยม

ฉันพยายามแปลงข้อมูลแรสเตอร์เป็นข้อมูลรูปหลายเหลี่ยม ก่อนหน้านี้ฉันใช้ shapefile.py เพื่อสร้างข้อมูลจุดบน lat / lon เดี่ยว ๆ บน raster แต่นั่นจะสร้างข้อมูลขนาดใหญ่และการสืบค้นมันจะเป็นปัญหาใหญ่

ดังนั้นจะดีกว่าถ้าฉันสามารถแปลงข้อมูลแรสเตอร์เป็นข้อมูลรูปหลายเหลี่ยม ฉันสามารถค้นหาSAGAได้ แต่ดูเหมือนว่าจะสามารถแปลงไฟล์ . grdได้เท่านั้น มีวิธีแปลงรูปแบบ GeoTIFF หรือชุดข้อมูลแรสเตอร์โดยทั่วไปเป็นรูปหลายเหลี่ยมหรือไม่?

อัปเดต: การใช้ gdal_polygonize.py

ฉันกำลังเรียกใช้โค้ดต่อไปนี้

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

แต่ฉันได้รับข้อผิดพลาดดังต่อไปนี้ รุ่น gdal ของฉันคือ 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

ลองgdal_polygonize.py เช่นการสร้าง Shapefile mylayer.shpจากอินพุตmyraster.tif:

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

หรือส่งออกไปยังฐานข้อมูล PostGIS โดยตรง (ดูรายละเอียดไดรเวอร์ PostgreSQL / PostGIS สำหรับ GDAL ):

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

(ไม่สนใจข้อความ "ข้อผิดพลาด 1: ไม่พบคำจำกัดความของฟิลด์สำหรับ 'mylayer' เป็นตารางหรือไม่")

มีการผูก Python ไว้ใต้GDALPolygonizeหากคุณต้องการเครื่องมือนี้เอง