Страница 1 из 4

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

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

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

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

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

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

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

Добавлено: 03 мар 2015, 20:18
Эдуард Казаков
Здравствуйте, Илья!
В принципе, статью можно воспринимать как пошаговое руководство с ответвлениями, поэтому начинающим должно помочь :)

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

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

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

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

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

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

Добавлено: 03 мар 2015, 23:24
ericsson
Как-то странно выглядит пара "до и после" для IHS - синий нормализовался, а вроде как не должен.

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

Добавлено: 03 мар 2015, 23:39
voltusya
Чуть выше говорилось о наличие флага Rebalance, пример приведен именно с его использованием.

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

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

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

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

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

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

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

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

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

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

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

Добавлено: 04 мар 2015, 16:42
Эдуард Казаков
Добавил комментарий про rescale у паншарпенинга

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

Добавлено: 04 мар 2015, 19:40
darsvid
Спасибо, отличная статья :)

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

Добавлено: 05 мар 2015, 08:19
ericsson
Понятно, что так статья может разрастись, но еще есть http://grass.osgeo.org/grass64/manuals/ ... t.rgb.html

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

Добавлено: 10 мар 2015, 21:05
Александр Мурый
Спасибо за отличную статью!

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

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

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

Добавлено: 11 мар 2015, 00:07
Эдуард Казаков
С именем сцены глупая ошибка, конечно, везде должно быть LC80110462013131, именно она использовалась. Исправил.