Estou aprendendo a executar a correção atmosférica manualmente, portanto, não usando os métodos automatizados disponíveis no ENVI
os metadados (banda 1) para minha cena do Landsat 8 são os seguintes;
RADIANCE_MULT_BAND_1 = 1.2344E-02
RADIANCE_ADD_BAND_1 = -61.72016
REFLECTANCE_MULT_BAND_1 = 2.0000E-05
REFLECTANCE_ADD_BAND_1 = -0.100000
SUN_ELEVATION = 48.12237303
Então, quando eu uso a matemática da banda para calcular o TOA Radiance, digito o seguinte:
(0,012344 * B1) -61,72016
E produz uma estatística espectral que parece certa
Quando uso qualquer uma das seguintes fórmulas matemáticas da banda;
Refletância TOA (excluindo o ajuste da elevação do sol)
(B1 * 0,00002) -0,1
Refletância TOA (incluindo ajuste da elevação do sol)
((B1 * 0,00002) -0,1) / 0,74457226676389733207607359928648
Ou calcule diretamente do DN bruto para a refletância TOA (incluindo o ajuste da elevação do sol)
((((0,012344 * B1) -61,72016) * 0,00002) -0,1) / 0,74457226676389733207607359928648
Termino com algo assim:
o que parece incorreto.
Eu pesquisei no Google e meus termos de pesquisa são lixo ou ninguém usa matemática da banda manual para fazer isso.
[EDITAR]
A resposta de Mikkel abaixo com o cálculo matemático da banda revisado (omitindo as etapas intermediárias) produziu as seguintes estatísticas:
que parecem um pouco mais realistas.
Pensando se devo adicionar 13,4305 a cada banda ou ignorar o mínimo negativo?
[FIM EDIT]
[EDIT # 2]
Mikkel estava no local novamente, ao sugerir que os valores de -13.4305 foram isolados nas áreas "sem dados" ao redor da imagem;
como você pode ver, isolado apenas na extremidade norte da cena.
[END EDIT # 2]
Aqui está uma boa leitura sobre a física da correção atmosférica e pode ajudar com o problema ... Sou um cara que pressiona o botão, sinto sua dor, pelo menos você pode simplesmente clicar no lindo botão QUAC ou FLAASH quando estiver fora da academia.
http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf
fonte