Як точно працює зональна статистика?


10

Я трохи заплутаний у тому, як працюють інструменти зональної статистики в QGIS (або інструмент "Зональна статистика" під інструментами Растра з інструментарію інструментів QGIS з геоалгоритмами, або "Растрова статистика з полігонів" під векторним інструментом SAGA <-> Растр ).

В основному у мене є шар растрової карти та деякі невеликі багатокутники в окремому шарі. Я хочу призначити кожному багатокутнику значення растру нижче. У випадку, якщо багатокутник перекриває кілька растрових комірок, я просто хочу мати максимум.

Тим не менш, результат не має сенсу, і я не можу відстати від міркувань цих результатів (які, очевидно, не є максимумом основної растрової комірки):

введіть тут опис зображення

Будь-які ідеї, що я роблю неправильно, або те, що я ще не вважав? Результати з API python в qgis.analysis.QgsZonalStatistics()( відповідно до цього питання тут ) набагато кращі, але він пропонує лише запропонувати count, meanі sum... але чомусь немає maxабо min.


1
Щойно спробував це в ArcMap, перетворивши мій Raster в Points, і схоже, що інструмент QGIS Zonal Stats робить те ж саме всередині і використовує лише ті випадки, коли точка перетинає полігон. Однак, схоже, не існує відповідного варіанту відповідності у QGIS (щось порівнянне з тими, що пропонуються в інструменті просторового з'єднання ArcMap, наприклад, INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki

Як підказку для з'ясування цього, я обчислив "підрахунок" для грубого растрового шару (приблизно один піксель знаходиться в полігоні) із зоновою статистикою QGIS і отримав значення 0,45 - на диво не ціле число. Це говорить про те, що QGIS може якось зважувати значення пікселів. Це може бути причиною того, що макс / хв не відповідає тому, що ви очікуєте - можливо, погана реалізація цього зважування при обчисленні макс / хв.
Джон

Див. Gis.stackexchange.com/a/281753/18189 для опису того, як працює алгоритм зональної статистики QGIS.
дбастон

Відповіді:


4

Оновлення у квітні 2018 року та QGIS 2.18

Поточна версія плагіну зональної статистики може містити:

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

Оригінальна відповідь

Коротше кажучи, такий підхід із QGIS, наскільки я знаю, не в змозі переглянути максимальний растровий піксель. Тільки підрахунок, середня сума та сума. Я надав деякі ресурси, які охоплюють це.

Ось сторінка doc для плагіна Zonal Statistics Plugin

За допомогою плагіна icon_zonal_statistics Zonal statistics ви можете проаналізувати результати тематичної класифікації. Це дозволяє обчислити кілька значень пікселів растрового шару за допомогою полігонального векторного шару (див. Рисунок_zonal_statistics). Ви можете обчислити суму, середнє значення та загальну кількість пікселів, які знаходяться у полігоні. Плагін генерує вихідні стовпці у векторному шарі з визначеним користувачем префіксом.

Також цей підручник QGIS відбору проб растрових даних за допомогою точок або полігонів надає вибіркові дані для досягнення зонального повернення статистики. Я включив скріншот, який вони надають як приклад їх робочого процесу та результатів.

Увімкніть плагіни зональної статистики. Це основний плагін, тому він уже встановлений. Див. Розділ Використання плагінів, щоб дізнатися, як увімкнути основний плагін

введіть тут опис зображення

введіть тут опис зображення

введіть тут опис зображення

Після завершення обробки виберіть шар графств. Скористайтеся інструментом «Ідентифікувати» та натисніть на будь-який багатокутник округу. Ви побачите три нові атрибути, додані до шару: ZS_count, ZS_mean та ZS_sum. Ці атрибути містять кількість растрових пікселів, середнє значення піксельних растрових пікселів та суму значень пікселів растрових відповідно. Оскільки нас цікавить середня температура, поле ZS_mean буде саме таким.

введіть тут опис зображення


2

Щодо QGIS 3.0 (можливо, раніше, я не впевнений), доступна наступна статистика:

  • Рахувати
  • Сума
  • Середній
  • Середня
  • Стандартне відхилення
  • Мінімум
  • Максимум
  • Дальність
  • Меншість
  • Більшість (режим)
  • Різноманітність
  • Варіантність
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.