Konwertowanie danych rastrowych na dane wielokąta

Próbuję przekonwertować dane rastrowe na dane wielokąta. Wcześniej używałem shapefile.py do tworzenia danych punktowych na poszczególnych lat / lonach na rastrze. Ale to stworzy ogromne dane, a zapytania będą dużym problemem.

Byłoby lepiej, gdybym mógł przekonwertować dane rastrowe na dane wielokąta. Udało mi się znaleźć SAGA , ale wygląda na to, że można konwertować tylko pliki .grd . Czy istnieje sposób na konwersję formatu GeoTIFF lub ogólnie zestawów danych rastrowych na Polygon?

AKTUALIZACJA: Korzystanie z gdal_polygonize.py

Wykonuję następujący wiersz kodu,

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

Ale pojawia się następujący błąd. Moja wersja gdal to 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

Spróbuj gdal_polygonize.py . Na przykład, aby utworzyć plik Shapefile na mylayer.shppodstawie danych wejściowych myraster.tif:

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

lub do wyjścia bezpośrednio do bazy danych PostGIS (zobacz szczegóły sterownika PostgreSQL / PostGIS dla GDAL ):

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

(zignoruj ​​komunikat „BŁĄD 1: Nie znaleziono definicji pól dla„ mylayer ”, czy to tabela?)

Istnieją powiązania Python poniżej GDALPolygonize , jeśli potrzebujesz tego narzędzia dostosowanego.