ラスターデータをポリゴンデータに変換する

ラスタデータをポリゴンデータに変換しようとしています。以前は、shapefile.pyを使用して、ラスター上の個々の緯度/経度にポイントデータを作成していました。しかし、それは巨大なデータを作成し、それをクエリすることは大きな問題になります。

したがって、ラスターデータをポリゴンデータに変換できればさらに良いでしょう。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を試してください。たとえば、mylayer.shp入力からシェープファイルを作成するにはmyraster.tif：

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

または、PostGISデータベースに直接出力します（GDALのPostgreSQL / PostGISドライバーの詳細を参照）。

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

（「エラー1： 'mylayer'のフィールド定義が見つかりません。テーブルですか？」というメッセージは無視してください）

このツールをカスタマイズする必要がある場合は、GDALPolygonizeの下にPythonバインディングがあります。