Ecco un approccio che potrebbe essere utile. Lo strumento GRASS r.report
è in grado di calcolare una statistica di base per ciascun valore univoco in un livello raster:
L'output può essere scritto come un semplice file txt simile a questo (a sinistra: valori univoci, a destra: area coperta da ciascun valore):
+-----------------------------------------------------------------------------+
| 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|
+-----------------------------------------------------------------------------+
Il seguente script di geoprocessing prende questo file di testo come input e imposta un renderer adatto (con colori casuali per comodità) per il livello raster (Ispirato da LAWHEAD, J .: QGIS Python Programming Cookbook, p. 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)
Grazie alle leggendarie capacità di elaborazione di QGIS è possibile unire questi due elementi nel modellatore grafico e aggiungere un parametro raster per scegliere semplicemente un livello raster caricato:
Dopo aver salvato il modello, facendo doppio clic su di esso nella casella degli strumenti di elaborazione può essere utilizzato come un normale strumento:
Con un risultato simile al seguente (affioramento di un set di dati natura2000 rasterizzato):
O nelle proprietà del livello:
Penso che questo abbia bisogno di un certo perfezionamento nell'esperienza utente, ma è un primo passo.