Can ggplot use two coordinate systems to add a rose diagram to ggmap basemap?

ggplot2

#1

Hello All,
I am trying to add a rose diagram of wind speed on top of an aerial base map. I am stumped thus far because the method for making a rose diagram is to use coord_polar() on a geom_histogram(). However, the use of a polar coordinate system is incompatible with using a geographically referenced base map as is done with ggmap. Does any one have as solution to either: 1) use two coordinate systems on one plot, 2) have a different method for making a rose diagram that does not require coord_polar(), or 3) some awesome idea I have not tried yet?

Thanks!

map with wind speed and direction

rose diagram of wind speed histogram

Reproducible example:

 library("sf")
library("ggplot2")
library("ggmap")
library("tidyr")
library("dplyr")

####### functions #########
# functions to convert speed and direction to geographic coordinates
speed2u<-function(s,d){
  u <- -s * sin(d * pi/180)
  return (u) 
}
speed2v<-function(s,d){
  v <- -s * cos(d * pi/180)
  return (v)
}
########################

####### Data Simulation and Cleaning ###
# create fake data set
dataset <- data.frame(id = seq(1,300,1),
                      wind_speed = rnbinom(300,5,0.25),
                      wind_dir   = rnorm(300,0,20)-108)

# location of wind sensor (lat, long)
sensor_loc <- c(39.9524, -75.1636)
# transform data
dat <- dataset %>% 
  mutate(wind_u      = speed2u(wind_speed, wind_dir),
         wind_v      = speed2v(wind_speed, wind_dir),
         lng         = sensor_loc[2],
         lat         = sensor_loc[1],
         xend        = lng + wind_u * 0.00015,
         yend        = lat + wind_v * 0.00015) %>% 
  mutate(wind_to     = case_when(wind_dir <= 180 ~ wind_dir+180,
                                 wind_dir > 180  ~ wind_dir-180))
# create sf object from sensor location
sensor_sf <- data.frame(lng = sensor_loc[2], lat = sensor_loc[1]) %>% 
  st_as_sf(., coords = c("lng","lat"), crs = 4326)
# retrieve basemap with ggmap
base_map = get_map(location = unname(st_bbox(st_buffer(sensor_sf, 0.005))),
                   source = "google", maptype = "hybrid")

###### Make plots
# plot individual wind speed observations on map
ggmap(base_map) +
  geom_segment(data = dat, aes(x=lng, y=lat, xend=xend, yend=yend, color = wind_speed)) +
  geom_sf(data = sensor_sf, inherit.aes = FALSE, color = "black", size = 3) +
  scale_color_viridis_c(option = "A", name = "wind speed (mph)") +
  labs(x = NULL, y = NULL) +
  theme(rect = element_blank(), 
        axis.text = element_blank(), 
        axis.title = element_blank(), 
        axis.ticks.length = unit(0, "pt")
  )
# Plot rose diagram
ggplot(dat, aes(x = wind_to)) +
  coord_polar(theta = "x", start = 0, direction = 1) +
  geom_histogram(fill = "red", color = "gray10", bins = 30) +
  scale_x_continuous(breaks = seq(0, 360, 10), limits = c(0, 360)) +
  theme_minimal() +
  theme(
    axis.text.y = element_blank(),
    axis.title = element_blank())

#2

Hi Matt,

I think you could try plotting the rose diagram as a grob and then use annotation_custom to put the grob on top of the map. Maybe. I've never actually done it.

This might help:

Ron.


#3

cowplot can also draw independent plots on top of each other. The intro vignette has lots of examples:
https://cran.r-project.org/web/packages/cowplot/vignettes/introduction.html

Were you hoping for the rose diagram to be scaled to the map's scale in any way? Or just sized based on legibility?