Landsat 8, атмосферная коррекция при помощи grass i.landsat.

[Rumato]

Добрый день! Подскажите, пожалуйста, как сделать следующее:

Есть снимки Landsat 8, нужно сделать их атмосферную коррекцию. Я использую Qgis и плагин к Grass, добавляю в набор Grass-a растры LAndsat и потом делаю команды

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

i.landsat.toar input_prefix=L8. \
  output_prefix=l8_toar \
  metfile=L5121060_06020060714_MTL.txt 

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

Мне кажется, что нужно ещё какие-то манипуляции произвести или это мало применимо. Вопрос: как правильно при помощи данного ПО?

[bolotoved]

как-то не очень похоже, что качество улучшелось

i.landsat.toar переводит Digital Numbers (безразмерные еденицы в которых хранятся сырые данные Landsat) в Top-of-Atmosphere Radiances (излучение попадающее на сенсор, которое имеет размерность). Подробнее см. тут: http://gis-lab.info/qa/dn2radiance.html

Как видно из формул пересчета, никакого отношения к "улучшению" данная операция не имеет. Просто преобразование из одних единиц в другие.

Чтобы в QGIS адекватно открыть одноканальный снимок из GRASS, нужно в свойствах отображения поставить "одноканальное серое", затем растянуть гистограмму по вкусу. Если преобразование i.landsat.toar прошли верно, то визуальной разницы между исходным и конечным снимками не должно быть.

[Rumato]

А вообще для улучшения, качества снимков (убрать дымку, часть облаков, улучшить детализацию) как-то можно использовать GRASS\QGIS\GDAL? В них если какие-либо модули для этих процедур?

[bolotoved]

А вообще для улучшения, качества снимков (убрать дымку, часть облаков, улучшить детализацию) как-то можно использовать GRASS\QGIS\GDAL? В них если какие-либо модули для этих процедур?

Чтобы убрать облака, нужно, чтобы у вас был не один снимок: облачные участки должны быть продублироаны безоблачными участками на другом снимке, тогда да, можно.
Улучшить детализацию можно при наличии у снимка более детального панхроматического канала, у ландсата 7 и 8 такой есть. Посмотрите:

http://courses.neteler.org/processing-l ... ass-gis-7/
http://courses.neteler.org/processing-l ... harpening/

[Rumato]

Большое спасибо за помощь, на счёт добавления панхроматического канала Я пытаюсь его добавить вот так:
http://grass.osgeo.org/grass70/manuals/ ... arpen.html

Использую Qgis и плвгин для работы с GRASS, но вот на шаге

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

i.pansharpen ms3=b_b ms2=b_g ms1=b_r pan=b_h  output_prefix=ihs321 sharpen=ihs

Выдаёт ошибку "команда не найдена"

Мне кажется, что нужно до установить этот модуль, так?

Я правильно мыслю? Но если нужно доустановить, не можете ли подсказать как это сделать?

[bolotoved]

Использую Qgis и плвгин для работы с GRASS

Лучше установить GRASS-7 (это его модуль): http://grass.osgeo.org/news/37/117/GRAS ... 0-0-beta3/
Народ сейчас на него массово переходит и возможностей у него больше.