rMaps使用自定义地图/ shp文件的ichoropleth

46
在R中,我想尝试复制这里的教程,使我的自定义SHP(Shapefile)文件或地图成为交互式区域填充地图......
该地图是北爱尔兰小地区的地图。可以在此处找到。
以下是我迄今为止所采取的步骤...
我认为问题在于geographyConfig数据的设置...
非常感谢任何帮助...
# 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)

加载rMaps.html文件不能产生相关的地图,因为它只显示底部的切割而不是地图。


加载 rMaps.html 并不会生成相关地图,而只会生成底部部分,因为它会触发 JavaScript 错误,这可以在 Chrome 的 JavaScript 控制台或 Safari 的 Web Inspector 中看到。我不知道是什么导致了这些 JavaScript 错误。 - WhiteViking
2个回答

1

请在你的回答中添加更多细节。 - NSNoob
Leaflet很棒。这是一个开始研究leaflet的好地方。https://dev59.com/FYnca4cB1Zd3GeqP5xNz - GlennFriesen

0

你看过http://www.r-bloggers.com/rmaps-mexico-map/吗? 我之前也遇到了类似的问题,后来发现对于小范围区域,需要在下面的代码中调整比例参数。这样可以解决问题。

    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)

网页内容由stack overflow 提供, 点击上面的
可以查看英文原文,
原文链接