examples/github_examples.R

In harveybarnhard/smoov: Simple Mapping of our Vicinities

# Set libraries
library(smoov)
library(ggplot2)
library(data.table)

# Set outpath for figures
outpath = file.path(.github, "smoov/examples")

# Plot all US counties =========================================================
actcom = function(pop, stayhome, bcycle, walk){
  # Calculate proportion of active commuters
  commuters = pop
  commuters[commuters<0] = NA
  active = 100*(bcycle + walk+1)/(commuters+1)
  
  # Winsorize
  quants = quantile(active, probs=c(0.05,0.95), na.rm=TRUE)
  active[active<quants[1]] = quants[1]
  active[active>quants[2]] = quants[2]
  return(
    active
  )
}

data(county_commute)
dtc = data.table(county_commute)
dtc = dtc[, active := actcom(pop_workers, commute_stayhome,
                             commute_bcycle, commute_walked)]

# Create fips code and plot
dtc[, fips := create_fips(state, county)]
usa = smoov("counties", data=dtc, value="active")
ggsave(file.path(outpath, "county_example.png"), plot=usa)

# Plot Boston surroundings =====================================================
data(tract_commute)
dtt = data.table(tract_commute)
dtt = dtt[, active := actcom(pop_workers, commute_stayhome,
                             commute_bcycle, commute_walked)]
dtt[, fips := create_fips(state, county, tract)]

# Plot Boston surroundings
east_MA =  c(1,5,17,21,23,25,027)
bos1 = smoov("tracts", data=dtt, value="active", states=25, counties=east_MA)
ggsave(file.path(outpath, "tract_example1.png"), plot=bos1)

# Plot Boston surroundings with titles and formatted legend ====================
bos2 = bos1 +
  labs(title="Active Commuters in Eastern Massachusetts",
       subtitle="% of Commuters who Walk or Cycle by Census Tract") +
  theme(plot.title = element_text(size=20, face="bold", hjust = 0.5),
        plot.subtitle = element_text(size=15, face="bold", hjust = 0.5),
        legend.position=c(0.8,0.9),
        legend.justification=c(0.8, 0.9))
ggsave(file.path(outpath, "tract_example2.png"), plot=bos2)

# Add county borders ===========================================================
bos3 = bos2 +
  smoov_border(
    geo="counties", states=25, counties=east_MA, line_size=1, line_color="black"
  )
ggsave(file.path(outpath, "tract_example3.png"), plot=bos3)

# Add county highlights for Suffolk County =====================================
not_suffolk = east_MA[east_MA!=25]
bos4 = bos2 +
  smoov_border(
    geo="counties", states=25, counties=not_suffolk,
    line_alpha=0, fill_color="white", fill_alpha=0.5
  )
ggsave(file.path(outpath, "tract_example4.png"), plot=bos4)


# osm examples =================================================================
med_streets <- sf::st_bbox(chi$data)%>%
  osmdata::opq()%>%
  osmdata::add_osm_feature(key = "place", value = c("city")) %>%
  osmdata::osmdata_sf()


labellayer = ggrepel::geom_text_repel(
  data = med_streets$osm_multipolygons,
  aes(label=name, geometry=geometry),
  stat="sf_coordinates",
  min.segment.length=0,
  fontface="bold",
  segment.color=NA
)

chi = chi+labellayer

# Interactive examples  ========================================================

bos5 = plotly::ggplotly(bos3)

widget_file_size = function(p) {
  d = tempdir()
  withr::with_dir(d, htmlwidgets::saveWidget(p, "index.html"))
  f = file.path(d, "index.html")
  mb = round(file.info(f)$size/1e6, 3)
  message("File is: ", mb, " MB")
}

widget_file_size(plotly::partial_bundle(bostest))