Added first project.

Configured renv.

.Rprofile

@@ -0,0 +1 @@
+source("renv/activate.R")

.gitignore

@@ -0,0 +1,4 @@
+.Rproj.user
+.Rhistory
+.RData
+.Ruserdata

README.md

@@ -0,0 +1,3 @@
+# Proyectos: Sistemas de Información Geográfica
+
+Este repositorio contiene el código para cada uno de los poryectos desarrollados en la materia de Sistemas de Información Geográfica de la Universidad del Caribe en Cancún, Q. Roo, México.

incidencia_delictiva_nac/destdv250k_2gw/destdv250k_2gw.prj

@@ -0,0 +1 @@
+GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

incidencia_delictiva_nac/destdv250k_2gw/destdv250k_2gw.xml

@@ -0,0 +1 @@
+<?xml version="1.0" encoding="UTF-8"?><metadata xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://www.conabio.gob.mx/informacion/gis/schemas/fgdc/fgdc-std-001-1998.xsd"><idinfo><citation><citeinfo><origin>Comisi&#x00F3;n Nacional para el Conocimiento y Uso de la Biodiversidad (2008). &apos;Divisi&#x00F3;n Pol&#x00ED;tica Estatal&apos;. Versi&#x00F3;n 2. Escala 1:250000. Modificado de Conjunto de Datos vectoriales y toponimia de la carta topogr&#x00E1;fica. Serie III. Instituto Nacional de Estad&#x00ED;stica, Geograf&#x00ED;a e Inform&#x00E1;tica (2003-2004). Marco Geoestad&#x00ED;stico Municipal, Instituto Nacional de Estad&#x00ED;stica, Geograf&#x00ED;a e Inform&#x00E1;tica (2005). Escala 1:250000. M&#x00E9;xico.</origin><pubdate>20081028</pubdate><title>Divisi&#x00F3;n pol&#x00ED;tica estatal de M&#x00E9;xico 1:250000. Versi&#x00F3;n 2</title><geoform>Shapefile. Formato vectorial compuesto por 4 archivos (shp, shx, dbf,prj)</geoform><othercit>250000</othercit></citeinfo></citation><descript><abstract>El mapa representa los l&#x00ED;mites de los Estados y el contorno de la Rep&#x00FA;blica Mexicana,. Estos l&#x00ED;mites son la uni&#x00F3;n del material cartogr&#x00E1;fico indicado en la cita de la informaci&#x00F3;n, el resultado es un mapa con el datalle de la escala 1:250000 manejada por el INEGI.</abstract><purpose>No conocido</purpose><supplinf>Este mapa fue elaborado en CONABIO con informaci&#x00F3;n cartogr&#x00E1;fica proporcionada por el INEGI. El mapa contiene informaci&#x00F3;n sobre nombres de los Estados y Capitales, seg&#x00FA;n INEGI.</supplinf></descript><timeperd><timeinfo><sngdate><caldate>00000000</caldate></sngdate></timeinfo><current>Octubre de 2008</current></timeperd><status><progress>Complete</progress><update>No planeado</update></status><spdom><bounding><westbc>-118.366485</westbc><eastbc>-86.710742</eastbc><northbc>32.718729</northbc><southbc>14.534964</southbc></bounding></spdom><keywords><theme><themekt>ESTRUCTURA</themekt><themekey>dipol</themekey><themekey>estata</themekey><themekey>destdv250k_2gw</themekey></theme><theme><themekt>CONABIO</themekt><themekey>Divisi&#x00F3;n pol&#x00ED;tica Estatal 1:250000 Divisi&#x00F3;n pol&#x00ED;tica estatal Entidades federativas Estados Rep&#x00FA;blica Mexicana </themekey></theme><place><placekt>CONABIO</placekt><placekey>Aguascalientes Baja California Baja California Sur Campeche Chiapas Chihuahua Coahuila Colima Distrito Federal Durango Estado de M&#x00E9;xico Estatal Guanajuato Guerrero Hidalgo Jalisco Michoac&#x00E1;n de Ocampo Morelos Nayarit Nuevo Le&#x00F3;n Oaxaca Puebla Quer&#x00E9;taro de Arteaga Quintana Roo Rep&#x00FA;blica Mexicana San Luis Potos&#x00ED; Sinaloa Sonora Tabasco Tamaulipas Tlaxcala Veracruz Ignacio de la Llave Yucat&#x00E1;n Zacatecas </placekey></place></keywords><accconst>Sin restricciones</accconst><useconst>Sin restricciones</useconst><browse><browsen>http://www.conabio.gob.mx/informacion/gis/layouts/destdv250k_2gw.png</browsen><browsed>Mapa ilustrativo de Divisi&#x00F3;n pol&#x00ED;tica estatal de M&#x00E9;xico 1:250000. Versi&#x00F3;n 2</browsed><browset>png</browset></browse><native>SIG ARC-INFO, versi&#x00F3;n 9.2, UNIX</native></idinfo><dataqual><logic>Gabinete</logic><complete>El mapa se obtuvo a partir de la uni&#x00F3;n de dos mapas: 1) divisi&#x00F3;n pol&#x00ED;tica estatal y 2) l&#x00ED;mite nacional en Arc/Info.</complete><lineage><srcinfo><srccite><citeinfo><origin>No conocido</origin><pubdate>00000000</pubdate><title>No conocido</title></citeinfo></srccite><srcscale>250000</srcscale><typesrc>Digital</typesrc><srctime><timeinfo><sngdate><caldate>00000000</caldate></sngdate></timeinfo><srccurr>No conocido</srccurr></srctime><srccitea>No conocido</srccitea><srccontr>Instituto Nacional de Estad&#x00ED;stica Geograf&#x00ED;a e Inform&#x00E1;tica. (2005). Marco Geoestad&#x00ED;stico Municipal, Escala 1:250000. D.F. M&#x00E9;xico.</srccontr></srcinfo><srcinfo><srccite><citeinfo><origin>No conocido</origin><pubdate>00000000</pubdate><title>No conocido</title></citeinfo></srccite><srcscale>250000</srcscale><typesrc>Digital</typesrc><srctime><timeinfo><sngdate><caldate>00000000</caldate></sngdate></timeinfo><srccurr>No conocido</srccurr></srctime><srccitea>No conocido</srccitea><srccontr>Instituto Nacional de Estad&#x00ED;stica, Geograf&#x00ED;a e Inform&#x00E1;tica (2003-2004). Conjunto de Datos vectoriales y toponimia de la carta topogr&#x00E1;fica. Serie III. Escala 1:250 000. D.F. M&#x00E9;xico.</srccontr></srcinfo><procstep><procdesc>La divisi&#x00F3;n estatal se obtuvo del MGM, 2005 a trav&#x00E9;s de un dissolve; para el l&#x00ED;mite nacional se realizo un proceso de conversi&#x00F3;n de formato para trabajar todo en cobertura. El trabajo de edici&#x00F3;n y conversi&#x00F3;n se hizo en Arc/info, tomando como base para la revisi&#x00F3;n material cartogr&#x00E1;fico de la misma escala y de la misma Instituci&#x00F3;n; por &#x00FA;ltimo la revisi&#x00F3;n y edici&#x00F3;n de etiquetas (nombres de los estados y capitales) se realiz&#x00F3; en Arcview.</procdesc><procdate>00000000</procdate></procstep></lineage></dataqual><spdoinfo><direct>Vector</direct><ptvctinf><sdtsterm ><sdtstype>G-polygon</sdtstype><ptvctcnt>1213</ptvctcnt></sdtsterm></ptvctinf></spdoinfo><spref><horizsys><geograph><latres>0.0174532925199433</latres><longres>0.0174532925199433</longres><geogunit>Decimal degrees</geogunit></geograph><geodetic><horizdn>WGS84</horizdn><ellips>WGS84</ellips><semiaxis>6378137.000000</semiaxis><denflat>298.257224</denflat></geodetic></horizsys></spref><eainfo><detailed><enttyp><enttypl>DESTDV250K_2GW.DBF</enttypl><enttypd>Divisi&#x00F3;n pol&#x00ED;tica estatal extra&#x00ED;do de datos vectoriales. Escala 1:250000.</enttypd><enttypds>No conocido</enttypds></enttyp><attr><attrlabl>Area</attrlabl><attrdef>Num&#x00E9;rico</attrdef><attrdefs>Calculado por Arc-Info</attrdefs><attrdomv><udom>&#x00C1;rea del pol&#x00ED;gono asociado a una Entidad Federativa</udom></attrdomv><attrvai><attrva>1</attrva><attrvae>Grados decimales al cuadrado</attrvae></attrvai></attr><attr><attrlabl>Perimeter</attrlabl><attrdef>Num&#x00E9;rico</attrdef><attrdefs>Calculado por Arc-Info</attrdefs><attrdomv><udom>Per&#x00ED;metro del pol&#x00ED;gono asociada a una Entidad Fedrativa</udom></attrdomv><attrvai><attrva>1</attrva><attrvae>Grados decimales</attrvae></attrvai></attr><attr><attrlabl>Num_edo</attrlabl><attrdef>Car&#x00E1;cter</attrdef><attrdefs>INEGI</attrdefs><attrdomv><udom>Identificador &#x00FA;nico relacionado con la Entidad Federativa</udom></attrdomv></attr><attr><attrlabl>Entidad</attrlabl><attrdef>Car&#x00E1;cter</attrdef><attrdefs>INEGI</attrdefs><attrdomv><udom>Nombre propio de la Entidad Federativa</udom></attrdomv></attr><attr><attrlabl>Capital</attrlabl><attrdef>Car&#x00E1;cter</attrdef><attrdefs>INEGI</attrdefs><attrdomv><udom>Nombre propio de la capital de la Entidad Federativa</udom></attrdomv></attr><attr><attrlabl>Rasgo_geog</attrlabl><attrdef>Car&#x00E1;cter</attrdef><attrdefs>INEGI</attrdefs><attrdomv><udom>Presencia de una isla</udom></attrdomv></attr></detailed></eainfo><distinfo><distrib><cntinfo><cntorgp><cntorg>CONABIO</cntorg><cntper>Jose Manuel Davila Rosas</cntper></cntorgp><cntpos>Subdirector de Sistemas de Informaci&#x00F3;n Geogr&#x00E1;fica</cntpos><cntaddr><addrtype>mailing and physical address</addrtype><address>Liga Perif&#x00E9;rico Insurgentes Sur Num. 4903 1er piso. Col. Parques del Pedregal. Delg. Tlalpan</address><city>Ciudad de M&#x00E9;xico</city><state>D.F.</state><postal>14010</postal><country>M&#x00E9;xico</country></cntaddr><cntvoice>50045012</cntvoice><cntfax>50044931</cntfax><cntemail>[email protected]</cntemail><hours>Lunes a Viernes de las 8 hasta 18 hrs.</hours></cntinfo></distrib><resdesc>destdv250k_2gw</resdesc><distliab>No se permite utilizar estos datos con fines lucrativos y se debe citar la fuente del mapa y a CONABIO. La Subdirecci&#x00F3;n de Sistemas de Informaci&#x00F3;n Geogr&#x00E1;fica de la CONABIO solo se hace responsable de la informaci&#x00F3;n presentada en el mapa y se deslinda del uso que se le de a la misma.</distliab><stdorder><digform><digtinfo><formname>ESRI Shapefile(SHP)</formname><formverd>199807</formverd><filedec>compactado en ZIP</filedec><transize>10.17</transize></digtinfo><digtopt><onlinopt><computer><networka><networkr> http://www.conabio.gob.mx/informacion/gis/maps/geo/destdv250k_2gw.zip</networkr></networka></computer><accinstr>Coordenadas geogr&#x00E1;ficas WGS84. (EPSG:4326). Formato vectorial de almacenamiento de elementos geogr&#x00E1;ficos y sus atributos asociados. </accinstr></onlinopt><onlinopt><computer><networka><networkr> http://www.conabio.gob.mx/informacion/gis/maps/ccl/destdv250k_2gw_c.zip</networkr></networka></computer><accinstr>Proyecci&#x00F3;n C&#x00F3;nica Conforme de Lambert WGS84. (EPSG:48402). Formato vectorial de almacenamiento de elementos geogr&#x00E1;ficos y sus atributos asociados. </accinstr></onlinopt></digtopt></digform><digform><digtinfo><formname>Google Earth KML</formname><formvern>2.2</formvern><filedec>compactado como KMZ</filedec></digtinfo><digtopt><onlinopt><computer><networka><networkr>  http://www.conabio.gob.mx/informacion/explorer/wms/kml?layers=destdv250k_2gw&amp;mode=download&amp;legend=true</networkr></networka></computer><accinstr>Coordenadas geogr&#x00E1;ficas WGS84. (EPSG:4326). Los archivos KML representan datos geogr&#x00E1;ficos para Google Earth, a menudo suelen distribuirse comprimidos como KMZ.</accinstr></onlinopt></digtopt></digform><digform><digtinfo><formname>Servicio de im&#x00E1;genes (WMS)</formname><formvern>1.1.1</formvern></digtinfo><digtopt><onlinopt><computer><networka><networkr>http://www.conabio.gob.mx/informacion/explorer/wms?Format=image.png&amp;request=GetMap&amp;version=1.1.1&amp;width=800&amp;height=424&amp;srs=EPSG:4326&amp;bbox=-120.0,13.8,-85.1,33.5&amp;layers=destdv250k_2gw</networkr></networka></computer><accinstr>Coordenadas geogr&#x00E1;ficas WGS84. (EPSG:4326). Servicio Web Map Service. Los mapas de datos referenciados espacialmente producidos se generan normalmente en un formato de imagen como PNG,.png o JPEG, pueden se vistos en un SIG </accinstr></onlinopt></digtopt></digform><fees>Acceso libre sin ning&#x00FA;n cargo al p&#x00FA;blico en general.</fees></stdorder><techpreq>Tener Arc-Info, Arcview o sistemas compatibles, en caso que existan se pueden visualizar en Google Earth KML y/o conexi&#x00F3;n al servicio WMS.</techpreq></distinfo><metainfo><metd>20081028</metd><metc><cntinfo><cntorgp><cntorg>CONABIO</cntorg></cntorgp><cntpos>Norma Moreno</cntpos><cntaddr><addrtype>mailing and physical address</addrtype><address>Liga Perif&#x00E9;rico Insurgentes Sur Num. 4903 1er piso. Col. Parques del Pedregal. Delg. Tlalpan</address><city>Ciudad de M&#x00E9;xico</city><state>D.F.</state><postal>14010</postal><country>M&#x00E9;xico</country></cntaddr><cntvoice>50045018 50044938</cntvoice><cntfax>50044931</cntfax><cntemail>[email protected]</cntemail><hours>Lunes a Viernes de 8 a 18 hrs.</hours></cntinfo></metc><metstdn>FGDC Content Standards for Digital Geospatial Metadata, 1998</metstdn><metstdv>FGDC-STD-001-1998 version 2.0</metstdv></metainfo></metadata>

incidencia_delictiva_nac/practica_1.R

@@ -0,0 +1,135 @@
+# Práctica: Incidencia Delictiva Nacional
+library(RODBC)
+library(readxl)
+library(tidyr)
+library(dplyr)
+library(stringr)
+library(sf)
+library(cartography)
+setwd("incidencia_delictiva_nac/")
+
+# Incidencia delictiva del fuero común
+# Fuente: Secretariado Ejecutivo del Sistema Nacional de Seguridad Pública
+seg_pub_db <- odbcConnectExcel2007("Estatal-Delitos - diciembre 2019.xlsb")
+
+seg_pub <- as_tibble(sqlFetch(seg_pub_db, "Hoja1$")) %>%
+  mutate_at(vars(Clave_Ent), list(~as.character(.))) %>%
+  mutate_if(is.numeric, as.integer) %>%
+  mutate_at(vars(Clave_Ent), list(~str_pad(., width = 2, side = "left", pad = "0")))
+
+odbcClose(seg_pub_db)
+rm(seg_pub_db)
+
+# Tamaño de población por estado
+# Fuente: INEGI Encuesta Intercensal 2015
+inegi <- read_excel("01_poblacion.xls", sheet = 4, range = "A7:E3573")
+
+# Polígonos estatales
+# Fuente: CONABIO
+polig_est <- read_sf("destdv250k_2gw/destdv250k_2gw.shp") %>%
+  select(-one_of("AREA", "PERIMETER", "COV_", "COV_ID", "CAPITAL", "RASGO_GEOG")) %>%
+  rename(Entidad = ENTIDAD, Clave_Ent = NUM_EDO) %>%
+  group_by(Clave_Ent, Entidad) %>%
+  summarise() %>%
+  mutate_at(vars(Entidad), list(~str_to_title(., locale = "es")))
+
+# Filtro: Homocidios dolosos por arma de fuego
+# Datos agrupados por año y por entidad
+hom_dol_arma_fuego <- seg_pub %>%
+  # Obtención de los homocidios dolosos por arma de fuego
+  filter(`Subtipo de delito` == "Homicidio doloso", Modalidad == "Con arma de fuego") %>%
+  select(-one_of("Bien jurídico afectado", "Tipo de delito", "Subtipo de delito", "Modalidad")) %>%
+  pivot_longer(
+    -one_of("Año", "Clave_Ent", "Entidad"),
+    names_to = "Mes",
+    values_to = "Total"
+  ) %>%
+  # Agrupación de los datos por año y por entidad
+  group_by(`Año`, Clave_Ent, Entidad) %>%
+  summarise_at(vars(Total), list(~sum(.))) %>%
+  rename(Hom_Dol_Arma_Fuego = Total)
+
+# Filtro: Tamaño de población por entidad
+tam_pob <- inegi %>%
+  filter(
+    `Entidad federativa` != "Estados Unidos Mexicanos",
+    Sexo == "Total",
+    `Grandes grupos de edad` == "Total",
+    Estimador == "Valor"
+  ) %>%
+  separate(
+    `Entidad federativa`, 
+    into = c("Clave_Ent", "Entidad"),
+    sep = 3
+  ) %>%
+  mutate_at(vars(Clave_Ent), list(~str_trim(.))) %>%
+  rename(Poblacion = `Población total`) %>%
+  select(-one_of("Sexo", "Grandes grupos de edad", "Estimador"))
+
+# Cálculo de la tasa de homicidios dolosos con arma de fuego por entidad por año
+tasa_hom <- hom_dol_arma_fuego %>%
+  full_join(tam_pob, by = "Clave_Ent") %>%
+  mutate(Tasa_Hom = Hom_Dol_Arma_Fuego / Poblacion * 100000) %>%
+  select(`Año`, Clave_Ent, Tasa_Hom)
+
+# Unión con los poligonos estatales
+sf_tasa_hom <- polig_est %>%
+  ungroup() %>%
+  full_join(tasa_hom, by = "Clave_Ent") %>%
+  select(`Año`, Clave_Ent, Entidad, Tasa_Hom)
+
+### Creación de los mapas temáticos ###
+
+# Mapa base de OpenStreetMap
+mex.osm <- getTiles(
+  x = polig_est,
+  type = "osm",
+  zoom = 5,
+  crop = TRUE
+)
+
+# División de la información por años
+tasa_hom_anio <- list(
+  "2015" = filter(sf_tasa_hom, `Año` == 2015),
+  "2016" = filter(sf_tasa_hom, `Año` == 2016),
+  "2017" = filter(sf_tasa_hom, `Año` == 2017),
+  "2018" = filter(sf_tasa_hom, `Año` == 2018),
+  "2019" = filter(sf_tasa_hom, `Año` == 2019)
+)
+
+# Cortes y paleta de colores para el Choropleth
+pal <- carto.pal(pal1 = "green.pal", n1 = 1, pal2 = "orange.pal", n2 = 9, transparency = FALSE)
+breaks <- getBreaks(seq(0, 100, by = 10), n = 10, method = "equal")
+
+# Creación del PDF con los 5 mapas temáticos
+pdf("inc_delic_anio_mex.pdf", width = 7, height = 5.3)
+for (info_anio in tasa_hom_anio) {
+  # Creación del mapa temático
+  tilesLayer(x = mex.osm)
+  choroLayer(
+    x = info_anio,
+    var = "Tasa_Hom",
+    breaks = breaks,
+    col = pal,
+    border = "black",
+    lwd = 0.2,
+    legend.pos = "n",
+    add = TRUE
+  )
+  legendChoro(
+    pos = "left",
+    title.txt = "Tasa de\nhomicidios",
+    breaks = breaks,
+    col = pal,
+    nodata = FALSE
+  )
+  layoutLayer(
+    title = paste0("Tasa de homicidios dolosos por arma de fuego (", unique(info_anio$`Año`), ")"),
+    sources = "Fuentes: INEGI (2016) y Secretariado Ejecutivo del SNSP (2020)\n© OpenStreetMap contributors.\nTiles style under CC BY-SA, www.openstreetmap.org/copyright.",
+    author = paste0("Por: Fernando Gómez Perera, usando cartography ", packageVersion("cartography")),
+    frame = TRUE, north = TRUE, tabtitle = TRUE
+  )
+}
+dev.off()
+
+setwd("../")

proyectos-sig.Rproj

@@ -0,0 +1,13 @@
+Version: 1.0
+
+RestoreWorkspace: Default
+SaveWorkspace: Default
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 2
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX

renv.lock

@@ -0,0 +1,19 @@
+{
+  "R": {
+    "Version": "3.6.2",
+    "Repositories": [
+      {
+        "Name": "CRAN",
+        "URL": "https://cran.rstudio.com"
+      }
+    ]
+  },
+  "Packages": {
+    "renv": {
+      "Package": "renv",
+      "Version": "0.9.3",
+      "Source": "Repository",
+      "Repository": "CRAN"
+    }
+  }
+}

renv/.gitignore

@@ -0,0 +1,3 @@
+library/
+python/
+staging/