Manual de Pre Procesamiento de Imágenes Satelitales Landsat 8 A Través Del Software Erdas 2015
Manual de Pre Procesamiento de Imágenes Satelitales Landsat 8 A Través Del Software Erdas 2015
Manual de Pre Procesamiento de Imágenes Satelitales Landsat 8 A Través Del Software Erdas 2015
Anexo 2: Manual de pre procesamiento de imágenes satelitales Landsat 8, mediante el uso del
software ERDAS 2015.
NOVIEMBRE 2017.
2
CONTENIDO
3. Unión de bandas…………………………………………………………………………….......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.
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
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
Bandas que
componen la
imagen L8.
Archivo
MTL.
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
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.
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.
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.
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.
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
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.
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.
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
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.
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
Banda 8.
Pixel Imagen de salida.
15mx15m.
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.
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.
Además se puede hacer uso de la función Transitión Extent, para visualizar de mejor
manera los cambios efectuados en la imagen.
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.
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.
Shapefile cargado.
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).
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.
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
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.
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.
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.