Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 13 vistas

    Practica02 - José Paulo Muñoz Catalán

    Cargado por

    JOSÉ PAULO MUÑOZ CATALÁN
    Este documento describe los pasos para estimar la temperatura de la superficie a partir de imágenes Landsat utilizando dos algoritmos (ETR y SW). Se obtienen datos de corrección atmosférica de MODTRAN y se aplican ganancias y offsets a las bandas. Se calcula el NDVI para estimar la emisividad y luego la radiación de superficie para calcular la temperatura con ETR. El mismo proceso se aplica a SW usando dos bandas térmicas. Las temperaturas de ambos métodos se comparan con un producto Landsat corregido y

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Practica02 - José Paulo Muñoz Catalán

    Cargado por

    JOSÉ PAULO MUÑOZ CATALÁN
    13 vistas23 páginas
    Este documento describe los pasos para estimar la temperatura de la superficie a partir de imágenes Landsat utilizando dos algoritmos (ETR y SW). Se obtienen datos de corrección atmosférica de MODTRAN y se aplican ganancias y offsets a las bandas. Se calcula el NDVI para estimar la emisividad y luego la radiación de superficie para calcular la temperatura con ETR. El mismo proceso se aplica a SW usando dos bandas térmicas. Las temperaturas de ambos métodos se comparan con un producto Landsat corregido y

    Título original

    Practica02_José Paulo Muñoz Catalán

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe los pasos para estimar la temperatura de la superficie a partir de imágenes Landsat utilizando dos algoritmos (ETR y SW). Se obtienen datos de corrección atmosférica de MODTRAN y se aplican ganancias y offsets a las bandas. Se calcula el NDVI para estimar la emisividad y luego la radiación de superficie para calcular la temperatura con ETR. El mismo proceso se aplica a SW usando dos bandas térmicas. Las temperaturas de ambos métodos se comparan con un producto Landsat corregido y

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    13 vistas23 páginas

    Practica02 - José Paulo Muñoz Catalán

    Cargado por

    JOSÉ PAULO MUÑOZ CATALÁN
    Este documento describe los pasos para estimar la temperatura de la superficie a partir de imágenes Landsat utilizando dos algoritmos (ETR y SW). Se obtienen datos de corrección atmosférica de MODTRAN y se aplican ganancias y offsets a las bandas. Se calcula el NDVI para estimar la emisividad y luego la radiación de superficie para calcular la temperatura con ETR. El mismo proceso se aplica a SW usando dos bandas térmicas. Las temperaturas de ambos métodos se comparan con un producto Landsat corregido y

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    Está en la página 1de 23

    PRACTICA 02

    Con los datos de corrección atmosférica obtenidos desde MODTRAN de la practica 01 para banda 10 y
    banda 11 realizaremos procesos para estimar la temperatura a partir de dos algoritmos (ETR y SW)
    con el objetivo final de comparar los resultados con un producto corregido de manera oficial por
    Landsat.

    Utilizamores imagenes Landsat:

    Las imágenes corresponden a dos niveles L1(radiancias capturadas por sensor) y L2 (corregidas
    atmosféricamente)
    Importamos las imágenes:

    y aplicamos un recorte para reducir la información y focalizar el análisis a:


    A lo anterior debemos aplicarle un descodificador para cada una de las bandas, para ello usamos el
    MTL de cada fichero para obtener valores asociados a Gain(multiplicador) y Offset(suma o resta) en
    relación a esta fórmula base :

    gain*BANDA_LANDSAT+offset

    Las fórmulas corresponderían a:

    L1 B10=3.8000E-04*B10+0,1
    L1 B11=3.4900E-04*B11+0,1
    L2 B4=2.75e-05*B4-0,2
    L2 B5=2.75e-05*B5-0,2
    L2 ST B10=0.00341802*B10+149

    Para B4 y B5 obtendremos reflectividad a nivel de superficie. Para B10 y B11 radiancias


    detectadas por sensor. Para ST B10 temperatura a nivel de superficie.
    Con la herramienta Pixel Value podemos evaluar la magnitud de las radiancias, estas se encuentran
    entre 5 y 15 por lo tanto validamos la descodificación, adicionalmente los valores de B10 son mayores
    a B11.

    Estimación de emisividad:

    Existen variadas formas de estimar la emisividad y una de ellas es partir del estudio de la vegetación.
    El NDVI nos hace inferir que proporciones altas de vegetación en una zona se relaciona directamente
    con valores alto de emisividad. Con los valores de NDVI podemos obtener el proporción o fracción de
    vegetación de la zona de interés :

    PV: NDVI-NDVIs/NDVIv-NDVIs

    Por lo tanto, con las bandas B4 y B5 se efectuara el calculo de un incide normalizado para determinar
    estado de la vegetación y así la proporción de estas en la zona de interés a partir de los valores
    extremos asociados a pixeles puros de suelo desnudo y pixeles puros de vegetación.

    NDVI: B4-B5/B4+B5
    Cálculo de NDVI y validación de rangos de reflectividad(valores entre 0 y 1):

    Para obtener los valores extremos (NDVIs y NDVIv) crearemos un roi sobre una zona de alto contraste
    que reflejaría los valores bajos y altos.
    La frecuencia más alta de valores cercano al 0 me indica suelo desnudo y los valores de frecuencias
    más altos cercanos al 1 me indican vegetación por lo tanto estos valores quedan indicados como:

    NDVIs : 0,1733 ; NDVIv: 0,874

    Con los valores obtenidos se hace posible calcular la proporción de vegetación, para ello aplicamos la
    siguiente formula :

    PV: (NDVI-0,17)/float(0,87-0,17)
    Los valores obtenidos de PV se deben normalizar para agrupar los valores en un rango de 0 y 1.Por lo
    tanto todos los valores negativos se agruparan en 0 y los valores sobre 1 quedan agrupados en dicho
    valor.

    PV <1 y PV >0

    Las imágenes obtenidas nos permiten aplicar las ecuaciones para obtener valores de emisividad,para
    ello aplicamos las siguientes condiciones por cada banda:

    Primera condición: NDVI<NDVIs


    b10 = 0.979-0.046 b4

    Formula a software: (NDVI lt 0.1733)*(0.979-0.046*B4)

    Los valores obtenidos deberían estar entre un rango de 0,9 y 1.


    Segunda Condición: NDVIs<NDVI<NDVIv :
    b10 = 0.987 Pv + 0.971 (1-Pv)
    Formula a software: (NDVI ge 0.1733)*(0.987* PV + 0.971* (1-PV))
    estas bandas son complementarias por lo tanto se suman:

    La emisividad total de la Banda10 nos permite aplicar la formula para obtener la energía real que emite
    la superficie (radiancia de la superficie):

    Con los valores obtenidos de partir de MODTRAN, aplicaremos la formula:


    Transmisividad 0,811
    Lup 1,56772568664247
    Ldown 2,37023593143593

    Fórmula para software: (b1-1.56)-(1-b2)*2.37*0.811/0.811*b2

    Donde b1=Lsen(radiancia medida por el sensor);b2=emisividad total.

    Con lo anterior podremos calcular Ts:

    Fórmula para software : (14388/10.9)*(alog((119110000/(b1*10.9^5))+1))^(-1)


    Donde b1= B(T).
    Comparamos con la L2:

    T(B)10 L2SP B10


    Para B11 realizaremos el mismo proceso aplicando las condiciones para obtener la emisividad total que
    nos permitirá calcular B(T)11 y Ts.

    NDVI<NDVIs
    b11 = 0.982-0.027 b4
    Fórmula para software: (b1 lt 0.1733)*(0.982-0.027*b2)

    Segunda condición: NDVIs<NDVI<NDVIv :

    Fórmula para software: (b1 ge 0.1733)*(0.989* b2 + 0.977* (1-b2) )


    Donde b1=NDVI;b2=PV
    ahora sumamos ambas condiciones:
    Con los valores:

    Transmisividad 0,807
    Lup 1,71690463874642
    Ldown 2,33069849294152

    calculamos B(T) para B11


    Fórmula para software: (b1-1.71)-(1-b2)*2.33*0.807/0.807*b2

    Luego aplicamos : (14388/12,0)*(alog((119110000/(b1*12,0^5))+1))^(-1) para obtener Ts


    T(B)10 T(B)11
    Comparamos la imagen Ts10 y L2SP 10:

    Los datos del Scatter Plot los llevaremos a Excel para obtener las diferencias en grados Celsius y
    valores estadísticos:

    Promedio 1.45730892
    Desviación 1.11095136
    RMS 1.83247434

    Los valores de promedio y desviación son bajos (menores a 2) por lo tanto el modelo es preciso exacto
    en relación al valor verdadero.
    Ts 10 v/s L2SP10
    55

    50

    45

    40

    35

    30

    25

    20
    0 10 20 30 40 50 60

    Cálculo de SW

    Para calcular temperaturas a partir de Split Window debemos tener dos bandas.Para ellos obtendremos
    temperaturas de brillo para las bandas B10 y B11 aplicando las fórmulas:

    TbrilloB10: (14387.7/12.001)*(alog((119110000/(float(b1)*12.001^5))+1))^(-1)

    TbrilloB11: (14387.7/10.9)*(alog((119110000/(float(b1)*10.9^5))+1))^(-1)
    TbrilloB10 TbrilloB11
    Se cumple que TsB10 es mayor TsB11,lo anterior ya que la longitud de onda en 10 um recibe más
    energía.
    Luego split window:

    Ts=b1+(1.378*(b1-b2))+(0.183*((b1-b2)^2))−0.268 +((54.30−2.238*2.22)*(1-
    b3))+((−129.20+16.40*2.2)*b4)

    Donde:
    b1=Banda10
    b2=banda11
    W:2,22
    b3= promedio de emisividad
    b4: diferencia de emisividad de b10-b11
    Ahora comparamos Ts SW y L2SP10
    Promedio -0.19392432
    Desviación 1.78606512
    RMS 1.79656206

    Los valores de promedio y desviación son bajos (menores a 2) por lo tanto el modelo es preciso exacto
    en relación al valor verdadero.

    Ts SW v/s L2SP10
    55

    50

    45

    40

    35

    30

    25

    20
    0 10 20 30 40 50 60
    Comparación ETS y SW:

    Promedio 1.53733606
    Desviación 1.65794197
    RMS 2.26101166

    Los valores de promedio y desviación son bajos (menores a 2) por lo tanto el modelo es preciso exacto
    en relación al valor verdadero.

    ETS vs SW
    55

    50

    45

    40

    35

    30

    25

    20
    0 10 20 30 40 50 60

    También podría gustarte