Индекс NDVI

[olechka9585]

Добрый день!
Стоит задача рассчитать индекс NDVI по снимкам с разных спутников Landsat. Чтобы иметь возможность сравнивать между собой индексные карты от разных спутников с разными камерами (калибруются по-разному) необходимо произвести конвертацию данных в показатели излучения на сенсоре (radiance), а затем конвертацию radiance в неатмосферное отражение (reflectance). Расчеты произвожу в ArcGIS с помощью калькулятора растра. Подскажите, я правильно понимаю, что после пересчетов индекс NDVI также должен получится от -1 до +1? Сколько не пытаюсь посчитать такого результата не получается. Может есть какой-то секрет?
Пока считаю для спутника Landsat 5 TM.
Расчеты произвожу по формулам:
L_λ=(L_maxλ-L_minλ)/Q_(cal max) ∙Q_cal+L_minλ (Здесь у меня был вопрос чему равно Q_(cal min), т.к. в метафайле написано, что равно 1, а во всех источниках говориться, что для Landsat 5 TM равно 0. Взяла 0)
L_maxλ и L_minλ беру из метафайла.
Пересчитав для 3 и 4 каналов произвожу конвертацию в reflectance:
ρ=(π∙L_λ∙d^2)/(〖ESUN〗_λ∙cosθ_s )
ESUN беру из таблицы. для 3 канала=1554, для 4=1036
d тоже из таблицы в соответствии с датой съемки.
θ_s =90°-SUN_ELEVATION (SUN_ELEVATION из метафайла), перевожу в радианы и считаю косинус.
получив значения в reflectance для 3 и 4 каналов считаю NDVI
NDVI=(NIR-RED)/(NIR+RED)

[evreyyy]

а какой получается результат?

[olechka9585]

от -7,2489 до 40,3626

[gamm]

от -7,2489 до 40,3626

у вас явные проблемы с арифметикой. Имеем RED>0, NIR>0, пусть RED=a*NIR, a>=0. Тогда (NIR-RED)/(NIR+RED) = (NIR-a*NIR)/(NIR+a*NIR) = (1-a)/(1+a). Это всегда в диапазоне от 1 (для а=0) до -1 (для а стремится к бесконечности)

[olechka9585]

При конвертации данных в показатели излучения на сенсоре тоже должно быть RED>0, NIR>0?