commit 7bc5ac27f8f6007795a175ae46a65f1bb500300e
Author: Matt <bedrocks@gmail.com>
Date:   Sun Jan 3 17:09:28 2021 -0600

    Init.

diff --git a/.gitignore b/.gitignore
new file mode 100644
index 0000000..5b6a065
--- /dev/null
+++ b/.gitignore
@@ -0,0 +1,4 @@
+.Rproj.user
+.Rhistory
+.RData
+.Ruserdata
diff --git a/arps-r.Rproj b/arps-r.Rproj
new file mode 100644
index 0000000..066341e
--- /dev/null
+++ b/arps-r.Rproj
@@ -0,0 +1,13 @@
+Version: 1.0
+
+RestoreWorkspace: Default
+SaveWorkspace: Default
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 4
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX
diff --git a/bin/parse_aprs.R b/bin/parse_aprs.R
new file mode 100644
index 0000000..264af15
--- /dev/null
+++ b/bin/parse_aprs.R
@@ -0,0 +1,45 @@
+
+library(dplyr)
+library(ggplot2)
+library("leaflet")
+
+# Read data
+d_raw <- readLines("data/APRS data from Pembina Gorge 2020-09-18.txt") %>% 
+            as.data.frame() %>% 
+            setNames("packet")
+
+# Grab lines with location
+d <- d_raw %>% 
+        # Keep only rows with locations
+        filter(grepl(".*:=(.*)>(.*)", packet)) %>% 
+        # Pull out raw location parts
+        mutate(location_ns     = gsub(".*:=(.*)[NS]/(.*)", "\\1", packet) %>% as.numeric(),
+               location_ns_dir = gsub(".*[0-9]([NS])/(.*)", "\\1", packet),
+               location_ew     = gsub(".*/(.*)[EW]>(.*)", "\\1", packet) %>% as.numeric(),
+               location_ew_dir = gsub(".*[0-9]([EW])>(.*)", "\\1", packet)) %>% 
+        # Set directions and fix to long/lat
+        mutate(lat = if_else(location_ns_dir == "N",
+                             location_ns / 100,
+                             -location_ns / 100),
+               lon = if_else(location_ew_dir == "E",
+                             location_ew / 100,
+                             -location_ew / 100))
+
+# Simple map
+ggplot() +
+    geom_point(data = d,
+               aes(x = lon,
+                   y = lat))
+
+# Interactive map
+# https://towardsdatascience.com/making-interactive-maps-in-r-with-less-than-15-lines-of-code-bfd81f587e12
+d %>% 
+    leaflet() %>%
+    addProviderTiles(providers$Esri.WorldImagery, group = "World Imagery") %>%
+    addProviderTiles(providers$Stamen.TonerLite, group = "Toner Lite") %>%
+    addLayersControl(baseGroups = c("Toner Lite", "World Imagery")) %>%
+    addMarkers(label = NA, 
+               popup = d$packet) # Val False
+    
+
+