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

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

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

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

Обещанные ранее замечания.
  • Может быть, стоит добавить в текст прямую ссылку на набор данных, используемый в статье? http://earthexplorer.usgs.gov/download/ ... 13131LGN01
    Чтобы читателям не искать самостоятельно. Зарегистрированные пользователи смогут залогиниться и быстро скачать архив.
  • Может быть, есть смысл дать ссылки на похожие статьи по этой теме от самих разработчиков GRASS?
    http://courses.neteler.org/processing-l ... ass-gis-7/
    http://courses.neteler.org/processing-l ... harpening/
  • В тексте при первом упоминании каждого модуля было не неплохо давать ссылку на документацию, причём для именно для версии GRASS 7 (http://grass.gis-lab.info/grass70/manuals/index.html)
  • Работоспособность самой актуальной сборки GRASS GIS 7 на день публикации статьи является адекватной только в операционных системах семейства Linux и в Mac OS. Поведение пакета в Windows непредсказуемо и нестабильно, большая часть из описанных в статье операций в ней не работает (например, генерируются пустые файлы или программа завершает работу с ошибкой).
    Я специально провёл эксперимент: установил с нуля "stand-alone"-версию GRASS 7 и проделал бОльшую часть действий из статьи. Результат: ни одной ошибки не выдалось, никаких пустых файлов.
    Предлагаю подробнее обсудить это, т.к. основная часть пользователей всё-таки под вёндами сидит, и напрасно отпугивать их от GRASS тоже не хотелось бы :)
  • При установке GRASS 7 на Windows в первый раз она запускается с указанием установленной вместе с программой т.н. "пробной" области ("Demolocation"), которая находится по умолчанию в уже созданной папке "grassdata" в директории типа "Мои документы" текущего пользователя. В статье же описана ситуация, когда у пользователя не создана директория с данными GRASS (т.н. "база данных") и нет выбора области, т.е. вообще пусто. Это может запутать начинающего пользователя (опытные разберутся и так).
  • GRASS GIS всегда показывает текстовую команду, соответствующую выполняемым в интерфейсе действиям
    Это не совсем так. Не надо путать с окном MapBasic в MapInfo :) Не все действия, производимые в интерфейсе, отображаются во вкладке "Command console".
  • При импорте данных можно сразу активировать флаг модуля <r.in.gdal> "Extend region extents based
    on new dataset" , и не надо будет задавать вычисл. регион потом отдельно.
  • Есть установленный термин для "region" - "вычислительный регион". Предлагаю использовать его, либо просто "регион", а не "region". Всё-таки статья на русском языке :)
  • В сообществе пользователей GRASS есть старая-добрая традиция всегда приводить точный текст выполняемой команды помимо скриншотов, причём именно в виде текста, чтобы можно было применить волшебный copy/paste :) И здесь это тоже очень просится. В вики лучше сделать в виде листинга (<pre>).
  • Как верно отметили коллеги выше, стоит упомянуть спец. модуль <i.landsat.rgb>, который в GRASS 7 переименован в <i.colors.enhance> (в меню Imagery - Manage image colors)
  • Может быть, имеет смысл дать прямые ссылки на данные SRTM, используемые в статье?
    http://gis-lab.info/data/srtm-tif/srtm_21_08.zip
    http://gis-lab.info/data/srtm-tif/srtm_21_09.zip
    Ну или хотя бы можно облегчить поиск, дав номера CGIAR-мозаик (21_08 и 21_09).
  • Поскольку индексное изображение одноканально, по умолчанию GRASS отрисует его в оттенках серого
    Откуда такой вывод? GRASS не умеет хранить многоканальные растры - они всегда разбиваются по отдельным каналам. Откуда информация, что по умолчанию для растров серая цветовая шкала? У меня после расчёта NDVI цветная получилась ("rainbow").
Редактор материалов, модератор форума

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

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

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

Спасибо, Александр!

Прямые ссылки на использованные данные мы не вставляли, т.к. не предполагали, что пользователю будет интересно работать с той же территорией, всё же это просто пример. Но, наверное, вы правы - новичку в первый раз захочется делать тоже самое и видеть тоже самое, что и в примере, чтобы убедиться, что он всё правильно делает. Добавлю.

Ссылки на Нетелера добавлю в доп. литературу, хотя там ничего нового относительно статьи вроде бы нет.

Статья писалась ещё до стабильного релиза (работаем с beta3). Видимо, в этом причина неработоспособности под Windows, хотя на наших операционках эта Beta работает безупречно. И Demo-location там тоже нет. Исправлю.

По поводу команды - в нижней части окна модуля всегда пишется текст команды, аналогичной настроенной в GUI, то имелось ввиду. Но и правда, лучше сделать. Правда, возможны сложности с операциями в интерактивном режиме.

Маленький раздел про i.colors.enhance добавлю.

В общем, работаем :) Глядишь, и правда адекватная статья получится.

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

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

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

Спасибо за ответы :)
Со своей стороны, могу предложить наделать новых скриншотов в PNG. Только вот вопрос в чём делать - в рабочем Дебиане или в виртуальной винде? Наверное, в Win всё же лучше, т.к. привычнее глазу потенциального читателя и лучше унифицировано будет (а то "один скрин в Lin, другой в Win, третий в Mac" таки мешает восприятию.
Редактор материалов, модератор форума

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

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

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

Честно говоря, сомневаюсь, что это занятие стоит вашего времени :) Сделаем и скриншоты! Уж если идти, то до конца. Заодно опробуем всё в Windows на стабильной версии.

KolesovDmitry
Гуру
Сообщения: 806
Зарегистрирован: 22 авг 2007, 14:58
Статьи: 18
Проекты: 4
Репутация: 116
Откуда: Казань

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

Сообщение KolesovDmitry » 15 мар 2015, 20:16

Спасибо за статью! Очень большая и качественная работа.

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

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

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

Новая итерация.
Исправления:
  • Всё переделано в новой стабильной версии Grass 7 в ОС Windows.
  • Добавлены прямые ссылки на использовавшиеся данные: Landsat, SRTM.
  • Переделаны скриншоты в .png.
  • Добавлены ссылки на документацию для каждого модуля grass при первом его (модуля) упоминании.
  • В разделе про Pan Sharpening дополнительно рассказывается про i.colors.enhance. Синтез в естественных цветах (обычный и Pan Sharpening) визуализирован после этой операции.
  • Почти каждая операция сопровождается полноценным текстом консольной команды (кроме выполняющихся в интерактивном режиме).
  • Небольшие корректировки по всему тексту. Убраны необоснованные заявления. Добавлены новые примечания и особенности работы.

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

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

Сообщение Александр Мурый » 31 мар 2015, 09:11

Спасибо за работу. Пара небольших замечаний:
-
Создадим векторный файл с помощью меню Vector - Develop Vector maps - Create new vector map.
Создавать пустой новый слой — лишнее действие. Проще просто включить оцифровку через "Digitize" в правом углу окна, окно с диалогом всё равно появится.

- Для оформления карт можно (и логичнее, кстати) открывать данные прямо в QGIS через модуль GRASS.
Редактор материалов, модератор форума

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

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

Сообщение voltusya » 31 мар 2015, 10:53

Здравствуйте!

- Создание нового слоя предполагалось для ознакомления с тем, как это делается в принципе. Возможно, правда, стоит упомянуть и о другой возможности создания файла сразу при векторизации, спасибо

-Такая же история, показать как экспортировать файлы, если их будут использовать для других целей.

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

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

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

Добавил пояснения про создание нового слоя и про поддержку GRASS в qgis и скриншоты.

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

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

Сообщение SS_Rebelious » 04 апр 2015, 03:09

Спасибо за статью!

В глаза бросается следующая фраза:
Следует отметить, что при создании "улучшенного" изображения отдельно взятого канала происходит, по большому счету, изобретение информации, которой в исходном наборе данных нет. Производить классификацию или какой-либо спектральный анализ по таким изображениям бессмысленно.
Подобные утверждения необходимо обосновывать. Вот, например, статья GeoEye-1 and WorldView-2 рan-sharpened imagery for object-based classification in urban environments (PDF) - улучшенные изображения успешно применяются для классификации урбанизированных территорий. Сохранению спектральной информации (или оценке её потерь) при паншарпенинге тоже посвящён ряд статей.
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
Гуру
Сообщения: 3154
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 664
Ваше звание: Вредитель полей

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

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

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

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

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

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

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

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

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

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

Не думаю, что имеет смысл вести подробное обсуждение в данной статье.
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
Гуру
Сообщения: 3154
Зарегистрирован: 27 июл 2009, 19:26
Репутация: 664
Ваше звание: Вредитель полей

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

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

Я бы просто поменял спорную фразу на что-то вроде "pansharpening существенно изменяет свойства изображения с точки зрения классификации и спектральных характеристик снимка".

Аватара пользователя
tsatsara
Интересующийся
Сообщения: 28
Зарегистрирован: 30 авг 2011, 12:22
Репутация: 3
Откуда: Красная Поляна

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

Сообщение tsatsara » 04 апр 2015, 16:14

Все очень здорово, вы молодцы!
не могу удержаться, простите: сельскохозяйственные угодья в легенде
а не сельско-хозяйственные
но это, конечно, ерунда и к теме не относится

Ответить

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