Función Zonal Histogram de ArcGIS en R
En este tutorial explicare una aplicación de la función “ zonal histogram” de ArcGIS en R, el objetivo, será calcular el área de las coberturas para los municipios del Valle del Cauca a partir de la capa de cobertura de la NASA y la división político administrativa del SIGOT. Los datos a utilizar, en coordenadas geográficas WGS 1984, son los siguientes (disponibles aquí ): 1. Cobertura de la tierra reclasificada de la NASA. 2. Municipios del Valle del Cauca A continuación los pasos: 1. Cargamos las librerías a utilizar, en caso que usted no las tenga instalada podría utilizar install.packages(“nameLibrary”) # Cargamos librerias require(raster) require(rgdal) require(foreign) require(ggplot2) require(dplyr) 2. Cargamos los archivos, que son el raster, la tabla de la leyenda del raster, y el shapefile de municipios. path <- "E:/R/_tutoriales/_tabulateArea/_blog/_data" cob <- raster(paste0(path, "/_raste