Обработка и интерпретация данных Landsat 8 в GRASS 7
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Обработка и интерпретация данных Landsat 8 в GRASS 7
Подробная статья, посвященная вопросам базовой обработки и интерпретации данных Landsat 8 (сенсор OLI) в среде GRASS GIS 7. Подготовлена вместе с моими студентами Натальей Волгушевой и Александром Семёновым.
Прошу вас высказать все замечания (по поводу орфографии, пунктуации и стилистики - в личные сообщения, остальное здесь, на всеобщее обсуждение) и пожелания.
===============
P.S. Статья опубликована.
Прошу вас высказать все замечания (по поводу орфографии, пунктуации и стилистики - в личные сообщения, остальное здесь, на всеобщее обсуждение) и пожелания.
===============
P.S. Статья опубликована.
-
- Активный участник
- Сообщения: 158
- Зарегистрирован: 11 сен 2011, 14:11
- Репутация: 96
- Откуда: Москва
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Спасибо, очень большая работа. Непонятно, может ли она служить инструкцией для начинающих пользователей GRASS — сам я не пользуюсь, — но явно поможет ищущим подсказок.
Мне, как непрофессионалу, не хватает картинок «до и после» в главах про коррекции и преобразования. Ещё странно, что на такой картинке в разделе про паншарпенинг разные цвета. В остальных разделах есть картинки «после», но не видя исходного куска карты, сложно их оценить.
Мне, как непрофессионалу, не хватает картинок «до и после» в главах про коррекции и преобразования. Ещё странно, что на такой картинке в разделе про паншарпенинг разные цвета. В остальных разделах есть картинки «после», но не видя исходного куска карты, сложно их оценить.
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Здравствуйте, Илья!
В принципе, статью можно воспринимать как пошаговое руководство с ответвлениями, поэтому начинающим должно помочь
Что касается картинок. В случае радиометрической калибровки изображения "до и после" визуально никак не отличаются. Мы пересчитываем одни единицы (нормализованные) в другие (физические), при этом отображение самой картинки на экране не меняется, смысла их вставлять нет. Наверное, есть смысл добавить табличное сравнение значений одних и тех же пикселей в трех типах изображений: исходных DN, radiance и reflectance. Если ещё кто-нибудь поддержит эту идею, добавлю.
У атмосферной коррекции ситуация получше, хотя визуально картинка и не очень сильно меняется. Просто для этой процедуры принципиально не изменение изображения внешне (малозаметное, особенно на небольших скриншотах), а корректировка значений пикселей для их дальнейшей алгоритмической обработки.
Ну и вообще внешний вид снимка, сами понимаете - покрутили яркость/контраст, потаскали гистограмму, и всё уже выглядит иначе.
Исходное изображение в естественных цветах точно нужно будет добавить, вечером сделаю.
Внешний вид результатов паншарпенинга не настраивался, то есть что склеилось модулем, то и вставлялось. Причем разные методы паншарпенинга дают разные оттенки.
В принципе, статью можно воспринимать как пошаговое руководство с ответвлениями, поэтому начинающим должно помочь
Что касается картинок. В случае радиометрической калибровки изображения "до и после" визуально никак не отличаются. Мы пересчитываем одни единицы (нормализованные) в другие (физические), при этом отображение самой картинки на экране не меняется, смысла их вставлять нет. Наверное, есть смысл добавить табличное сравнение значений одних и тех же пикселей в трех типах изображений: исходных DN, radiance и reflectance. Если ещё кто-нибудь поддержит эту идею, добавлю.
У атмосферной коррекции ситуация получше, хотя визуально картинка и не очень сильно меняется. Просто для этой процедуры принципиально не изменение изображения внешне (малозаметное, особенно на небольших скриншотах), а корректировка значений пикселей для их дальнейшей алгоритмической обработки.
Ну и вообще внешний вид снимка, сами понимаете - покрутили яркость/контраст, потаскали гистограмму, и всё уже выглядит иначе.
Исходное изображение в естественных цветах точно нужно будет добавить, вечером сделаю.
Внешний вид результатов паншарпенинга не настраивался, то есть что склеилось модулем, то и вставлялось. Причем разные методы паншарпенинга дают разные оттенки.
-
- Гуру
- Сообщения: 3321
- Зарегистрирован: 27 июл 2009, 19:26
- Репутация: 748
- Ваше звание: Вредитель полей
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Как-то странно выглядит пара "до и после" для IHS - синий нормализовался, а вроде как не должен.
-
- Участник
- Сообщения: 78
- Зарегистрирован: 04 ноя 2014, 18:40
- Репутация: 48
- Откуда: Санкт-Петербург
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Чуть выше говорилось о наличие флага Rebalance, пример приведен именно с его использованием.
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
В раздел Pan Sharpening, где впервые возникает осмысленная потребность визуально сравнивать что-то с исходным изображением, добавил полноразмерный "сырой" RGB-композит, с которым предполагается сопоставление полученных продуктов.
-
- Гуру
- Сообщения: 3321
- Зарегистрирован: 27 июл 2009, 19:26
- Репутация: 748
- Ваше звание: Вредитель полей
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Про rebalance - понятно.
Еще маленькое замечание по скриншотам. Скриншоты интерфейса стоит всегда сохранять в PNG и никогда в JPEG, а то они "грязными" выглядят.
Еще маленькое замечание по скриншотам. Скриншоты интерфейса стоит всегда сохранять в PNG и никогда в JPEG, а то они "грязными" выглядят.
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Да уж, с jpeg я уже осознал провал... Очевидные вещи ведь. Надеюсь в данном случае это не критично, очень уж много скриншотов.
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Хорошо бы аргументировать этот шаг.
После того, как мы сменили регион, нужно изменить диапазон значений у всех предполагаемых к использованию файлов, т.е. у стандартных каналов (например, для синтеза в естественных цветах, B2,B3,B4) и у B8 на единый диапазон [0;255].
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Участник
- Сообщения: 78
- Зарегистрирован: 04 ноя 2014, 18:40
- Репутация: 48
- Откуда: Санкт-Петербург
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
В документации к модулю i.pansharpen сказано, что он работает только с 8-битными изображениями. Мы проверяли, действительно не работает без изменеия диапазона всех 4 каналов. Может это стоит отметить в статье, спасибо
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Добавил комментарий про rescale у паншарпенинга
- darsvid
- Активный участник
- Сообщения: 163
- Зарегистрирован: 29 июн 2012, 12:40
- Репутация: 115
- Откуда: Kyїv, Ukraine
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Спасибо, отличная статья
-
- Гуру
- Сообщения: 3321
- Зарегистрирован: 27 июл 2009, 19:26
- Репутация: 748
- Ваше звание: Вредитель полей
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Понятно, что так статья может разрастись, но еще есть http://grass.osgeo.org/grass64/manuals/ ... t.rgb.html
-
- Гуру
- Сообщения: 5173
- Зарегистрирован: 26 сен 2009, 16:26
- Репутация: 793
- Ваше звание: званий не имею
- Откуда: Москва
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
Спасибо за отличную статью!
В статье сказано: "В примере будет использован набор данных LC81740292013131, охватывающий юго-восточное побережье Кубы, поэтому используется имя Cuba. "
Указанный район имеет номера Path=11 и Row=45-47. В статье Path=174 и Row=29. Название набора данных перепутано.
Мелких замечаний, в основном по GRASS, довольно много, но это чуть позже, отдельно.
В статье сказано: "В примере будет использован набор данных LC81740292013131, охватывающий юго-восточное побережье Кубы, поэтому используется имя Cuba. "
Указанный район имеет номера Path=11 и Row=45-47. В статье Path=174 и Row=29. Название набора данных перепутано.
Мелких замечаний, в основном по GRASS, довольно много, но это чуть позже, отдельно.
Редактор материалов, модератор форума
- Эдуард Казаков
- Гуру
- Сообщения: 549
- Зарегистрирован: 23 апр 2014, 17:11
- Репутация: 534
- Откуда: Planet Earth
- Контактная информация:
Re: Обработка и интерпретация данных Landsat 8 в GRASS 7
С именем сцены глупая ошибка, конечно, везде должно быть LC80110462013131, именно она использовалась. Исправил.
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя