1

I use the tmap package. I plot a map with the qtm() function but the legend is not displayed. Indeed, I understand that the legend is automatically display, unless you specify to hide it.

I also try to use a work-around with the tm_add_legend() but it doesn't work.

MyPlot

Here is my code :

library(sp) ; library(rgdal) ; library(rgeos) ;
library(tmap) ; library(maptools)

Directory <- ""

MapFrance <- readOGR(dsn = paste0(Directory, "Departements"),  layer = "DEPARTEMENT")

MapFrance@data$Score <- round(runif(96, min = 1, max = 10), 0)

brk <- c(1, 2.5, 5, 7.5, 10)
ClassScore <- classIntervals(MapFrance@data$Score, style = "fixed", fixedBreaks = brk)
Palette <- brewer.pal(length(brk) - 1, "Greens")

MapFrance@data$Categorie <- as.character(cut(MapFrance@data$Score,
                                             breaks = ClassScore$brks,
                                             labels = Palette,
                                             include.lowest = TRUE))

ExampleData <- MapFrance[MapFrance$CODE_REG == "11", ]

qtm(ExampleData, fill = "Categorie", text = "CODE_DEPT", text.size = "AREA", style = "gray",
    text.root = 5, fill.title = "Opinion scale")

And a link to the geographic data : https://drive.google.com/file/d/1wksMnkufQPdOfSbhWoZbPKMqlcLcll_j

4

1 回答 1

1

具有三个优点的解决方案:

  • 解决图例问题
  • 解决标题问题
  • 避免classIntervals在绘图之前调用离散化连续变量

编码 :

library(sp) ; library(rgdal) ; library(rgeos) ;
library(tmap) ; library(maptools)

Directory <- ""

MapFrance <- readOGR(dsn = paste0(Directory, "Departements"),  layer = "DEPARTEMENT")

MapFrance@data$Score <- round(runif(96, min = 1, max = 10), 0)

ExampleData <- MapFrance[MapFrance$CODE_REG == "11", ]

brk <- c(1, 2.5, 5, 7.5, 10)

tm_shape(ExampleData) +
  tm_fill("Score",
          title = "Satisfaction score",
          style = "fixed",
          breaks = brk,
          palette = brewer.pal(length(brk) - 1, "Greens"),
          auto.palette.mapping = FALSE) +
  tm_borders() +
  tm_text("CODE_DEPT",
          size = "AREA",
          root = 10) +
tm_layout(main.title = "Satisfaction score in France",
          legend.position = c("left", "bottom"))

剧情 :

剧情

于 2017-11-17T10:29:58.297 回答