Error applying facet_wrap() in geom_sf() incomplete polygon in ggplot R-CodePudding

The

OBS EDIT

I tried @stefan's answer which partly worked, but generated a facet called "NA" bad.

new code:

p_total_especies_shp <- ggplot(total_especies_shp) geom_sf(data=regiao) geom_sf(aes(fill=factor(Total_esp))) scale_fill_brewer( palette = "Spectral", na.value = "grey", direction = -1, "Total de\nSepertens Regatadas") facet_wrap(~nome_popular) p_total_especies_shp

CodePudding user response：

The issue is that with faceting the data is splitted in groups and only the polygons contained in the splitted data will show up.

If you want all regions to be shown in each facet then one option would be to add a base map via second geom_sf layer. In your case geom_sf(regiao) geom_sf() should do the job.

As an example I make use of the default example from ?geom_sf:

library(ggplot2) set.seed(42) nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE) base <- nc nc$facet <- sample(c("a", "b", "c", "d"), size = nrow(nc), replace = TRUE) ggplot(nc) geom_sf(data = base) geom_sf(aes(fill = AREA)) facet_wrap(~facet)

CodePudding user response：

Working on @stefan answer, if you want to get rid of the NA panel, you need to provide the data in the second geom with the na.omit, leaving the ggplot call empty:

p_total_especies_shp <- ggplot() geom_sf(data = regiao) geom_sf(aes(fill = factor(Total_esp)), data = na.omit(total_especies_shp)) scale_fill_brewer( palette = "Spectral", na.value = "grey", direction = -1, "Total de\nSepertens Regatadas" ) facet_wrap(~nome_popular, drop = TRUE) p_total_especies_shp

Which gives the result you want: