From 66392afe4556fb6beb757351e527ebea194ab236 Mon Sep 17 00:00:00 2001 From: Ben Varick Date: Mon, 16 Oct 2023 15:24:53 -0500 Subject: [PATCH] changed figure margins --- 02_process_data.R | 22 ++++++++++++++ 03_make_figures.R | 74 +++++++++++++++++++++++------------------------ 2 files changed, 59 insertions(+), 37 deletions(-) diff --git a/02_process_data.R b/02_process_data.R index f0bf72a..672ac62 100644 --- a/02_process_data.R +++ b/02_process_data.R @@ -1,6 +1,28 @@ # load libraries library(sf) +# reload extent data so you don't need to run the download data script again +data_dir <- "data" +nhdplusTools_data_dir(dir = paste0(data_dir, "/hydrologic")) + +extent <- read.csv(file = "extent.csv") +crs <- 4269 + +extent_poly <- st_polygon( + x = list( + cbind( + extent$longitude[c(1,2,2,1,1)], + extent$latitude[c(1,1,2,2,1)]) + ) +) + +extent_poly <- st_sfc(extent_poly, crs=4326) +extent_bbox <- st_bbox(st_transform(x = extent_poly, crs)) +extent <- list(longitude_max = max(extent$longitude), + longitude_min = min(extent$longitude), + latitude_max = max(extent$latitude), + latitude_min = min(extent$latitude)) + # load and crop hydrologic data data <- list(NULL) sf_use_s2(FALSE) diff --git a/03_make_figures.R b/03_make_figures.R index e9f5ddc..9ab5308 100644 --- a/03_make_figures.R +++ b/03_make_figures.R @@ -7,54 +7,54 @@ ifelse(!dir.exists(file.path(getwd(), figure_dir)), dir.create(file.path(getwd() # set colors colors <- list(darkblue = "#062e57", - lightblue = "#b1dcf3") + lightblue = "#b1dcf3") + +# set dimensions of figures +width <- 8.5 +height <- 11 +units <- "in" +res <- 600 +margin <- 0 + # plot map plot_map <- function () { - plot(sf::st_geometry(extent_poly), - col = colors$darkblue, - extent = extent_bbox, - xlim = c(extent$longitude_min, extent$longitude_max), - ylim = c(extent$latitude_min, extent$latitude_max), - border = "black") - plot(sf::st_geometry(data$NHDWaterbody), - col = colors$lightblue, - border = NA, - extent = extent_bbox, - xlim = c(extent$longitude_min, extent$longitude_max), - ylim = c(extent$latitude_min, extent$latitude_max), - add = TRUE) - plot(sf::st_geometry(data$NHDArea), - col = colors$lightblue, - border = NA, - extent = extent_bbox, - xlim = c(extent$longitude_min, extent$longitude_max), - ylim = c(extent$latitude_min, extent$latitude_max), - add = TRUE) - plot(sf::st_geometry(data$NHDFlowline), - col = colors$lightblue, - lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446, - border = NA, - extent = extent_bbox, - xlim = c(extent$longitude_min, extent$longitude_max), - ylim = c(extent$latitude_min, extent$latitude_max), - add = TRUE) +plot(sf::st_geometry(extent_poly), + col = colors$darkblue, + border = "black") +plot(sf::st_geometry(data$NHDWaterbody), + col = colors$lightblue, + border = NA, + add = TRUE) +plot(sf::st_geometry(data$NHDArea), + col = colors$lightblue, + border = NA, + add = TRUE) +plot(sf::st_geometry(data$NHDFlowline), + col = colors$lightblue, + lwd = data$NHDFlowline$TotDASqKM^0.3204*0.0446, + border = NA, + add = TRUE) } # save figure tiff(filename = paste0(figure_dir,"/map.tiff"), - width = 11, - height = 17, - units = "in", - res = 600, + width = width, + height = height, + units = units, + res = res, compression = "lzw") +par(mai=c(margin, margin, margin, margin)) +par(mar=c(1,1,1,1)) plot_map() dev.off() png(filename = paste0(figure_dir, "/map.png"), - width = 11, - height = 17, - units = "in", - res = 600) + width = width, + height = height, + units = units, + res = res) +par(mai=c(margin, margin, margin, margin)) +par(mar=c(1,1,1,1)) plot_map() dev.off()