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