Meteorite Falls and Finds in the UK

[This article was first published on R, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

(Tidy Tuesday is a project to supply weekly data sets for R users to practice their coding skills on. You can find full details here.)

The Tidy Tuesday dataset for this week, compiling meteorite strikes across the world, grabbed my interest so I decided to finally dive in and give it a go. I ended up looking at strikes in the UK, noted by mass.

I exported the original code image as a PDF and then applied some effects in Affinity Photo to give it an old, weathered look. This was my first time experimenting with something like this, so it’s a little sketchy, but I enjoyed seeing how I could pull a plot apart in Photo to apply image effects (although sorry, this does mean the code for this plot doesn’t get you all the way to the finished produced).

Code:

library(tidyverse)
library(maps)

meteorites <- readr::read_csv("https://raw.githubusercontent.com/rfordatascience/tidytuesday/master/data/2019/2019-06-11/meteorites.csv")

ukFalls <- meteorites %>% 
  filter(lat >= 50 & lat <= 60 & long <= 2 & long >= -11) %>% 
  filter(!name %in% c("Limerick", "Dundrum", "Mooresfort", "Leighlinbridge", "Pettiswood")) %>%
  mutate(massKG = mass/1000) %>% 
  mutate(nameYear = paste(name, year, sep = ", ")) %>% 
  arrange(year)

ukMap <- map_data(regions = "uk", map = "world") %>% filter(lat <=58.73)

ggplot(ukMap, aes(x = long, y = lat)) + 
  geom_polygon(aes(group = group), colour = "black", fill = "#222831") +
  geom_point(data = ukFalls, aes(x = long, y = lat, size = massKG), colour = "#d65a31") +
  geom_text(data = subset(ukFalls, massKG >=10), aes(label = nameYear), hjust = 0, nudge_x = 0.2, size = 3) +
  scale_size(range = c(0, 5), breaks = c(5,15,25,50), labels = c("5 kg", "15 kg", "25 kg", "50 kg")) +
  theme_void() +
  theme(plot.background = element_rect(fill = "#5588a3"),
        plot.margin = margin(.5, .5, .5, .5, "cm")) +
  labs(size = "Mass",
       alpha = "Mass")

ggsave("map.pdf")

To leave a comment for the author, please follow the link and comment on their blog: R.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)