datamarindo

(click para agrandar)    En el post anterior hicimos un scrapping de las concesiones mineras de la página del SIAM; sin  embargo, después de revisar, el código sólo rascaba las geometrías, por lo cual tuvimos que minar en otra parte del REST query que nos daba un json por cada concesión minera de las 32,000 cambiando la URL por esta: foreach(i=1:32000) [...] sprintf("http://187.191.71.167/ArcGIS/rest/services/DGM/ConcesionesMineras/MapServer/0/%s?f=json&pretty=true", i)    Esta URL no deja mandar consultas en bloques de cientos/miles y además no permite más de 1 consulta por segundo (otras dejan más de 10 consultas por segundo); así que con paciencia (relativa) en unas cuantas horas podemos obtener un bello geojson con atributos de las concesiones mineras, que trae información del nombre del titular (compañía minera), superficie y sustancias que van a extraer. Descargar y descomprimir     Me preguntaran: ¿Por qué geojson? Porque estamos en 2018; además, c

El sistema SIAM (http://www.cartografia.economia.gob.mx/cartografia/#) contiene información de las concesiones mineras; sin embargo, el acceso a la información geográfica que tiene la página hace difícil la descarga y obliga a hacerlo concesión por concesión; por esto minamos la página con R y sacamos este archivo de 32000 concesiones mineras (el que contiene atributos necesitó un poco más de código): Concesiones_con_atributos Concesiones mineras (solo_geometrias) El código necesario en R para obtener todos los polígonos: library (foreach) library (doMC) library (sf) registerDoMC(12)                      # candela baba <- seq(1, 32000, by=100) # para ir sacando de cien en cien foreach (i=baba, .combine=rbind, .inorder=F) %dopar% {   tet <- paste(as.character(i:(i+99)), collapse = "%2C")  #para crear la serie de cien en la URL   url <- sprintf ("http://187.191.71.167/ArcGIS/rest/services/DGM/ConcesionesMineras/MapServer/0/query?text=&geom