GEOG208 Final

This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code.

Try executing this chunk by clicking the Run button within the chunk or by placing your cursor inside it and pressing Cmd+Shift+Enter.

library("sf")

Linking to GEOS 3.6.1, GDAL 2.1.3, PROJ 4.9.3

Map1

head(fire_csv)

fire_csv <- read.csv("/Users/krgr.df/Desktop/GEOG208/Final Assignment/Fire Data/NASA VIIRs CSV/fire_nrt_V1_32986.csv")
str(fire_csv)

'data.frame':   1845 obs. of  14 variables:
 $ latitude  : num  34.2 34.1 34.2 34.2 34.2 ...
 $ longitude : num  -118 -118 -118 -118 -118 ...
 $ bright_ti4: num  303 298 345 299 297 ...
 $ scan      : num  0.41 0.41 0.47 0.42 0.39 0.39 0.56 0.34 0.47 0.4 ...
 $ track     : num  0.37 0.37 0.4 0.37 0.36 0.36 0.43 0.56 0.4 0.6 ...
 $ acq_date  : Factor w/ 14 levels "2018-11-08","2018-11-09",..: 2 2 2 3 8 8 9 10 13 14 ...
 $ acq_time  : int  948 948 2106 930 936 936 918 900 924 1024 ...
 $ satellite : Factor w/ 1 level "N": 1 1 1 1 1 1 1 1 1 1 ...
 $ instrument: Factor w/ 1 level "VIIRS": 1 1 1 1 1 1 1 1 1 1 ...
 $ confidence: Factor w/ 3 levels "h","l","n": 3 3 3 3 3 3 3 3 3 3 ...
 $ version   : Factor w/ 1 level "1.0NRT": 1 1 1 1 1 1 1 1 1 1 ...
 $ bright_ti5: num  286 285 313 282 277 ...
 $ frp       : num  0.5 0.3 6.4 0.6 0.3 0.7 0.8 0.6 0.7 0.8 ...
 $ daynight  : Factor w/ 2 levels "D","N": 2 2 1 2 2 2 2 2 2 2 ...

fire_csv$acq_date <- as.Date(fire_csv$acq_date)
hist(fire_csv$acq_date, "days", 
     freq = TRUE,
     breaks = 15,
     xlab = "Date Acquired",
     col = "firebrick3",
     border = "orange"
     )

##################Animated Map##################
install.packages('devtools')
devtools::install_github('thomasp85/gganimate')
library(ggplot2)
library(gganimate)


fire_csv$acq_date <- as.Date(fire_csv$acq_date)
str(fire_csv)

animated_map <- ggplot(fire_csv, (aes(longitude, latitude, color = bright_ti4))) +
  borders("county", ".*ventura",colour = "gray85", fill = "gray80", ) +
  theme_map() +
  borders("county", ".*los angeles",colour = "gray85", fill = "gray80", ) +
  theme_map() +
  geom_point(alpha = .5) +
  scale_color_gradient(low="yellow", high="red", name = "Brightness Value") + 
  labs(title = 'Date: {frame_time}', size = 3) +
  transition_time(acq_date) +
  enter_fade() + 
  ease_aes("linear")

animated_map
animate(animated_map, fps=1.75)
ease_aes('linear')

https://photos.google.com/share/AF1QipNZN2aIUorKUYKMkM5VbLYoWQUNbYvYlYubbeNgeh5SEIdfSBAkNcCuuEAu__RiHg/photo/AF1QipM8bCC5syWDq5n3_9v2GmFVQO7jzpSIIFAcZCIM?key=NDFycU16eExNSl81dUpRLXAwVEtWZmtFMkhuX3V3

##################Satellite Imagery GIF##################
#https://cran.r-project.org/web/packages/magick/vignettes/intro.html

before_fire <- image_scale(image_read("/Users/krgr.df/Desktop/Screen Shot 2018-12-04 at 2.20.06 AM.png"), "x450")
before_points <- image_scale(image_read("/Users/krgr.df/Desktop/Screen Shot 2018-12-04 at 2.07.35 AM.png"), "x450")
with_points <- image_scale(image_read("/Users/krgr.df/Desktop/Screen Shot 2018-12-04 at 2.07.52 AM.png"), "x450")
frames <- image_morph(c( before_fire, with_points, before_points, before_points, before_points, before_fire), frames = 20)
image_animate(frames)

https://photos.google.com/share/AF1QipNZN2aIUorKUYKMkM5VbLYoWQUNbYvYlYubbeNgeh5SEIdfSBAkNcCuuEAu__RiHg/photo/AF1QipNkvlR-hjBYCTTb4qZYSz-cO5hbNM0GNOsES2rP?key=NDFycU16eExNSl81dUpRLXAwVEtWZmtFMkhuX3V3