使用 sf 对象在 ggplot2 中自定义图例

dja*_*jas 7 r ggplot2 r-sf

我正在sf使用 绘制(映射)对象ggplot2。我的理解是,由于版本 2.2.1ggplot2包含 geom geom_sf,用于简单的特征对象。

我可以通过执行以下操作来生成我想要的精确地图:

library(sf)
library(ggplot2)

# some points to start with
a <- st_as_sf(data.frame(lon = c(1,4,6), lat = c(0,0,-3)), coords = c('lon', 'lat'))
b <- st_as_sf(data.frame(lon = c(2.5,4), lat = c(-4.5,-5)), coords = c('lon', 'lat'))

# circles around those points
buf.a <- st_buffer(a, 1)
buf.b <- st_buffer(b, 1)

# colors to mark the points
sol.a = rgb(1,0,0) 
sol.b = rgb(0,0,1) 

# colors to fill the circles
fil.a = adjustcolor(sol.a, alpha.f = .25)
fil.b = adjustcolor(sol.b, alpha.f = .25)

# the plot I want
g = ggplot() +
    geom_sf(data = buf.a, fill =  fil.a, color = NA) +
    geom_sf(data = buf.b, fill =  fil.b, color = NA) +
    geom_sf(data = a,     color = sol.a, shape = 20, size = 3) +
    geom_sf(data = b,     color = sol.b, shape = 20, size = 3)
g
Run Code Online (Sandbox Code Playgroud)

产生

在此输入图像描述

这就是我想要的,只是它缺少一个图例。为此,我正在做

cols.fill = c("GROUP A" = fil.a, "GROUP B" = fil.b)
cols.sol = c("GROUP A" = sol.a, "GROUP B" = sol.b)

g = ggplot() +
    geom_sf(data = buf.a, color = NA, aes(fill = 'GROUP A')) +
    geom_sf(data = buf.b, color = NA, aes(fill = 'GROUP B')) +
    geom_sf(data = a,     shape = 20, size = 3, aes(color = 'GROUP A')) + 
    geom_sf(data = b,     shape = 20, size = 3, aes(color = 'GROUP B')) +
    scale_fill_manual(name = "circles", values = cols.fill) +
    scale_color_manual(name = "points", values = cols.sol)
g
Run Code Online (Sandbox Code Playgroud)

这使

在此输入图像描述

这不是我想要的,因为传说中:

  1. “点”应该是点(而不是正方形);和
  2. “圆圈”应该是圆圈(同样,不是正方形)

如果图例能够尊重我的颜色的透明度(在本例中就是这样做的),那就太好了。

我尝试将上面的最后几行更改为类似的内容

scale_fill_manual(name = "circles", values = cols.fill,
                  guide = guide_legend(override.aes = list(shape = c(19, 19)))) +
scale_color_manual(name = "points", values = cols.sol,
                   guide = guide_legend(override.aes = list(shape = c(20, 20))))
Run Code Online (Sandbox Code Playgroud)

但这对我的情节没有任何影响。

有想法吗?

注意:如果绘图最终变得更简单,我可以更改数据的结构,例如,通过将对象a和组合b在同一简单特征对象中,并添加指示组的列(与buf.a和相同buf.b)。

dja*_*jas 5

这是我设法达到的程度。

g = ggplot() +
    geom_sf(data = buf.a, color = NA, aes(fill = 'GROUP A'), show.legend = "point") +
    geom_sf(data = buf.b, color = NA, aes(fill = 'GROUP B'), show.legend = "point") +
    geom_sf(data = a,     shape = 20, size = 3, aes(color = 'GROUP A'), show.legend = "point") + 
    geom_sf(data = b,     shape = 20, size = 3, aes(color = 'GROUP B'), show.legend = "point") +
    scale_color_manual(name = "points", values = cols.sol,
                       guide = guide_legend(override.aes = list(shape = c(20, 20)))) +
    scale_fill_manual(name = "circles", values = cols.fill,
                      guide = guide_legend(override.aes = list(shape = c(20, 20), color = cols.fill, size = 8)))
g
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

要消除图例符号中的灰色背景,

g + theme(legend.key = element_rect(fill = "white"))
Run Code Online (Sandbox Code Playgroud)

在此输入图像描述

这里唯一的问题是圆圈没有我想要的透明度。这很奇怪。