Как "сшить" Landsat?

[Nule]

Здравствуйте!
Существует ли инструмент для "сшивки" сцен landsat в мозаику. Или только обрезать вручную?
Спасибо!

[rhot]

Существует, и не один. Вот, например, статья.

[Nule]

Спасибо! А в ENVI есть такое?

[rhot]

есть, конечно

[nadiopt]

basic tools - mosaiking, чего уж проще

[Nule]

basic tools - mosaiking, чего уж проще

Может, я совсем отупел, но для меня это сложно:
По вертикали склейка не происходит:



По горизонтали еще интереснее, инструмент сам решает, что ему склеивать а чего так оставить:



А вот шов:



Так сшить и в gdal можно. А точнее не сшить, а наложить, как здесь, собственно, и сделано.

Подскажите, пожалуйста, как правильно использовать этот инструмент?

[Nule]

Насчет сшивает\не сшивает разобрался. Надо последовательно импортировать каждый файл, тогда всё нормально. Если импортировать пакетно, получается ерунда.
Можно как-нибудь настроить, чтобы боковые артефакты сцен обрезались, а не просто одна накладывалась на другую?

[Saimon]

Создавайте линии сшивки/порезов

[Nule]

Создавайте линии сшивки/порезов

В инструменте basic tools - mosaiking не могу найти ничего похожего

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

*Nule
Стоит попробовать ещё GRASS с модулем r.patch для сшивки. Сшивать можно в разном порядке (как укажете). Для более тонкой сшивки нужно будет иметь дело уже с растровыми масками и растровой алгеброй.

[Nule]

*Nule
Стоит попробовать ещё GRASS с модулем r.patch для сшивки. Сшивать можно в разном порядке (как укажете). Для более тонкой сшивки нужно будет иметь дело уже с растровыми масками и растровой алгеброй.

Я надеялся, что для такой массовой вещи как landsat всё уже давно создано и оптимизировано.

[KolesovDmitry]

*Nule
Стоит попробовать ещё GRASS с модулем r.patch для сшивки. Сшивать можно в разном порядке (как укажете). Для более тонкой сшивки нужно будет иметь дело уже с растровыми масками и растровой алгеброй.

Если использовать чистый r.patch, то после сшивки будет виден шов. А с помощью растровой алгебры можно получить бесшовное соединение. Нужно найти пересечение сшиваемых снимков P1 и P2 (условно первый и второй снимки), потом расчитать растр

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

b*P1+(1-b)P2

где величина b монотонно меняется так, что она принимает значение 0 на "дальней" от шва границе первого снимка, 1 на дальней границе второго снимка и 0.5 на самом шве.
Основная сложность была в вычислении b (это даже не очень сложно -- я делал пару лет назад, но подробности сейчас не вспомню).

[nadiopt]

насчет убирания краевых артефактов могу посоветовать только обрезать растр по маске (создайте roi) размером чуть поменьше снимка, чтоб артефакты остались за бортом

[Nule]

В конечном итоге решил использовать MosaicPro из Erdas Imagine 2011. Поканально мозаики получаются. А можно как-нибудь сразу и сборку по каналам сделать? А то приходится 3 раза одно и тоже делать. Если просто открывать несколько каналов, то они накладываются, но не смешиваются.

[Ремесленник]

Уважаемые Господа!
А каким образом (кроме модуля r.patch) в GRASS можно сшить многоканальные или одноканальные растры?