в чем рассчитать матрицу ошибок и статистику каппа

[nadiopt]

Коллеги, нужна помощь
Есть результаты классификации снимков на одну и ту же территорию, но разными алгоритмами и в разных программах. Результаты можно представить как угодно (вектор, грид...). Классы, естественно, одинаковые
вопрос: в чем бы рассчитать матрицу ошибок и статистику каппа (меру согласованности) между всеми классификациями? Эрдас и Энви не едят результаты классификации из других программ:( в распоряжении есть еще Сканекс Имидж Процессор, Грасс с кугисом.
Напрашивается R, можно ли там решить такую задачу "в лоб"? нет ли инструкции, как это сделать?

[KolesovDmitry]

Есть результаты классификации снимков на одну и ту же территорию, но разными алгоритмами и в разных программах.<..> в чем бы рассчитать матрицу ошибок и статистику каппа (меру согласованности) между всеми классификациями?<..> в распоряжении есть еще Сканекс Имидж Процессор, Грасс с кугисом.

В GRASS есть специальный модуль r.kappa поэтому, если в у вас есть хотя бы базовое знакомство с GRASS, проблем возникнуть не должно.

[nadiopt]

совсем поверхностное, но разберусь. Спасибо!

[gamm]

Напрашивается R, можно ли там решить такую задачу "в лоб"?

там все можно

первое, что попало под руку - asbio::Kappa (функция в пакете asbio). Данные вытащить в Сканексе сэмплером по векторной маске, в текстовый файл. Там еще десятка два пакетов с каппой, но думаю разницы нет.

[Максим Дубинин]

вот тут еще
http://gis-lab.info/qa/geoproc-python-ag.html

[Александр Мурый]

Нашёл каппа-индекс в SEXTANTE (у меня в gvSIG OADE Edition 2010, но есть, наверное, и под Kosmo, и под OpenJUMP).

SEXTANTE_Kappa_index.png (22.33 КБ) 12647 просмотров

А вот в SAGA, как ни странно, не оказалось..

[nadiopt]

а для гвсиг растр надо в какой-нибудь хитрый формат преобразовывать, или грид сойдет?

[Александр Мурый]

Что есть "грид"? Аркгисовский "ASCII Grid", что ли? Мне кажется, gvSIG должен открывать все GDAL-совместимые форматы растров. Пробовать надо, короче

[nadiopt]

арковский, да. попробую, о результатах отпишусь

[nadiopt]

вот, добралась. Гвсиг (в секстанте) делает таблицу, которую кладу сюда, и некий растр. ВОПРОС: а где значения статистики-то самой?
да, то что там в одном гриде -99999, а в другом -9999, я заметила, это влияет?

[bolotoved]

Гвсиг (в секстанте) делает таблицу, которую кладу сюда, и некий растр.

Это что-то не то.
Растра по идее вообще не должно быть. Какова последовательность ваших действий?

На входе должен быть проверяемый классифицированный растр и проверочный растр (карта) или точки с известными значениями классов.
На выходе матрица ошибок по каждому классу и различные показатели точности, например, каппа.

[bolotoved]

Нашёл каппа-индекс в SEXTANTE (у меня в gvSIG OADE Edition 2010, но есть, наверное, и под Kosmo, и под OpenJUMP).

Судя по тому, что стоит в компании с kappa индексом - разные ландшафтные метрики из FragStat, то это совсем не то, что надо.

[bolotoved]

У меня есть инструкция как расчитать каппу в GRASS - делал для студней

[nadiopt]

спасибо
есть такой модуль в секстанте - кросс чек (каппа статистикс). скормила ему asc, сконвертированное из двух гридов? на выходе он дает дбф под названием кросс-классификайшн и растр

[Александр Мурый]

У меня есть инструкция как расчитать каппу в GRASS - делал для студней

Однозначно в статьи