[doc sps="1.8" acron="rge" jtitle="Revista Geográfica de América Central" stitle="Rev. Geog. Amér. Central" issn="1011-484X" pissn="1011-484X" eissn="2215-2563" pubname="Universidad Nacional, Costa Rica" license="http://creativecommons.org/licenses/by-nc-sa/4.0" volid="65" issueno="2" dateiso="20201200" season="Jul/Dec" order="01" fpage="17" lpage="32" pagcount="17" doctopic="oa" language="es"][doi]10.15359/rgac.65-2.1[/doi]
[toctitle]Artículo[/toctitle]


[doctitle language="es"]Estimación de la temperatura superficial en imágenes históricas Landsat 5 mediante corrección atmosférica monocanal en el térmico para la cuenca del río Tempisque-Bebedero (Costa Rica)[/doctitle]
[doctitle language="en"]Estimation of the surface temperature in Landsat 5 historical images by means of single-channel atmospheric correction in the thermal for the Tempisque-Bebedero river basin (Costa Rica)[/doctitle]

[author role="nd"][surname]Martínez-Barbáchano[/surname] [fname]RUBÉN[/fname][xref ref-type="aff" rid="aff1"]1[/xref][authorid authidtp="orcid"]0000-0002-9176-3112[/authorid][/author]

[normaff id="aff1" ncountry="Costa Rica" norgname="Universidad de Costa Rica" icountry="CR"][label][sup]1[/sup][/label][role]Máster[/role]. [orgname]Universidad de Costa Rica[/orgname], [country]Costa Rica[/country]. Correo electrónico: [email]ruben.martinezbarbachano@ucr.ac.cr[/email][/normaff] 

[xmlabstr language="es"][sectitle]Resumen:[/sectitle]
[p]En las últimas décadas la utilización de imágenes satelitales se ha incrementado al punto de constituir una herramienta fundamental en el monitoreo de procesos ambientales y de cambio global. En tal sentido, la importancia histórica de la familia de satélites Landsat marca un hito en 1985 con la introducción del sensor Tematic Mapper, que ofrecía la posibilidad a los usuarios de adquirir información de una base histórica de más de 25 años de imágenes térmicas a 120 metros de resolución espacial. El presente trabajo se presenta como una metodología de obtención de temperatura superficial a partir de una serie temporal de imágenes Landsat-5 TM adquiridas en el año 2001 sobre el Pacífico norte de Costa Rica. Se presentan los distintos algoritmos que pueden utilizarse para la estimación de temperatura superficial a partir de datos Landsat con el método monocanal por la función inversa de Planck. Aspectos como la transmisividad media de la atmósfera, la radiancia ascendente y descendente, o la emisividad de la superficie terrestre fueron tomados en cuenta, previamente al cálculo de la temperatura de brillo de las imágenes.[/p][/xmlabstr]
[kwdgrp language="es"][sectitle]Palabras clave:[/sectitle] [kwd]Teledetección[/kwd], [kwd]infrarrojo térmico[/kwd], [kwd]emisividad[/kwd], [kwd]temperatura superficial[/kwd][/kwdgrp].

[xmlabstr language="en"][sectitle]Abstract:[/sectitle]
[p]In recent decades, the use of satellite imagery has increased to the point of being a fundamental tool in the monitoring of environmental processes and global change.  In this sense, the historical importance of the Landsat family of satellites marks a milestone in 1985 with the introduction of the Tematic Mapper sensor, which offered users the possibility of acquiring information from an over 25-year historical database of thermal imagery at 120 meters of spatial resolution.  This paper is presented as a methodology for obtaining surface temperature data from a time series of Landsat-5 TM images acquired in 2001 on the North Pacific Region of Costa Rica.  The different algorithms that can be used for the estimation of surface temperature from Landsat data via the single-channel method by way of the inverse Planck function are presented.  Prior to the calculation of the brightness temperature of the images, aspects such as average atmospheric transmissivity, ascending and descending radiance or earth´s surface emissivity were taken into account.[/p][/xmlabstr]
[kwdgrp language="en"][sectitle]Keywords:[/sectitle] [kwd]Remote sensing[/kwd]; [kwd]Thermal infrared[/kwd]; [kwd]Emissivity[/kwd]; [kwd]Surface temperature[/kwd][/kwdgrp].

[hist]Recibido: [received dateiso="20190406"]6/4/2019[/received]
Aceptado: [accepted dateiso="20191029"]29/10/2019[/accepted][/hist]


[xmlbody][sec sec-type="intro"][sectitle]Introducción[/sectitle]


[p]La teledetección espacial ha experimentado un renovado impulso en las últimas décadas, con nuevos sensores y plataformas que han permitido la obtención de magnitudes o variables geofísicas de enorme valor geográfico. Dentro del amplio abanico de dichas variables, la temperatura terrestre es el principal objeto de estudio de la teledetección térmica. De hecho, la temperatura superficial es una fuente fundamental de información, tanto cualitativa como cuantitativa, acerca de los procesos que ocurren en la superficie terrestre, permitiendo por tanto su caracterización, análisis y modelización ([xref  ref-type="bibr" rid="r15"]Quattrochi y Luvall, 2000[/xref]).[/p]
[p]Las imágenes en infrarrojo térmico de media resolución permiten, además, estimar otras variables como la evapotranspiración a partir de la ecuación de balance energético, de gran utilidad en la gestión del recurso hídrico.[/p]
[p]La puesta en órbita en 1982 del sensor Thematic Mapper (TM) a bordo de la plataforma Landsat-4 y posteriormente (en 1984) con la plataforma Landsat-5, ha permitido el acceso a una base histórica de más de 25 años de imágenes térmicas. El lanzamiento en 1999 de la plataforma Landsat-7 con el sensor Enhanced Thematic Mapper plus (ETM+) permitió el acceso a imágenes térmicas con una resolución de 60 metros, aunque un fallo técnico del sensor limitó severamente su uso, por lo que fue sustituido en 2013 por los instrumentos OLI y TIRS de la plataforma Landsat 8.[/p]


[/sec][sec][sectitle]Marco conceptual[/sectitle]


[p]Teledetección térmica es un término que hace referencia a las medidas obtenidas por un sensor ubicado a bordo de una plataforma (un avión o satélite) que capta la radiación electromagnética emitida por la superficie terrestre en la región del infrarrojo térmico, típicamente entre los 8 y 14 µm ([xref  ref-type="bibr" rid="r16"]Sobrino, et al. 2000[/xref]).[/p]
[p]En el intervalo espectral comprendido entre 8 y 14 µm (infrarrojo térmico) se produce el 80 % de la emisividad de la superficie terrestre. En dicho intervalo es donde, además, la transmisividad atmosférica es más elevada y por tanto, donde se ubican las medidas que diferentes sensores toman a bordo de las misiones de los satélites ([xref  ref-type="bibr" rid="r13"]Pérez y Muñoz, 2002[/xref]). [/p]
[p]La atenuación de la señal en dicha región del espectro se debe a la absorción de aerosoles, pero se trata de una cantidad despreciable debido a que las concentraciones de vapor de agua (H2O) ozono (O3) y dióxido de carbono (CO2) se encuentran en menor proporción ([xref ref-type="bibr" rid="r14"]Prata et al., 1994[/xref]).[/p]
[p]Sin embargo, la transmisión atmosférica de la radiación se ve afectada por el vapor de agua cinco veces más que por el resto de los gases juntos, absorbiendo la radiación de onda larga del infrarrojo térmico y onda corta en la región de las microondas (desde los 22 µm a 1 mm). En cuanto al dióxido de carbono, este absorbe fuertemente la radiación en la región del infrarrojo lejano del espectro electromagnético, mientras que el ozono lo hace para longitudes de ondas más cortas y energéticas ([xref  ref-type="bibr" rid="r1"]Atkins & Jones, 2006[/xref]).[/p]
[p]En consecuencia, la estimación de la temperatura superficial está condicionada, principalmente, por el contenido de vapor de agua de la atmósfera (especialmente sensible en la zona intertropical) y la emisividad de la superficie estudiada. Por todo ello, una correcta estimación de la temperatura superficial implica necesariamente la consideración y corrección de estos factores.[/p]
[p]La radiación emitida por una superficie a una determinada longitud de onda, dentro del rango infrarrojo térmico, es función no sólo de su temperatura, sino que depende también de su emisividad. Es por ello, un parámetro crítico para el estudio del balance radiativo terrestre. La ecuación de la emisividad está representada como la energía emitida por una superficie eficiente ([xref  ref-type="bibr" rid="r18"]Valor y Caselles, 1996[/xref]) y viene dada por:[/p]

[equation id="e1"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 1. Ecuación de la emisividad ([xref  ref-type="bibr" rid="r18"]Valor y Caselles, 1996[/xref])[/p]

[p]Siendo ei la emisividad de una superficie (considerando la transmisividad e irradiancia de la atmósfera), Ri(Ts) la radiancia espectral emitida por un cuerpo a una determinada longitud de onda y temperatura, Bi(Ts) la radiancia emitida por un cuerpo negro a una temperatura dada, medida en campo, y por último Li 8 la irradiancia atmosférica.[/p]
[p]Se suele definir un cuerpo negro como un emisor perfecto, es decir, con una emisividad de 1; sin embargo, la vegetación tiene valores de emisividad menores, entre 0,97 - 0,985, debido a que parte de la energía incidente es reflejada, absorbida y transmitida en proporciones diferentes, según la longitud de onda (Rubio et al., 2003). Su valor es altamente dependiente de la composición de la superficie, sobre todo con su contenido de humedad, rugosidad y tamaño de sus partículas ([xref  ref-type="bibr" rid="r3"]Barreto, et al. 2009[/xref]).[/p]
[p]También varía con la longitud de onda en la cual es medida, la resolución y ángulo de observación, dando una idea de la complejidad que representa su estimación. Algunas de las técnicas que se han utilizado para estimar la emisividad de la superficie son: [/p]

- [list listtype="bullet"][li]Método de emisividad normalizada (NEM) propuesto por [xref  ref-type="bibr" rid="r10"]Gillespie et al., (1985[/xref]).[/li]
- [li]El índice Independiente de la Temperatura (TISI) propuesta por [xref  ref-type="bibr" rid="r4"]Becker y Li, (1990[/xref]).[/li]
- [li]La separación de emisividad (TES) propuesta por [xref  ref-type="bibr" rid="r11"]Gillespie et al. (1999[/xref]).[/li]
- [li]El método basado en índice de vegetación (NDVI) propuesto por [xref  ref-type="bibr" rid="r18"]Valor y Caselles, (1996[/xref]) que fue el utilizado en el presente estudio.[/li][/list]

[p]El canal térmico (6 para Landsat TM-ETM+) convertido a radiancia se puede transformar a temperatura de brillo efectiva por medio de la ecuación de cálculo de la temperatura superficial a partir de los valores de brillo de la imagen térmica, que requiere de dos constantes de calibración como muestra el [xref ref-type="table" rid="t1"]cuadro 1[/xref]:[/p]

[tabwrap id="t1"][label]Cuadro 1[/label]. [caption]Constantes de calibración de la banda térmica del satélite Landsat[/caption]
[table][tbody][tr][td align="center"]Satélites[/td]
[td align="center"]Bandas[/td]
[td align="center"]K1[/td]
[td align="center"]K2[/td]
[/tr][tr][td align="center" colspan="2"] 
 [/td]
[td align="center"](W/m-2/sr-1/µm-1)[/td]
[td align="center"]Kelvin[/td]
[/tr][tr][td align="center"]Landsat 4[/td]
[td align="center"]6[/td]
[td align="center"]671,62[/td]
[td align="center"]1284,30[/td]
[/tr][tr][td align="center"]Landsat 5[/td]
[td align="center"]6[/td]
[td align="center"]607,76[/td]
[td align="center"]1260,56[/td]
[/tr][tr][td align="center"]Landsat 7[/td]
[td align="center"]6[/td]
[td align="center"]666,09[/td]
[td align="center"]1282,71[/td]
[/tr][tr][td align="center"]Landsat 8[/td]
[td align="center"]10[/td]
[td align="center"]774,89[/td]
[td align="center"]1321,08[/td]
[/tr][tr][td align="center"] [/td]
[td align="center"]11[/td]
[td align="center"]480,89[/td]
[td align="center"]1201,14[/td][/tr][/tbody][/table]
[fntable id="TFN1"]Fuente: [xref  ref-type="bibr" rid="r7"]Chander et al., 2009[/xref].[/fntable][/tabwrap]

[p]Por otra parte, el uso de la banda térmica en el estudio de la temperatura superficial de la Tierra requiere de la eliminación de los efectos atmosféricos en dicha longitud de onda. Además, la emisividad de un objeto en el terreno es igualmente atenuada y aumentada por la atmósfera. Por todo ello, resulta imprescindible un conocimiento apropiado de la atmósfera para la estimación de la transmitancia, la radiancia ascendente y la radiancia descendente. Cuando se conocen los valores de estas 3 variables, se puede convertir la radiancia medida en el espacio a la radiancia emitida por una superficie ([xref ref-type="bibr" rid="r2"]Barsi et al., 2003[/xref]).[/p]


[/sec][sec][sectitle]Área geográfica de estudio[/sectitle]


[p]La cuenca del río Tempisque-Bebedero, en el noroeste de Costa Rica, cubre en total 5455 km² en la provincia de Guanacaste. La cuenca se encuentra situada 10º 58' 25'' de latitud norte, 10º 06' 41'' de latitud sur, 85º 46' 59'' oeste y 84º 49' 40'' este. El clima predominante en la cuenca es de tipo tropical con una prolongada estación seca, que va de finales de noviembre a principios de mayo. La estación de lluvias comienza en mayo y se prolonga hasta noviembre. Existe por tanto una fuerte estacionalidad en las precipitaciones, con un valor promedio anual de 1900 mm. En cuanto a las temperaturas, no existe la marcada estacionalidad que se da en las precipitaciones, con un promedio anual de 28 ºC. Sin embargo, es habitual en plena estación seca que la temperatura promedio ascienda hasta los 34 ºC de promedio, con máximos puntuales de hasta 40 ºC. Desde un punto de vista geomorfológico se trata de una zona de origen predominantemente sedimentario, y por ello reciente. El complejo Nicoya configura la red de drenaje en virtud de colinas calizas que constituyen un espacio formado fundamentalmente a partir del Mioceno ([xref  ref-type="bibr" rid="r9"]Denyer, Montero & Flores, 2014[/xref]). La existencia de espacios abiertos, llanuras inundables, cultivos y pastizales facilita la obtención de temperatura superficial del suelo, tanto por la homogeneidad de los pixeles obtenidos por el satélite como por la existencia de suelo sin cobertura forestal densa.[/p]


[/sec][sec sec-type="materials|methods"][sectitle]Materiales y métodos[/sectitle]


[p]Para el presente trabajo se utilizaron exclusivamente recursos libres online y software de código abierto. Todo el procesamiento de imágenes que sustenta la presente metodología se realizó con la herramienta "calculadora de mapas" de GvSIG 2.4 ([xref ref-type="fig" rid="f1"]figura 1[/xref]).[/p]

[figgrp id="f1"][graphic href="?a01n65"][/graphic]
[label]Figura 1[/label]. [caption]Calculadora de mapas de GvSIG[/caption][/figgrp]

[p]El problema fundamental con la estimación de la temperatura superficial de la Tierra a partir de datos de satélite es la necesidad de corregir en los mismos, los efectos debidos a la absorción atmosférica y a la emisividad de la superficie medida. Para la producción del mapa de emisividad, por tanto, se parte del modelo presentado por ([xref  ref-type="bibr" rid="r18"]Valor y Caselles, 1996[/xref]). Pero previamente, es conveniente la realización de una corrección DOS1 a las bandas roja e infrarroja de Landsat5, para remover el efecto de la dispersión atmosférica.[/p]
[p]Seguidamente, calculamos el índice de vegetación diferencial normalizado, o NDVI (como se muestra en la [xref ref-type="fig" rid="f2"]figura 2[/xref]) a partir de las bandas 3 y 4, ya corregidas atmosféricamente. La fórmula del NDVI es:[/p]

[equation id="e2"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 2. Índice de vegetación diferencial normalizado[/p]

[figgrp id="f2"][graphic href="?a01n65"][/graphic]
[label]Figura 2[/label]. [caption]Índice de Vegetación Diferencial Normalizado o NDVI de la cuenca del Río Tempisque-Bebedero, Costa Rica[/caption][/figgrp]

[p]Una vez obtenido el NDVI, se puede utilizar este directamente para calcular la emisividad de los pixeles, o bien se puede generar un PV (Proportion of Vegetation). La ventaja de esta segunda opción es que permite aislar la emisividad de las diferentes cubiertas, tomando como referencia la de la vegetación que puede ser calculada a partir de este índice ([xref  ref-type="bibr" rid="r6"]Carlson & Ripley, 1997[/xref]) mostrado en la ecuación 3.[/p]

[equation id="e3"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 3. Proporción de vegetación[/p]
[p]Fuente: [xref  ref-type="bibr" rid="r6"]Carlson & Ripley, 1997[/xref].[/p]

[figgrp id="f3"][graphic href="?a01n65"][/graphic]
[label]Figura 3[/label]. [caption]Índice de Proporción de Vegetación o PV, de la cuenca del Río Tempisque-Bebedero, Costa Rica[/caption][/figgrp]

[p]En la [xref ref-type="fig" rid="f3"]figura 3[/xref] se puede constatar que la cantidad de superficie cubierta por vegetación es sensiblemente menor en el PVI que en el NDVI.[/p]
[p]Seguidamente, se procede a mapear la emisividad de los pixeles a partir de la vegetación, como se muestra en la [xref ref-type="fig" rid="f4"]figura 4[/xref]. Asumiendo que el valor típico de la emisividad en el térmico es de 0.99 para la vegetación, escoger un valor típico de emisividad para el suelo puede ser algo más complejo. Una posibilidad es tomar como referencia el valor propuesto a partir de las 49 muestras de la biblioteca espectral de Aster (http://asterweb.jpl.nasa.gov). En tal caso el valor medio es de 0.973, con una desviación estándar de 0.004 y la expresión final para calcular la emisividad en el monocanal térmico sería la correspondiente con la ecuación 4:[/p]

[equation id="e4"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 4. Cálculo de emisividad a partir de la proporción de vegetación o Pv[/p]
[p]Fuente: [xref  ref-type="bibr" rid="r17"]Sobrino, Jiménez-Muñoz & Paolini, 2004[/xref].[/p]

[figgrp id="f4"][graphic href="?a01n65"][/graphic]
[label]Figura 4[/label]. [caption]Mapa de emisividad de la cuenca del Río Tempisque-Bebedero, Costa Rica[/caption][/figgrp]

[p]Se va a tomar como referencia el valor de emisividad del pixel correspondiente a la estación meteorológica de Liberia ([xref ref-type="fig" rid="f5"]figura 5[/xref]) con coordenadas UTM 16 N 1172043 X / 660453 Y[/p]

[figgrp id="f5"][graphic href="?a01n65"][/graphic]
[label]Figura 5[/label]. [caption]El valor de emisividad para el pixel de la estación de Liberia es de 0.987321[/caption][/figgrp]

[p]Se calculó la temperatura superficial a partir del método monocanal, basado en la ecuación de transferencia radiativa de [xref  ref-type="bibr" rid="r8"]Coll, Galve, Sánchez y Caselles (2010[/xref]). Para ello, se tuvieron en cuenta los aportes de [xref  ref-type="bibr" rid="r18"]Valor y Caselles (1996[/xref]) y [xref  ref-type="bibr" rid="r7"]Chander et al. (2009[/xref]). [/p]
[p]En primer lugar, se realizó una calibración radiométrica para llevar los valores digitales a radiancia espectral (W/m-2/sr-1/µm-1) como muestra la [xref ref-type="fig" rid="f6"]figura 6[/xref]. La calibración radiométrica o transformación a valores de radiancia es un paso esencial para poder comparar datos o información de múltiples sensores y plataformas en una escala radiométrica común.[/p]
[p]La radiancia se calcula de acuerdo a los metadatos del archivo MTL de la banda 6, a través de la ecuación 5:[/p]

[equation id="e5"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 5. Cálculo de radiancia[/p]

[p]Donde L? es la radiancia espectral de la banda térmica que va a ser calculada; LMAX y LMIN son los valores de radiancia máxima y mínima para la banda 6; QCALMAX y QCALMIN son los valores máximo y mínimo de la radiancia calibrada, es decir, modificados a escala en unidades de números digitales (en función de la resolución radiométrica del sensor) y DN son los niveles digitales (la banda 6 en este caso).[/p]

[figgrp id="f6"][graphic href="?a01n65"][/graphic]
[label]Figura 6[/label]. [caption]Radiancia espectral de la cuenca del Río Tempisque-Bebedero, Costa Rica[/caption][/figgrp]

[p]A diferencia de imágenes anteriores, los valores de radiancia más elevados (en verde) se dan en zonas de escasa cobertura vegetal.[/p]
[p]Para calcular la temperatura superficial se consideró la siguiente ecuación de transferencia radiativa de C. [xref  ref-type="bibr" rid="r8"]Coll, J. M. Galve, J. M. Sánchez y V. Caselles (2010[/xref]):[/p]

[equation id="e6"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 6. Ecuación de transferencia radiativa[/p]
[p]Fuente: [xref  ref-type="bibr" rid="r8"]Coll, Galve, Sánchez y Caselles, 2010[/xref].[/p]

[p]Donde; Lsen es la radiancia medida por el sensor (W/m-2/sr-1/µm-1), e es la emisividad de la superficie, B es el parámetro derivado de la Ley de Planck calculado a partir de la ecuación de cálculo de temperatura de la superficie terrestre (°K), L?: radiancia atmosférica descendente (W/m-2/sr-1/µm-1), t: transmisividad atmosférica y L?: radiancia atmosférica ascendente (W/m-2/sr-1/µm-1).[/p]
[p]La transmisividad de la atmósfera, y la radiancia ascendente y descendente no son parámetros presentes en los metadatos de una imagen satelital. Para su obtención se recurrió a una herramienta operativa en línea de corrección atmosférica que permite calcular la temperatura superficial con un sesgo inferior a 0.5 ± 0.8ºK ([xref ref-type="bibr" rid="r2"]Barsi y otros, 2003[/xref]). La herramienta se denomina Atmospheric Correction Parameter Calculator y está disponible en el sitio web; https://atmcorr.gsfc.nasa.gov/. [/p]

[figgrp id="f7"][graphic href="?a01n65"][/graphic]
[label]Figura 7[/label]. [caption]Para la fecha de la escena Landsat (6 de febrero de 2001) se obtuvieron unos valores de transmisividad atmosférica de 0,54, una radiancia ascendente de 3,66 (W/m-2/sr-1/µm-1) y una radiancia descendente de 5,50 (W/m-2/sr-1/µm-1).[/caption][/figgrp]
[p]Fuente: https://atmcorr.gsfc.nasa.gov/[/p]

[p]Barsi y otros (2003) afirman que los valores de temperatura superficial pueden ser estimados con un error de ±2ºK cuando la atmósfera está relativamente despejada y se conoce la emisividad de la superficie. Esta herramienta utiliza perfiles atmosféricos del National Center for Environmental Prediction (NCEP) interpolados a una posición, fecha y hora dadas y el modelo MODTRAN 4 ([xref  ref-type="bibr" rid="r5"]Berk et al. 1999[/xref]).[/p]
[p]Seguidamente, se aplicó la ecuación 7 de transferencia radiativa propuesta por [xref  ref-type="bibr" rid="r8"]Coll, Galve, Sánchez y Caselles (2010[/xref]). El error que introduce esta ecuación se incrementa con el ancho de longitud de onda, la humedad atmosférica y las variaciones de emisividad de la superficie ([xref ref-type="bibr" rid="r8"]Coll et al, 2010[/xref]).[/p]

[p]El siguiente paso consistió en aplicar la ecuación 7 de cálculo de temperatura de brillo:[/p]

[equation id="e7"][graphic href="?a01n65"][/graphic][/equation]
[p]Ecuación 7. Cálculo de la temperatura a partir de los valores de brillo de la imagen ([xref  ref-type="bibr" rid="r7"]Chander, 2009[/xref])[/p]

[p]Donde k2 y k1 son los coeficientes de calibración, con valores de 1260,56 y 607,76 (W/m-2/sr-1/µm-1), respectivamente, ([xref  ref-type="bibr" rid="r7"]Chander et. al. 2009[/xref]) mientras que Lsen es el ráster de radiancia corregida obtenido en el paso anterior. El resultado se obtiene en grados kelvin, que fueron transformados en Celsius.[/p]

[figgrp id="f8"][graphic href="?a01n65"][/graphic]
[label]Figura 8[/label]. [caption]Temperatura superficial en grados Celsius de la cuenca del río Tempisque-Bebedero Costa Rica[/caption][/figgrp]


[/sec][sec sec-type="conclusions"][sectitle]Conclusiones[/sectitle]


[p]La temperatura finalmente obtenida para el aeropuerto de Liberia fue de 44,379 ºC a las 15:40 del 6 de febrero de 2001. De acuerdo con la información proporcionada por el Instituto Meteorológico Nacional de Costa Rica, la temperatura máxima en el aeropuerto de Liberia ese día fue de 34,1ºC. La obtención de la temperatura sin tomar en cuenta parámetros como la emisividad, la transmisividad y la radiancia bidireccional (es decir, aplicando la ecuación 8 directamente sobre la radiancia) genera valores de temperatura más bajos, y por tanto más parecidos a la temperatura del aire a la sombra y a cierta distancia del suelo. Concretamente, en el caso del pixel correspondiente a la estación meteorológica del aeropuerto de Liberia, la temperatura obtenida con este método fue de 32ºC. Considerando la presencia de suelo asfaltado y el hecho de que la temperatura superficial puede superar ampliamente a la temperatura medida a la sombra y a un metro del suelo, la diferencia de 12ºC parece más que plausible. Por otra parte, cuanto menor es la irradiancia solar y mayor es la velocidad del viento, menores diferencias existen entre estas temperaturas, y por ello más próxima es la temperatura del aire a la temperatura superficial medida desde satélite. Además, estas diferencias también disminuyen a mayor proporción de cobertura vegetal y cuando no existe estrés hídrico ([xref  ref-type="bibr" rid="r12"]Niclós, et al, 2010[/xref]). [/p]
[p]La existencia de correlaciones entre temperatura del aire y temperatura superficial en diferentes tipos de coberturas de la zona intertropical y en zonas con diferentes relieves debe ser estudiada con mayor profundidad. También se requiere analizar y validar la introducción de nuevos parámetros en las ecuaciones propuestas, con objeto de obtener un algoritmo general que permita estimar la temperatura del aire a partir de la observación desde satélite.[/p]
[p]La principal desventaja que presenta esta metodología deriva del hecho de que el valor de emisividad que se incorpora a la ecuación 5 proviene de un pixel del terreno (no de un valor promedio de la imagen) por lo que la estimación del valor de temperatura superficial de la imagen sólo tendría una precisión de ±2ºK en el pixel del cual se tomó el respectivo valor de emisividad.[/p]
[p]Finalmente, la utilidad de esta metodología monocanal de obtención de temperatura superficial permite aproximaciones históricas (como la presente, del año 2001) muy precisas acerca de cómo contribuían las coberturas artificiales en la formación de islas de calor urbano en el pasado.[/p][/sec][/xmlbody]


[refs][sectitle]Referencias[/sectitle]


[ref id="r1" reftype="book"][authors role="nd"][pauthor][surname]Atkins,[/surname] [fname]P.[/fname][/pauthor] & [pauthor][surname]Jones,[/surname], [fname]L.[/fname][/pauthor][/authors] ([date dateiso="20060000" specyear="2006"]2006[/date]). [source]Principios de Química. Los caminos del descubrimiento[/source]. [edition]Ed. 5ta[/edition]). [publoc]Madrid. España[/publoc]: [pubname]Médica Panamericana[/pubname]. ISBN: [isbn]978-950-06-0080-4[/isbn].[/ref]
[ref id="r2" reftype="journal"][authors role="nd"][pauthor][surname]Barsi J.[/surname], [fname]Barker J. L.[/fname][/pauthor] & [pauthor][surname]Schott[/surname], [fname]J. R.[/fname][/pauthor][/authors] ([date dateiso="20030000" specyear="2003"]2003[/date]). [arttitle]An Atmospheric Correction Parameter Calculator for a Single Thermal Band Earth-Sensing Instrument[/arttitle]. [source]International Geoscience and Remote Sensing Symposium (IGARSS)[/source]. [volid]5[/volid]. [pages]3014 - 3016[/pages] vol.5. [moreinfo]10.1109/IGARSS.2003.1294665.[/moreinfo][/ref]
[ref id="r3" reftype="journal"][authors role="nd"][pauthor][surname]Barreto[/surname], [fname]Á.[/fname][/pauthor], [pauthor][surname]Arbelo[/surname], [fname]M.[/fname][/pauthor], [pauthor][surname]Núñez-Casillas[/surname], [fname]L.[/fname][/pauthor], [pauthor][surname]Hernández-Leal[/surname], [fname]P.[/fname][/pauthor], [pauthor][surname]González-Calvo[/surname], [fname]A.[/fname][/pauthor] & [pauthor][surname]Alonso-Benito[/surname], [fname]A.[/fname][/pauthor][/authors] ([date dateiso="20090000" specyear="2009"]2009[/date]). [arttitle]Estimación de la emisividad espectral con datos térmicos ASTER de la Isla de Tenerife[/arttitle]. [source]IEEE International Geoscience & Remote Sensing Symposium, IGARSS 2009[/source], [publoc]University of Cape Town[/publoc], [pubname]South Africa[/pubname][/ref]
[ref id="r4" reftype="journal"][authors role="nd"][pauthor][surname]Becker[/surname], [fname]F.[/fname][/pauthor] & [pauthor][surname]Z.[/surname], [fname]Li.[/fname][/pauthor][/authors] ([date dateiso="19900000" specyear="1990"]1990[/date]). [arttitle]Temperature-independent spectral indices in thermal infrared bands, Remote Sens[/arttitle]. [source]Environ[/source]., [volid]32[/volid], [pages]17- 33[/pages].[/ref]
[ref id="r5" reftype="book"][authors role="nd"][pauthor][surname]Berk[/surname], [fname]A.[/fname][/pauthor], [pauthor][surname]Anderson[/surname], [fname]G. P.[/fname][/pauthor], [pauthor][surname]Acharya[/surname], [fname]P. K.[/fname][/pauthor], [pauthor][surname]Chetwynd[/surname], [fname]J. H.[/fname][/pauthor], [pauthor][surname]Bernstein[/surname], [fname]L. S.[/fname][/pauthor], [pauthor][surname]Shettle[/surname], [fname]E.P.[/fname][/pauthor], [pauthor][surname]Matthew[/surname], [fname]M.W.[/fname][/pauthor] & [pauthor][surname]Adler-Golden[/surname], [fname]S. M.[/fname][/pauthor][/authors] ([date dateiso="19990000" specyear="1999"]1999[/date]). [source]MODTRAN 4 user's manual[/source]. [pubname]Air Force Research Laboratory[/pubname], [publoc]Space Vehicles Directorate, Air Force Materiel Command, Hascom AFB, MA.[/publoc][/ref]
[ref id="r6" reftype="journal"][authors role="nd"][pauthor][surname]Carlson[/surname], [fname]T.[/fname][/pauthor] & [pauthor][surname]Ripley[/surname], [fname]D.[/fname][/pauthor][/authors] ([date dateiso="19970000" specyear="1997"]1997[/date]). [arttitle]On the Relation between NDVI, Fractional Vegetation Cover, and Leaf Area Index[/arttitle]. [source]Remote Sensing of Environment[/source]. [volid]62[/volid]. [pages]241-252[/pages]. [moreinfo]10.1016/S0034-4257(97)00104-1[/moreinfo].[/ref]
[ref id="r7" reftype="journal"][authors role="nd"][pauthor][surname]Chander[/surname], [fname]G.[/fname][/pauthor], [pauthor][surname]Markham[/surname], [fname]B. L.[/fname][/pauthor], & [pauthor][surname]Helder[/surname], [fname]D. L.[/fname][/pauthor][/authors] ([date dateiso="20090000" specyear="2009"]2009[/date]). [arttitle]Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors[/arttitle]. [source]Remote Sensing of Environment[/source] [volid]113[/volid]: [pages]893-903[/pages].[/ref]
[ref id="r8" reftype="journal"][authors role="nd"][pauthor][surname]Coll[/surname], [fname]C.[/fname][/pauthor], [pauthor][surname]Galve[/surname], [fname]J.[/fname][/pauthor], [pauthor][surname]Sánchez-Tomás[/surname], [fname]J.[/fname][/pauthor] & [pauthor][surname]Caselles[/surname], [fname]V.[/fname][/pauthor][/authors] ([date dateiso="20100000" specyear="2010"]2010[/date]). [arttitle]Validation of Landsat-7/ETM+ Thermal-Band Calibration and Atmospheric Correction with Ground-Based Measurements[/arttitle]. [source]Geoscience and Remote Sensing, IEEE Transactions on[/source]. [volid]48[/volid]. [pages]547 - 555[/pages]. [moreinfo]10.1109/TGRS.2009.2024934[/moreinfo].[/ref]
[ref id="r9" reftype="book"][authors role="nd"][pauthor][surname]Denyer[/surname], [fname]P.[/fname][/pauthor], [pauthor][surname]Aguilar[/surname], [fname]T.[/fname][/pauthor] & [pauthor][surname]Montero[/surname], [fname]W.[/fname][/pauthor][/authors] ([date dateiso="20140000" specyear="2014"]2014[/date]). [source]Cartografía geológica de la península de Nicoya[/source]. [pubname]Universidad de Costa Rica[/pubname]. [/ref]
[ref id="r10" reftype="journal"][authors role="nd"][pauthor][surname]Gillespie[/surname], [fname]A. R.[/fname][/pauthor][/authors], ([date dateiso="19850000" specyear="1985"]1985[/date]). [arttitle]Lithologic mapping of silicate rocks using TIMS[/arttitle], in [source]The TIMS Data Users' Workshop[/source], [pubname]JPL Publ.[/pubname], [pages]86-38, 29-44[/pages][/ref]
[ref id="r11" reftype="book"][authors role="nd"][pauthor][surname]Gillespie[/surname], [fname]A. R.[/fname][/pauthor], [pauthor][surname]Rokugawa[/surname], [fname]S.[/fname][/pauthor], [pauthor][surname]Hook[/surname], [fname]S. J.[/fname][/pauthor], [pauthor][surname]Matsunaga[/surname], [fname]T.[/fname][/pauthor] & [pauthor][surname]Kahle[/surname], [fname]B.[/fname][/pauthor][/authors] ([date dateiso="19990000" specyear="1999"]1999[/date]). [source]Temperature/emissivity separation algorithm theoretical basis document, version 2.4.[/source] [contract]Contract NAS5-31372[/contract], [pubname]Natl. Aeronaut. and Space Admin.[/pubname], [publoc]Washington, D. C.[/publoc][/ref]
[ref id="r12" reftype="journal"][authors role="nd"][pauthor][surname]Niclòs[/surname], [fname]R.[/fname][/pauthor], [pauthor][surname]Estrela[/surname], [fname]M. J.[/fname][/pauthor], [pauthor][surname]Valiente[/surname], [fname]J. A.[/fname][/pauthor] & [pauthor][surname]Barberà[/surname], [fname]M. J.[/fname][/pauthor][/authors] ([date dateiso="20100000" specyear="2010"]2010[/date]). [arttitle]Análisis de correlaciones entre la temperatura del aire y la temperatura de las superficies vegetadas medida con radiometría térmica[/arttitle]. [source]Revista de Teledetección[/source], [volid]34[/volid]: [pages]36-43[/pages][/ref]
[ref id="r13" reftype="book"][authors role="nd"][pauthor][surname]Pérez[/surname], [fname]C.[/fname][/pauthor] y [pauthor][surname]Muñoz[/surname], [fname]A.L.[/fname][/pauthor][/authors] ([date dateiso="20020000" specyear="2002"]2002[/date]). [source]Teledetección: nociones y aplicaciones[/source]. [pubname]Universidad de Salamanca[/pubname], [publoc]España[/publoc].[/ref]
[ref id="r14" reftype="journal"][authors role="nd"][pauthor][surname]Prata[/surname], [fname]J. A.[/fname][/pauthor][/authors] ([date dateiso="19940000" specyear="1994"]1994[/date]). [arttitle]Land surface temperatures derived from the advanced very high resolution radiometer and the along-track scanning radiometer: 2. Experimental results and validation of AVHRR algorithms[/arttitle]. [source]Journal of Geophysical Research[/source], [volid]99[/volid], [issueno]NO. D6[/issueno], [pages]P. 13,025-13,058[/pages].[/ref]
[ref id="r15" reftype="book"][authors role="nd"][pauthor][surname]Quattrochi[/surname], [fname]D. A.[/fname][/pauthor], & [pauthor][surname]Luvall[/surname], [fname]J. C.[/fname][/pauthor][/authors] ([date dateiso="20000000" specyear="2000"]2000[/date]). [source]Thermal remote sensing in land surface processes[/source]. [publoc]Boca Raton, FL, USA[/publoc], [pubname]CRC Press.[/pubname][/ref]
[ref id="r16" reftype="book"][authors role="nd"][pauthor][surname]Sobrino[/surname], [fname]J. A.[/fname][/pauthor], [pauthor][surname]Sobrino[/surname], [fname]J.[/fname][/pauthor], [pauthor][surname]Raissouni[/surname], [fname]N.[/fname][/pauthor], [pauthor][surname]Kerr[/surname], [fname]Olioso[/fname][/pauthor], [pauthor][surname]López[/surname], [fname]M. J.[/fname][/pauthor], [pauthor][surname]Belaid[/surname], [fname]El-Kharraz[/fname][/pauthor], [pauthor][surname]J.[/surname], [fname]Cuenca[/fname][/pauthor], [pauthor][surname]J.[/surname] & [fname]Dempere[/fname][/pauthor][/authors] ([date dateiso="20000000" specyear="2000"]2000[/date]). [pubname]Teledetección[/pubname]. [moreinfo]10.13140/2.1.2883.3283[/moreinfo].[/ref]
[ref id="r17" reftype="journal"][authors role="nd"][pauthor][surname]Sobrino[/surname], [fname]J.[/fname][/pauthor] & [pauthor][surname]Jiménez-Muñoz[/surname], [fname]J. C.[/fname][/pauthor] & [pauthor][surname]Paolini[/surname], [fname]L.[/fname][/pauthor][/authors] ([date dateiso="20040000" specyear="2004"]2004[/date]). [arttitle]Land surface temperature retrieval from LANDSAT TM 5[/arttitle]. [source]Remote Sensing of Environment[/source] . [moreinfo]90. 434-440. 10.1016/j.rse.2004.02.003[/moreinfo].[/ref]
[ref id="r18" reftype="journal"][authors role="nd"][pauthor][surname]Valor[/surname], [fname]E.[/fname][/pauthor] y [pauthor][surname]Caselles[/surname], [fname]V.[/fname][/pauthor][/authors] ([date dateiso="19960000" specyear="1996"]1996[/date]). [arttitle]Mapping Land Surface Emissivity from NDVI: Application to European, African, and South American Areas[/arttitle]. [source]Remote Sensing of Environment[/source], [volid]57[/volid], [pages]167-184[/pages][/ref][/refs][/doc].