Example showing Blue/Orange Diverging Color Palette

blue_orange_diverging_palette.R

library(sf)
library(ggplot2)

nc <- st_read(system.file("shapes/sids.shp", package="spData")[1], quiet=TRUE)

nc %>% 
  ggplot(aes(fill = AREA)) + 
  geom_sf(size = 0.05, color = 'white') + 
  scale_fill_gradientn(
    colors = c("#DA8907", "#ECCD00", "white", "#62AEDC", "#203856"))

nc %>% 
  ggplot(aes(fill = AREA)) + 
  geom_sf(size = 0.05, color = 'white') + 
  scale_fill_stepsn(
    colors = c("#DA8907", "#ECCD00", "white", "#62AEDC", "#203856"))

# note that the breaks and range of values can be specified manually by using the values argument