Plot with earthquake data

This package allows you to clean the data then plot the data with either timeline or map fashion.

Clean the data

The eq_clean_data function provide necessary functionalities to transform the data into a form that can be plotted with ggplot2 or leafleat.

library(noaatk)
library(ggplot2)
library(leaflet)
library(dplyr,quietly = T)

cleaned <- eq_clean_data(noaa)

Draw Timeline plot

There are 2 ggplot2 goems in the package that allow you to draw Timeline graph. One without label, while the other provide location labels.

plot_dat <- cleaned %>% filter(DATE< lubridate::ymd("2009-1-1") & DATE > lubridate::ymd("2005-1-1")) %>% filter(COUNTRY %in% c("CHINA", "JAPAN", "USA"))
ggplot(plot_dat, aes(x = DATE))+ geom_timeline(aes(fill = DEATHS, size = EQ_MAG_MS)) + theme_bw()
ggplot(plot_dat, aes(x = DATE,y = COUNTRY))+ geom_timelinelabel(aes(fill = DEATHS, size = EQ_MAG_MS, loc = LOCATION_NAME_)) + theme_bw()

Draw Map with earthquke annotations

eq_map function allow you to use the power of leaflet to produce an interactive map with earthquakes ploted at their 'quake-location'. There is also a very useful helper function eq_create_label that can be used to generate labels.

map_dat <- plot_dat %>% filter(COUNTRY == "CHINA")
labeled_map_dat <- map_dat %>% mutate(pop_text = eq_create_label(.))
labeled_map_dat %>% eq_map(DATE)
labeled_map_dat <- map_dat %>% mutate(pop_text = eq_create_label(.))


esppk/noaatk documentation built on May 25, 2019, 5:21 p.m.