scrapping de un rest query de concesiones Mineras
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=&geometry=&geometryType=esriGeometryPolygon&inSR=&spatialRel=esriSpatialRelIntersects&relationParam=&objectIds=%s&where=&time=&returnIdsOnly=false&returnGeometry=true&maxAllowableOffset=&outSR=&outFields=&f=pjson", tet)
read_sf(url) # leemos el archivo geoJson
} -> concesiones # nuestro archivo
plot(concesiones) # crea la imagen de arriba
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=&geometry=&geometryType=esriGeometryPolygon&inSR=&spatialRel=esriSpatialRelIntersects&relationParam=&objectIds=%s&where=&time=&returnIdsOnly=false&returnGeometry=true&maxAllowableOffset=&outSR=&outFields=&f=pjson", tet)
read_sf(url) # leemos el archivo geoJson
} -> concesiones # nuestro archivo
plot(concesiones) # crea la imagen de arriba
hola tengo problemas al descargar el paquete doMC ya que al instalarlo me sale que mi version de R (3.4.2 bajo windows 10) no es compatible con dicho paquete, me podrias recomendar otro paquete o alguna funcion para poder suplir esto o algun otro ejemplo? te lo agradeceria
ResponderBorrarMira, existe el DoSNOW, pero de verdad si quieres andar haciendo este tipo de cosas lo mejor es que dejes windows
Borrar