Форумы GIS-Lab.info

Подскажите,пожалуйста, правильно и я делаю коррекцию?
Открываю модуль 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)

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

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

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

Перед тем, как использовать модуль, очень рекомендую вдумчиво читать мануал.

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

проштудировала мануал, пробую делать коррекцию и опять не получается...что-то,видимо,делаю не так.
Мои действия
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кб

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

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

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

Действительно в ГРАСС 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

Вопрос к Юлии, а какие вы выставляете параметры для пересчета?

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

Значит я таки был неправ, пардоньте. Просто разработчики в свежей версии изменили значение параметра с "ot8" на "oli8".

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

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

.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)

Неясно, почему проходит коррекция только первых двух каналов?

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

Похоже, у вас только два растра, начинающихся на "LC81130252013190LGN00_B" давайте проверим: выполните команду:

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

g.mlist rast pattern="LC81130252013190LGN00_B*"

и сообщите, что она выводит.

пока не добралась до рабочего компьютера, но в папке файлы все называются одинаково

Озарян Юлия писал(а):пока не добралась до рабочего компьютера, но в папке файлы все называются одинаково

Всё-таки хочется взглянуть на вывод команды GRASS, а не на содержимое каталога, из которого вы импортировали данные --- вопрос то был как раз о том, всё ли импортировалось

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

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

Если вы пользуетесь интерфейсом wxGUI, то см. скриншот:

Форумы GIS-Lab.info

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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