спектральные характеристики

[Anton]

Прошу прощения за беспокойство.
Помогите, если можете, разобраться с таким вопросом.
Во многих публикациях по дистанционному анализу приводятся графики т.н. спектральных характеристик объектов. Но нигде не сказано в какой программе их делают! Авторам, наверное, кажется, что это настолько элементарно, что об этом и не стоит говорить. Но я никак не могу понять, как, например, в эрдасе сделать графики спектральных характеристик, допустим, полигона сплошного покрытия фукусов на литорали? Чтобы по оси Х отображались нанометры, а по оси Y коэффициент отражения или обилие? Такие графики постоянно мелькают в публикациях, может быть Вы знаете алгоритм их построения по многозональным космоснимкам?
Спасибо.

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

http://gis-lab.info/qa/zonalstats.html
или
http://gis-lab.info/qa/extract-values-grids.html

в зависимости от того сразу на всю площадь нужен график или по точкам/пикселям.
А вообще это просто график, нанометры берутся из обозначений каналов тех данных которые используются.

[gamm]

http://gis-lab.info/qa/zonalstats.html
или
http://gis-lab.info/qa/extract-values-grids.html

в зависимости от того сразу на всю площадь нужен график или по точкам/пикселям.
А вообще это просто график, нанометры берутся из обозначений каналов тех данных которые используются.

ответ неправильный - по горизонтали получатся никак не нанометры

для получения спектрального профиля в нанометрах нужна модель (кривая) отклика каждого сенсора (канала), пересчет значений в reflectance + атмосферная коррекция, и синтез спектрального профиля. Сделать это технически несложно на любом языке программирования, если есть данные по сенсорам и программа для атмосферной коррекции (коэффициенты пересчета из попугаев DN в "натуральные показатели" лежат в метаданных снимка). Наверняка во всяких ердасах это есть.

В силу небольшого числа каналов и их перекрытия, такой спектральный профиль будет весьма приближенным. Для построения нормального спектрального профили используются гиперспектральные системы с сотнями узких каналов, именно такие картинки обычно и приводятся в статьях.

Эталонные спектральные портреты "фикусов на литорали" получают на земле, с использованием спектрометров.

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

хм, почему не нанометры (диапазон соотв. каналу), я не понял:



мне кажется глубина вопроса автора была на два порядка ниже ответа, но это пусть он сам решает. Мне показалось ему просто нужен поканальный профиль по DN.

[alexandr cherepanov]

ответ неправильный - по горизонтали получатся никак не нанометры

по оси Х может быть что угодно в том числе и длина волны в нанометрах

для получения спектрального профиля в нанометрах нужна модель (кривая) отклика каждого сенсора (канала)

не нужна, все данные по Solar Irradiance давно доступны

В силу небольшого числа каналов и их перекрытия

И опять не угадали, у современных систем перекрытие почти отсутствует (обычно только синий с зеленым), по крайне мере в области FWHM и на уровне 33% и ниже, а вот у гипперспектральных оно как раз есть.

такой спектральный профиль будет весьма приближенным

Любой спектральный профиль будет приближенным, даже самый детальный.

Эталонные спектральные портреты "фикусов на литорали" получают на земле, с использованием спектрометров.

на дерево сами полезете или с крана будете спектрометрировать?

PS
Максим извини, не удержался, впредь обещаю не троллить

[gamm]

ответ неправильный - по горизонтали получатся никак не нанометры

по оси Х может быть что угодно в том числе и длина волны в нанометрах

видимо, я неверно понял те инструменты, которые привел Максим - мне показалось, что по оси X они дают либо значения (т.е. получим гистограмму) , либо порядковый номер точки в профиле. Тогда был не прав. Остальное - из опыта работы со спектрометром с 2700 каналами

Кстати, бегать по лесу не нужно - самолетные гиперспектральные системы никто не отменял. А листочки у нас в лабораторию носили.

[gamm]

для получения спектрального профиля в нанометрах нужна модель (кривая) отклика каждого сенсора (канала)

не нужна, все данные по Solar Irradiance давно доступны

речь была о модели отклика в каналах (spectral curves), см. http://landsat.usgs.gov/tools_spectralViewer.php
Solar Irradiance ни причем

[alexandr cherepanov]

Посмотрите для чего нужны Spectral Response Curves и как раз увидите то, что они используются для определения Solar Spectral Irradiance в каждом из каналов. Именно об этом и идет речь

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

alexandr cherepanov, по-моему все конструктивно и у тебя и у gamm, информация полезная.

Остается только выяснить, что именно хотел узнать автор исходного вопроса.

[SS_Rebelious]

У меня, видимо, такая же проблема, как у автора поста. Хочу на основе снимка WorldView-2 построить график для определённых объектов в котором ось Y - коэффициент спектральной яркости, ось - X - длина волны в нанометрах.

Как это осуществить в QGIS, Optics или в другом свободно ПО?

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

КСЯ и DN это не обязательно одно и тоже, но на графике показывать можно и DN, подписи с длинами волн это эстетика.

удобно смотреть профили в кугисе в расширении Spectral Profile (без эстетики)

[SS_Rebelious]

КСЯ и DN это не обязательно одно и тоже, но на графике показывать можно и DN, подписи с длинами волн это эстетика.

К сожалению, мне нужен именно КСЯ, так как у меня есть КСЯ для похожего типа объектов и их надо бы сравнить. Мне собственно не до эстетики: меня интересует не инструмент, который построит мне график, а данные для построения - эти самые КСЯ для различных длин волн. Как их рассчитать для WorldView2?

удобно смотреть профили в кугисе в расширении Spectral Profile (без эстетики)

Я думал о Spectral Profile, но там же указано, что он работает с 8-битными файлами, а WorldView2 - 16-битный...

[syroezhka]

WorldView2 - 16-битный...

Это по паспорту, а в реале в 16ти битном файле самое большое значение 2047. т.е. он 11-битный.

К сожалению, мне нужен именно КСЯ, так как у меня есть КСЯ для похожего типа объектов и их надо бы сравнить.

Что бы сравнивать значения, с уже имеющимися, надо их привести к одному стандарту. Например альбедо или коефициенту отражения.
Для ВВ2 есть описание как перевести в эти коефициенты. Сейчас я попробую откопать

[syroezhka]

Такс...
Перевод из DN в Radiance
http://www.digitalglobe.com/downloads/p ... use_v1.pdf

Перевод из Radiance в альбедо или в коефициент отражения с. 293
http://biblio.ugent.be/input/download?f ... OId=469003

Вообще-то проблема хорошого спектрального визуализатора и для меня актуальна, т.к. сама искала недавно. Закончилось тем, что вытгивала аркгисом значения из снимка, а потом скармливала их матлабу. И там уже строила профиль. Правда не интерактивно это получилось. Матлаб, конечно не относится к бесплатным программам, но можно в R попробовать или экселе. Было дело пользовались последним для этих целей.
А че! - Дешево и сердито
Вообще я это все мутила, для того что бы сравнить данные ВВ2 с эталонами с спектральной библиотеки гислаба.

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

spectral profie (тот что мы правили) работает с 16бит и достаточно производителен, для себя делал