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


12

ラスタデータをポリゴンデータに変換しようとしています。以前は、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

1
TIFFを使用して土地分類を行う場合、教師なし最尤分類を使用してTIFFをポリゴンに変換する方法があるリンク(gis.stackexchange.com/q/29570/8104)を確認してください。
アーロン

@Aaronに感謝しますが、このプロジェクトのオープンソースで作業する必要があります
-Sam007

一度に200 TIFを変換する方法はありますか?
ベンジョー

回答:


21

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バインディングがあります。


マイクありがとう。試してみましょう。また、シェープファイルとして出力するのではなく、これを直接PostGISに割り当て、Geoserverからアクセスできます。
Sam007

@ Sam007はい、PostGIS出力も動作します。更新を参照してください
マイクT

シェープファイルを作成する理由はこのエラーです。postgisデータベースに追加しませんでした。dpaste.de/L8tSK GeoTIFFが別のディレクトリにある場合、問題になりますか?
-Sam007

これは、Pythonで使用されるSWIGバインディングに関連する気を散らすエラーのように見えます。使用しているGDALのバージョンとインストール方法は?
マイクT

1
@ Sam007試してみてくださいsudo apt-get install python-gdal。それでもうまくいかない場合は、新しい質問を投稿してください。OpenGeoでの動作の有無をテストすることはまだできていません。
マイクT
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.