Причина ошибки при расчете NDVI по Landsat-8 (завышенные значения)

[_taras_]

По снимку Landsat-8 считаю NDVI, предварительно сделав коррекцию влияния атмосферы
В R

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

B4 = 0.0000275*(raster ("LC08_L2SP_181028_20150217_20200910_02_T1_SR_B4.TIF"))- 0.2
B5 = 0.0000275*(raster ("LC08_L2SP_181028_20150217_20200910_02_T1_SR_B5.TIF"))- 0.2
NDVI = (B5 - B4)/(B5 + B4)

Получаю зля зимнего времени значения в р-не 0,2-0,3, что много для февраля.
Ежели делать без коррекции, то значения будут близки к 0,1, что соответствует истине.

Подскажите, где ошибаюсь? Может для Landsat-8 уже ненужна коррекция?

[gamm]

По снимку Landsat-8 считаю NDVI, предварительно сделав коррекцию влияния атмосферы

кто, где и как делал атмосферную коррекцию? насколько я помню, у Ландсата Level 2 все лежит уже пересчитанное в наземный reflectance, с учетом атмосферы

[_taras_]

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

Насколько я понял из объяснений antonv

Несмотря на то, что это продукт, уже пересчитанный в surface reflectance, в пикселях записаны не эти самые reflectances в долях единицы, а отмасштбированные целочисленные значения в диапазоне от 0 до 65 тысяч.

И приведенная выше формула перерасчета вроде подходит для всех Ландсатов (для 4-7 так точно). Но 8 Ландсат для зимы дает завышенные значения.

[gamm]

вы не ответили, кто и как делал атмосферную коррекцию.

Если это файлы, скачанные с амариканских серверов, то SR говорит нам о Surface reflectance, т.е. атмосферная коррекция у них уже сделана. На что умножать и что добавлять, смотрите в описании конкретного продукта, у них константы могут меняться.

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

Reflectance bands have new scaling factors.
Collection 1 used a 0.0001 scale factor. 
Collection 2 uses a 2.75e-05 scale factor and -0.2 offset.

ваши файлы, похоже, Collection 2, см. стр.2RTFM