ben/pretty_rivers
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

edited figure script

Browse Source
This commit is contained in:
Ben Varick 2023-10-16 16:53:19 -05:00
parent 9cf982c483
commit 529d3ae923
1 changed files with 22 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

44
03_make_figures.R
View File

@ -7,7 +7,7 @@ ifelse(!dir.exists(file.path(getwd(), figure_dir)), dir.create(file.path(getwd()
# set colors # set colors
colors <- list(darkblue = "#062e57", colors <- list(darkblue = "#062e57",
lightblue = "#b1dcf3") lightblue = "#b1dcf3")
# set dimensions of figures # set dimensions of figures
width <- 8.5 width <- 8.5
@ -20,40 +20,40 @@ margin <- 0
# plot map # plot map
plot_map <- function () { plot_map <- function () {
plot(sf::st_geometry(extent_poly), plot(sf::st_geometry(extent_poly),
col = colors$darkblue, col = colors$darkblue,
border = "black") border = "black")
plot(sf::st_geometry(data$NHDWaterbody), plot(sf::st_geometry(data$NHDWaterbody),
col = colors$lightblue, col = colors$lightblue,
border = NA, border = NA,
add = TRUE) add = TRUE)
plot(sf::st_geometry(data$NHDArea), plot(sf::st_geometry(data$NHDArea),
col = colors$lightblue, col = colors$lightblue,
border = NA, border = NA,
add = TRUE) add = TRUE)
plot(sf::st_geometry(data$NHDFlowline), plot(sf::st_geometry(data$NHDFlowline),
col = colors$lightblue, col = colors$lightblue,
lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446, lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446,
border = NA, border = NA,
add = TRUE) add = TRUE)
} }
# save figure # save figure
tiff(filename = paste0(figure_dir,"/map.tiff"), tiff(filename = paste0(figure_dir,"/map.tiff"),
width = width, width = width,
height = height, height = height,
units = units, units = units,
res = res, res = res,
compression = "lzw") compression = "lzw")
par(mai=c(margin, margin, margin, margin)) par(mai=c(margin, margin, margin, margin))
par(mar=c(1,1,1,1)) par(mar=c(1,1,1,1))
plot_map() plot_map()
dev.off() dev.off()
png(filename = paste0(figure_dir, "/map.png"), png(filename = paste0(figure_dir, "/map.png"),
width = width, width = width,
height = height, height = height,
units = units, units = units,
res = res) res = res)
par(mai=c(margin, margin, margin, margin)) par(mai=c(margin, margin, margin, margin))
par(mar=c(1,1,1,1)) par(mar=c(1,1,1,1))
plot_map() plot_map()