рассчёт произвольной зональной статистики

[jerry-maori]

Ситуация следующая:
есть исходный растр и сеть полигонов
Задача:
посчитать внутри каждого полигона энтропию (соотношение пикселей разного цвета)
С ходу задача не решилась. Пришлось тащить матрицу в виде ASCII в Matlab и считать там.
А есть ли какие решения внтури самого ArcGIS ?

Заранее спасибо.

[lam]

посчитать внутри каждого полигона энтропию (соотношение пикселей разного цвета)

Уточните по какой формуле Вы определяете "соотношение пикселей разного цвета внутри каждого полигона".

[jerry-maori]

Шеннон, -sum(p.^log(p)), где p - отношение числа пикселей данного цвета к общему числу пикселей внутри полигона

[lam]

Шеннон, -sum(p.^log(p)), где p - отношение числа пикселей данного цвета к общему числу пикселей внутри полигона

Что такое p. и чем оно отличается от p? p=quantity pixel red/quantity pixel all?

[jerry-maori]

пусть в квадрте N пикселей разного цвета. тогда pi=ni/N доля пикселя цвета i.
в итоге получим вектор p
и считаем для него индекс Шеннона

[lam]

Вот эта тема должна Вам помочь - viewtopic.php?f=16&t=12306
Чтобы работать со значениями каждого пикселя в гриде, предлагаю перейти к узлам этого грида в виде точечного шейпа.