Jak dokładnie działa statystyki strefowe?

10

Jestem trochę zdezorientowany, jak działają narzędzia statystyk strefowych w QGIS (narzędzie „Statystyka strefowa” w narzędziach Raster z zestawu narzędzi geoalgorytmowych QGIS lub narzędzie „Statystyki rastrowe z wielokątów” w narzędziach SAGA Vector <-> Raster ).

Zasadniczo mam warstwę mapy rastrowej i kilka małych wielokątów w osobnej warstwie. Chcę przypisać każdemu wielokątowi wartość rastra poniżej. W przypadku gdy wielokąt nachodzi na wiele komórek rastrowych, chcę tylko mieć maksimum.

Jednak dane wyjściowe nie mają żadnego sensu i nie mogę się oprzeć za uzasadnieniem tych wyników (które oczywiście nie są maksimum podrzędnej komórki rastrowej):

wprowadź opis zdjęcia tutaj

Jakieś pomysły, co robię źle lub czego jeszcze nie rozważałem? Wyniki z API Pythona qgis.analysis.QgsZonalStatistics()( zgodnie z tym pytaniem tutaj ) są znacznie lepsze, ale nie tylko count, meana sum... ale z jakiegoś powodu nie maxlub min.

GeoEki
źródło
1
Właśnie wypróbowałem to w ArcMap, konwertując mój Raster na punkty. Wygląda na to, że narzędzie QGIS Zonal Stats robi to samo wewnętrznie i wykorzystuje tylko te przypadki, w których punkt przecina wielokąt. Wydaje się jednak, że nie ma odpowiedniej opcji dopasowania w QGIS (coś porównywalnego z tymi oferowanymi w narzędziu ArcMap's Spatial Join, np. INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki,
Jako wskazówkę, aby to rozgryźć, obliczyłem „liczbę” grubej warstwy rastrowej (w przybliżeniu jeden piksel znajduje się w wielokącie) za pomocą statystyk strefowych QGIS i otrzymałem wartość 0,45 - zaskakująco nie jest liczbą całkowitą. To sugeruje, że QGIS może w jakiś sposób ważyć wartości pikseli. Może to być przyczyną, dla której maks./min nie odpowiada oczekiwaniom - być może słaba implementacja tego ważenia przy obliczaniu maks./min.
Jon
Zobacz gis.stackexchange.com/a/281753/18189 dla opisu jak QGIS strefowego statystyki algorytm działa.
dbaston

Odpowiedzi:

4

Aktualizacja w kwietniu 2018 r. I QGIS 2.18

Obecna wersja wtyczki Zonal Statistics może zapewniać:

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

Oryginalna odpowiedź

Krótko mówiąc, takie podejście z QGIS, o ile mi wiadomo, nie jest możliwe wyświetlenie maksymalnej liczby pikseli rastrowych. Tylko liczba, średnia i suma. Udostępniłem niektóre zasoby, które to obejmują.

Oto strona dokumentacji wtyczki Zonal Statistics Plugin

Dzięki wtyczce icon_zonal_statistics Zonal Statistics możesz analizować wyniki klasyfikacji tematycznej. Pozwala obliczyć kilka wartości pikseli warstwy rastrowej za pomocą wielobocznej warstwy wektorowej (patrz rysunek_zonal_statistics). Możesz obliczyć sumę, średnią wartość i całkowitą liczbę pikseli w obrębie wielokąta. Wtyczka generuje kolumny wyjściowe w warstwie wektorowej ze zdefiniowanym przez użytkownika prefiksem.

Ponadto w tym samouczku QGIS Próbkowanie danych rastrowych za pomocą punktów lub wielokątów podano przykładowe dane pozwalające uzyskać zwrot statystyk strefowych. Dołączyłem zrzuty ekranu, które podają jako przykład ich przepływu pracy i wyników.

Włącz wtyczki statystyk strefowych. Jest to podstawowa wtyczka, więc jest już zainstalowana. Zobacz Korzystanie z wtyczek, aby dowiedzieć się, jak włączyć podstawową wtyczkę

wprowadź opis zdjęcia tutaj

wprowadź opis zdjęcia tutaj

wprowadź opis zdjęcia tutaj

Po zakończeniu przetwarzania wybierz warstwę powiatów. Użyj narzędzia Identyfikuj i kliknij dowolny wielokąt hrabstwa. Zobaczysz trzy nowe atrybuty dodane do warstwy: ZS_count, ZS_mean i ZS_sum. Te atrybuty zawierają odpowiednio liczbę pikseli rastrowych, średnią wartości pikseli rastrowych i sumę wartości pikseli rastrowych. Ponieważ interesuje nas średnia temperatura, skorzystamy z pola ZS_mean.

wprowadź opis zdjęcia tutaj

whyzar
źródło
2

Od QGIS 3.0 (być może wcześniej, nie jestem pewien) dostępne są następujące statystyki:

  • Liczyć
  • Suma
  • Oznaczać
  • Mediana
  • Odchylenie standardowe
  • Minimum
  • Maksymalny
  • Zasięg
  • Mniejszość
  • Większość (tryb)
  • Różnorodność
  • Zmienność
dbaston
źródło