Снимок Landsat - естественные цвета

[Воронин С.В.]

Господа.
Скачал снимки Landsat с сайта USGS.
Но не могу получить из снимков композит с естественными цветами.
У меня такое ощущение, что где-то я такую возможность видел.
То ли на самом сайте USGS, то ли ещё где.
А сейчас не могу найти. Может кто подскажет.

И ещё.
Я понимаю, что каждый снимок сделан в некотором диапазоне.
Но какой в каком?
Подскажите диапазон снимков.

[gimran]

Скачайте полную версию снимка (там много файлов в архиве, несколько сот мегабайт).
И укажите версию ландсата

[xen87]

Статья (перевод) на хабре, посвященная диапазонам Ландсата 8: http://habrahabr.ru/post/183416/

[Воронин С.В.]

Господа.
Если кому интересно.
Получать композиционный растр лучше с использованием программ ENVI.
Композиты получаемые в ArcGIS значительно уступают по качеству.

[ericsson]

http://www.geosage.com/highview/features_landsat8.html
А использовать ENVI - это межконтинентальной ракетой по воробьям.

[ericsson]

Еще есть чудесное средство IJFusion http://www.ijfusion.es/ , но странные испанцы, которые его разработали, по не менее странной причине не дают его просто скачать, а выдают по запросу. Чем они думают - непонятно. Но там есть куча разных алгоритмов, я даже не про все слышал.

[ericsson]

А тут и в следующей статье "пошагово" учат делать это в GRASS (я не пользуюсь, но все понятно) http://courses.neteler.org/processing-l ... ass-gis-7/

[bolotoved]

Сравнение каналов landsat 7 и 8: http://landsat.usgs.gov/L8_band_combos.php
Т.о. для "естественных" цветов нужна комбинация 4:3:2

Композиты получаемые в ArcGIS значительно уступают по качеству

Как такое возможно? Входящие в композит слои при сборке не должны изменяться, по этому без разницы, в чем вы их собираете. А вот отображаются на мониторе композиты в разном ПО по разному - с пом разных алгоритмов, нужно разбираться с настройками.

P.S. Фантастические возможности регуляции отображения композитов присутствуют в QGIS.

[ericsson]

Композиты получаемые в ArcGIS значительно уступают по качеству.

Скорее всего, это означает, что приведение 16-битной картинки к 8-битной сделано неоптимальным способом, и наверняка это можно сделать иначе.

[ericsson]

Я тут попробовал для обработки Landsat 8 такую чудную штуку, как ImageJ.
В ней реально очень удобно превращать 16-битное изображение в 8-битное: и гистограмму для выбранного участка покажет, и потом найденные по ней границы применить даст, и из 16 в 8 бит одним щелчком переводит. Единственное - пока не нашел для него плагин для чтения/записи GeoTIFF, потому надо метаданные дампить listgeo/geotifcp или gdal.
Вообще, она даже сама композиты делать может, при том очень шустро. Нужно только под 64-битной Java ее запускать, а то памяти не хватит.

[Александр Мурый]

Я тут попробовал для обработки Landsat 8 такую чудную штуку, как ImageJ.
В ней реально очень удобно превращать 16-битное изображение в 8-битное: и гистограмму для выбранного участка покажет, и потом найденные по ней границы применить даст, и из 16 в 8 бит одним щелчком переводит. Единственное - пока не нашел для него плагин для чтения/записи GeoTIFF, потому надо метаданные дампить listgeo/geotifcp или gdal.
Вообще, она даже сама композиты делать может, при том очень шустро. Нужно только под 64-битной Java ее запускать, а то памяти не хватит.

ericsson, а вы не могли бы написать заметку об использовании ImageJ для обработки Landsat 8?

[Воронин С.В.]

А как использовать панхроматический 8-й диапазон из снимка Landsat 8 для большей детализации композита.
В статье http://habrahabr.ru/post/183416/ показан только результат.

[rhot]

потому надо метаданные дампить listgeo/geotifcp или gdal

Расскажите как.

[dime1]

А как использовать панхроматический 8-й диапазон из снимка Landsat 8 для большей детализации композита.
В статье http://habrahabr.ru/post/183416/ показан только результат.

Выбрав композит и панхроматический канал через Ctrl, нажать Pan-sharpening http://resources.arcgis.com/en/help/mai ... 00m7000000

[ericsson]

Да ну, какая там статья про ImageJ, там действий - раз, два, и готово.

Сначала чтобы дамп сделать:

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

listgeo yourfile.tiff >>yourfile.meta

Чтобы в конце вернуть:

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

geotifcp -g yourfile.meta yourfile_without_meta.tiff yourfile_with_meta.tiff

А в самом ImageJ все предельно просто:
Открыть файлы каналов в том порядке, в каком они пойдут в RGB.
Далее их нужно объединить в так называемый Stack, для чего Images->Stacks->Images to Stack.
Далее есть два варианта - либо превратить это в обычное RGB изображение, либо сделать отдельный тип данных Composite.
Для первого - Image->Color->Stack to RGB
Для второго - Image->Color->Make Composite

Статистика работает элементарно:
Analyze->Histogram.
Если нажать кнопку Live, то будет отображать все изменения. Если на изображении выбрать ROI (первые пять кнопок слева на панели) - будет показывать гистограмму оттуда. При возюкании курсором по гистограмме, показывает число пикселей каждой яркости. Так что подобрать ширину окна вручную - не проблема.

Для каждого изображения (будть то изначальный отдельный канал, готовое RGB-изображение или один из Slice внутри Stack/Composite ) работают любые фильтры.
Process->Enhance Contrast, например (самый простой и полезный для пробы) - там простой клиппинг и нормализация гистограммы.
Image->Adjust->Brightness/Contrast - это и так понятно что.
Image->Adjust->Window/Level - почти то же самое, только задает диапазон входных значений яркости и диапазон их отображения после коррекции.

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

Единственное, что жаль - это то, что испанская приблуда для pansharpening/fusion только по запросу, хотя и бесплатная. А то там возможностей очень много.