Zonal Statistics ツール (ArcGIS Spatial Analyst) のオープン ソースの代替手段はありますか? ラスターファイルの統計をカウントするための最良のツール(スクリプトで使用できる)は何ですか?
4 に答える
Rでラスターパッケージを使用できます
library(raster)
v <- raster('raster filename')
z <- raster('zones raster filename')
zv <- zonal(v, z, fun=mean)
間違っていたら訂正してください、RobertH ですが、ゾーンが何らかの意味で既に「ラスタライズ」されている必要があると思いますが、多くの場合、ポリゴン内にあるラスタ セルの統計が必要になります。これには sp パッケージ内の R のさまざまなオーバーレイ メソッド (参照: ?"overlay-methods" ) が必要ですが、私が間違っている場合は喜んで聞いてくれます。私は、SpatialGridsDataFrames を使用するよりもラスター パッケージをかなり好みますが、ポリゴンとグリッド データを混在させるには sp クラスに依存する必要があると思います。これは問題ありませんが、ラスター パッケージの優れたメモリ管理が不足しているため、大きなラスターで R でポイント イン ポリゴン スタイルの操作を実行するのが非常に難しくなるため、問題が発生します。
また、これはGRASS内で、および/またはQGISの次のリリース(1.7)で何らかの組み込みのゾーン統計機能を備えたQGISを介して実行できると信じていますが、試していません。
Rasterstats パッケージは、私にとってはうまく機能する素晴らしいオープン ソース ツールです: http://blog.perrygeo.net/2013/09/24/python-raster-stats/
arcpy の ZonalStatistics メソッドが、ラスターを配列に変換しようとすると奇妙なエラーにつながる問題のあるラスターを生成していたため、回避策として使用し始めました ( https://gis.stackexchange.com/questions/110274/save- fails-on-raster-object-created-from-numpyarraytoraster )。Rasterstats はうまく機能し、私の問題に対する効率的な解決策を提供してくれました。