Myujth

I am plotting (mapping) sf objects with ggplot2. My understanding is that since version 2.2.1 ggplot2 contains the geom geom_sf, for simple feature objects.

I can produce the exact map that I want by doing the following:

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

which produces

This is what I want except that it is missing a legend. For that, I am doing

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

which gives

That's not what I want, because in the legend:

1. 'points' should be points (not squares); and


2. 'circles' should be, well, circles (again, not squares)

Would be nice if the legend could respect the transparency of my colors (which it did in this example).

I tried to change the last couple of lines of the above to something like

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))))

but that didn't do anything to my plot.

Ideas?

Note: If it ends up being simpler for the plot, I could change the structure of the data, e.g., by combining objects a and b in the same simple feature object and add a column indicating the group (same for buf.a and buf.b).

Here's how far I managed to get to.

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

To get rid of the gray background in the legend symbols,

g + theme(legend.key = element_rect(fill = "white"))

The only issue here is that the circles do not have the transparency I wanted. This is odd.