Olá, estou fazendo um trabalho e preciso analisar uma série temporal de imagens landsat. Para isso estou usando 3 imagens do landsat 5 (para 3 anos diferentes) nas bandas 3 e 4, e uma mais atual do landsat 8 bandas 4 e 5. Acontece que depois de convertido os ND para Reflectância de todas as cenas, eu faço o cálculo do NDVI mas o resultado correspondente ao landsat 8 apresenta grande diferença para o landsat 5. Principalmente depois que eu aplico uma renderização no NDVI.
Estou usando Qgis 2.14.17 em Windows 10 64bits.
imagens landsat 5 b3 e b4 / landsat 8 b4 e b5
Já fiz tudo pela calculadora raster usando as formulas abaixo e também já fiz por plugins (Geosud e Semi-Automatic Classification Plugin), mas o resultado é o mesmo e não satisfatório!

Fórmula NDVI usado na calculadora raster após transformação de ND para Reflectancia: (B4 - B3) / (B4 + B3) e (B5 - B4) / (B5 + B4)
Eu percebi que os valores de NDVI que o LANDSAT 8 gera é sempre inferior, e isso atrapalha na hora que estou renderizando e atribuindo cores

percebam que ele sempre fica entre variando entre -0,6.... à 0,8... , enquanto os NDVI gerado a partir das bandas 3 e 4 do sensor TM ficaram entre (-0,9...à 0,9...) (-0,8... à -0,7...) (-0,8... à 0,8..) . Isso atrapalha quando vou abrir e colocar todos nessas "classes" de valores, pois como pode ser observado ao lado, a área fica muito diferente.

Estou a disposição para demais informações!