GIS-LAB

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


Форумы GIS-Lab.info

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


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

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

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

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

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

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

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

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

Zverik
Активный участник
 
Зарегистрирован: 11 сен 2011
Откуда: Москва (на карте)
 
Сообщения: 136
Репутация: 79
Статьи: 2

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

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

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

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

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

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

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

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

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

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

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

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

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

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

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

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

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

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

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

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

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

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

Хорошо бы аргументировать этот шаг.
После того, как мы сменили регион, нужно изменить диапазон значений у всех предполагаемых к использованию файлов, т.е. у стандартных каналов (например, для синтеза в естественных цветах, B2,B3,B4) и у B8 на единый диапазон [0;255].
  • 0

___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
Аватара пользователя
rhot
ВольныйИсследователь
 
Зарегистрирован: 25 янв 2011
Откуда: Архангельск (на карте)
 
Сообщения: 1521
Репутация: 112
Статьи: 1

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

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

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

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

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

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

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

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

Аватара пользователя
darsvid
Активный участник
 
Зарегистрирован: 29 июн 2012
Откуда: Kyїv, Ukraine (на карте)
 
Сообщения: 130
Репутация: 85
Статьи: 5
Проекты: 1

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

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

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

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

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

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

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

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

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

Редактор материалов, модератор форума
Александр Мурый
званий не имею
 
Зарегистрирован: 26 сен 2009
Откуда: Москва (на карте)
 
Сообщения: 4884
Репутация: 668
Статьи: 3
Проекты: 5/1

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

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

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

Аватара пользователя
Эдуард Казаков
Завсегдатай
 
Зарегистрирован: 23 апр 2014
Откуда: Санкт-Петербург (на карте)
 
Сообщения: 282
Репутация: 237
Статьи: 3

След.

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

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

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


(Геокруг)

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