Преобразование растровых данных в полигональные

Я пытаюсь преобразовать растровые данные в полигональные. Ранее я использовал shapefile.py для создания точечных данных по отдельным широтам / долготам на растре. Но это создаст огромные данные, и их запрос станет большой проблемой.

Таким образом, было бы лучше, если бы я мог преобразовать растровые данные в данные многоугольника. Мне удалось найти SAGA , но, похоже, он может конвертировать только файлы .grd . Есть ли способ преобразовать формат GeoTIFF или наборы растровых данных в целом в Polygon?

ОБНОВЛЕНИЕ: Использование 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 », это таблица?»)

Под GDALPolygonize есть привязки Python , если вам нужно настроить этот инструмент.