1

Anexo 2: Manual de pre procesamiento de imágenes satelitales Landsat 8, mediante el uso del
software ERDAS 2015.

MANUAL DE PRE PROCESAMIENTO DE

IMÁGENES SATELITALES LANDSAT 8
MEDIANTE EL USO DEL SOFTWARE ERDAS
2015.

Eliecer Armando Portilla

Iván Andrés Revelo

NOVIEMBRE 2017.
 2

CONTENIDO

1. Adquisición de imágenes satelitales Landsat 8…………………………………………………3

2. Importación de imagen Landsat 8.………………………………………………………...........6

3. Unión de bandas…………………………………………………………………………….......8

4. Fusión de imagen multibanda con imagen Banda 8 (Pancromática).…………………............13

5. Recorte de área de interés mediante Shapefile………………………………………………..20

 3

1. Adquisición de imágenes satelitales Landsat 8.

Este primer paso se realiza a través de la plataforma Earth Explorer del Servicio Geológico
de los Estados Unidos (USGS), la cual se encuentra al ingresar en el siguiente enlace
http://earthexplorer.usgs.gov. Se usa esta plataforma, puesto que allí se encuentra un amplio
catálogo de imágenes Landsat 8 las cuales pueden ser descargadas de manera gratuita.

Para poder realizar la descarga de un clic en Login o Regístrese y posteriormente introduzca

una cuenta de usuario.

Luego ingrese a la pestaña Address/Place, escriba el nombre de su área de interés (Para este
caso será el municipio de Guaitarilla, Nariño, Colombia) y finalmente de click en el botón Show.
De esta manera el programa mostrará un punto con la ubicación de su área de interés al dar click
sobre el resultado sugerido.

Figura 1. Visualización inicial del área de interés (Municipio de Guaitarilla, Nariño, Colombia).

Además, se puede seleccionar en el rango de datos las fechas de las imágenes satelitales, en
este caso se buscarán imágenes actualizadas, por lo tanto el rango establecido ha sido entre el
 4

01/01/2015 y 27/10/2017 y se seleccionaran todos los meses. Finalmente seleccione la opción

Data Sets para buscar dentro del catálogo.

Figura 2. Despliegue de catálogo de imágenes, para este caso Landsat.

Selecciona imágenes Landsat y a su vez, Landsat Collection 1 Level-1. Para el caso de

estudio Landsat 8 OLI/TIRS C1 Level-1. (8 OLI hace referencia a las últimas imágenes lanzadas
por el satélite Landsat) y por último se da click en el botón Results.

Una vez indique la pestaña los resultados selecciona la imagen que cumpla con los requisitos
de baja cobertura de nubes y cubrimiento total del área de estudio, en este caso se seleccionará la
escena correspondiente al paquete de imágenes
ID:LC08_L1TP_010059_20170920_20171012_01_T1, del 20 de Septiembre de 2017 las cuales
cubren el área total de estudio y además no presentan nubosidad en el área de interés.
 5

Opciones de descarga

Figura 3. Selección de la escena Landsat 8.

Figura 4. Selección para descarga de imágenes Geotiff.

Para iniciar la descarga se da click en el icono Downolad y como muestra la figura 4

selecciona la opción Level-1 GeoTIFF Data Product. De esta manera iniciará el proceso de
descarga de las imágenes satelitales Landsat 8.
 6

2. Importación de imagen Landsat 8.

Al descargar la imagen satelital Landsat 8 la misma vendrá compuesta por un número de

bandas (11 en total) y un archivo MTL, el cual contiene los datos de cada banda (Metadata).
Además como se puede observar las correspondientes bandas que componen la imagen se
encontrarán en formato .TIFF.

Bandas que
componen la
imagen L8.

Archivo
MTL.

Figura 5. Archivo MTL y bandas que componen la imagen Landsat 8.

Por lo tanto, el proceso consiste en convertir cada una de las bandas al formato de ERDAS
Imagine 2015, es decir con extensión “.img”. Para iniciar el proceso, abra el programa y seleccione
Manage Data | Import Data. Diríjase a la ventana Import, en el campo Format debe seleccionar
la opción TIFF (Aunque esta suele venir por defecto), a continuación vaya a la ventana en la cual
se encuentre la imagen .TIFF a importar, que en este caso corresponde a la imagen de la banda 1
de la escena de satélite Landsat 8 y en el botón Output File seleccione el destino de salida de la
imagen con un nombre adecuado, en este caso “banda_1” puesto que se sigue el orden de la banda
que se está importando. Esta misma acción debe repetirse con las bandas 1 a 8, (se omiten las
demás bandas puesto que el software ERDAS 2015 no reconocería la posterior unión de bandas
como una imagen Landsat 8).
 7

Archivo destino en salida

formato .img.

Figura 6. Importación de bandas en formato .img.

Posteriormente, tras haber seleccionado la imagen .TIFF de entrada y nombrado la imagen

.img de salida se da click en OK y automáticamente el software iniciará el proceso.

Figura 7. Importación TIFF a .img.

 8

Figura 8. Automatización de la importación.

3. Unión de bandas

Una vez haya concluido la importación de las respectivas bandas que componen la imagen
Landsat 8, se procedera a realizar la unión de las mismas con el fin de poder trabajar toda la escena
desde un solo archivo en formato .img.

Primero, se debe abrir el programa ERDAS 2015, dirigirse a la pestaña Raster | Spectral |
Layer Stack.

Figura 9. Selección de la herramienta Layer Stack.

 9

Inmediatamente se despliega una ventana llamada Layer Selection and Stacking, allí se
cargará cada una de las bandas importadas que componen la imagen, desde la opción Input File.

Figura 10. Cuadro dialogo Layer Selection and Stacking.

Cuando se encuentre la respectiva banda en el Input File, se debe dar click en el botón Add
e inmediatamente se adicionara a una lista, será necesario tener en cuenta que la banda inicialmente
este seleccionada en Layer 1 y el Output Unsigned 16 bit, además en Output Options la opción
Union.

Listado de bandas agregadas.

Figura 11. Adicionamiento de bandas en Layer Seletion and Stacking.

 10

Como se mencionó anteriormente el proceso debe repetirse para cada una de las bandas que
componen la imagen exceptuando la Banda 8 (Pancromática), la cual al presentar una resolución
espacial distinta al resto de bandas (15 metros), genera conflictos al momento de la unión. Por lo
tanto deben adicionarse las bandas 1 a 7 en su respectivo orden.

Figura 12. Compilamiento de bandas 1 a 7.

Una vez hayan cargado las respectivas bandas se dirige a la pestaña Output File (*.img) y
selecciona en el botón Open el destino del archivo de salida con su respectivo nombre, para este
caso se llamará “Union_2017” el cual se enviará a una carpeta que ha sido titulada con el mismo
nombre. (Resulta importante en cada proceso ir introduciendo subcarpetas para mantener un orden
adecuado y evitar posteriores confusiones).

A continuación se verifica que en la pestaña Layer este seleccionado All, y de igual manera
este llena la casilla Ignore Zero in Stats.
 11

Figura 13. Generación de archivo de salida.

Figura 14. Automatización del proceso de unión de bandas.

Una vez se culmine el proceso se puede visualizar la imagen final dirigiéndose a la pestaña
File y automáticamente en documentos recientes encontrara el nombre del último archivo que el
software ha generado.
 12

Último archivo
generado.

Figura 15. Pestaña File para carga rápida de la imagen más reciente.

Cuando haya seleccionado la imagen diríjase a la pestaña Home y seleccione la opción Fit to
Frame e inmediatamente el programa cargará el último archivo generado.

Visualización
de la imagen.

Figura 16. Visualización de la imagen una vez realizada la unión de bandas.

Por último, para visualizar la imagen con la respectiva combinación de bandas diríjase a la
pestaña Multispectral y en la opción Bands seleccione Landsat 8 MS, a continuación realice la
combinación adecuada para su objeto de análisis (Para este caso al desear identificar coberturas de
la tierra se seleccionará RGB 5,6,4 idónea para análisis de vegetación).
 13

Combinación RGB
Selección
del sensor.

Figura 17. Combinación de bandas en Landsat 8.

Figura 18. Escena Landsat 8 aplicada combinación de bandas 5,6,4.

4. Fusión de imagen multibanda con imagen Banda 8 (Pancromática).

La fusión de imágenes es un proceso que tiene como finalidad mejorar la resolución espacial
de una imagen multibanda, (en este caso la última imagen generada a partir de la unión de 7 bandas)
con el fin de lograr una más alta interpretación visual de las imágenes.
 14

En el caso de la imágenes satelitales Landsat 8 se va a combinar una imagen pancromática

(Banda 8), la cual tiene como resolución espacial un tamaño por pixel de 15m x 15m, con la imagen
obtenida a partir de la unión de bandas (imagen multibanda), la cual tiene una resolución espacial
de 30m x 30m. El resultado esperado será una imagen en la cual se ha conservado la calidad
espectral de la imagen multibanda y se ha ganado la resolución espacial de la imagen pancromática.

Para ello, en primer lugar una vez abierto el programa ERDAS 2015 y cargada la imagen
obtenida a partir de la unión de bandas se dirige a la pestaña Raster y dentro de la misma se ubica
en el botón Pan Sharpen, posteriormente se selecciona la opción HPF Resolution Merge.

Figura 19. Selección de la opción HPF Resolution Merge.

 15

Figura 20. Cuadro de dialogo HPF Reolution Merge.

Una vez haya abierto el cuadro de dialogo en la opción High Resolution Input File cargue la
imagen importada de la banda 8 (Pancromática) y en la opción Multispectral Input File cargue
la imagen multibanda obtenida a partir de la unión de bandas (1 a 7), seleccione el fichero y nombre
de archivo de salida en Output File, seleccione la opción Ignore Zero in St y deje las demás
opciones por defecto por último presione el botón OK. Inmediatamente el software de manera
automática iniciará el proceso (puede tardar varios minutos).
 16

Imagen unión, multibanda (7 Bandas).

Pixel 30m x 30m.

Banda 8.
Pixel Imagen de salida.
15mx15m.

Figura 21. Cuadro de dialogo HPF Resoution Merge.

Figura 22. Automatización del proceso de fusión de imágenes.

Una vez finalice el proceso se puede visualizar la imagen final dirigiéndose a la pestaña File
y automáticamente en documentos recientes encontrará el nombre del último archivo que el
software ha generado, para este caso “fusion_2017”. Al cargar la imagen de la fusión han quedado
abiertas ambas imágenes, (Unión y fusión) en la ventana principal del programa.
 17

Imágenes
producto de los
procesos de
fusión y unión.

Figura 23. Fusión de la imagen cargada en el programa.

Para comprobar si ha mejorado la resolución espacial de la última imagen (para este caso
“fusion_2017”) se dirige al metadato de la imagen, para ello da click derecho sobre el nombre de
la imagen cargada en 2D View #1 y selecciona la opción Metadata.

Inmediatamente aparecerá un cuadro de dialogo con la información de la imagen, allí se

puede comprobar que el tamaño de sus pixeles es ahora de 15 metros x 15 metros, es decir se ha
aumentado su resolución espacial y por lo tanto el nivel de detalle.

Figura 24. Metadata de la imagen.

 18

Tamaño de cada pixel

que compone la imagen.

Figura 25. Cuadro de dialogo Image Metadata.

El siguiente paso es comprobar en una determinada zona de la imagen dicho aumento en la

resolución espacial, para ello se dirige nuevamente a la pestaña principal del programa (Home) y
selecciona la opción Swipe, esta herramienta sirve para realizar una comparación entre las dos
imágenes que se encuentran abiertas dentro del programa. Al seleccionarla el programa carga de
manera simultánea ambas imágenes, además brinda la opción de mostrar la transición (para mayor
detalle se recomienda hacer Zoom a una zona determinada de la imagen).
 19

Figura 26. Selección de la herramienta Swipe.

Además se puede hacer uso de la función Transitión Extent, para visualizar de mejor
manera los cambios efectuados en la imagen.

Figura 27. Herramienta Swipe haciendo uso de la opción "Transition Extent".

 20

De esta manera se concluye la fusión de imágenes logrando así una imagen final multibanda
con una resolución espacial de 15m x 15m que ayudará a la mejor identificación de coberturas
presentes en el área de interés.

5. Recorte de área de interés mediante Shapefile.

Este proceso tiene como finalidad realizar el recorte de una determinada zona de interés
dentro de una imagen, teniendo como base un “Shapefile”, en este caso se tomará como base el
shapefile del municipio de Guaitarilla en el departamento de Nariño, lo primero que se debe tener
en cuenta es que ambos Layer (Raster y Vectorial), es decir la imagen y el shapefile estén en el
mismo sistema de referencia geográfica. Para esto se abre tanto la imagen en este caso, la imagen
obtenida en el paso anterior “fusión_2017.img” y el Shapefile que para este caso se ha
denominado “silueta_guaitarilla”, este último se abre dando click derecho en 2D View #1, luego
en Open Vector Layer y seleccionando cargando el shape desde el lugar donde este guardado.

Desde este campo abre el shapefile del área de interés.

Figura 28. Open vector Layer, para abrir el Shapefile.

Se carga el shapefile, en este caso denominado “SILUETA.shp” e inmediatamente

aparecerá dentro de la imagen el área de interés, seguramente si tienen el mismo sistema de
referencia aparecerá en un lugar dentro de la imagen.
 21

Figura 29. Cuadro de dialogo Select Layer To Add.

El shapefile en este caso “silueta.shp”, aparecerá con un color distinto.

Shapefile cargado.

Figura 30. Carga del shapefile del área de interés.

El paso siguiente es crear un nuevo AOI Layer, para ello se da primero un click en el
shapefile ubicado sobre la imagen (este tomará un color amarillo) a continuación se dirige a la
 22

pestaña Drawing y nuevamente da click sobre el icono Paste From Selected Object.
Automáticamente se creara un nuevo AOI en base al Shapefile, (se comprueba porque el área
donde se encuentra el shapefile tomará una textura punteada).

Figura 31. Herramienta Paste From Selected Object.

Figura 32. Aoi, municipio de Guaitarilla.

 23

Una vez se ha creado el AOI, se procede a desactivar la capa del shapefile y quedara
únicamente el área punteada que contiene la información del área de interes a recortar.

Unicamente activadas las capas de AOI

y Shapefile.

Figura 33. Selección del AOI del municipio de Guaitarilla.

Ahora se debe guardar el AOI creado para ello deber hacer click sobre el mismo y
seleccionar la opción Save Layer As , darle un nombre y un fichero de destino, para este caso será
nombrado AOI_guaitarilla.

Guardado AOI

Figura 34. Guardado AOI, municipio de Guaitarilla.

 24

A continuación se dirige a la pestaña Raster, en Subset & Chip selecciona la opción Create
Subset Image. Donde se desplegará un cuadro de dialogo. Allí automáticamente el programa
pondrá como archivo de entrada la imagen con la cual se viene trabajando (es decir para este caso
fusión_2017), por lo tanto selecciona el destino de la imagen de salida (output_file), el cual se
nombrará “recorte_2017”, se dejan los demás valores por defecto, se da  sobre Ignore Zero in
Outputs y al dar click sobre el botón AOI se debe cargar el AOI creado anteriormente y dar Ok.

Figura 35. Herramienta Create Subset Image.

Figura 36. Cuadro de dialogo Create Subset Image.

Inmediatamente el programa iniciará el proceso de recorte automático.
 25

Figura 37. Recorte automático del área de interés.

Finalmente se obtendrá el recorte del área de estudio, para trabajarlo se puede apagar o
eliminar el resto de capas del visor y cargar únicamente el recorte para este caso “recorte_2017”.
Se aplica a la imagen la combinación de bandas 5,6,4 y se obtiene el recorte final.

Figura 38. Recorte final, municipio de Guaitarilla.

 26

Finalmente al dirigirse a la ventana Multispectral se podrán aplicar algunos filtros a la

imagen para mejorar calidad visual, como brillo, contraste nitidez entre otros, los cuales se dejan
al criterio de la persona que use el software.

Radiometría Brillo y Nitidez

Contraste

Figura 39. Mapa final del municipio de Guaitarilla, realizado el preprocesamiento.

 27

BIBLIOGRAFIA CONSULTADA

ERDAS®, Inc. (Julio de 2015). ERDAS Field Guide. Obtenido de The RS/GIS Laboratory

Jensen, J. R. (1996). Introductory Digital Image Processing. A Remote Sensing Perspective. USA:
USA: Prentice-Hall.

Unidad SIG. Centro de Ciencias Humanas y Sociales. (2011). Introducción a la teledetección con
ERDAS. Madrid: CISC.

USGS. (Septiembre de 2017). Earth Explorer. Obtenido de Visualización de imágenes Landsat:

http://earthexplorer.usgs.gov