Chuyển đổi dữ liệu raster thành dữ liệu đa giác

Tôi đang cố gắng chuyển đổi dữ liệu raster thành dữ liệu đa giác. Trước đây tôi đã sử dụng shapefile.py để tạo dữ liệu điểm trên lat / lon riêng lẻ trên raster. Nhưng điều đó sẽ tạo ra một dữ liệu khổng lồ và truy vấn nó sẽ là một vấn đề lớn.

Vì vậy, sẽ tốt hơn nếu tôi có thể chuyển đổi dữ liệu raster thành dữ liệu đa giác. Tôi đã có thể tìm thấy SAGA , nhưng dường như nó chỉ có thể chuyển đổi các tệp .grd . Có cách nào để chuyển đổi định dạng GeoTIFF hoặc bộ dữ liệu raster nói chung sang Đa giác không?

CẬP NHẬT: Sử dụng gdal_polygonize.py

Tôi đang thực thi dòng mã sau đây,

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

Nhưng tôi nhận được lỗi sau. Phiên bản gdal của tôi là 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

Hãy thử gdal_polygonize.py . Ví dụ: để tạo Shapefile mylayer.shptừ đầu vào myraster.tif:

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

hoặc để xuất trực tiếp vào cơ sở dữ liệu PostGIS (xem chi tiết trình điều khiển PostgreSQL / PostGIS cho GDAL ):

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

(bỏ qua thông báo "LRI 1: Không tìm thấy định nghĩa trường nào cho thông báo 'mylayer', đây có phải là bảng không?"

Có các ràng buộc Python bên dưới GDALPolygonize , nếu bạn cần công cụ này tùy chỉnh.