Mapa de Cobertura Forestal para la República de Guatemala 2010 y Dinámica de Cobertura Forestal Periodo 2006-2010 Alejandro Suarez A Para la elaboración del mapa nacional de cobertura forestal 2010 toma como base las imágenes satelitales obtenidas por los satélites LANDSAT 5 y 7 dependiendo la región del país( Path 21 Row 50, Path 20 Row 50 de la serie LANDSAT 7 y Path 19 Row 50 de la serie LANDSAT 5), para el periodo de estudio. Luego de haber realizado procesos de, verificación de la Precisión geométrica y reproyección, Conversión a reflectancia y normalización radiométrica, se procedió a realizar la clasificación de las imágenes para posteriores ediciones manuales. Para poder realizar esta, se tomaron en cuenta factores de ISODATAS, puntos de control y verificación y por último se evaluó la precisión de la exactitud de la clasificación para cada imagen. Imagen P21R50 clasificada P21R50 primera clasificación P21R50 imagen con ediciones manuales y AOIs Figura No.1 imagen P21R50 recortada y clasificada. Obtención de puntos de control y verificación Este procedimiento permite evaluar, mediante el uso de ortofotos de alta definición, la cobertura “real” en un determinado punto geográfico. A través de un proceso de selección aleatoria, se distribuyó un total de 30,000 puntos sobre el territorio nacional. Estos puntos sirvieron como medio de verificación del tipo de cobertura del suelo al evaluar cada uno de manera individual sobre las ortofotos. Con cada punto se desplegaron en pantalla tres circunferencias con distintos diámetros: la mayor con 350 m, la intermedia de 60 m y la central y más pequeña con 30 m. Las circunferencias de 60 m y 30 m (también llamados “buffers”) fueron utilizadas para determinar si un punto tenía homogeneidad de cobertura del suelo. De existir homogeneidad en los buffers, el punto se consideró “válido” para ser utilizado y se procedió a clasificar su tipo de cobertura. Corte imagen P21R50 clasificada para nubes P21R50 corte primera clasificación P21R50 corte imagen ediciones manuales y AOIs Figura No.2 Corte de imagen P21R50 para edición de nubes al ser del departamento. Clasificación para edición manual La clasificación de las imágenes satelitales se hizo con un procedimiento supervisado utilizando para ello los puntos de control clasificados sobre las ortofotos que sirven como puntos de entrenamiento para que la computadora pueda clasificar toda la imagen. Estos puntos de entrenamiento en efecto definen áreas donde el uso del suelo es conocido y de donde se extrae información espectral para definir las llamadas firmas espectrales. Estas firmas determinan de forma única el espectro de luz reflejada de los usos del suelo a clasificar Se definieron firmas espectrales para cada uno de los usos del suelo y con esa información se procedió a realizar la clasificación de todas las imágenes. Este proceso de clasificación automatizado tiene sus limitantes principalmente porque tiende a confundir áreas con firmas espectrales muy parecidas. Esto hace necesaria una edición manual del producto clasificado para reclasificar áreas con error. Esta edición manual se hizo convirtiendo a formato vectorial las imágenes ya clasificadas para luego revisar individualmente los vectores que tenían posibilidad de haber sido clasificados erróneamente. Imagen P20R50 clasificada P20R50 primera clasificación P20R50 imagen con ediciones manuales y AOIs Figura No.3 Imagen P20R50 sin y con ediciones manuales. Para el presente mapa se utilizaron las siguientes clases para las posteriores ediciones: a) b) c) d) e) f) Clase 1 Bosque Clase 2 Matorrales Clase 3 Plantaciones Forestales Clase 4 Plantaciones Herbáceas Clase 5 Plantaciones Bajo Sombra (café) Clase 6 No Bosque Como se puede ver en las distintas imágenes, el proceso de clasificación no es del todo certero, es por eso que se utiliza otro proceso para darle mayor precisión a la clasificación. En este proceso se re clasifican las áreas que no fueron clasificadas correctamente por distintas razones. Estos procesos se realizan por medio de AOI´s una herramienta del programa ERDAS. Un claro ejemplo de esto son las imágenes dentro de este documento. Los AOI´s se relazando sobre la imagen RASTER y no sobre la clasificada, ya que en la imagen RASTER por medio de las distintas combinaciones de bandas y firmas espectrales, se pueden detectar y diferenciar entre las distintas clases previamente descritas. Imagen P19R50 clasificada P19R50 primera clasificación P19R50 imagen con ediciones manuales y AOIs Figura No.4 Imagen P19R50 sin y con ediciones manuales. A la primera clasificación que salió del proceso, se le aplico un “recode”, que es una herramienta de ERDAS, en donde se agruparon las clases más parecidas y problemáticas para poder tener mejor visualización de las mismas. Al final de este proceso se obtuvieron solo cuatro clases en vez de seis, ya que de esta manera es más fácil trabajar las ediciones manuales y poder editar lo que está confundido como bosque y viceversa. Clase 1: Bosque y plantación forestal Clase 2: No bosque Clase 3: Agua Clase 4: Matorral Al tener ya las imágenes con sus cuatro clases principales, se procedió a hacer los cortes de cada imagen con un “buffer” de 500m para cada una. Luego de tener ya les imágenes recortadas para el país, se procedió a editar manualmente y cambiar lo que estaba confundido en la primera versión de la clasificación. A partir de esto se fueron haciendo ediciones manuales en donde se cambio el matorral (color amarillo Figura No.4) ya se bosque o no bosque (color verde “bosque” o amarillo pálido “no bosque”) dependiendo lo que diga la firma espectral de la imagen RASTER, luego de varios procesos a los cuales fue sometida con anterioridad. Evaluación de exactitud de la clasificación Se hizo la evaluación de la exactitud de la clasificación de las imágenes usando el segundo grupo de puntos de control que fueron clasificados sobre las ortofotos independientemente por un grupo externo. Comparando estos puntos contra el mapa clasificado se generaron estadísticas para estimar el error del mapa, incluyendo matrices de contingencia con porcentaje de coincidencia entre la clasificación de las imágenes y la clasificación independiente de las ortofotos. Figura No.5 Ejemplo de evaluación de exactitud para la imagen P20R50 A todas las imágenes que se clasificaron se les realizo este procesos de evaluación de exactitud, dependiendo si la imagen clasificada no llegaba a un 80% de precisión con respecto a los puntos de control, se tuvo que revisar que proceso no se realizo correctamente o si enserio la imagen había tenido mucha dinámica (cambio en la cobertura donde se encuentran los puntos de control) y de ser este el caso, se reevaluaron los puntos y se procedió a realizar la clasificación y verificación correspondiente. Se trato del que el porcentaje de la primera clasificación cruda estuviese arriba de 80%. Ya que con eso se tiene mayor certeza para empezar con las ediciones manuales. Esta evaluación de exactitud también se realiza al final de las ediciones manuales tratando de que este arriba de un 95%. Bibliografía Gonzalez, Rafael C. and Paul Wintz. 1977. Digital Image Processing. Reading, Massachusetts: Addison-Wesley Publishing Company. Hodgson, M. and B. Shelley. 1994. “Removing the Topographic Effect in Remotely Sensed Imagery.” The ERDAS Monitor, Vol. 6, No. 1: 4-6. Pratt, William K. 1991. Digital Image Processing. New York: John Wiley & Sons, Inc. Sabins, Floyd F., Jr. 1987. Remote Sensing Principles and Interpretation. New York: W. H. Freeman and Co. Schowengerdt, Robert A. 1983. Techniques for Image Processing and Classification in Remote Sensing. New York: Academic Press Universidad del Valle de Guatemala, Instituto Nacional de Bosques, Consejo Nacional de Áreas Protegias, Universidad Rafael Landivar. 2011. Mapa de Cobertura Forestal de Guatemala 2006 y Dinámica de la Cobertura Forestal 2001-2006. Guatemala, 97pp.
© Copyright 2024