descarga_procesa_NOAA

## NOAA
Visita el sitio de la NOAA, https://ngdc.noaa.gov/eog/viirs/
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.  
$ gdalwarp -q -cutline ../marco_geo/raw/areas_geoestadisticas_estatales.shp \
  -crop_to_cutline --config GDALWARP_IGNORE_BAD_CUTLINE YES \
  "../luminosidad_tile1_2015.tif" \
  "../luminosidad_mexico_2015.tif"

  
# Procesar datos de NOAA.

Los datos de la NOAA se procesan aparte en QGIS. 
Topamos a 175 debido a zonas con luminosidad extrema.  (Ver ficha técnica)

0. Empezar nuevo proyecto en QGIS como: ../data/qgis/imco_magda.qgs
1. Abrir capa de ráster del .tif generado en: ../data/viirs/luminosidad_mexico.tif
    - Renombrar como viirs
2. Aplicar raster-calculator:  viirs * ( viirs <= 175), 
    - guardar en ../data/viirs/processed/tope_175.tif
3. Abrir capas vectoriales de .shp en ../data/inegi/marco_geo/raw/
    - poligonos_localidades_urbanas_y_rurales.shp
    - areas_geoestadisticas_municipales.shp
4. 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
5. 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