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

Обсуждение материалов сайта: вопросы, замечания, предложения
Ответить
Аватара пользователя
SS_Rebelious
Гуру
Сообщения: 1304
Зарегистрирован: 24 фев 2009, 16:51
Статьи: 10
Проекты: 3/1
Репутация: 97
Ваше звание: GIS pro-fan
Откуда: Lahti / Газ-ПУТИНбург
Контактная информация:

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

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

Очередная статья.
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!

ericsson
Гуру
Сообщения: 3166
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 672
Ваше звание: Вредитель полей

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

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

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

ericsson
Гуру
Сообщения: 3166
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 672
Ваше звание: Вредитель полей

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

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

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

Аватара пользователя
Denis Rykov
Гуру
Сообщения: 3315
Зарегистрирован: 11 апр 2008, 21:09
Статьи: 33
Проекты: 9
Репутация: 492
Ваше звание: Author
Контактная информация:

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

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

А не пробовали реализовать алгоритм расчета из OTB на R? Интересно было бы сравнить результаты.
Spatial is now, more than ever, just another column- The Geometry Column.

alexandr cherepanov
Завсегдатай
Сообщения: 492
Зарегистрирован: 30 ноя 2006, 13:31
Статьи: 3
Проекты: 1
Репутация: 72
Откуда: Moscow

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

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

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

Аватара пользователя
nadiopt
Гуру
Сообщения: 4606
Зарегистрирован: 29 янв 2009, 13:27
Статьи: 4
Репутация: 434
Ваше звание: хрюкало
Откуда: Москва-Питер
Контактная информация:

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

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

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

Аватара пользователя
SS_Rebelious
Гуру
Сообщения: 1304
Зарегистрирован: 24 фев 2009, 16:51
Статьи: 10
Проекты: 3/1
Репутация: 97
Ваше звание: GIS pro-fan
Откуда: Lahti / Газ-ПУТИНбург
Контактная информация:

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

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

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

Добавил оригинальные гистограммы. Те гистограммы которые вы пытаетесь строить по продемонстрированным картинкам - бесполезны, так как они - всего-лишь скриншоты оригинальных изображений (к тому же подвергнутых операциям по улучшению визуализации).
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
Гуру
Сообщения: 1304
Зарегистрирован: 24 фев 2009, 16:51
Статьи: 10
Проекты: 3/1
Репутация: 97
Ваше звание: GIS pro-fan
Откуда: Lahti / Газ-ПУТИНбург
Контактная информация:

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

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

nadiopt писал(а):а мне вот что-то показалось, что картинка после орфео получше для визуального дешифрирования. они не перепутаны?
Картинки не перепутаны. Что касается визуального восприятия, то тут всё индивидуально.
alexandr cherepanov писал(а):А для чего нужна реализация именно через R? Тем более такого примитивного алгоритма не нуждающегося в функционале R.
Чтобы показать, как это делается в R. Статья адресована в первую очередь новичкам в R.
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!

ericsson
Гуру
Сообщения: 3166
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 672
Ваше звание: Вредитель полей

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

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

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

Аватара пользователя
SS_Rebelious
Гуру
Сообщения: 1304
Зарегистрирован: 24 фев 2009, 16:51
Статьи: 10
Проекты: 3/1
Репутация: 97
Ваше звание: GIS pro-fan
Откуда: Lahti / Газ-ПУТИНбург
Контактная информация:

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

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

ericsson писал(а):Э... то есть вы для иллюстрации результата приводите то, что еще непонятно каким преобразованиям подвергалось? Хмм...
Разъясняю. Спутниковые снимки поставляются не в виде RGB-композитов (хотя, и такие опции иногда присутствуют), а в виде наборов каналов, где значения пикселей соответствуют отражательной способности объектов местности в данном диапазоне длин волн. Любой растр при визуализации подвергается тем или иным преобразованиям которые контролируются настройками ПО. В статье паншарпенингу подвергался восьмиканальный растр, то есть я мог использовать для иллюстрации резултата либо 1 из каналов (для оттенков серого), либо 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!

ericsson
Гуру
Сообщения: 3166
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 672
Ваше звание: Вредитель полей

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

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

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

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

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

Аватара пользователя
SS_Rebelious
Гуру
Сообщения: 1304
Зарегистрирован: 24 фев 2009, 16:51
Статьи: 10
Проекты: 3/1
Репутация: 97
Ваше звание: GIS pro-fan
Откуда: Lahti / Газ-ПУТИНбург
Контактная информация:

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

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

ericsson писал(а):исходный композит по некой причине измучен настройками до того, что в нем начинает быть заметен клиппинг (не совсем верно сделанная компрессия диапазона)
Не понимаю, на каком основании вы сделали подобные выводы, разъясните поподробнее, пожалуйста. Сцена мной относительно значения пикселей никак не модифицировалась и не "настраивалась".
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!

ericsson
Гуру
Сообщения: 3166
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 672
Ваше звание: Вредитель полей

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

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

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

Ответить

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