タグ付けされた質問 「zonal-statistics」

3
ArcGIS Zonal統計の主要なバグですか?
更新:バグはArcGIS 10.4リリースで修正されました ArcGIS 10.2.2を使用して、いくつかのゾーンのゾーン統計を決定しています。値ラスターにNoDataがある場合、正確にツールの説明で通知されているように、ゾーンの結果が「NoData」になるようにします。このツールの説明の状態: DATA —特定のゾーン内では、入力値ラスターに値を持つセルのみが、そのゾーンの出力値の決定に使用されます。値ラスターのNoDataセルは、統計計算では無視されます。 NODATA —特定のゾーン内で、値ラスタにNoDataセルが存在する場合、そのゾーン内のすべてのセルの統計計算を実行するには情報が不十分であると見なされます。したがって、ゾーン全体が出力ラスタでNoData値を受け取ります。 この写真で私のセットアップを見てください: 1つのNoDataピクセルを持つ値ラスターでNODATAオプションを使用しているため、結果のゾーン値(ゾーン61154)が「NoData」になると予想されます。代わりに、12.74の値(画像では13に四捨五入)を取得します。これは2つのレベルで混乱させます。1つ目は「NoData」、2つ目は、結果が12.74になることです。値ラスターの最大値(この場合は10)よりも小さい。 DATAオプションを使用している場合、約9.1の値が得られますが、これは理にかなっています。さまざまなデータセット、コンピューター、およびArcGISバージョンでこれをテストしました。 ここで何が欠けていますか? 編集/追加コメント:「Count」属性もその特定のゾーンに対して間違っていることに気付きました。そのゾーンには実際に421個のセルがありますが、このツールは297個しかカウントしません。421から297個を計算すると124個になります。ゾーン内で。ツールがセル数を間違えている(低すぎる)可能性があり、これが平均の増加を説明している可能性があります。 編集:これは私が使用しているデータへのリンクです。 編集:ダン・パターソンと私はさらにいくつかのデバッグをした、ここで ESRIのフォーラムで。

2
ゾーン統計はどのように正確に機能しますか?
ゾーン統計ツールがQGISでどのように機能するかについて少し混乱しています(QGISジオアルゴリズムツールボックスのラスターツールの「ゾーン統計」ツールまたはSAGAベクトル<->ラスターツールの「ポリゴンのラスター統計」ツール) )。 基本的に、ラスターマップレイヤーといくつかの小さなポリゴンが別のレイヤーにあります。各ポリゴンに下のラスターの値を割り当てます。ポリゴンが複数のラスターセルに重なっている場合は、最大値が欲しいだけです。 しかし、出力は意味をなさず、それらの結果の推論(下にあるラスターセルの最大値ではないことは明らかです)を後回しにすることはできません。 私が間違っていることや、まだ考慮していないことについてのアイデアはありますか?でPythonのAPIとの結果qgis.analysis.QgsZonalStatistics()(ここでは、この質問に応じては)あまり良くありませんが、それだけの申し出を行いcount、meanそしてsum...しかし、何らかの理由のためのmaxかmin。

3
PythonコンソールからQGIS Zonal Statsプラグインを使用していますか?
QGISのゾーン統計プラグインを使用して、オーバーレイポリゴンからラスター統計を抽出しています(種の分布のシェープファイルがあり、各種の範囲内から環境データを抽出したい)。データを取得する必要のあるファイルが300個ほどあるので、Pythonコンソールで実行するスクリプトを記述したいのですが、Pythonの初心者であり、これを行う方法がわかりません。

3
Rを使用してポリゴン内の非NAラスターセルの数を取得する方法
私はArcGIS ZonalStatsを使用してあらゆる種類の問題に遭遇しており、Rが優れた方法であると考えました。私はRにかなり慣れていると言いますが、コーディングの背景を得ました。 状況としては、いくつかのラスターと、さまざまなサイズの多くのフィーチャを持つポリゴンシェープファイルがあります(すべてのフィーチャはラスターセルよりも大きく、ポリゴンフィーチャはラスターに位置合わせされています)。抽出付きのラスターライブラリを使用して、各ポリゴンフィーチャの平均値を取得する方法を理解しました。 #load packages required require(rgdal) require(sp) require(raster) # ---Set the working directory------- datdir <- "/test_data/" #Read in grid of water depth ras <- raster("test_data/raster/pl_sm_rp1000/w001001.adf") #read in polygon shape file proxNA <- shapefile("test_data/proxy/PL_proxy_WD_NA_test") #calc mean depth per polygon feature #unweighted - only assigns grid to district if centroid is …
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.