Форумы GIS-Lab.info

У меня большая проблема, может быть, вы могли бы посоветовать, в чем тут может быть дело?
Мне нужно посчитать NDVI для Landsa5 TM: LT52190762011111CUB00
Я уже использовала коррекцию QUAC для ландсатов 5, потому, хотела и сейчас так сделать. Но где-то у меня ошибка, поскольку не считается правильно (белый экран).

Мои шаги
1. Открываю канал 2, 3 и 4 Landsa5
2. Через Layer Stacking склеиваю в один файл
3. Открываю его через Edit ENVI HEADER.
4. В нем через Edit Attributes вношу следующие параметры:
1. Wavelength
band2 = 0.60
band3 = 0.69
band4 = 0.90

5. Units - ставлю Micrometers
6. Sensor Type - Landsat TM
7. Pixel Sizes - X Y = 30.0 meter
8. Провожу атм коррекцию QUAC
9. В Band math считаю NDVI по формуле
10. (B1-B2)/(B1+B2)
11. "B1 = band4
12. B2 = band3"

Итогом – ничего не получается (белый экран, все значения по всему экрану Data = 0)

Что я сделала не так? Может быть какой-то параметр упущен в настройке на стадии Edit Envi Header


Может быть кто-нибудь мог бы скачать тот же снимок - получается ли у него коррекция и расчет NDVI?

сейчас Гловис не качает
мне помнится, что в Эрдасе это дело еще как-то нормализовалось. так же не могу взять в толк, а второй канал-то зачем?

nadiopt писал(а):сейчас Гловис не качает
мне помнится, что в Эрдасе это дело еще как-то нормализовалось. так же не могу взять в толк, а второй канал-то зачем?

Для QUAC нужно 3 канала, иначе он не работает. Прибавила просто 2-й

nadiopt писал(а):сейчас Гловис не качает
мне помнится, что в Эрдасе это дело еще как-то нормализовалось. так же не могу взять в толк, а второй канал-то зачем?

А вносимые параметры на стадии Envi Header - полностью правильны? Может быть там какой-нибудь integer - float нужно где-то ставить, или еще что-то?


Вроде бы делала точно по этому алгоритму. А сейчас что-то не выходит.

как говорит Александр Черепанов, float лишним не бывает. Пропишите его к каждому бэнду

nadiopt писал(а):как говорит Александр Черепанов, float лишним не бывает. Пропишите его к каждому бэнду

А это, оказывается, невозможно
Data Type - может быть только Byte Другое почему то не возможно

При всем ином значении Data Type сам QUAC не считается. А раньше сам QUAC считался. но вот NDVI по нему - нет

раз одна сцена, может ну ее, атмосферную коррекцию?

Это какой то бред. ГДЕ там ошибка???

band4 +band3 - считает, выдает растр
band4 - band3- считает, выдает растр

Делю второе на первое - что в ENVI что ArcGIs - ОШИБКА. Либо ноль либо -1

может тупо скобочки поставить?

nadiopt писал(а):может тупо скобочки поставить?

Я вычисляла просто по отдельности, и со скобками.

Кто сможет помочь - прошу скачать этот снимок и попробовать повторить те шаги по расчету NDVI Может там в чем-то другом ошибка?
Или в FLAASH как сделать расскажите - я уже много раз пробовала освоить, не могу. Шаги описываются в хелпе упираясь в непонятное сразу же. Разбираюсь с одним - опять лабиринт шагов. На форуме - не отвечают. Напишите раз по шагам, как сделать эту коррекцию во Flaash. Не знаю почему всем наплевать на начинающих пользователей. Почему хелпы пишутся так, что без высшего специального образования (или долго пользования продуктом и прохождения мучений поиска самостоятельно) их не поймешь?

По FLAASH выложили же ссылку на полную инструкцию по обработке Landsat 7, SPOT 4 и MODIS там есть ошибки, но не настолько критичные чтобы по ней нельзя было сделать Landsat 5 уж точно. С Landsat 8 чуть сложнее, но тоже уже с полтора года ссылки на форуме на пошаговые инструкции есть.
http://blog.sina.com.cn/s/blog_764b1e9d01018tc5.html
http://blog.sina.com.cn/s/blog_764b1e9d01018ud7.html
Если ENVI совсем старая нужно только найти или сделать ldcm_oli.sli.

8. Провожу атм коррекцию QUAC
9. В Band math считаю NDVI по формуле
10. (B1-B2)/(B1+B2)

QUAC/FLAASH выдают вроде UINT на выходе, но точно не FLOAT.
Запишите формулу по другому уже много раз было на форуме.
(float(b1)-b2)/(float(b1)+b2)

alexandr cherepanov писал(а): QUAC/FLAASH выдают вроде UINT на выходе, но точно не FLOAT.
Запишите формулу по другому уже много раз было на форуме.
(float(b1)-b2)/(float(b1)+b2)

Теперь вспомнила. Но всё равно не получилось. И вы не напутали формуле, как ее читать? Почему float(b1) - а потом - b2 не вкоючено в скобки?

(float(b1)-b2)/(float(b1)+b2)[/quote]

Но я пробовала считать и прямо по вашей формуле.
И так, как мне показалось логичным

float((B1-B2)/(B1+B2))

Итог - всё тот же. Не работает.

alexandr cherepanov писал(а): QUAC/FLAASH выдают вроде UINT на выходе, но точно не FLOAT.

За ссылки про Flaash - спасибо. Может быть там понятно написано, посмотрю.

Но в моем примере - нет ли возможности разобраться? Там точно какая-то ошибка на стадии Edit Header.

Я открыла сейчас откорректированный файл и увидела, что Band4 в нем обнулился (все значения = 0)! То есть в это и есть причина ошибки. Но сам исходный канал - нормальный. То есть что-то такое произошло после Edit Header и QUAC от чего канал обнулился. Нет ли у вас идей, в чем причина?

А если посмотреть на сам header-файл? Или, на крайний случай, поделиться им с общественностью?
Потому что я более чем уверен, что ошибка относится к очевидным, а потому у других людей она, естественно, не воспроизведется.