Conversión de datos ráster a datos poligonales

Estoy tratando de convertir datos ráster en datos poligonales. Anteriormente estaba usando shapefile.py para crear datos de puntos en lat / lon individuales en el ráster. Pero eso creará una gran cantidad de datos y su consulta será un gran problema.

Por lo tanto, sería mejor si pudiera convertir los datos ráster en datos de polígono. Pude encontrar SAGA , pero parece que solo puede convertir los archivos .grd . ¿Hay alguna forma de convertir el formato GeoTIFF o datasets ráster en general a Polygon?

ACTUALIZACIÓN: Uso de gdal_polygonize.py

Estoy ejecutando la siguiente línea de código,

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

Pero recibo el siguiente error. Mi versión gdal es 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

Prueba gdal_polygonize.py . Por ejemplo, para crear un Shapefile mylayer.shpdesde la entrada myraster.tif:

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

o para enviar directamente a una base de datos PostGIS (vea los detalles del controlador PostgreSQL / PostGIS para GDAL ):

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

(ignore el mensaje "ERROR 1: No se encontraron definiciones de campo para 'mylayer', ¿es una tabla?")

Hay enlaces de Python debajo de GDALPolygonize , si necesita esta herramienta personalizada.