Interactive visualization of VMS tracking data of fishing vessels

p <- mpa_tracks %>% 
        ggplot()+
        geom_polygon(data = mex, aes(x = long, y = lat, group=group), colour = "grey", fill = "grey", alpha = .3)+
        geom_polygon(data = eez, aes(x = long, y = lat, group=group), fill = "#68829E", alpha=.3)+
        geom_polygon(data = mpa, aes(x = long, y = lat, group=group), colour = "red", fill = NA, alpha = .3)+
        geom_path(data = filter(mpa_tracks, overall_EM_beh == "not_hauling"), aes(x=lon_bin, y=lat_bin, key=name, labels=date), col = "gray80", size = 1)+
        geom_point(data = filter(mpa_tracks, overall_EM_beh == "not_hauling"), aes(x=lon_bin, y=lat_bin, key=name, labels=date), col = "white", size = 1)+
        geom_point(data = filter(mpa_tracks, overall_EM_beh == "hauling"), aes(x=lon_bin, y=lat_bin, key=name, labels=date), colour = "red", size = 2, alpha = .3)+
        coord_quickmap()+
        coord_sf(xlim = c(-116,-109), ylim = c(17,21))+
        facet_grid(~year)+
        theme_seaecolab()

plotly::ggplotly(p)

Red dots are suspicious activities, white dots are logging position (each hour), grey lines are inferred track of the fishing vessel.

The Revillagigedo protection polygon is a no fishing zone and all fishing activities are prohibited.