Странный результат при расчете NDVI по Landsat 8

[Ostran2004]

День добрый, столкнулся со странным результатом при расчете NDVI по Landsat 8.
Использовались снимки 2 уровня обработки (L2) с сайта NASA, расчет в ArcGIS 10.8 через растровый калькулятор.
формула:
Float("LC08_L2SP_111026_20210818_20210827_02_T1_B5.TIF" - "LC08_L2SP_111026_20210818_20210827_02_T1_B4.TIF")
/ Float("LC08_L2SP_111026_20210818_20210827_02_T1_B5.TIF" + "LC08_L2SP_111026_20210818_20210827_02_T1_B4.TIF")

В результате получились "перевёрнутые" цифры для растительности: максимальные значения имеют гари и вырубки, а коренные темнохвойные и смешанные леса значения такие же или ниже.
Столкнулся с таким впервые, в чем может быть причина?

[gamm]

я бы для начала сравнил с NDVI на эту же дату и территорию, полученным с других спутников, например MOD13/MYD13 с Terra/Aqua

[Константин Силкин]

Добрый день! Радиометрическую коррекцию проводили?

[gamm]

Радиометрическую коррекцию проводили?

level2 - это уже с коррекцией, surface reflectance.

[Константин Силкин]

По-моему нормальные значения. Низковатые, но не перевёрнутые

[Ostran2004]

Да нет, перевернутые, см картинку

[gamm]

Да нет, перевернутые, см картинку

судя по картинке, на свежей гари 0. Остальное смотреть нужно за прошлые годы на разных спутниках. Этот фрагмент снимка выложите в GeoTIFF, чтобы посмотреть можно было, если снимки под рукой найду.

[Ostran2004]

Спасибо, не надо.
Посмотрел другие снимки Landsat - 5, 7, 8, для этой и соседних территорий (Юг Хабаровского края, Сихотэ-Алинь). Оказалось елово-пихтовые и другие темнохвойные леса высокой сомкнутости (полнота 0.7 и выше) из-за сильного затенения на нашей территории имеют низкие значения NDVI, так как практически отсутствует подлесок и травяной покров, а так же лиственные деревья в первом ярусе. Как результат - значения NDVI 0.2 - 0.35, то есть для детектирования рубок и оценки лесовосстановления нужно или использовать иной сезон, или подбирать иные вегетационные индексы.