Атмосферная коррекция снимка Landsat8Oli с помощью модуля i.

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Атмосферная коррекция снимка Landsat8Oli с помощью модуля i.

Сообщение Озарян Юлия » 08 май 2014, 08:06

Подскажите,пожалуйста, правильно и я делаю коррекцию?
Открываю модуль i.atcorr, в поле "Имя исходной растровой карты" вставляю слой с необходимым каналом снимка, в поле Name of input text file вставляю файл с метаданными, забиваю имя выходного слоя, нажимаю "Выполнить", получаю такой ответ :?:
(Thu May 08 15:53:27 2014)
i.atcorr iimg=LC81130252013190LGN00_B3@PERMANENT icnd=C:\ya102330\OZARYAN\GLOVIS\Landsat8 Oli\LC81130252013190LGN00.tar\LC81130252013190LGN00\LC81130252013190LGN00_MTL.txt oimg=L8B3atcorr
wavelength less than 0.25 micron:
let's take s(l)=s(0.25)
WARNING: Unknown atmospheric model!
WARNING: Unknown aerosol model!
(Thu May 08 15:54:08 2014) Операция завершена (41 sec)

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Александр Мурый » 08 май 2014, 08:44

Сенсор Landsat 8 не поддерживается модулем <i.atcorr>. Для атмосферной коррекции Landsat 8 используйте i.landsat.toar ("sensor=oli8"). Подробнее почитайте на спец. вики-странице.
Редактор материалов, модератор форума

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 08 май 2014, 08:53

так и подумала

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 08 май 2014, 08:56

а если делать коррекцию более ранних ландсатов, всё так вбивается?

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Александр Мурый » 08 май 2014, 09:38

Перед тем, как использовать модуль, очень рекомендую вдумчиво читать мануал. :) Там ответы почти на все возникающие вопросы.
В параметре "icnd" указывается не файл с метаданными (*_MTL.txt), а создаваемый вручную текстовый файл с параметрами атм.коррекции. Для разных сенсоров параметры, ест-но, будут разными. Полное описание там же в мануале, пересказывать здесь смысла нет.
Если у вас что-то не будет получаться с конкретными данными — пишите, будем вместе разбираться.
Редактор материалов, модератор форума

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 27 май 2014, 04:12

проштудировала мануал, пробую делать коррекцию и опять не получается...что-то,видимо,делаю не так.
Мои действия
1.с помощью r.in.gdal импортирую файлы с каналами 1-11
в модуле i.landsat.tour вбиваю имя входного и выходного растра
2. прикрепляю фай с метаданными
3. выбираю сенсор Oli8 (кстати? у меня в выпадающем списке он значится как ot8)
4. кликаю "выполнить" , получаю
Calculating...
Writing reflectance of <LC81130252013190LGN00_B1> to <B.toar1>...
Writing reflectance of <LC81130252013190LGN00_B2> to <B.toar2>...

5. в результате, созданы два файла B.toar1 и B.toar2 оба весят по 0кб
Вложения
screenshot_2014_05_27_12_06_42.png
screenshot_2014_05_27_12_06_42.png (697.58 КБ) 6779 просмотров

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

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

Обнаружил, что у меня в GRASS 6.4.3svn в любом интерфейсе к модулю <i.landsat.toar> вообще нет варианта "oli8" для параметра "sensor", т.е. Landsat 8 в этой версии ещё не поддерживался. Обновился — появилось.

Так что очень рекомендую обновить GRASS до свежей стабильной версии 6.4.4 или (что ещё лучше) поставить активно разрабатываемую версию GRASS 7.
Никаких "ot8" в списках быть не может, вы, вероятно, спутали с "tm.."

Вопрос к Юлии: интересно, а как вы оцениваете "вес" выходных файлов? Мы же где-то раньше уже обсуждали, что подглядывание за папками/файлами, которые создаёт GRASS в своей базе данных, обычно ничем не помогает в работе, а наоборот, путает начинающего пользователя и мешает ему.
Редактор материалов, модератор форума

xen87
Активный участник
Сообщения: 168
Зарегистрирован: 08 апр 2009, 14:44
Проекты: 1
Репутация: 18
Откуда: Кострома

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение xen87 » 27 май 2014, 14:11

Действительно в ГРАСС 6.4.3 в модуле i.landsat.toar для снимков Landsat 8 в качестве сенсора указан "ot8".
В качестве подтверждения выкладываю скриншот.
Сейчас проверил у себя - все сработало - без ошибок.
В мое варианте выставил такие параметры:

Код: Выделить всё

i.landsat.toar input_prefix=LC81760202013183LGN00_B output_prefix=LC81760202013183LGN00_B_toar_ metfile=H:\Arhiv_snimkov\2013\LC81760202013183LGN00\LC81760202013183LGN00_MTL.txt sensor=ot8
Вопрос к Юлии, а какие вы выставляете параметры для пересчета?
Вложения
landsat_8_oli.png
landsat_8_oli.png (48.85 КБ) 6726 просмотров

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Александр Мурый » 27 май 2014, 14:20

xen87 писал(а):Действительно в ГРАСС 6.4.3 в модуле i.landsat.toar для снимков Landsat 8 в качестве сенсора указан "ot8".
Значит я таки был неправ, пардоньте. Просто разработчики в свежей версии изменили значение параметра с "ot8" на "oli8".
Редактор материалов, модератор форума

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 28 май 2014, 03:59

Версию программу обновлю, спасибо.
Мои параметры и результат такие :

Код: Выделить всё

.landsat.toar input_prefix=LC81130252013190LGN00_B output_prefix=LC81130252013190LGN00_B_toar_ metfile=C:\ya102330\OZARYAN\GLOVIS\Landsat8 Oli\LC81130252013190LGN00.tar\LC81130252013190LGN00\LC81130252013190LGN00_MTL.txt sensor=ot8
Calculating...
WARNING: 'cell/LC81130252013190LGN00_B1' был найден в нескольких наборах (также найден в <PERMANENT>)
WARNING: Используется<LC81130252013190LGN00_B1@NDVI1>
Writing reflectance of <LC81130252013190LGN00_B1> to <LC81130252013190LGN00_B_toar_1>...
WARNING: 'cell/LC81130252013190LGN00_B2' был найден в нескольких наборах (также найден в <PERMANENT>)
WARNING: Используется<LC81130252013190LGN00_B2@NDVI1>
Writing reflectance of <LC81130252013190LGN00_B2> to <LC81130252013190LGN00_B_toar_2>...
(Wed May 28 11:57:52 2014) Операция завершена (2 min 3 sec) 
Неясно, почему проходит коррекция только первых двух каналов?

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение KolesovDmitry » 28 май 2014, 09:08

Озарян Юлия писал(а):Неясно, почему проходит коррекция только первых двух каналов?
Похоже, у вас только два растра, начинающихся на "LC81130252013190LGN00_B" давайте проверим: выполните команду:

Код: Выделить всё

g.mlist rast pattern="LC81130252013190LGN00_B*"
и сообщите, что она выводит.

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 29 май 2014, 15:47

пока не добралась до рабочего компьютера, но в папке файлы все называются одинаково :roll:
Вложения
Новый рисунок (1).JPG
Новый рисунок (1).JPG (21.79 КБ) 6625 просмотров

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

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение KolesovDmitry » 29 май 2014, 17:35

Озарян Юлия писал(а):пока не добралась до рабочего компьютера, но в папке файлы все называются одинаково :roll:
Всё-таки хочется взглянуть на вывод команды GRASS, а не на содержимое каталога, из которого вы импортировали данные --- вопрос то был как раз о том, всё ли импортировалось

Озарян Юлия
Активный участник
Сообщения: 131
Зарегистрирован: 22 май 2013, 01:35
Репутация: 11

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение Озарян Юлия » 30 май 2014, 08:14

а куда команду вводить? ....в командную строку?

xen87
Активный участник
Сообщения: 168
Зарегистрирован: 08 апр 2009, 14:44
Проекты: 1
Репутация: 18
Откуда: Кострома

Re: Атмосферная коррекция снимка Landsat8Oli с помощью модул

Сообщение xen87 » 30 май 2014, 08:59

а куда команду вводить? ....в командную строку?
Если вы пользуетесь интерфейсом wxGUI, то см. скриншот:
Вложения
screenshot_2014_05_30_09_55_37.png
screenshot_2014_05_30_09_55_37.png (56.85 КБ) 6571 просмотр

Ответить

Вернуться в «GRASS»

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

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