Repositorio con código de magda
Se utiliza el software QGIS y R con RStudio. También el paquete GDAL por terminal. En los sitios correspondientes se indica cómo instalarlos.
Algunos paquetes de R se cargan automáticamente al abrir el proyecto. Verificar la instalación previa de stringr, data.table, lubridate, readr, dtplyr, tidyr, ggplot2, dplyr, magrittr.
En particular, el paquete readxl tiene que estar actualizado, pues algunas funciones rompen el comportamiento.
Además encontramos funciones que se cruzan entre paquetes. Por ejemplo, filter se confunde de paquetes, por lo que la asignamos en el archivo .Rprofile: filter <- dplyr::filter contains <- dplyr::contains
Descargamos datos de los siguientes sitios.
CNBV, INEGI, NOAA, BIE, y datos de referencia.
Vínculo: http://portafolioinfo.cnbv.gob.mx/PUBLICACIONES/IO/Paginas/bm.aspx
Archivos desde abril 2011, hasta el más reciente.
Se puede correr el archivo de R 0_descargar_cnbv.R que los guarda en ../data/cnbv/raw/
Vínculo: http://www.inegi.org.mx/sistemas/bie/ Camino: - Indicadores económicos de coyuntura > ITAEE > Series desestacionalizadas > todas ... Desde 2010, como CSV, y demás opciones tal cual. Guárdala en: ../data/bie/raw/itaee.csv Una vez guardada, borra la selección anterior.
Camino: - Cuentas nacionales > PIBE > Act. Eco. y Ent. Federativa > Precios corrientes 2008 > Total de Act. Eco. Desde 2010, como CSV, y demás opciones tal cual. Guardar en: ../data/bie/raw/pibe.csv
Marco Geoestadístico:
https://www.inegi.org.mx/temas/mg/default.html#Descargas
http://www.inegi.org.mx/geo/contenidos/geoestadistica/m_geoestadistico.aspx Seleccionar Datos vectoriales > descarga. Descomprimir en: ../data/inegi/ Y ejecutar el siguiente comando para copiar los datos $ mv ../data/inegi/conjunto_de_datos/* ../data/inegi/marco_geo/raw/
Visita el sitio de la NOAA, https://ngdc.noaa.gov/eog/viirs/ *** Cambió dirección a https://eogdata.mines.edu/download_dnb_composites.html
VIIRS DNB Nighttime Lights > 2014 > Monthly > 201412 > Tile1_75N180W > VCMSLCFG Guardar en: ../data/viirs/raw/ (Puede tardar varios minutos)
Descomprimir el archivo y renombrar el que termina con "avg_rade0.tif" como: ../data/viirs/raw/luminosidad_tile1.tif
Este archivo tiene datos de casi medio hemisferio norte, el siguiente comando recorta la información de la República Mexicana en sistema Linux.
En una misma carpeta poner archivos vectoriales de laás reas geoestadísticas estatales
y el archivo a procesar
$ gdalwarp -q -cutline areas_geoestadisticas_estatales.shp -crop_to_cutline --config GDALWARP_IGNORE_BAD_CUTLINE YES "luminosidad_tile1_2018.tif" "luminosidad_mexico.tif"
Los datos de la NOAA se procesan aparte en QGIS. Topamos a 175 debido a zonas con luminosidad extrema. (Ver ficha técnica)
- Empezar nuevo proyecto en QGIS como: ../data/qgis/imco_magda.qgs
- Abrir capa de ráster del .tif generado en: ../data/viirs/luminosidad_mexico.tif
- Renombrar como viirs
- Aplicar raster-calculator: viirs * ( viirs <= 175),
- guardar en ../data/viirs/processed/tope_175.tif
- Abrir capas vectoriales de .shp en ../data/inegi/marco_geo/raw/
- poligonos_localidades_urbanas_y_rurales.shp
- areas_geoestadisticas_municipales.shp
- Guardar capa de localidades y municipios con la misma proyección del ráster del INEGI.
- Esto equivale a cambiar el sistema de coordenadas de referencia CRS (en inglés): CRS = EPSG:4326 - WGS 84,
- botón derecho sobre cada capa vectorial > guardar como > cambiar CRS Con el mismo formato .shp en ../data/inegi/marco_geo/processed/datos_localidades.shp ../data/inegi/marco_geo/processed/datos_municipios.shp
- Habilitar plugin, de zonal-statistics y aplicar función en el menú Raster.
- capa ráster: tope_175
- capa de polígonos: datos_localidades, (después lo mismo con datos_municipios)
- columna de salida: 175_
- estadísticos: únicamente la función Suma
Los datos de referencia vienen incluidos en el repo de Git.