ESTHER QUIÑONES LUNA
Aquí se muestra el código para elaborar mapas temáticos en R-STUDIOS, puedes descargar la base de datos de producción bruta total del Censos Económicos 2019 de Instituto Nacional de Estadística y Geografía (INEGI) y la capa de Entidades de México. La actividad económica es la edificación de vivienda unifamiliar (La Clasificación Industrial de América del Norte 2018 (SCIAN) es 236111).
CODIGO
###Mapas Censos Económicos
###Esther Quiñones Luna
##Paqueteria
library(udunits2)
library(units)
library(ggspatial)
library(ggplot2)
library(maps)
library(sp)
library(rgdal)
library(tibble)
library(tidyr)
library(readr)
library(purrr)
library(dplyr)
library(raster)
library(rgeos)
library(gpclib)
library(maptools)
library(dplyr)
library(stringr)
##Base de datos
setwd("C:\\")
entidad <- readOGR(dsn = "C:\\", layer="ENTIDAD_WG84")
enti <- read.csv('DMS.csv', encoding = "UTF-8")
enti <- rename(enti, CVEGEO=X.U.FEFF.CVENT)
enti$CVEGEO<- str_pad(string = enti$CVEGEO, width = 2, side = "left", pad = 0)
entidad <- merge(x=entidad, y=enti, by="CVEGEO", all.x=TRUE, incomparables=FALSE)
##Mapa opción 1
pal <- brewer.pal(5,"PuBuGn")
min(entidad$c236111)
max(entidad$c236111)
breaks_qt <- classIntervals(entidad$c236111, n=5, style ="quantile")
br <- breaks_qt$brks
offs <- 0.000001
br[1] <- br[1]-offs
br[length(br)] <- br[length(br)]+offs
entidad$c236111_bracket <- cut_interval(entidad$c236111, n=5)
entidad$c236111_bracket
spplot(entidad,"c236111_bracket", scales = list(draw=TRUE), main="Producción Brura Total, Clase 236111", col.regions=pal)
require(grid)
grid.text(expression("Millones de pesos"), x=unit(0.90, "npc"), y=unit(0.65, "npc"))
#Mapa opción 2
pl <- hsv(seq(0,1 - 1/20,length.out = 20), 0.5 , 1)
spplot(entidad, "c236111", sp.layout = spl, scales = list(draw=TRUE), main="Producción por Estado", col.regions=pl, cex.legend=6, cex.axis=2)
lbls <- as.character(paste(entidad$CVE_ENT))
spl <- list('sp.text', coordinates(entidad), lbls, cex=0.95)
spplire(grid)
grid.text(expression("Millones de pesos"), x=unit(0.90, "npc"), y=unit(0.5, "npc"), rot=90, gp=gpar(fontsize=8))
BASES DE DATOS
Comments