Convertendo dados de varredura em dados de polígono

Estou tentando converter dados rasterizados em dados poligonais. Anteriormente, eu estava usando o shapefile.py para criar dados de pontos em lat / lon individuais na varredura. Mas isso criará uma enorme quantidade de dados e a consulta será um grande problema.

Portanto, seria melhor se eu pudesse converter os dados rasterizados em dados poligonais. Consegui encontrar o SAGA , mas parece que ele só pode converter os arquivos .grd . Existe uma maneira de converter o formato GeoTIFF ou conjuntos de dados raster em geral para Polygon?

ATUALIZAÇÃO: Usando gdal_polygonize.py

Estou executando a seguinte linha de código,

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

Mas estou recebendo o seguinte erro. Minha versão 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

Tente gdal_polygonize.py . Por exemplo, para criar um Shapefile a mylayer.shppartir da entrada myraster.tif:

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

ou para enviar diretamente para um banco de dados PostGIS (consulte os detalhes do driver PostgreSQL / PostGIS para GDAL ):

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

(ignore a mensagem "ERRO 1: Nenhuma definição de campo encontrada para 'mylayer', é uma tabela?")

Existem ligações Python abaixo do GDALPolygonize , se você precisar dessa ferramenta personalizada.