##Raster processing=group           # задаём имя группы скриптов SEXTANTE, где будет находиться скрипт 
##input=raster                      # задаём поле выбора растра
##round_values_to_ndigits=string    # задаём поле ввода значения для округления
##output_file=output html           # задаём вывод результатов в виде html-файла (путь к файлу)
 
# импортируем необходимые библиотеки
from osgeo import gdal
import sys
import math
 
# загружаем растр
gdalData = gdal.Open(str(input))
 
# получаем размеры растра в пикселях
xsize = gdalData.RasterXSize
ysize = gdalData.RasterYSize
 
# получаем количество слоёв растра
bands = gdalData.RasterCount
 
# начинаем запись результирующего html-файла
f = open(output_file, 'a')
f.write('<TABLE>\n<TH>Band Number</TH> <TH>Cell Value</TH> <TH>Count</TH>\n')
 
# подсчёт пикселей и запись результата
for i in xrange(1, bands + 1):
  band_i = gdalData.GetRasterBand(i)
  raster = band_i.ReadAsArray()
 
  # создаём словарь в который будем записывать результат
  count = {}
 
  # подсчитываем уникальные значения пикселей для данного слоя
  for col in range(xsize):
    for row in range(ysize):
      cell_value = raster[row, col]
 
      # проверяем, не соответствует ли cell_value тип NaN
      if math.isnan(cell_value):
        cell_value = 'Null' # если да, то считаем этот пиксель как 'Null' 
 
      # округляем значения пикселя если это необходимо
      elif round_values_to_ndigits:
        try:
          cell_value = round(cell_value, int(round_values_to_ndigits))
        except:
          # если в поле 'round_values_to_ndigits' было введено нечто непонятное, 
          # будем считать, что пользователь хотел округлить до целых чисел
          cell_value = round(cell_value)
 
      # наконец, добавляем / обновляем счётчик для данного значения пикселя
      try:
        count[cell_value] += 1
      except:
        count[cell_value] = 1
 
  # записываем отсортированные результаты обработки данного слоя растра в таблицу
  for key in sorted(count.iterkeys()):
    line = "<TD>%s</TD> <TD>%s</TD> <TD>%s</TD>" %(i, key, count[key])
    f.write('<TR>'+ line + '</TR>' + '\n')
 
# заканчиваем запись таблицы и закрываем файл  
f.write('</TABLE>')
f.close