diff --git a/R/plot_scripts/Fig_autotrophy_map_only.R b/R/plot_scripts/Fig_autotrophy_map_only.R index f5fd611..9699977 100644 --- a/R/plot_scripts/Fig_autotrophy_map_only.R +++ b/R/plot_scripts/Fig_autotrophy_map_only.R @@ -6,7 +6,7 @@ lapply(c("lubridate","cowplot", "tidyverse", "reshape2", "ggmap","maps","mapdata", "mapproj", - "ggsn", "sf", "ggspatial"), require, character.only=T) + "viridis", "sf", "ggspatial"), require, character.only=T) ## Import Data GPP_info <- read_csv("data_working/across_sites_model_data.csv") %>% @@ -98,12 +98,12 @@ my_breaks <- c(0.05, 0.15, 0.4, 1) fill = "white", color = "black") + # map of states geom_point(data = AI_sf, aes(x = lon, y = lat, color = mean_PR_annual, - size = mean_GPP_annual), - alpha = 0.95) + # map of sites + size = mean_GPP_annual)) + # map of sites labs(x = "Longitude", y = "Latitude") + scale_color_gradientn("Mean Annual P:R", - colors = c("#F0C9C0", "#EDB48E", - "#E8C32E", "#E6E600", "#00A600"), + colors = c("#440154FF", "#482677FF", + "#404788FF", "#2D708EFF", + "#73D055FF", "#95D840FF"), trans = "log", breaks = my_breaks, labels = my_breaks) + # based on terrain.colors(n = 10) # custom colors @@ -154,7 +154,7 @@ my_breaks <- c(0.05, 0.15, 0.4, 1) # export exploratory figures # ggsave(sitemap, -# filename = "figures/Annual_GPP_PR_USmap_043024.png", +# filename = "figures/Annual_GPP_PR_USmap_071924.png", # width = 22, # height = 11, # units = "cm" diff --git a/figures/Annual_GPP_PR_USmap_071924.png b/figures/Annual_GPP_PR_USmap_071924.png new file mode 100644 index 0000000..b714550 Binary files /dev/null and b/figures/Annual_GPP_PR_USmap_071924.png differ