GRASS, проверка качества спектральных эталонов
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
GRASS, проверка качества спектральных эталонов
Есть ли в GRASS возможность визуально или количественно оценить качество спектральных эталонов для классификации с обучением?
Думаю, задачу не сложно решить в R, но м.б. этот алгоритм уже реализован в самом GRASS?
Алгоритм для визуальной оценки качества эталонов мне видится таковым:
1. Собрать статистику по растру и построить график плотности спектральных яркостей по каким-либо 2-м каналам, т.ч. диапазоны спектральных яркостей которые втречаются чаще, отображались ярче.
2. Собрать спектральную статистику по эталонам и построить границы эталонов разными цветами в пространстве спектральных яркостей тех же 2-х каналов.
3. Совместить два графика.
На полученной диаграмме будет видно пересекаются ли эталоны и достаточно ли равномерно все спектральное пространство покрыто эталонами.
Думаю, задачу не сложно решить в R, но м.б. этот алгоритм уже реализован в самом GRASS?
Алгоритм для визуальной оценки качества эталонов мне видится таковым:
1. Собрать статистику по растру и построить график плотности спектральных яркостей по каким-либо 2-м каналам, т.ч. диапазоны спектральных яркостей которые втречаются чаще, отображались ярче.
2. Собрать спектральную статистику по эталонам и построить границы эталонов разными цветами в пространстве спектральных яркостей тех же 2-х каналов.
3. Совместить два графика.
На полученной диаграмме будет видно пересекаются ли эталоны и достаточно ли равномерно все спектральное пространство покрыто эталонами.
-
- Гуру
- Сообщения: 5173
- Зарегистрирован: 26 сен 2009, 16:26
- Репутация: 793
- Ваше звание: званий не имею
- Откуда: Москва
Re: GRASS, проверка качества спектральных эталонов
Не видел ничего подобного. Статистика по растру считается в <r.statistics> и <r.stats>, далее строим графики в R или gnuplot и т.д.
Редактор материалов, модератор форума
-
- Гуру
- Сообщения: 810
- Зарегистрирован: 22 авг 2007, 14:58
- Репутация: 123
- Откуда: Казань
Re: GRASS, проверка качества спектральных эталонов
Откровенного говоря, существуют алгоритмы "более лучшие" -- эта тема довольно неплоха проработана, и для оценки качества эталонов разработаны специальные методы и метрики, которые учитывают не два канала, а их произвольный набор. Лично мне нравится книга Дэвида Ландгребе, в которой все это довольно подробно разжевано.bolotoved писал(а):Есть ли в GRASS возможность визуально или количественно оценить качество спектральных эталонов для классификации с обучением?
...
Алгоритм для визуальной оценки качества эталонов мне видится таковым
...
Что касается скрипта, то пару лет назад я писал модуль для рассчета расстояния Джефриса-Мацуситы, которое используется для оценки качества выбранного набора признаков.
-
- Гуру
- Сообщения: 4067
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1062
- Ваше звание: программист
- Откуда: Казань
Re: GRASS, проверка качества спектральных эталонов
там в формуле, похоже, опечатка - в правой части, внутри ln(), не хватает нормы для полусуммы матриц классов ... но вообще-то статистику лучше делать в специально предназначенном для этого месте, в R, тем более, что gdal там подключается.KolesovDmitry писал(а):Что касается скрипта, то пару лет назад я писал модуль для рассчета расстояния Джефриса-Мацуситы,
-
- Гуру
- Сообщения: 810
- Зарегистрирован: 22 авг 2007, 14:58
- Репутация: 123
- Откуда: Казань
Re: GRASS, проверка качества спектральных эталонов
Да там вообще кошмар творится -- сейчас посмотрел свежим взглядом, за голову схватился: там помимо всего прочего еще и корней не хватает внутри формул. И если норма для полусуммы хотя бы в коде модуля была, то корней нет нигде.gamm писал(а):там в формуле, похоже, опечатка - в правой части, внутри ln(), не хватает нормы для полусуммы матриц классов
Вечером переделаю - перезалью новую версию модуля.
Так получилось после того, как в очередной раз мне не хватило памяти прочитать из R снимок Landsat, я на R обиделся и создал свой велосипед (как сейчас выясняется, с квадратными колесами).но вообще-то статистику лучше делать в специально предназначенном для этого месте, в R, тем более, что gdal там подключается.
А вообще, конечно, для статистики лучше R. Тем более, что судя по новостям, в 3-м R проблему с памятью решили.
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: GRASS, проверка качества спектральных эталонов
Не сразу, но нашлось: http://grass.osgeo.org/grass70/manuals/ ... scatt.htmlАлександр Мурый писал(а):Не видел ничего подобного.
Выглядит сыровато, поддерживает работу только с 8-битными снимками, т.е. очевидно, не будет работать с атмосферно скорректированными. Зато есть некоторые плюшки, типа подсветка пикселов которые будут отнесены к данному классу при заданной сигнатуре, скатерплоты и др.
-
- Гуру
- Сообщения: 920
- Зарегистрирован: 30 дек 2008, 14:11
- Репутация: 236
- Откуда: Ханты-Мансийск
- Контактная информация:
Re: GRASS, проверка качества спектральных эталонов
Судя по оглавлению, книжка интересная, спасибо. Но пока что мне нужно более простая штука, выяснилось, что в 7-м GRASSe ее почти реализовали:KolesovDmitry писал(а):Откровенного говоря, существуют алгоритмы "более лучшие"
http://grasswiki.osgeo.org/wiki/WxIClass
http://grass.osgeo.org/grass70/manuals/ ... scatt.html
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей