Estoy aprendiendo cómo realizar la corrección atmosférica manualmente, por lo tanto, no estoy usando los métodos automatizados disponibles en ENVI
los metadatos (banda 1) para mi escena Landsat 8 son los siguientes;
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
SEL_ELEVATION = 48.12237303
Entonces, cuando uso matemática de banda para calcular TOA Radiance, escribo esto:
(0.012344 * B1) -61.72016
Y produce estadísticas espectrales que parecen correctas
Cuando uso cualquiera de las siguientes fórmulas matemáticas de banda;
Reflectancia TOA (excluido el ajuste de elevación del sol)
(B1 * 0.00002) -0.1
Reflectancia TOA (incluido el ajuste de elevación del sol)
((B1 * 0.00002) -0.1) /0.74457226676389733207607359928648
O calcule directamente desde el DN bruto hasta la reflectancia TOA (incluido el ajuste de elevación del sol)
((((0.012344 * B1) -61.72016) * 0.00002) -0.1) /0.74457226676389733207607359928648
Termino con algo como esto:
que se ve incorrecto
Busqué en Google y mis términos de búsqueda son basura o nadie usa matemática de banda manual para hacer esto.
[EDITAR]
La respuesta de Mikkel a continuación con el cálculo matemático de banda revisado (omitiendo los pasos intermedios) arrojó las siguientes estadísticas:
que parecen algo más realistas.
¿Se pregunta si debería agregar 13.4305 a cada banda o ignorar el mínimo negativo?
[FIN EDITAR]
[EDITAR # 2]
Mikkel volvió a acertar al sugerir que los valores de -13.4305 se aislaron en las áreas "sin datos" alrededor de la imagen;
como puede ver, aislado solo en el borde norte de la escena.
[EDICIÓN FINAL # 2]
Aquí hay una buena lectura sobre la física de la corrección atmosférica, y podría ayudar con el problema ... Soy un tipo que presiona un botón, siento tu dolor, al menos puedes hacer clic en ese hermoso botón QUAC o FLAASH una vez que estás afuera de la academia
http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf
fuente