mapeo - mapas de densidad en r

rMaps icoropleta con mapa personalizado/shp

Question

En R, me gustaría ayudar a tratar de replicar el tutorial here para que mi archivo o mapa SHP (Shapefile) personalizado sea un mapa interactivo de coropletas ...

El mapa es de las Pequeñas Áreas de Irlanda del Norte. Que se puede encontrar here .

Debajo están los pasos que he tomado hasta ahora ...

Creo que el problema es la configuración de los datos de geographyConfig ...

Cualquier ayuda sería muy apreciada ...

# Download and unzip the data
system(''wget http://www.nisra.gov.uk/archive/geography/digital_products/SA2011_Esri_Shapefile.zip'')
system(''unzip SA2011_Esri_Shapefile.zip'')
# Load libraries
library(rgdal)
library(rgeos)
library(rMaps)
shp.file <- ''SA2011.shp''
# Convert projection
system(paste0(''ogr2ogr tmp.shp '',
shp.file,
'' -t_srs "+proj=longlat +ellps=WGS84 +no_defs +towgs84=0,0,0"''))
# Read in the data
xx <- readOGR(dsn=getwd(),layer=''tmp'')
mm <- xx@data
head(mm)
n <- nrow(mm)
dat.val <- mm$Hectares
# Add extra year data
mm <- mm[rep(seq(n),3),]
mm$Hectares <- c(dat.val,rev(dat.val),dat.val/2)
mm$year <- rep(c(2000:2002),each=n)
colnames(mm)[1] <- ''ID''
id.var <- ''SA2011''
# Convert to json
system(paste0(''topojson -o tmp.json -s 1e-7 -q 1e5 tmp.shp -p ID='',
id.var,
'' --id-property '',
id.var))
d1 <- ichoropleth(Hectares ~ ID, data = mm, ncuts = 9, pal = ''YlOrRd'',
animate = ''year'', map = ''states''
)
d1$set(
geographyConfig = list(
dataUrl = "tmp.json"
),
scope = ''states'',
setProjection = ''#! function( element, options ) {
var projection, path;
projection = d3.geo.mercator()
.center([-7, 55]).scale(element.offsetWidth)
.translate([element.offsetWidth / 2, element.offsetHeight / 2]);
path = d3.geo.path().projection( projection );
return {path: path, projection: projection};
} !#''
)
d1$save(''rMaps.html'', cdn = TRUE)

La carga de rMaps.html no produce el mapa relevante, ya que solo muestra los cortes en la parte inferior pero no en el mapa.

Answer 1

¿Has visto http://www.r-bloggers.com/rmaps-mexico-map/ ? Estaba teniendo un problema similar y descubrí que para las áreas de pequeña escala, necesita manipular el argumento de escala en el fragmento de código siguiente. Funciona.

 d1$set(
geographyConfig = list(
dataUrl = "tmp.json"
),
scope = ''tmp'',
setProjection = ''#! function( element, options ) {
var projection, path;
projection = d3.geo.mercator()
.center([-5.832087, 54.605035]).scale(element.offsetWidth)
.translate([element.offsetWidth / 2, element.offsetHeight / 2]);
path = d3.geo.path().projection( projection );
return {path: path, projection: projection};
} !#''
)
d1$save(''rMaps.html'', cdn = TRUE)

Answer 2

Varias personas en recordnotfound.com le pidieron al autor que respondiera a este hilo (en recordnotfound.com). Me puse en contacto con el propietario del proyecto rMaps, Ramnath. Sugirió que hay un nuevo proyecto que proporciona un conjunto mejorado de características: http://github.com/rstudio/leaflet

Aquí hay más detalles sobre el proyecto del folleto: https://recordnotfound.com/leaflet-rstudio-35205