Advanced Plotting

# Add REI to stations2
rei <- REI(dets_with_stations, Rxdeploy)
stations2 <- left_join(stations2, rei) %>% select(-latitude, -longitude)
stations2 <- stations2 %>% rename(REI=rei)

head(stations2)

lon_range <- range(stationFishID$lon) + c(-10, 10)
lat_range <- range(stationFishID$lat) + c(-10, 10)

library(mapdata)
w <- map_data("lakes", ylim = lat_range, xlim = lon_range)

# new range to pass to plot for Nova Scotia region
lon_range <- range(stationFishID$lon) + c(-.5, .5)
lat_range <- range(stationFishID$lat) + c(-.5, .5)

library(plotly)

p <-  ggplot(stations2) +
  geom_polygon(data = w, aes(x = long, y = lat, group = group), fill = "darkblue") +
  theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
  #geom_sf(data = area, fill = 'white') +
  geom_point(data=stations2, aes(lon, lat, size=detections, col=REI))+
  scale_color_continuous(low="yellow", high="red")+
  coord_sf(xlim = lon_range, ylim = lat_range) +
  xlab("Longitude")+ ylab("Latitude")+
  theme(panel.background = element_rect(fill = 'lightgreen'),
        legend.position = 'none')

ggplotly(p)

library(sf)
library(gganimate)
library(rnaturalearth)
library(gifski)

dets_with_stations$year <- strftime(dets_with_stations$detection_timestamp_utc, format="%Y")
dets_with_stations$month <- strftime(dets_with_stations$detection_timestamp_utc, format="%Y-%m")

plot_data <- dets_with_stations %>% filter(year>"2017")

p <- ggplot(dets_with_stations) +
  geom_polygon(data = w, aes(x = long, y = lat, group = group), fill = "darkblue") +
  theme_bw() +
  theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())+
  #geom_sf(data = area, fill = 'white') +
  geom_point(data=dets_with_stations, aes(deploy_long, deploy_lat, col = animal_id), size = 3) +
  coord_sf(xlim = lon_range, ylim = lat_range) +
  labs(title = '',
       subtitle = 'Date: {format(frame_time, "%d %b %Y")}',
       x = "Longitude", y = "Latitude") +
  theme(panel.background = element_rect(fill = 'lightgreen'))+
  transition_time(as.POSIXct(detection_timestamp_utc, tz= "UTC", format = "%Y-%m-%d %H:%M"))+
  shadow_wake(wake_length = 0.5, alpha = FALSE)
  

pAnim <- animate(p, duration=24, nframes=96, height = 900, width = 900)
#watch it:
pAnim
#save it:
anim_save("pAnim.gif")