Aster in ERDAS

[Вадимчик]

Я тут новичок. Подскажите пожалуйста как ERDASом(других пакетов не имеем) в астеровских снимках можно избавиться от вертикальных полос во внировских и свировских каналах? Горизонтальные в тировских и некоторых свировских каналах я убрал с помощью Fourier Analysis - Periodic Noise Removal.

И в чем разница при импорте из ASTER HDF EOS формата в IMG с включенной галочкой Destripping и без него?

Премного благодарен.

[alexandr cherepanov]

версия erdas какая?
ASTER L1A надо полагать.

Чтоб избавиться от полосок нужна радиометрическая коррекция для vnir для swir в radiance для tir в С или К
радиометрическая коррекция по форулам:
VNIR
L =AV*V/GV+DV
где:
Dv, Av, Gv radiometric correction coefficients for VNIR from RadiometricCorrTable.
SWIR
L = AS*V/GS + DS
где:
Ds, As, Gs radiometric correction coefficients for SWIR from RadiometricCorrTable.
TIR
L = C0 + C1 * V + C2 * V ^ 2
где:
C0, C1, C2 radiometric correction coefficients for TIR from RadiometricCorrTable.

Destripping - радиометрическая коррекция + перевод в unsigned 8 bit - VNIR, SWIR и unsigned 16 bit - TIR.
полосы также исчезнут что и при просто conversion radiance

viewtopic.php?t=1409
http://gis-lab.info/qa/aster_radiocorr.html

в принципе при импорте полосы должны были убраться и без всяких дополнительных преобразований, но это если версия от 8.6 и выше, в 8.5 фигня вроде была полная, уже не помню что именно.

[Вадимчик]

Спасибо!
Да ASTER L1A
Erdas 8.7
conversion radiance - не избавил от вертикальных полос.
Destripping - тоже не избавил.
Но как с помощью только ердаса убрать полосы во внировских и свировских каналах?? Даже если таблицы извлечем, как их связать с пикселами??

[Максим Дубинин]

ERDAS всегда плохо делал радиометрическую коррекцию ASTER L1A, лучше всего с этим справлялась ENVI, правда было это давно, но не уверен, что что-то поменялось. Про то, как связать данные с пикселами и рассказывается в http://gis-lab.info/qa/aster_radiocorr.html

[alexandr cherepanov]

Erdas 8.7 должен делать точно, наверное где-то допустили ошибку, может забыли переключать каналы при обработке? (next image), радиометрическую коррекцию erdas сейчас делает нормально но не совсем удобно для пользователя. Если снимков Aster предстоит обрабатывать много, то лучше поискать ENVI 4.3 делает автоматом радиометричекую коррекцию для всех каналов при открытие снимка + можно выполнить геометрическую коррекцию чтобы решить вот эту проблему viewtopic.php?t=1409

[Вадимчик]

Частично проблему удалось решить импортом внировских и свировских каналов с "дестриппингом", с последующей обработкой Fourier Periodic Noise Removal.
Теорию по радиометрической корррекции я понял. А вот как это (" Радиометрическая коррекция VNIR данных ASTER ") осущесвить в Ердасе ?

[Максим Дубинин]

получить коэффициенты из HDF и применить к растру через modeller

[Вадимчик]

Так коэффициенты, насколько я понял - для каждого пиксела свои.

[Максим Дубинин]

коэффициенты свои для каждого элемента пзс линейки, то есть для пикселей одной колонки они одинаковые

[aspirant-1]

Доброго времени суток!

Не совсем понял про коэффициенты, поясните пожалуйста.

Открываю в Excel файл Band4.RadiometricCorrTable (Band4 - это просто к примеру), вижу 3 колонки значений, так какие из них необходимо использовать? какие из них Ds, As, Gs?

Данные в формате GeoTIFF

[alexandr cherepanov]

А зачем Вам они вообще?

В таком порядке и должны идти - D, A, G
A=slope
G=gain
D=оффсет

[aspirant-1]

Чтобы избавиться от полос в SWIR каналах по этому выражению:

L = AS*V/GS + DS

Выполять всё буду в ERDASе, так уж поручили.

Так каие же коэффициенты брать? самые верхние

"коэффициенты свои для каждого элемента пзс линейки, то есть для пикселей одной колонки они одинаковые"

Так какие же коэффициенты брать? где они находятся в этой таблице? в каком ряду?

[alexandr cherepanov]

Ну ERDAS должен уметь это делать и без дополнительных манипуляций руками, или у вас совсем древний ERDAS?

Каждая группа коэффициентов для одного столбца изображения, количество строк с коэффициентами=количеству столбцов на изображение.

[aspirant-1]

ERDAS 9.2. Автоматически он делает destripe, если данные в HDF формате, а у меня в GeoTIFF.

Ну подскажите, пожалуйста, пример выражения, которое вводится в spatial modeller.) А то я уже весь лоб разбил.

[Saimon]

Добрый день,

возможно ли получить Ваши данные для эксперимента, а то у меня только в HDF.
Данные в GeoTIFF хранятся отдельными каналами или по типам каналов?

С уважением,