changed figure margins

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)

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()