Automatycznie stylizować raster przy użyciu unikalnych wartości w QGIS?

W ArcMap możliwe jest automatyczne obliczenie unikalnych wartości rastra, a następnie zastosowanie innego stylu do każdej unikalnej wartości (patrz 1. zdjęcie).

Jednak w QGIS muszę ręcznie dodawać wartości do stylizacji podczas korzystania ze stylu „Jednopasmowy pseudokolor” (patrz drugi obraz). Czy istnieje sposób automatycznego zapełniania unikalnych wartości, tak jak ArcMap?

ArcMap:

QGIS (jak automatycznie wypełnić listę wartości unikalnymi wartościami?):

Dla tych, którzy wciąż tego szukają. Unikalne wartości rastrowe zostały dodane do QGIS 3. ”Dodane przez Nyall Dawson około rok temu

[FUNKCJA] Zezwól na klasyfikowanie renderowanego palety przy użyciu unikalnych wartości z warstwy rastrowej

Dodaje łatwy sposób stylizowania dyskretnych rastrów, takich jak klasy użytkowania terenu, za pomocą renderera Paletted. Wystarczy wybrać renderer z paletą, wybrać zespół, a następnie nacisnąć przycisk „Dodaj unikalne wartości”. Unikalne wartości pikseli zostaną pobrane z warstwy, a kolor zostanie przypisany do każdej z nich przy użyciu aktualnie wybranej rampy kolorów. ”

Oto podejście, które może być pomocne. Narzędzie GRASS r.reportjest w stanie obliczyć podstawową statystykę dla każdej unikalnej wartości w warstwie rastrowej:

Dane wyjściowe można zapisać jako zwykły plik tekstowy wyglądający tak (po lewej: wartości unikalne, po prawej: obszar objęty każdą wartością):

+-----------------------------------------------------------------------------+
|                      Category Information                        |    square|
|   #|description                                                  |     miles|
|-----------------------------------------------------------------------------|
|1111| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  1.249845|
|1113| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  0.062666|
[...]
|3412| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  0.013926|
|4111| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  4.713902|
|4211| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  0.083555|
|4212| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . |  0.135777|
|-----------------------------------------------------------------------------|
|TOTAL                                                             | 17.048727|
+-----------------------------------------------------------------------------+

Poniższy skrypt geoprzetwarzania pobiera ten plik tekstowy jako dane wejściowe i konfiguruje odpowiedni renderer (dla wygody dla losowych kolorów) dla warstwy rastrowej (zainspirowany przez LAWHEAD, J .: QGIS Python Programming Cookbook, s. 135):

##giswg=group
##thexml=file
##theraster=raster
from random import randint
from PyQt4.QtGui import QColor
from qgis.core import *
myraster = processing.getObject(theraster)
infile = open(thexml, 'r')
univalues = [l.split('|')[1] for l in [f for f in infile.readlines()][4:-4]]
s = QgsRasterShader()
c = QgsColorRampShader()
c.setColorRampType(QgsColorRampShader.EXACT)
i = []
for u in univalues:
    i.append(QgsColorRampShader.ColorRampItem(float(u), QColor('#%06x' % randint(0, 2**24)), u))
c.setColorRampItemList(i)
s.setRasterShaderFunction(c)
ps = QgsSingleBandPseudoColorRenderer(myraster.dataProvider(), 1, s)
myraster.setRenderer(ps)

Dzięki legendarnym możliwościom przetwarzania QGIS można połączyć te dwa elementy w graficznym modelerze i dodać parametr rastrowy, aby po prostu wybrać załadowaną warstwę rastrową:

Po zapisaniu modelu, klikając go dwukrotnie w przyborniku przetwarzania, można go używać jak zwykłego narzędzia:

Z wynikiem wyglądającym tak (wychylenie zrasteryzowanego zestawu danych natura2000):

Lub we właściwościach warstwy:

Myślę, że wymaga to nieco udoskonalenia doświadczenia użytkownika, ale to pierwszy krok.