rutils/hex.R

In ibarraespinosa/eixport: Export Emissions to Atmospheric Models

library(vein)
library(sf)
library(hexSticker)
library(ggplot2)
data(net)

cellsize_backup <- 0.01
color_backup <- 3681

cs <- seq(0.000, 0.01 * 3, 0.002)
cs <- cs[2:length(cs)]

i <- 8

pngi <- ifelse(nchar(seq_along(cs)) < 2,
      paste0(0, seq_along(cs)),
      seq_along(cs)
)


g2 <- st_make_grid(st_as_sfc(st_bbox(st_as_sf(net))),
      cellsize = cs[i],
      square = FALSE
)
ge <- st_sf(id = 1:length(g2), geometry = st_sfc(g2))
hdv <- sf::st_as_sf(net)[, "hdv"]
netg <- emis_grid(spobj = hdv[as.numeric(hdv$hdv) > 600, ], g = ge)
hdv2 <- netg[as.numeric(netg$hdv) > 760, ]
hdv2$hdv[25] <- NA
# plot(hdv2["hdv"],
#      axes = F,
#      main = "",
#      key.pos=NULL,
#      pal = cptcity::cpt(3681, colorRampPalette = T, rev = T),
#      lty= 0)



sticker(~ plot(hdv2["hdv"],
      axes = F,
      main = "",
      key.pos = NULL,
      pal = cptcity::lucky(colorRampPalette = T, rev = T),
      lty = 0
),
package = "eixport",
s_x = 0.8,
s_y = 0.7,
s_width = 2,
s_height = 2,
p_x = 1,
p_y = 1.6,
p_color = "darkblue",
p_family = "sans",
p_size = 18,
h_size = 1.2,
h_fill = "white",
h_color = "black",
spotlight = FALSE,
l_x = 1,
l_y = 0.5,
l_width = 3,
l_height = 3,
l_alpha = 0.4,
url = "",
u_x = 1,
u_y = 0.08,
u_color = "black",
u_family = "Aller_Rg",
u_size = 1.5,
u_angle = 30,
white_around_sticker = FALSE,
filename = paste0("rutils/logo", pngi[i], ".png"),
asp = 1, dpi = 220
)
rm(list = ls())
gc()