GIS-LAB

Географические информационные системы и дистанционное зондирование


Форумы GIS-Lab.info

Геоинформационные системы (ГИС) и Дистанционное зондирование Земли


Паншарпенинг при помощи R

Обсуждение материалов сайта: вопросы, замечания, предложения

Паншарпенинг при помощи R

Сообщение SS_Rebelious » 19 апр 2015, 03:43

Очередная статья.
  • 3

Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!
Аватара пользователя
SS_Rebelious
GIS pro-fan
 
Зарегистрирован: 24 фев 2009
Откуда: Lahti / Газ-ПУТИНбург (на карте)
 
Сообщения: 1304
Репутация: 96
Статьи: 10
Проекты: 3/1

Re: Паншарпенинг при помощи R

Сообщение ericsson » 19 апр 2015, 04:59

Добавьте для наглядности
Orfeo.png
гистограммы результата из Orfeo Toolbox
Orfeo.png (6.48 КБ) Просмотров: 10431

R.png
гистограммы результата из R
R.png (6.48 КБ) Просмотров: 10431
  • 0

ericsson
Вредитель полей
 
Зарегистрирован: 27 июл 2009
 
Сообщения: 3103
Репутация: 640

Re: Паншарпенинг при помощи R

Сообщение ericsson » 19 апр 2015, 05:05

original_multispectral.png
оригинальное мультиспектральное изображение
original_multispectral.png (6.53 КБ) Просмотров: 10427

А вот это гистограмма исходника. Довольно плохо выглядит - в ней уже адский клиппинг в светах (это такое поставщик отдает или это все же последствие издевательств над исходником вроде неправильного перевода из 16 бит?) и ясно видно, что кто-то растягивал зеленый канал.
  • 0

ericsson
Вредитель полей
 
Зарегистрирован: 27 июл 2009
 
Сообщения: 3103
Репутация: 640

Re: Паншарпенинг при помощи R

Сообщение Denis Rykov » 19 апр 2015, 08:25

А не пробовали реализовать алгоритм расчета из OTB на R? Интересно было бы сравнить результаты.
  • 0

Spatial is now, more than ever, just another column- The Geometry Column.
Аватара пользователя
Denis Rykov
Author
 
Зарегистрирован: 11 апр 2008
 
Сообщения: 3278
Репутация: 482
Статьи: 33
Проекты: 9

Re: Паншарпенинг при помощи R

Сообщение alexandr cherepanov » 19 апр 2015, 11:38

А для чего нужна реализация именно через R? Тем более такого примитивного алгоритма не нуждающегося в функционале R. Давайте еще напишем статьи про реализацию fusion c использованием pilow, imagemagick, opencv и т.д.
  • 1

alexandr cherepanov
Завсегдатай
 
Зарегистрирован: 30 ноя 2006
Откуда: Moscow (на карте)
 
Сообщения: 481
Репутация: 64
Статьи: 3
Проекты: 1

Re: Паншарпенинг при помощи R

Сообщение nadiopt » 19 апр 2015, 11:59

ну мало ли, может надо срочно что-то пошарпить, а под рукой только R
а мне вот что-то показалось, что картинка после орфео получше для визуального дешифрирования. они не перепутаны?
  • 0

ин гроссен фамилен нихт клювен клац клац
Аватара пользователя
nadiopt
хрюкало
 
Зарегистрирован: 29 янв 2009
Откуда: Москва-Питер (на карте)
 
Сообщения: 4560
Репутация: 426
Статьи: 4

Re: Паншарпенинг при помощи R

Сообщение SS_Rebelious » 19 апр 2015, 12:08

Denis Rykov писал(а):А не пробовали реализовать алгоритм расчета из OTB на R? Интересно было бы сравнить результаты.

В документации к OTB я не нашёл конкретного описания используемого алгоритма.

Добавил оригинальные гистограммы. Те гистограммы которые вы пытаетесь строить по продемонстрированным картинкам - бесполезны, так как они - всего-лишь скриншоты оригинальных изображений (к тому же подвергнутых операциям по улучшению визуализации).
  • 0

Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!
Аватара пользователя
SS_Rebelious
GIS pro-fan
 
Зарегистрирован: 24 фев 2009
Откуда: Lahti / Газ-ПУТИНбург (на карте)
 
Сообщения: 1304
Репутация: 96
Статьи: 10
Проекты: 3/1

Re: Паншарпенинг при помощи R

Сообщение SS_Rebelious » 19 апр 2015, 12:15

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

Картинки не перепутаны. Что касается визуального восприятия, то тут всё индивидуально.

alexandr cherepanov писал(а):А для чего нужна реализация именно через R? Тем более такого примитивного алгоритма не нуждающегося в функционале R.

Чтобы показать, как это делается в R. Статья адресована в первую очередь новичкам в R.
  • 0

Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!
Аватара пользователя
SS_Rebelious
GIS pro-fan
 
Зарегистрирован: 24 фев 2009
Откуда: Lahti / Газ-ПУТИНбург (на карте)
 
Сообщения: 1304
Репутация: 96
Статьи: 10
Проекты: 3/1

Re: Паншарпенинг при помощи R

Сообщение ericsson » 19 апр 2015, 22:32

SS_Rebelious писал(а):Добавил оригинальные гистограммы. Те гистограммы которые вы пытаетесь строить по продемонстрированным картинкам - бесполезны, так как они - всего-лишь скриншоты оригинальных изображений (к тому же подвергнутых операциям по улучшению визуализации).

Э... то есть вы для иллюстрации результата приводите то, что еще непонятно каким преобразованиям подвергалось? Хмм...
  • 0

ericsson
Вредитель полей
 
Зарегистрирован: 27 июл 2009
 
Сообщения: 3103
Репутация: 640

Re: Паншарпенинг при помощи R

Сообщение SS_Rebelious » 20 апр 2015, 00:44

ericsson писал(а):Э... то есть вы для иллюстрации результата приводите то, что еще непонятно каким преобразованиям подвергалось? Хмм...


Разъясняю. Спутниковые снимки поставляются не в виде RGB-композитов (хотя, и такие опции иногда присутствуют), а в виде наборов каналов, где значения пикселей соответствуют отражательной способности объектов местности в данном диапазоне длин волн. Любой растр при визуализации подвергается тем или иным преобразованиям которые контролируются настройками ПО. В статье паншарпенингу подвергался восьмиканальный растр, то есть я мог использовать для иллюстрации резултата либо 1 из каналов (для оттенков серого), либо 3 (для цветного). Естественно, был выбран второй вариант - ибо а) это красивше, и б) видны артефакты паншарпенинга, связанные с искажением цветов (например, зелёный цвет может становится серым или коричневым). Для каждого восьмиканального растра, используемого в качестве иллюстрации использовались идентичные комбинации каналов и настройки отображения. Таким образом у нас есть возможность сравнивать результаты работы алгоритмов, а также растра до и после преобразований.
  • 0

Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!
Аватара пользователя
SS_Rebelious
GIS pro-fan
 
Зарегистрирован: 24 фев 2009
Откуда: Lahti / Газ-ПУТИНбург (на карте)
 
Сообщения: 1304
Репутация: 96
Статьи: 10
Проекты: 3/1

Re: Паншарпенинг при помощи R

Сообщение ericsson » 20 апр 2015, 03:08

Я прекрасно знаю, в форме чего поставляются спутниковые снимки и что с ними делают, так что приведенный ликбез - лишнее.

Наличие в приведенных вами гистограммах каналов, которые отсутствуют в композите, который служит для иллюстрации, не добавляет им наглядности. Я не против RGB-композита, он имеет смысл, как иллюстрация, безусловно.

Просто то, что и исходный композит по некой причине измучен настройками до того, что в нем начинает быть заметен клиппинг (не совсем верно сделанная компрессия диапазона), и результирующие композиты, подвергшиеся pansharpening-у, грешат ровно тем же, но (судя по анализу) в еще большей степени - все это снижает иллюстративную ценность картинок и вызывает ложное впечатление, что pansharpening улучшает разрешение, но портит все остальное.
  • 0

ericsson
Вредитель полей
 
Зарегистрирован: 27 июл 2009
 
Сообщения: 3103
Репутация: 640

Re: Паншарпенинг при помощи R

Сообщение SS_Rebelious » 20 апр 2015, 11:41

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


Не понимаю, на каком основании вы сделали подобные выводы, разъясните поподробнее, пожалуйста. Сцена мной относительно значения пикселей никак не модифицировалась и не "настраивалась".
  • 0

Look for something long enough, and you will find it. Look for something without understanding, and it will find you...
"All paid jobs absorb and degrade the mind." Aristotle
If you take 1 step towards freedom it'll take 2 steps towards you!
Аватара пользователя
SS_Rebelious
GIS pro-fan
 
Зарегистрирован: 24 фев 2009
Откуда: Lahti / Газ-ПУТИНбург (на карте)
 
Сообщения: 1304
Репутация: 96
Статьи: 10
Проекты: 3/1

Re: Паншарпенинг при помощи R

Сообщение ericsson » 21 апр 2015, 04:27

Так, у вас в статье есть картинка, подписанная "Мультиспектральный растр с пространственным разрешением 2,0 м" - это исходный RGB-композит низкого разрешения.
У него не все хорошо с уровнем черного (первое значение появляется на уровне 57 для красного, 63 для зеленого и 103 для синего).
А в светах картина обратная: в изображении для красного канала 1227 пикселей имеют значение 255, и только 83 - значение 254, 53 - значение 253 (далее в сторону теней идет "плато" на среднем уровне 25) - если это не клиппинг, то я уже не знаю, что можно назвать клиппингом. Картина с зеленым каналом всего в два раза лучше (там 466 пикселей с макс. насыщенностью), с синим - в четыре (там их 317).
  • 0

ericsson
Вредитель полей
 
Зарегистрирован: 27 июл 2009
 
Сообщения: 3103
Репутация: 640


Вернуться в Материалы сайта

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 1


(Геокруг)

© GIS-Lab и авторы, 2002-2013. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).