nrg_te_oil_1_Data <- read_excel("europa_exp_imp_petr/nrg_te_oil_1_Data.xlsx", 
                                col_types = c("text", "text", "text", 
                                              "text", "text", "numeric", "text"))
ue_exp <- nrg_te_oil_1_Data %>% select("TIME","GEO", "SIEC", "UNIT", "Value" ) %>%
  group_by(TIME) %>%
  arrange(TIME, desc(Value)) %>%
  filter(GEO != "European Union - 27 countries (from 2020)") %>%
  filter(GEO != "European Union - 28 countries (2013-2020)") %>%
  filter(GEO != "Euro area - 19 countries  (from 2015)") %>%
  mutate(ranking = row_number())  %>%
  filter(ranking <=15)

# Animacion UE Exp
animacion_1 <- ue_exp %>%
  ggplot() +
  geom_col(aes(ranking, Value, fill = GEO)) +
  geom_text(aes(ranking, Value, label = Value), hjust=-0.1) +
  geom_text(aes(ranking, y=0 , label = GEO), hjust=1.1) + 
  geom_text(aes(x=15, y=max(Value) , label = as.factor(TIME)), vjust = 0.2, alpha = 0.5,  col = "gray", size = 20) +
  coord_flip(clip = "off", expand = FALSE) +  scale_x_reverse() +
  theme_igray() + theme(
    panel.grid = element_blank(), 
    legend.position = "none",
    axis.ticks.y = element_blank(),
    axis.title.y = element_blank(),
    axis.text.y = element_blank(),
    plot.margin = margin(1, 4, 1, 4, "cm")
  ) +
  transition_states(TIME, state_length = 0, transition_length = 2) +
  enter_fade() +
  exit_fade() + 
  ease_aes('quadratic-in-out') + 
  labs(title = "Gráfico: Top 10 Países EU por exportaciones",
       subtitle = "(oil and petroleum products)",
       caption = "Datos provenientes de Eurostat",
       y = "Valor en miles de toneladas") 

animacion_1
animate(animacion_1, width = 700, height = 432, fps = 25, duration = 20, rewind = FALSE)

# Datos Europa Imp
nrg_ti_oil_1_Data <- read_excel("europa_exp_imp_petr/nrg_ti_oil_1_Data.xlsx", 
  col_types = c("text", "text", "text", 
  "text", "text", "numeric", "text"))

ue_imp <- nrg_ti_oil_1_Data %>% select("TIME","GEO", "SIEC", "UNIT", "Value" ) %>%
  group_by(TIME) %>%
  arrange(TIME, desc(Value)) %>%
  filter(GEO != "European Union - 27 countries (from 2020)") %>%
  filter(GEO != "European Union - 28 countries (2013-2020)") %>%
  filter(GEO != "Euro area - 19 countries  (from 2015)") %>%
  mutate(ranking = row_number())  %>%
  filter(ranking <=15)
# Animación UE Imp
animacion_2 <- ue_imp %>%
  ggplot() +
  geom_col(aes(ranking, Value, fill = GEO)) +
  geom_text(aes(ranking, Value, label = Value), hjust=-0.1) +
  geom_text(aes(ranking, y=0 , label = GEO), hjust=1.1) + 
  geom_text(aes(x=15, y=max(Value) , label = as.factor(TIME)), vjust = 0.2, alpha = 0.5,  col = "gray", size = 20) +
  coord_flip(clip = "off", expand = FALSE) + scale_x_reverse() +
  theme_igray() + theme(
    panel.grid = element_blank(), 
    legend.position = "none",
    axis.ticks.y = element_blank(),
    axis.title.y = element_blank(),
    axis.text.y = element_blank(),
    plot.margin = margin(1, 4, 1, 4, "cm")
  ) +
  transition_states(TIME, state_length = 0, transition_length = 2) +
  enter_fade() +
  exit_fade() + 
  ease_aes('quadratic-in-out') + 
  labs(title = "Gráfico: Top 10 Países EU por importaciones",
       subtitle = "(oil and petroleum products)",
       caption = "Datos provenientes de Eurostat",
       y = "Valor en miles de toneladas")

animacion_2
animate(animacion_2, width = 700, height = 432, fps = 25, duration = 20, rewind = FALSE)

# Datos Exp e Imp UE 

nrg_ti_te_oil_1_Data <- read_excel("europa_exp_imp_petr/nrg_ti_te_oil_1_Data.xlsx", 
  col_types = c("text", "numeric", "text", 
  "text", "text", "text", "numeric", "numeric"))

ue_exp_imp <- nrg_ti_te_oil_1_Data %>% select("TIME","GEO", "SIEC", "UNIT", "Valueexp", "Valueimp", "Code") %>%
  filter(GEO != "European Union - 27 countries (from 2020)") %>%
  filter(GEO != "European Union - 28 countries (2013-2020)") %>%
  filter(GEO != "Euro area - 19 countries  (from 2015)")
# Animacion Exp e Imp UE 


animacion_3 <- ggplot(
  ue_exp_imp, 
  aes(x = Valueexp , y= Valueimp, size = 0,7, colour = GEO, country = Code)
) +
  geom_point(show.legend = FALSE, alpha = 0.7) +
  scale_color_viridis_d() +
  scale_size(range = c(2, 12)) +
  labs(x ="Exportaciones", y = "Importaciones") + 
  theme_igray() + 
  transition_time(TIME) +
  labs(title ="Gráfico: Países EU por importaciones y exportaciones", 
       subtitle = "Year: {frame_time}") +
  geom_flag(size = 10)
  
animacion_3

nrg_ti_te_oil_1_Data <- read_excel("europa_exp_imp_petr/nrg_ti_te_oil_1_Data.xlsx", 
  col_types = c("text", "numeric", "text", 
  "text", "text", "text", "numeric", "numeric"))
# Datos Exp e Imp España 
spain_exp_imp <- nrg_ti_te_oil_1_Data %>% select("TIME","GEO", "SIEC", "UNIT", "Valueexp", "Valueimp" ) %>%
  filter(GEO == "Spain")
# Animacion Exp e Imp España 
spain_imp <- ggplot(spain_exp_imp, aes(TIME, Valueimp, fill = Valueimp)) +
  geom_col() +
  scale_fill_distiller(palette = "Reds", direction = 1) +
  theme_minimal() +
  theme(
    panel.grid = element_blank(),
    panel.grid.major.y = element_line(color = "white"),
    panel.ontop = TRUE
  ) + transition_states(TIME, wrap = FALSE) +
  shadow_mark()
spain_imp 

spain_exp <- ggplot(spain_exp_imp, aes(TIME, Valueexp, fill = Valueexp)) +
  geom_col() +
  scale_fill_distiller(palette = "Blues", direction = 1) +
  theme_minimal() +
  theme(
    panel.grid = element_blank(),
    panel.grid.major.y = element_line(color = "white"),
    panel.ontop = TRUE
  ) + transition_states(TIME, wrap = FALSE) +
  shadow_mark()
spain_exp 

# Estados Unidos
us_hex <- geojson_read("us_states_hexgrid.geojson",  what = "sp")
us_hex@data = us_hex@data %>%
  mutate(google_name = gsub(" \\(United States\\)", "", google_name))

us_hex@data = us_hex@data %>% mutate(google_name = gsub(" \\(United States\\)", "", google_name))
us_hex_fortified <- tidy(us_hex, region = "google_name")


centers <- cbind.data.frame(data.frame(gCentroid(us_hex, byid=TRUE), id=us_hex@data$iso3166_2))


us_oil <- read_excel("PET_CRD_CRPDN_ADC_MBBL_A.xls", 
  col_types = c("text", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric", "numeric", "numeric", 
                "numeric"))
us_hex_oil <- us_hex_fortified %>%
  left_join(. , us_oil, by=c("id"="state"))


# Preparación 1 
us_hex_oil$bin <- cut( us_hex_oil$y_2020 , breaks=c(1,1000, 5000, 10000, 50000, 100000, Inf), labels=c("1-1000", "1001-5000", "5001-10000", "10001-50000", "50001-100000", "100000+" ), include.lowest = TRUE )

# Preparación 2 
my_palette <- rev(magma(8))[c(-1,-8)]

# plot
ggplot() +
  geom_polygon(data = us_hex_oil, aes(fill = bin, x = long, y = lat, group = group) , size=0, alpha=0.9) +
  geom_text(data=centers, aes(x=x, y=y, label=id), color="white", size=3, alpha=0.6) +
  theme_void() +
  scale_fill_manual( 
    values=my_palette, 
    name="Unidad: miles de barriles     Nota: 1 barril = 159L aprox", 
    guide = guide_legend( keyheight = unit(3, units = "mm"), keywidth=unit(12, units = "mm"), label.position = "bottom", title.position = 'top', nrow=1) 
  ) +
  ggtitle( "Crude Oil production USA 2020" ) +
  theme(
    legend.position = c(0.5, 0.9),
    text = element_text(color = "#22211d"),
    plot.background = element_rect(fill = "#f5f5f2", color = NA), 
    panel.background = element_rect(fill = "#f5f5f2", color = NA), 
    legend.background = element_rect(fill = "#f5f5f2", color = NA),
    plot.title = element_text(size= 22, hjust=0.5, color = "#4e4d47", margin = margin(b = -0.1, t = 0.4, l = 2, unit = "cm")),
  )

us_oil <- us_oil %>% filter(state != "US") %>% filter(state != "Texas") %>% 
  mutate(expresion_log = log(y_2020))


us <- ggplot(us_oil, aes(state, expresion_log, fill = expresion_log)) +
  geom_col() + coord_polar() +
  scale_y_continuous(
    limits = c(0, 14),
    expand = c(0, 0),
    breaks = c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14)
  ) + 
  scale_fill_gradientn(
    "Escala logarítmica",
    colours = c( "#6C5B7B","#C06C84","#F67280","#F8B195")
  ) +
  guides(
    fill = guide_colorsteps(
      barwidth = 15, barheight = .5, title.position = "top", title.hjust = .5
    )
  ) +
  theme(
    axis.title = element_blank(),
    axis.ticks = element_blank(),
    axis.text.y = element_blank(),
    axis.text.x = element_text(color = "gray12", size = 6),
    legend.position = "bottom",
  )+labs(title = "Producción por Estados",
         subtitle = "Unidades originales: miles de barriles ")
us

knitr::include_graphics(here::here("imagenes", "aguacate-o-petroleoo.gif")  )