Landsat8 DN a TOA Reflectancia (menos valores)

8

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 SpectralStats

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:

TOA Reflectance

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:

ingrese la descripción de la imagen aquí

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;

ingrese la descripción de la imagen aquí

como puede ver, aislado solo en el borde norte de la escena.

[EDICIÓN FINAL # 2]

Kyuss
fuente

Respuestas:

7

Un defecto en tu enfoque. No necesita pasar por DN para brillar.
Puede ir directamente al DN a la reflectancia.
Solo quédate ((B1*0.00002)-0.1)/0.74457226676389733207607359928648.

Mikkel Lydholm Rasmussen
fuente
esto funcionó bastante bien! Supongo que tampoco se explicó que la luminosidad TOA / reflectancia TOA excl. la elevación del sol no fueron pasos necesarios en la secuencia de corrección.
Kyuss
según la publicación editada arriba, ¿es necesario agregar 13.4305 a cada banda para llevarlos a min = 0, o es innecesario?
Kyus
@kyuss, intente verificar dónde se encuentran esos -0.134305. Sospecho que son las áreas sin datos alrededor del borde de la imagen.
Mikkel Lydholm Rasmussen
genial, ni siquiera comprobé nada más que las fronteras del sur donde no había valores negativos. Gracias por toda tu ayuda.
kyuss
@kyuss ahora el siguiente paso en su proceso es "simplemente" hacer la corrección atmosférica.
Mikkel Lydholm Rasmussen
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

geodranico
fuente
gracias Steven He aprendido, como has deducido, mucho durante mi estudio, pero ese documento modis puede ayudar a llenar mis vacíos de comprensión con respecto a la transferencia radiativa.
Kyus