# NOT RUN {
# Workflow: create a dynamic tileset for dot-density mapping
library(tidycensus)
library(sf)
library(mapboxapi)
# Get population data for Census tracts in Vermont
vt_population <- get_decennial(
geography = "tract",
variables = "P001001",
state = "Vermont",
year = 2010,
geometry = TRUE
)
# Convert to representative dots - 1 per person
vt_dots <- st_sample(
vt_population,
size = vt_population$value
)
# Use tippecanoe to create dynamic tiles
tippecanoe(
input = vt_dots,
output = "vt_population.mbtiles",
layer_name = "vermont_population",
max_zoom = 18,
drop_rate = 1.5
)
# Upload to your Mapbox account for visualization
# A Mapbox secret access token must be set with mb_access_token()
# to upload data to your account
upload_tiles(
input = "vt_population.mbtiles",
username = "kwalkertcu",
tileset_id = "vt_population_dots",
multipart = TRUE
)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab