Вот подход, который может быть полезным. Инструмент GRASS r.report
способен вычислять базовую статистику для каждого уникального значения в растровом слое:
Вывод может быть записан как простой текстовый файл, который выглядит следующим образом (слева: уникальные значения, справа: область, охватываемая каждым значением):
+-----------------------------------------------------------------------------+
| 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|
+-----------------------------------------------------------------------------+
Следующий скрипт геообработки принимает этот текстовый файл в качестве входных данных и устанавливает подходящий рендерер (со случайными цветами для удобства) для растрового слоя (Вдохновленный LAWHEAD, J .: QGIS Python Programming Cookbook, стр. 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)
Благодаря легендарным возможностям обработки QGIS вы можете объединить эти два в графическом редакторе моделей и добавить параметр растра, чтобы просто выбрать загруженный растровый слой:
После сохранения модели, дважды щелкнув ее в панели инструментов обработки, ее можно использовать как обычный инструмент:
Результат выглядит примерно так (выход растеризованного набора данных natura2000):
Или в свойствах слоя:
Я думаю, что это требует некоторого уточнения в пользовательском опыте, но это первый шаг.