Перетворення растрових даних у дані багатокутника

Я намагаюся перетворити растрові дані в полігонові дані. Раніше я використовував shapefile.py для створення точкових даних про окремі lat / lon на растрі. Але це створить величезні дані і запити їх будуть великою проблемою.

Отже, було б краще, якщо я зможу перетворити растрові дані в дані багатокутника. Мені вдалося знайти 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 ", це таблиця?")

Під GDALPolygonize є прив'язки Python , якщо вам потрібен цей інструмент.