Улучшить разрешение с помощью панхрома

[tanyshka]

Доброе время суток!!! Нужна помощь!!! Как с помощью панхрома улучшить разрешение многозонального изображения? речь идет о ландсате!!! ведь панхом 15 м, а так 30 м... если есть ссылка, буду признательна! Спасибо!

[Анна]

какой софт есть под рукой?
процедура называется pansharpening.
В ERDAS - Resolution merge (Interpreter->Spatial Enchancement -> Resolution merge)
В ENVI - Image sharpening (Transform -> Image sharpening)

[Михаил Жиленев]

Здравствуйте, tanyshka.
Чтобы сделать синтез pansharpening ENVI или Resolution merge в ERDAS (делал по Lansat в обоих данных продуктах и еще в PCI Geomatica), необходимо чтобы пространственное разрешение панхрома было в 3-5 выше, чем у мультиспектра. Соответственно, для 28-30м Landsat-овских данных надо иметь панхром около 6-10м. Я делал синтез наших панхроматических 8-ми метровых данных с КА "Монитор-Э" на Московскую обл. и Москву, Ленинградскую обл. и Санкт-Петербург. Получилось очень даже впечатляющее. Попутно проявился еще один крайне интересный эффект.Данные Landsat и "Монитор-Э" были значительно разнесены во времени. Если заинтересует панхром с "Монитор-Э", то можно посмотреть возможности как Вам его получить.
Желаю успехов.
С уважением

[Дмитрий Барышников]

какой софт есть под рукой?
процедура называется pansharpening.
В ERDAS - Resolution merge (Interpreter->Spatial Enchancement -> Resolution merge)
В ENVI - Image sharpening (Transform -> Image sharpening)

Еще в ArcGIS на лету можно это же сделать. Но его модели менее качественные.

[KolesovDmitry]

Сделать это можно в любом ПО, предоставляющем возможность работы с растровой арифметикой (растровые калькуляторы). Идеи и разные варианты описаны здесь. Возьмите такую формулу (используется для создания RGB растров):

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

                                   DN_b1 
DN_fused =      -----------------------------------      * DN_pan
                         DN_b1 + DN_b2 + DN_b3

Где DN_fused - улучшенный растр DN_b... - соответствующие каналы, DN_pan - панхром.

Конкретно для LandSat получается:

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

DN_fused.r  = tm.5 / (tm.2 + tm.4 + tm.5) * etmpan
DN_fused.g = tm.4 / (tm.2 + tm.4 + tm.5) * etmpan
DN_fused.b = tm.2 / (tm.2 + tm.4 + tm.5) * etmpan