COMPOSICIONES LANDSAT

EN ARCGIS
Guía Básica

Herramientas en este ejercicio:

• El sensor Landsat, bandas y visualización on-
line
• Búsqueda y descarga de imágenes
• Selección de combinaciones
Rodolfo Franco
• Composición RGB
• Recortar según área de estudio Especialista SIG
© 2017
 Composiciones Landsat en ARCGIS

“Todo lo que le ocurra a la tierra,

le ocurrirá a los hijos de la tierra”

Jefe Indio Seattle, 1855

Imagen de portada
Escena Landsat 8, Path–Row 7-55,
Composición RGB 652, enero 16
de 2016, Sierra Nevada del Cocuy,
Boyacá, Colombia

FRANCO, Rodolfo, 2017. Composiciones

Landsat en ARCGIS. Guía Básica. Bogotá,
Colombia. 45p. Disponible en
http://wp.me/p2IwQU-1bh

http://mixdyr.wordpress.com
 CONTENIDO

Página
PRESENTACIÓN…………………………………………….....….. 4

1. EL SENSOR LANDSAT……………………………………….…. 5

1.1. El programa Landsat…………………………………….… 5

1.2. Landsat 8………………………………………………….… 5
1.3. La cuadrícula Landsat…………………………………..… 7
2. BANDAS LANDSAT……………………………….…………..… 10
2.1. Paralelo Bandas Landsat 7 y Landsat 8………..…......... 10
2.2. Perfil o respuesta espectral………………………….……. 10
2.3. Uso típico de cada banda....................................……. 12
2.4. Ejemplos de imagen por banda………………….......... 13
3. COMPOSICIONES LANDSAT…………………………………. 17
3.1. Principales combinaciones y su uso………..…………… 17
3.2. Landsat on AWS……………………………………………. 18
3.3. Ejemplos de composiciones……………………….......... 21
4. BÚSQUEDA Y DESCARGA DE ESCENAS LANDSAT…….......... 25
4.1. Sitios de descarga…………………………………………. 25
4.2. Ejemplo de descarga en Earth Explorer………….......... 26
5. CONFIGURAR UNA COMPOSICIÓN EN ARCGIS…………… 32
5.1. Preparación del proyecto y archivos……………………. 32
5.2. Carga de bandas……………………………………........ 33
5.3. Composición……………………………………………….. 35
5.4. Recorte al área de estudio……………………………….. 38
5.5. Otro método: Image Analysis……………………………. 40
6. COMENTARIOS FINALES………………………………………. 45
 PRESENTACIÓN

La presente mini-guía muestra un ejemplo de caso para la generación de

una composición satelital con bandas del sensor Landsat mediante el
software ArcGis de ESRI.

Aunque ArcGis es más conocido como un poderoso programa en Sistemas de

Información Geográfica también realiza algunos procesos de
procesamiento digital de imágenes. Ya su antecesor, ArcView, hace unos
15 años, permitía algunos tratamientos mediante una extensión
llamada Image Analysis. Hoy ArcGis incluye cada vez más herramientas
raster con múltiples opciones para la manipulación de escenas satelitales.

Se ilustran en esta guía algunos conceptos básicos sobre el sensor Landsat,

la forma en que se pueden buscar, visualizar y descargar imágenes On-
line, consejos sobre la selección de bandas, la generación de una
composición y finalmente el recorte una zona de estudio específica dentro
de la escena.

Para conocer más de ARCGIS se puede visitar mi web en

https://mixdyr.wordpress.com/

Rodolfo Franco
Ingeniero Forestal, Especialista SIG

http://mixdyr.wordpress.com
4
 1. EL SENSOR LANDSAT

1.1. El programa Landsat

La serie de satélites estadounidenses Landsat constituye el programa de

mayor trayectoria en adquisición de imágenes de nuestro planeta. El primer
instrumento se lanzó en 1972 y desde entonces la serie ha evolucionado como se
muestra en la siguiente tabla:

Satélite Lanzamiento Notas

Landsat 1 1972 Primero de la serie, se llamó originalmente Earth Resources
Technology Satellite 1. Estuvo activo durante casi tres años
y
fue operado por la NASA. Incorporaba un escáner
multiespectral (MSS) de cuatro canales.
Landsat 2 1975 De características similares a Landsat 1 y estuvo activo
por
cerca de 7 años. Operado por la NASA.
Landsat 3 1978 Tenía el mismo diseño de Landsat 2, estuvo activo por
unos 5 años y fue operado por la NASA. A diferencia de
sus antecesores incluyó además un sensor en una banda
térmica pero este instrumento falló al poco tiempo del
despliegue del
satélite.
Landsat 4 1982 Operado por Earth Observation Satellite Company (EOSAT)
Fue el primero en incluir un sensor Thematic Mapper (TM)
de 7 bandas, una de estas termal.
Landsat 5 1984 Idéntico a Landsat 4. Fue operado por EOSAT y posee
el record como el satélite de observación terrestre de
mayor tiempo de actividad en la historia con 29 años de
servicio.
Transmitió más de 2,5 millones de imágenes.
Landsat 6 1993 Fue lanzado el 5 de octubre de 1993 pero no alcanzó su
órbita.
Landsat 7 1999 Incluye una banda en pancromático con 15m de
resolución espacial. Aún está activo pero la calidad de sus
imágenes bajo al ocurrir una falla en el Scan Line Corrector
(SLC) en mayo de 2003 que provoca que se pierdan
líneas de escaneo que
suman un 22% de cada escena. Es operado por el
USGS.

Curso ARCGIS- Profesor Rodolfo Franco, conozca más en

http://mixdyr.wordpress.com
5
Landsat 8 2013 Activo desde febrero de 2013. Incorpora dos sensores:
por una parte el Operational Land Imager (OLI) con 9
bandas que incluyen una pancromática y, por otro, el
Thermal Infrared
Sensor (TIRS) con dos bandas. Es operado por el USGS
Landsat 9 2020 Lanzamiento tentativo para diciembre de 2020. Lo
construirá
la NASA y sus datos serán administrados por el USGS.

http://mixdyr.wordpress.com
6
 Línea de tiempo de la serie satelital Landsat

1.2. Landsat 8

El satélite Landsat 8 fue lanzado el 11 de febrero de 2013 en un proyecto

conjunto entre el Servicio Geológico de los Estados Unidos (USGS) y la NASA. Se
tiene proyectada una vida útil de entre 5 y 10 años pero puede ser mayor
considerando el tiempo de su antecesor, Landsat 7 que aún está activo en órbita
luego de más de 17 años muy por encima de sus expectativas.

Landsat 8 se encuentra orbitando a una distancia de entre 701 y 703 km sobre la

superficie terrestre, pasa por un mismo sitio del planeta cada 16 días y a los 8 días
del paso de Landsat 7. El satélite captura unas setecientas imágenes cada día y
cada escena tiene un cubrimiento de 185 x 180 km.

El producto Landsat es del tipo Level 1T, es decir, con calibración radiométrica y
con ortorectificación sustentada con puntos de control en terreno y un modelo
digital de elevación (DEM). El desempeño radiométrico es cuantificado sobre un
rango dinámico de 12-bits.

Los instrumentos principales de Landsat 8 son el OLI (Operational Land Imager) y el

TIRS (Thermal Infrared Sensor). OLI colecta imágenes en nueve bandas espectrales
con resolución espacial de 30 m a excepción de la imagen pancromática con
15 metros. OLI arroja un error circular de 12 metros en la precisión cartográfica
para un nivel de confianza del 90%. Las escenas capturadas con este sensor
permiten distinguir áreas urbanas, cultivos, bosques, usos del suelo, entre otros.
El sensor TIRS, por su parte, colecta escenas térmicas en dos bandas de infrarrojo
de onda larga. La resolución especial es de 100 metros pero se puede
remuestrear a 30 m por fusión con las imágenes de OLI. Respecto a la precision
cartográfica, TIRS arroja error circular de 41 metros para un nivel de confianza del
90%. Las escenas con este instrumento son útiles para mediciones de la
temperatura de la superficie, evapotranspiración, entre otros.

Landsat 8 en órbita. Representación NASA

1.3. La cuadrícula Landsat

El catálogo de escenas Landsat usa una cuadrícula mundial de referencia

también conocida como WRS (Worldwide Reference System). Los satélites Landsat
1,2 y 3 emplearon la cuadrícula WRS-1 mientras que Landsat 7 y 8 emplean la
WRS-2, la misma que usó Landsat 4 y 5. La cuadrícula viene conforme el recorrido
que hace el satélite sobre el globo así que tiene una configuración ascendente o
nocturna y la descendente, que es la que corresponde a las escenas capturadas
con la luz diurna. La cuadrícula cubre el globo totalmente en longitud y entre los
82,6° y -86,2° de latitud. Shapefiles de los sistemas WRS-1 y WRS-2 se
pueden descargar en https://landsat.usgs.gov/pathrow-shapefiles

En este sistema cada escena es denotada por un número específico de columna

y fila o Path-Row dentro de la cuadrícula.

Para el territorio colombiano las escenas Landsat se codifican entre el Path 3 al 14

y de Row 49 al 63. La siguientes imágenes muestran el cubrimiento Landsat para el
país y se puede consultar on-line en http://arcg.is/2bcGOxw. La versión en formato
Shapefile se puede descargar en: https://drive.google.com/open?
id=0B48yy4whz_7veTBzdGMyQkRDVEk
Cuadrícula WRS-2 de Landsat para Colombia
Mapa disponible en http://arcg.is/0vWm5K
 Cuadrícula WRS-2 de Landsat para
Colombia Shapefile disponible en
https://drive.google.com/open?id=0B48yy4whz_7vTEZOcklXdGlnb2M
 2. BANDAS LANDSAT

2.1. Paralelo bandas Landsat 7 y Landsat 8

La siguiente tabla presenta un paralelo y detalle de las bandas espectrales de los

sensores de Landsat 7 ETM+ y Landsat 8 OLI –TIRS. Se indica el nombre de banda,
rango de ancho de banda en µm y la resolución espacial en metros:

LANDSAT 7 ETM+ LANDSAT 8 OLI -TIRS

Banda Ancho (µm) Resolución Banda Ancho (µm) Resolución
Banda 1 Coastal 0.435 - 0.451 30 m
Banda 1 Blue 0.441 - 30 m Banda 2 Blue 0.452 - 0.512 30 m
0.514
Banda 2 Green 0.519 - 30 m Banda 3 Green 0.533 - 0.590 30 m
0.601
Banda 3 Red 0.631 - 30 m Banda 4 Red 0.636 - 0.673 30 m
0.692
Banda 4 NIR 0.772 - 30 m Banda 5 NIR 0.851 - 0.879 30 m
0.898
Banda 5 SWIR 1.547 - 30 m Banda 6 SWIR 1 1.566 - 1.651 30 m
1 1.749
Banda 5 SWIR 1.547 - 30 m Banda 7 SWIR 2 2.107 – 2.294 30 m
1 1.749
Banda 8 Pan 0.515 - 15 m Banda 8 Pan 0.503 - 0.676 15 m
0.896
Banda 9 Cirrus 1.363 - 1.384 30 m
Banda 10 TIRS 10.60 - 11.19 100 m
Banda 6 TIR 10.31 – 60 m 1
12.36
Banda 11 TIRS 11.50 - 12.51 100 m
2

2.2. Perfil o respuesta espectral

Cada banda según sus características ofrece ventajas para el reconocimiento de

determinados objetos de la superficie en función al comportamiento espectral de
los mismos.
Las diferentes coberturas terrestres responden de forma diferente en términos de
reflectividad de la luz solar según la longitud de onda. La siguiente gráfica muestra
la respuesta espectral típica de algunos elementos en relación con las
nueve bandas del sensor OLI:
 Respuesta espectral de algunas coberturas en las bandas del sensor OLI

La firma espectral obedece a la composición y estructura misma de los objetos.

La respuesta de la vegetación, por ejemplo, está asociada a la absorción por la
clorofila, el contenido de agua, entre otros. Los suelos, por su parte, reaccionan
distinto bien se trate de suelos arcillo-limosos, areniscas, suelos orgánicos, etc. En
la nieve, por otro lado, tiene mucha incidencia el tamaño granular y, en el en
caso del agua, su turbidez.
2.3. Uso típico de cada banda

Cada imagen capturada por los sensores en cada banda tendrá una
apariencia distintiva según las coberturas presentes. La siguiente tabla resume
las ventajas y uso típico de cada banda para el caso de Landsat 8:

Banda en Landsat Aplicación principal / Descripción

8
Banda 1 Para mapeo de costas y estudios de aerosol.
Coastal
Banda 2 Útil para mapeo batimétrico, delimitar costas,
Azul diferenciar suelo de vegetación, diferenciar
coníferas de latifoliadas, detección de rasgos
urbanos, vías y
construcciones.
Banda 3 Empleada para discriminar sedimentos en
Verde suspensión, evaluar vigor de las plantas por la alta
reflectancia de la vegetación verde y sana,
delinear aguas poco
profundas, rasgos urbanos y de infraestructura.
Banda 4 Permite un mayor contraste de áreas con y sin
Rojo vegetación, discriminar gradientes de vegetación,
delimitar áreas urbanas y áreas agrícolas.
Banda 5 Útil para el cálculo de biomasa de vegetación,
NIR delimitar costas, para diferenciación suelos-cultivos
Infrarrojo cercano y suelos- agua, para geomorfología, suelos y geología.
Banda 6 Se alcanzan a penetrar nubes delgadas. Es útil
SWIR 1 para
discriminar contenido de humedad en los suelos y la
vegetación, diferenciar entre nubes, nieve y hielo.
Banda 7 Útil para discriminar tipos de rocas, para estudios de
SWIR 2 suelos y mejora la determinación de contenidos de
humedad en suelos y vegetación.
Banda 9 Es una imagen sensible a todo el espectro visible y más
Pancromática afinada en tanto su resolución es de 15m
Banda 8 Ofrece una mejora en la detección de nubes cirrus.
Cirrus
Banda 10 Útil para mapeo termal y estimación de humedad del
TIRS 1 suelo.
Térmica
Banda 11 Mapeo termal mejorado y estimación de humedad del
TIRS 2 suelo.
Térmica
2.4. Ejemplos de imagen por banda

A continuación un ejemplo de la apariencia de cada banda. La escena de muestra

es adquirida del sitio Landsat on AWS y corresponde a una porción de la imagen del
24 de febrero de 2017, Path-Row 10-53, sector de San Bernardo del Viento y San
Antero, en el departamento de Córdoba, costa atlántica de Colombia:

Océano Atlántico

COLOMBIA

Banda 1 Coastal

Para mapeo de
costas y estudios de
aerosol.

Banda 2 Blue

Útil para mapeo

batimétrico,
delimitar costas,
diferenciar suelo de
vegetación, diferencia
coníferas de
latifoliadas,
detección de rasgos
urbanos, vías y
construcciones.
Banda 3 Green

Usada para
discriminar
sedimentos en
suspensión, evaluar
vigor de las
plantas por la alta
reflectancia de la
vegetación verde y
sana, para delinear
aguas poco
profundas, rasgos
urbanos y de
infraestructura.

Banda 4 Red
para contraste de
áreas con y sin
vegetación, gradientes
de
vegetación, sirve para
delimitar áreas
urbanas y áreas
agrícolas.

Banda 5 NIR

Útil para el cálculo

de biomasa de
vegetación, delimitar
costas, para
diferenciación suelos-
cultivos y suelos-
agua, para
geomorfología, suelos
y geología.
Banda 6 SWIR 1

Se alcanzan a penetrar
nubes delgadas. Es útil
para discriminar
contenido de
humedad en los suelos
y la vegetación,
diferenciar entre
nubes, nieve y hielo.

Banda 7 SWIR 2

Útil para discriminar

tipos de rocas, para
estudios de suelos y
mejora la
determinación de
contenidos de
humedad en suelos y
vegetación.

Banda 8 Pancromática

Es una imagen sensible

a todo el espectro
visible y más afinada
en tanto su
resolución es de 15m
Banda 9 Cirrus

Ofrece una
mejora en la
detección de nubes
cirrus.

Banda 10 TIRS 1

Útil para mapeo

termal y
estimación de
humedad del suelo.

Banda 11 TIRS 2

Mapeo termal
mejorado y
estimación de
humedad del suelo.
 3. COMPOSICIONES LANDSAT

3.1. Principales combinaciones y su uso

Una composición satelital Landsat es una imagen compuesta por la combinación

de tres bandas diferentes del sensor y dispuestas cada una en los tres canales de
proyección en pantalla: Rojo (R), Verde (G) y Azul (B). La escena en RGB resulta en
una imagen a color. La siguiente es una relación de algunas de
las combinaciones más comunes y el uso típico de estas:

Combinación RGB en Aplicación principal / Descripción

LANDSAT 8
431 Batimetría
Permite destacar objetos en cuerpos de agua
432 Color Natural
Es una composición que se aproxima al color real de las
coberturas ante la vista humana
Su equivalente en Landsat 7 es la combinación 321
543 Infrarrojo Color
Esta combinación da como resultado una imagen muy similar
a la tradicional fotografía aérea infrarrojo color.
Es útil para estudios de vegetación, patrones de suelos,
crecimiento de cultivos y monitoreo de drenajes.
La vegetación saludable tiende a una apariencia rojo
brillante
562 Combinación útil para distinguir vegetación saludable
564 Combinación útil para diferenciar tierra / agua y para
destacar
áreas de bosque
652 Agricultura
654 Análisis de vegetación
742 Geología
Permite destacar elementos geológicos
753 Natural con remoción atmosférica
754 Combinación de infrarrojos
764 Permite destacar formaciones rocosas. También útil para
para
distinguir áreas urbanas
765 Penetración Atmosférica
Debido a que no incorpora bandas del espectro visible
esta combinación reduce las influencias atmosféricas en la
imagen. Es útil para estudios geológicos así como para
análisis de
humedad en suelos
3.2. Landsat on AWS

Una buena estrategia para agilizar el trabajo con escenas Landsat multibanda es
probar combinaciones en sitios on-line para tal fin y luego sí descargar las
imágenes que se consideren apropiadas. Uno de los servicios web
más destacados con este propósito es Landsat on AWS que usa un visor
desarrollado en Web app Builder for ARCGIS, tecnología ESRI.

En este aplicativo se puede hacer búsqueda a un lugar de interés según palabras

clave y visualizar escenas en nueve modalidades pre-configuradas:

Agricultura (RGB 652)

Color Natural (RGB 432 fusionada con pancromática)

Infrarrojo Color (RGB 543)

Infrarrojo de onda corta SWIR (RGB 764)

Geología (RGB 764)

Batimetría (RGB 431)

Pancromática (banda 8)

Índice de vegetación de diferencia normalizada NDVI

[Banda 5 - Banda 4] / [Banda 5+ Banda 4]
Índice de humedad de diferencia normalizada NDMI
[Banda 5 – Banda 6] / [Banda 5+ Banda 6]
 Composiciones Landsat en ARCGIS

Además de las escenas predefinidas, el usuario puede también configurar su

propia combinación mediante la herramienta Build:

Una herramienta muy valiosa es Identify, con esta opción el aplicativo devuelve el
código, fecha, Patch-Row y otros datos de la escena Landsat actual.
Luego, usando el puntero, se puede hacer clic sobre cualquier píxel y el
sistema devolverá su perfil espectral y además señalará cuál es el perfil
típico de cobertura al cuál más se le aproxima. Esto se convierte en una gran
ayuda para la interpretación de la imagen.

Como ejemplo en la imagen inferior, se ha consultado para un pixel particular en

la combinación 652 y el perfil devuelto se acerca a la respuesta típica de un
bosque. Efectivamente se trata aquí de un manglar.
Otra de las ventajas del sitio es la posibilidad de consultar parte del histórico de
escenas Landsat para el área de interés mediante la herramienta Time. El
sistema mostrará una barra de deslizamiento a través del tiempo y el
usuario podrá retroceder a viejas escenas de los primeros satélites con
fechas incluso de los años 70’s u 80’s según disponibilidad.

Las siguientes imágenes son ejemplos de escenas multitemporales del sector de

la Laguna de Fúquene, al norte del departamento de Cundinamarca, Colombia
(Path – Row 8-56).

Diciembre 16 de 1987
Landsat 5, RGB 541

Febrero 22 de 2007
Landsat 7, RGB 541

Noviembre 22 de 2016
Landsat 8, RGB 652
 Composiciones Landsat en ARCGIS

3.3. Ejemplos de composiciones

Como ejemplos de las principales combinaciones se han capturado del

sitio Landsat on AWS las siguientes escenas del sector de manglar entre San
Bernardo del Viento y San Antero, departamento de Córdoba, Colombia:

RGB 431 Batimetría

Es una combinación que

permite distinguir
objetos en cuerpos de
agua

RGB 432 Color Natural

Es una composición que

da una apariencia que
se aproxima al color real
de las coberturas ante
la vista humana.

RGB 543 Infrarrojo color

Similar a la
fotografía aérea
infrarrojo color.
Útil para estudios de
vegetación, patrones de
suelos, crecimiento de
cultivos y monitoreo
de drenajes. La
vegetación saludable
tiende a una
apariencia rojo brillante.
RGB 562

Combinación útil
para distinguir
vegetación saludable.

RGB 564

Combinación útil
para diferenciar tierra
de agua y para destacar
áreas de bosque.

RGB 652 Agricultura

Combinación que
ayuda para distinguir
áreas agrícolas.
RGB 654
Con especial uso para
análisis de vegetación

RGB 742 Geología

Permite distinguir
elementos geológicos.

RGB 753 Natural con

remoción atmosférica.
RGB 754
Combinación de
infrarrojos.

RGB 764 SWIR

Permite destacar
formaciones rocosas.
También útil para para
distinguir áreas urbanas.

RGB 765 Penetración

atmosférica

Debido a que no
incorpora bandas del
espectro visible esta
combinación reduce las
influencias
atmosféricas en la
imagen.
Es útil para estudios
geológicos así
como para análisis
de humedad en suelos
 4. BÚSQUEDA Y DESCARGA DE ESCENAS LANDSAT

4.1. Sitios de descarga

Son varios y cada vez más numerosos los sitios web que permiten la descarga de
imágenes Landsat. Entre los principales figuran varios del USGS como Earth
Explorer, Glovis, Landsatlook Viewer y Global Data Explorer. Por supuesto hay varios
sitios de la NASA como REVERB y Earth Data Search. Se destacan otros muy
populares como Vertex (Portal de datos de sensores remotos del Alaska Satellite
Facility) y el GLCF (Global Land Cover Facility – Universidad de Maryland). También
hay unos sitios de interfaz muy ligera como Libra, Remote Pixel, Land
Viewer y Open Aerial Map. Los siguientes son enlaces a los sitios mencionados:

USGS: EARTH EXPLORER USGS: GLOVIS USGS: GLOBAL DATA EXPLORER

USGS: LANDSATLOOK VIEWER NASA: REVERB NASA: EARTH DATA SEARCH

GLCF LIBRA VERTEX

LAND VIEWER OPEN AERIAL MAP REMOTE PIXEL

Un amplio listado de webs para descargas satelitales se puede conocer en el

siguiente link: https://mixdyr.wordpress.com/p-remota/descargas-satelitales/
4.2. Ejemplo de búsqueda y descarga en Earth Explorer

Se accede a Earth Explorer desde la dirección https://earthexplorer.usgs.gov/

Es mejor acceder como usuario registrado si de desean descargar escenas
completas. El registro se puede realizar por la opción Register en la parte superior
derecha:

El proceso de registro comprende

cuatro pestañas:

una primera sección es para indicar

un usuario y contraseña. La siguiente
sección es para brindar información
demográfica entre la cual figura
el tipo de trabajo que tiene el
usuario y el tipo de uso que va a
brindar a la información. La tercera
pestaña es para aportar datos de
contacto y la última pestaña para
completar el registro.

Por supuesto, en la próxima ocasión

que se acceda a la página se puede
hacer como usuario registrado desde la
opción Login
Nota importante: después de 2 horas de
inactividad se cierra la sesión de forma
automática. Para evitar inconvenientes
se puede actualizar la sesión desde
la opción Refresh Timer en la parte
superior derecha.
Una vez abierta la sesión se procede a realizar la búsqueda de imágenes
iniciando por la selección según área de estudio:

En la primera pestaña de esta sección

(Search Criteria) se puede hacer la
búsqueda según nombres de
lugares, marcando un punto o trazando
un polígono sobre el área de interés.

También puede ser indicando un Path-

Row como se ve en la imagen de
ejemplo. Para Landsat 7 y 8 el Path-Row
a indicar es sobre la cuadrícula WRS2.

Entre otras alternativas para

selección figuran la opción de indicar
coordenadas o usar una capa de
referencia en shapefile o en KML.

La búsqueda se puede restringir a

un período de tiempo específico si se
desea.

Ejemplo de
selección
de área de
studio por
trazado de un
polígono
La siguiente pestaña (Data Sets) permite elegir el tipo de producto a encontrar.

Considerando interés en imágenes

del sensor Landsat se
selecciona Landsat Archive. En
esta sección se podrán buscar
entre imágenes de diferentes
series de Landsat y que han sido
reorganizadas por el USGS según
sensor.

Con atención a estudios

multitemporales es de reconocer
que aquí se pueden buscar desde
escenas muy antiguas (Landsat 1-5
MSS) hasta las más
recientes (Landsat 8 OLI/TIRS).

También hay disponibilidad de

imágenes de Landsat 7 ETM+
desde 1999 pero hay que
considerar que las escenas
obtenidas desde fines de mayo de
2003 tienen líneas incompletas
debido a una falla en el SCL (Scan
Line Corrector).

Como se aprecia en la imagen de

ejemplo a la izquierda, en este
caso se pretende indagar sobre
imágenes Landsat 8 (2013 – al
presente)

En la sección Landsat Legacy se

pueden buscar mosaicos y otros
productos diversos anteriores a
2003.
La tercera pestaña (Additional Criteria) es donde se puede filtrar aún más la
búsqueda con criterios adicionales entre los cuales figuran los siguientes:

Si se conoce el identificador de la
imagen se puede
digitar directamente dicho
código en espacio para tal fin.

Se puede especificar un Path-Row

específico o una serie de Path-Rows.
Se puede elegir el sensor
específico de la serie, en el caso de
Landsat 8 se tiene:

• OLI (Operational Land Imager)

• OLI_TIRS (Operational Land Imager)
y TIRS (Thermal Infrared Sensor)
• TIRS (Thermal Infrared Sensor)

Se puede especificar el nivel de

procesamiento:

Level 1GT: con

calibración radiométrica y
correciones geométricas sistemáticas
con uso de un DEM pero sin puntos
de control en terreno.

Level 1T: es el producto de mayor

calidad en tanto se trata de
imágenes con calibración
radiométrica y ortorectificadas
usando puntos de control en terreno y
un modelo digital de elevación (DEM)
para corregir desplazamiento por el
relieve.

Como se aprecia en la imagen de ejemplo arriba en este caso se

buscarán imágenes Landsat 8 correspondientes al Path-Row 7-56 con
cualquier porcentaje de nubes y cualquier nivel de procesamiento.
Para terminar la búsqueda se ingresa a la pestaña de resultados (Results) y se elige
una imagen o imágenes con cualidades satisfactorias. Hay que considerar que se
trata de grandes escenas de unos 185 x 180 km pero donde podría suceder que
el área de estudio solo ocupe una porción. En este caso hay que examinar que
sea una imagen donde haya escasa o nula nubosidad para el área de estudio
en la fecha o fechas disponibles.

Haciendo clic sobre la

imagen miniatura se puede acceder a
tres vistas preliminares y el
metadato o ficha descriptiva de la
escena.

La primera vista previa corresponde a

una composición de falso color
RGB 654 (sería RGB 543 si es Landsat
7), se trata de una escena con
apariencia “natural” donde se han
ajustado los valores de reflectancia
para enfatizar la vegetación:
 Otra imagen en el caso de Landsat 8 es una vista
previa de una escena termal (Banda 10)

Una descripción detallada de esta vista y otras

se puede ver en el siguiente
enlace: https://landsat.usgs.gov/landsatlook-images

La última imagen en vista previa es una relacionada

con la calidad de la misma. Allí se marcan algunos
pixels con determinados colores ya que que por sus
características pueden afectar el uso e
interpretación de la escena.

En la imagen de ejemplo a la izquierda aparecen

pixels amarillos (nubes cirrus) y pixels blancos (nubes
en general).

Opciones de descarga

Para descargar la imagen se accede

al icono de Download options.

Los productos a descargar

corresponen a las tres vistas
previas anteriormente
mencionadas, las mismas
escenas pero
georreferenciadas y al final lo más
importante: la imagen con su
conjunto completo de bandas
en formato GeoTIFF y que puede
tener un tamaño entre 800 Mb y 1 GB.

Esta descarga resulta en un archivo

que puede ser descomprimido
fácilmente con Winrar, WinZip, etc.
 5. CONFIGURAR UNA COMPOSICIÓN EN ARCGIS

5.1. Preparación del proyecto y archivos

Como es habitual en ARCGIS se recomienda crear una carpeta para

el proyecto donde se trabajará con la imagen. También crear
una subcarpeta para las capas que se tengan, como por ejemplo, un
Shapefile del área de estudio) y las imágenes que se van a ir generando
durante el proceso de configuración de la composición.

Se descomprime el archivo de la imagen descargada en la carpeta del

proyecto.

Se ingresa a ArcMap, se guarda el proyecto con un nombre apropiado y,

como es habitual, se le da la ruta relativa (Menú File-Map document
properties – Store relative pathnames…). Aunque las
imágenes descargadas vienen en WGS84 a la vista se le puede asignar
el sistema proyectado local si se desea.
5.2. Carga de bandas

Para cargar las imágenes se usa el botón habitual de adicionar capas, sin
embargo, para agilizar el proceso se recomienda que el orden en que se
vayan adicionado los archivos sea el inverso a como se desea
luego configurar la composición. Por ejemplo, si se desea más adelante
obtener la combinación RGB 5-6-4 entonces se inicia primero cargando la
banda 4, luego la 5 y finalmente la 6.

Luego de seleccionar el archivo y adicionarlo, muy probablemente saldrá

un aviso sobre layers piramidales. Se trata de submuestreos del ráster
original a medida que se hace zoom a la imagen y permiten acelerar el
despliegue de una imagen a diferentes resoluciones. Cuando se trabajan
rasters tanto en SIG, en procesamiento digital de imágenes o en
diseño gráfico, es conveniente aceptar la creación de estos
elementos para mayor rendimiento. Así que se acepta su creación.

Para cada imagen agregada se le crearán los layers piramidales y como

fruto de esto se formará un archivo .ovr para cada archivo original 1
además del respectivo archivo .xml

1 En el caso de imágenes de ERDAS IMAGINE se formará un archivo .rrd

En el caso que para la vista se le tengan asignadas coordenadas
proyectadas el sistema avisará que la imagen a cargar está en otro
sistema. No hay inconveniente púes de todos modos la imagen está
georreferenciada.

Una vez agregadas las bandas, estas aparecerán como capas rasters en
tonos de grises como se ve en la siguiente imagen:
5.3. Composición

Para crear una composición de bandas se puede proceder con la

herramienta Composite Bands (Pestaña Catalog>>Toolboxes>> System
Toolboxes>> Data Management Tools>>Raster>>Raster
Processing>> Composite Bands).

Se cargan las bandas

en el mismo orden
en que se
desea
configurar la
combinación.

Como ubicación de
salida se
selecciona la
subcarpeta
preparada para tal
fin. El nombre de
salida debe ser corto,
sin espacios en
blanco y en el mismo
formato.(En este caso
es .tif)
 Nueva imagen

Aparecerá entonces la escena resultante como un nuevo raster. Se

puede verificar en la carpeta de destino que la nueva imagen se
compone de varios archivos y ya esta lista para ser utilizada en este o en
otros proyectos. Ahora si se desea se pueden remover de la tabla
de contenido las imágenes originales de las bandas.

Si la imagen sale en formato .TIFF, el siguiente es el conjunto de archivos

resultantes:
TFW: es el worldfile o archivo con las
coordenadas de la imagen.
TIF: es la imagen generada
.tif.aux (XML): es un archivo de
auxiliar con estadísticas de
la imagen
.ovr: es el archivo de los
layers piramidales
.tif (XML): es el archivo de metadato

Es importante entonces considerar que si la composición generada se

quiere copiar o llevar a otra carpeta, se deben llevar TODO el conjunto de
archivos mencionados.
Desde la opción de propiedades de la imagen (Clic derecho al nombre de la
composición >> Properties) se puede consultar, entre otros, lo siguiente:

En Display se
puede ajustar el
contraste, brillo y
transparencia.

En Simbology se
puede cambiar el
orden de las bandas,
elegir un stretch
determinado ((ajuste
de rango de valores
de histograma) y revisar
los histogramas por
cada canal de
proyección (RGB). Un
histograma es un
gráfico con la
frecuencia de pixeles
por cada valor.

El stretch que está por

defecto es aquel
por percent clip donde
es aplicado el ajuste
pero excluyendo
un
porcentaje de los
valores extremos
máximos y mínimos en
la amplitud del
histograma.

También estará activa

la opción de un stretch
gamma que permite
controlar el excedente
de brillo en la imagen.
5.4. Recorte al área de estudio
Para recortar la imagen al área de estudio primero se carga la capa que
sirva como límite (1) y se procede por la herramienta Clip (2) (Pestaña
Catalog>>Toolboxes>> System Toolboxes>> Data Management
Tools>>Raster>>Raster Processing>> Clip).

La imagen a recortar

Polígono límite

Activar esta opción

para que la imagen se
recorte con la forma
del polígono. De
lo contrario el recorte
será al rectángulo
mínimo que lo
contiene.
Ubicación y nombre de
la imagen resultante
incluyendo la extensión
respectiva (en este caso
.tif )
En valor para los
datos nulos dejar en
blanco.
El sistema procesará y aparecerá la escena recortada resultante como un
nuevo raster. Si se desea, se podrá remover de la tabla de contenido
la composición original. La escena nueva tendrá valor por defecto en
las propiedades de display y Simbology que pueden ajustarse
como se consideren convenientes para mejorar la visualización.

Este es el zoom a la escena recortada:

Al igual que la escena original, la imagen recortada también
conocida como subset, tendrá sus respectivo conjunto de archivos que la
conforman y que deberán copiarse en su totalidad cuando se desee llevar
el producto a otro proyecto.

Tanto la escena original como la imagen recortada, ya pueden

ser empleadas para ejercicios de interpretación, clasificación y otros
ejercicios típicos de procesamiento digital de imágenes multibanda.

5.5. Otro método: Image Analysis

A continuación se explica otro método para la creación de composición de

bandas y es por la vía de Image Analysis.

Se trata de una opción interesante donde se puede generar una composición

temporal para luego convertirla en un archivo definitivo y finalmente, si
es del caso, hacer el recorte respectivo a un área de estudio determinada.

(1) Primero se cargan las imágenes de las bandas de

interés de la misma forma que en el primer método y
nuevamente considerando abrirlas en el orden inverso a
como se desea la configuración.

(2) Para crear la composición se inicia el proceso

accediendo al menú Windows – Image Analysis.
 (3) Se seleccionan las bandas del listado con ayuda de la tecla Ctrl.

(4) Luego en la sección

Processing, se ejecuta
el botón de Composite
Bands

1 3

4
Luego de unos instantes aparecerá la composición como una nueva capa
raster. Sin embargo, esta escena generada solo es temporal, es decir, aún
no es un archivo independiente y solo estará visible para este proyecto por
ahora. Se puede cerrar ya la ventana de Image Analysis.

Conversión de la composición a un archivo

(1) Para convertir la escena temporal en un archivo se exporta la imagen con

clic derecho al nombre de la misma >>Data >> Export Data…
 Composiciones Landsat en ARCGIS

(2) Aparecerá la ficha para elegir parámetros de la conversión. Se

recomienda:
Dejar la misma
referencia espacial y el
mismo tamaño de
pixel.
Los datos nulos
(Nodata as) en blanco.
Como ubicación de
salida se selecciona la
subcarpeta preparada
para tal fin.
El nombre de
salida corto, sin
espacios en blanco y
en el mismo formato.

(3) El sistema procesará, aparecerá en la parte inferior una barra con el

porcentaje de progreso y al terminar preguntará si se desea agregar a la vista
y se acepta.

1
2

3
 Composiciones Landsat en ARCGIS

Aparecerá entonces la escena exportada como un nuevo raster pero en esta

oportunidad ya es un archivo como tal. Se puede verificar en la carpeta
de destino que la nueva imagen se compone de varios archivos y ya esta lista
para ser utilizada en este o en otros proyectos. Ahora la composición
inicial y las imágenes originales de las bandas se pueden quitar de la tabla de
contenido.

Para recortar la composición al área de estudio se puede usar la misma vía que
se ilustró con el primer método, es decir, por la herramienta Clip de Raster
Processing o mediante la opción Clip que también aparece en Image Analysis y
por supuesto, tendiendo una capa poligonal de referencia. El subset resultante
será una escena temporal que podrá convertirse a archivo definitivo si se desea y
por la misma rutina de la escena original (Data-Export data)

Clip
6. COMENTARIOS FINALES

Generar composiciones RGB es tan solo uno de los primeros insumos para
adelantar luego otros procesos como clasificaciones supervisadas y
no supervisadas e interpretaciones visuales de coberturas.

En futuras guías se mostrarán otros ejercicios fundamentados en escenas

multibanda y que son típicos del procesamiento digital de imágenes
tal como son los índices de vegetación, índices de humedad,
pan- sharpening, mosaicos, entre otros.

Rodolfo Franco
Especialista SIG
© 2017
https://mixdyr.wordpress.com