Обработка и интерпретация данных Landsat 8 в GRASS 7

Обсуждение материалов сайта: вопросы, замечания, предложения
Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 03 мар 2015, 12:42

Подробная статья, посвященная вопросам базовой обработки и интерпретации данных Landsat 8 (сенсор OLI) в среде GRASS GIS 7. Подготовлена вместе с моими студентами Натальей Волгушевой и Александром Семёновым.
Прошу вас высказать все замечания (по поводу орфографии, пунктуации и стилистики - в личные сообщения, остальное здесь, на всеобщее обсуждение) и пожелания.

===============
P.S. Статья опубликована.

Zverik
Активный участник
Сообщения: 141
Зарегистрирован: 11 сен 2011, 14:11
Статьи: 2
Репутация: 86
Откуда: Москва
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Zverik » 03 мар 2015, 19:10

Спасибо, очень большая работа. Непонятно, может ли она служить инструкцией для начинающих пользователей GRASS — сам я не пользуюсь, — но явно поможет ищущим подсказок.

Мне, как непрофессионалу, не хватает картинок «до и после» в главах про коррекции и преобразования. Ещё странно, что на такой картинке в разделе про паншарпенинг разные цвета. В остальных разделах есть картинки «после», но не видя исходного куска карты, сложно их оценить.

Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 03 мар 2015, 20:18

Здравствуйте, Илья!
В принципе, статью можно воспринимать как пошаговое руководство с ответвлениями, поэтому начинающим должно помочь :)

Что касается картинок. В случае радиометрической калибровки изображения "до и после" визуально никак не отличаются. Мы пересчитываем одни единицы (нормализованные) в другие (физические), при этом отображение самой картинки на экране не меняется, смысла их вставлять нет. Наверное, есть смысл добавить табличное сравнение значений одних и тех же пикселей в трех типах изображений: исходных DN, radiance и reflectance. Если ещё кто-нибудь поддержит эту идею, добавлю.

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

Ну и вообще внешний вид снимка, сами понимаете - покрутили яркость/контраст, потаскали гистограмму, и всё уже выглядит иначе.

Исходное изображение в естественных цветах точно нужно будет добавить, вечером сделаю.

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

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

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение ericsson » 03 мар 2015, 23:24

Как-то странно выглядит пара "до и после" для IHS - синий нормализовался, а вроде как не должен.

voltusya
Интересующийся
Сообщения: 23
Зарегистрирован: 04 ноя 2014, 18:40
Статьи: 1
Репутация: 3
Откуда: СПб

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение voltusya » 03 мар 2015, 23:39

Чуть выше говорилось о наличие флага Rebalance, пример приведен именно с его использованием.

Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 04 мар 2015, 00:32

В раздел Pan Sharpening, где впервые возникает осмысленная потребность визуально сравнивать что-то с исходным изображением, добавил полноразмерный "сырой" RGB-композит, с которым предполагается сопоставление полученных продуктов.

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

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

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

Про rebalance - понятно.
Еще маленькое замечание по скриншотам. Скриншоты интерфейса стоит всегда сохранять в PNG и никогда в JPEG, а то они "грязными" выглядят.

Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 04 мар 2015, 10:08

Да уж, с jpeg я уже осознал провал... Очевидные вещи ведь. Надеюсь в данном случае это не критично, очень уж много скриншотов.

Аватара пользователя
rhot
Гуру
Сообщения: 1521
Зарегистрирован: 25 янв 2011, 17:50
Статьи: 1
Репутация: 112
Ваше звание: ВольныйИсследователь
Откуда: Архангельск

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение rhot » 04 мар 2015, 11:47

Хорошо бы аргументировать этот шаг.
После того, как мы сменили регион, нужно изменить диапазон значений у всех предполагаемых к использованию файлов, т.е. у стандартных каналов (например, для синтеза в естественных цветах, B2,B3,B4) и у B8 на единый диапазон [0;255].
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

voltusya
Интересующийся
Сообщения: 23
Зарегистрирован: 04 ноя 2014, 18:40
Статьи: 1
Репутация: 3
Откуда: СПб

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение voltusya » 04 мар 2015, 12:33

В документации к модулю i.pansharpen сказано, что он работает только с 8-битными изображениями. Мы проверяли, действительно не работает без изменеия диапазона всех 4 каналов. Может это стоит отметить в статье, спасибо

Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 04 мар 2015, 16:42

Добавил комментарий про rescale у паншарпенинга

Аватара пользователя
darsvid
Активный участник
Сообщения: 135
Зарегистрирован: 29 июн 2012, 12:40
Статьи: 5
Проекты: 1
Репутация: 87
Откуда: Kyїv, Ukraine
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение darsvid » 04 мар 2015, 19:40

Спасибо, отличная статья :)

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

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

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

Понятно, что так статья может разрастись, но еще есть http://grass.osgeo.org/grass64/manuals/ ... t.rgb.html

Александр Мурый
Гуру
Сообщения: 4979
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 698
Ваше звание: званий не имею
Откуда: Москва

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Александр Мурый » 10 мар 2015, 21:05

Спасибо за отличную статью!

В статье сказано: "В примере будет использован набор данных LC81740292013131, охватывающий юго-восточное побережье Кубы, поэтому используется имя Cuba. "
Указанный район имеет номера Path=11 и Row=45-47. В статье Path=174 и Row=29. Название набора данных перепутано.

Мелких замечаний, в основном по GRASS, довольно много, но это чуть позже, отдельно.
Редактор материалов, модератор форума

Аватара пользователя
Эдуард Казаков
Завсегдатай
Сообщения: 307
Зарегистрирован: 23 апр 2014, 17:11
Статьи: 3
Репутация: 267
Откуда: Санкт-Петербург
Контактная информация:

Re: Обработка и интерпретация данных Landsat 8 в GRASS 7

Сообщение Эдуард Казаков » 11 мар 2015, 00:07

С именем сцены глупая ошибка, конечно, везде должно быть LC80110462013131, именно она использовалась. Исправил.

Ответить

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