Bagaimana cara kerja Statistik Zonal?

10

Saya agak bingung tentang cara kerja alat statistik zonal di QGIS (baik alat 'Statistik Zonal' di bawah alat Raster dari kotak geoalgoritma QGIS atau alat 'Statistik raster dari poligon' di bawah SAGA Vector <-> Alat raster ).

Pada dasarnya saya memiliki layer peta raster dan beberapa poligon kecil di layer terpisah. Saya ingin memberi setiap poligon nilai raster di bawah ini. Jika poligon tumpang tindih dengan beberapa sel raster, saya hanya ingin memiliki maksimum.

Namun hasilnya tidak masuk akal dan saya tidak bisa mendapatkan alasan di balik hasil tersebut (yang jelas tidak maksimal dari sel raster underlain):

masukkan deskripsi gambar di sini

Adakah ide yang saya lakukan salah atau yang belum saya pertimbangkan? Hasil dengan API python di qgis.analysis.QgsZonalStatistics()( sesuai dengan pertanyaan ini di sini ) jauh lebih baik, tetapi hanya menawarkan count, meandan sum... tetapi karena alasan tertentu tidak ada maxatau min.

qgis raster layers saga zonal-statistics GeoEki
sumber
1
Baru saja mencoba ini di ArcMap dengan mengonversi Raster saya ke Poin dan sepertinya alat Statistik QGIS Zonal melakukan hal yang sama secara internal dan hanya menggunakan kasus di mana titik memotong poligon. Namun, sepertinya tidak ada opsi pencocokan yang sesuai di QGIS (sesuatu yang sebanding dengan yang ditawarkan dalam alat Spatial Join ArcMap, misalnya INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki
Sebagai petunjuk untuk mencari tahu ini, saya menghitung "hitungan" untuk lapisan raster kasar (kira-kira satu piksel dalam poligon) dengan Statistik Zonal QGIS dan mendapat nilai 0,45 - yang mengejutkan bukan bilangan bulat. Ini menunjukkan bahwa QGIS entah bagaimana mungkin menimbang nilai piksel. Ini mungkin mengapa max / min tidak sesuai dengan apa yang Anda harapkan - mungkin implementasi yang buruk dari pembobotan ini ketika menghitung maks / min.
Jon
Lihat gis.stackexchange.com/a/281753/18189 untuk deskripsi tentang bagaimana algoritma statistik zona QGIS bekerja.
dbaston

Jawaban:

4

Pembaruan sebagai April 2018 dan QGIS 2.18

Versi Plugin Statistik Zonal saat ini dapat memberikan:

Count: to count the number of pixels
Sum: to sum the pixel values
Mean: to get the mean of pixel values
Median: to get the median of pixel values
StDev: to get the standard deviation of pixel values
Min: to get the minimum of pixel values
Max: to get the maximum of pixel values
Range: to get the range (max - min) of pixel values
Minority: to get the less represented pixel value
Majority: to get the most represented pixel value
Variety: to count the number of distinct pixel values

Jawaban asli

Singkatnya, pendekatan ini dengan QGIS, sejauh yang saya ketahui tidak mungkin untuk melihat maksimum pixel raster sebagai hasilnya. Hanya hitungan, berarti dan jumlah. Saya telah menyediakan beberapa sumber daya yang mencakup hal ini.

Ini adalah halaman dokumen untuk plugin Zonal Statistics Plugin

Dengan plugin statistik icon_zonal_statistics Zonal, Anda dapat menganalisis hasil klasifikasi tematik. Ini memungkinkan Anda untuk menghitung beberapa nilai piksel dari layer raster dengan bantuan layer vektor poligonal (lihat gambar_zonal_statistics). Anda dapat menghitung jumlah, nilai rata-rata dan jumlah total piksel yang ada dalam poligon. Plugin menghasilkan kolom output di lapisan vektor dengan awalan yang ditentukan pengguna.

Juga, Tutorial QGIS ini Mengambil Sampel Data Raster menggunakan Poin atau Poligon menyediakan data sampel untuk mencapai pengembalian statistik zona. Saya telah menyertakan tangkapan layar yang mereka berikan sebagai contoh alur kerja dan hasilnya.

Aktifkan Plugin Statistik Zonal. Ini adalah plugin inti sehingga sudah diinstal. Lihat Menggunakan Plugin untuk mengetahui cara mengaktifkan plugin inti

masukkan deskripsi gambar di sini

masukkan deskripsi gambar di sini

masukkan deskripsi gambar di sini

Setelah pemrosesan selesai, pilih layer county. Gunakan alat Identifikasi dan klik pada poligon county. Anda akan melihat tiga atribut baru ditambahkan ke layer: ZS_count, ZS_mean dan ZS_sum. Atribut ini berisi jumlah piksel raster, rata-rata nilai piksel raster, dan jumlah nilai piksel raster. Karena kami tertarik pada suhu rata-rata, bidang ZS_mean akan menjadi yang digunakan.

masukkan deskripsi gambar di sini

Whyzar
sumber
2

Pada QGIS 3.0 (mungkin lebih awal, saya tidak yakin), statistik berikut tersedia:

  • Menghitung
  • Jumlah
  • Berarti
  • Median
  • Simpangan baku
  • Minimum
  • Maksimum
  • Jarak
  • Minoritas
  • Mayoritas (mode)
  • Variasi
  • Perbedaan
dbaston
sumber