Allan Just allanjust

following on from https://bsky.app/profile/jerry-shannon.bsky.social/post/3lcisatdppk2h

gdalinfo /vsigs/gcp-public-data-arco-era5/ar/full_37-1h-0p25deg-chunk-1.zarr-v3
ERROR 15: No valid GCS credentials found. For authenticated requests, you need to set GS_SECRET_ACCESS_KEY, GS_ACCESS_KEY_ID, 
GS_OAUTH2_REFRESH_TOKEN, GOOGLE_APPLICATION_CREDENTIALS, or other configuration options, or create a /home/ubuntu/.boto file. Consult
https://gdal.org/en/stable/user/virtual_file_systems.html#vsigs-google-cloud-storage-files for more details. For unauthenticated requests
on public resources, set the GS_NO_SIGN_REQUEST configuration option to YES.

	library(mapgl)

	mapboxgl(
	style = mapbox_style("streets"),
	center = c(-87.61694, 41.86625),
	zoom = 15.99,
	pitch = 40,
	bearing = 20,
	antialias = TRUE
	) \|>

	# UPDATE: Instead of using this script, you should probably
	# use http://web.archive.org/web/2021id_/https://modis.ornl.gov/files/modis_sin.kmz

	#!/usr/bin/env python3

	"""
	Draw the shapes of tiles used by Moderate Resolution Imaging
	Spectroradiometer (MODIS) satellite products. The output is a TSV file
	with the shapes as well-known text (WKT) that can be read into QGIS.
	"""

	library(tidyterra)
	library(tidyverse)
	library(terra)
	library(tigris)
	library(sf)
	library(mapboxapi)

	# Get data from https://www.mrlc.gov/data/nlcd-2019-land-cover-conus
	# Warning - it's big!
	nlcd <- rast("~/data/nlcd_2019_land_cover_l48_20210604/nlcd_2019_land_cover_l48_20210604.img")

	# Context: https://twitter.com/grant_mcdermott/status/1487528757418102787

	library(data.table)
	library(fixest)

	bboot =
	function(object, reps = 100L, cluster = NULL, ...) {

	fixest_obj = inherits(object, c('fixest', 'fixest_multi'))


	library(tidyverse)
	library(sf)
	library(spData)
	library(gganimate)
	data(world)

	# see also: https://proj.org/operations/projections/index.html
	projs <- list(
	"Mercator" = "+proj=merc",

	## Visualizing authorship networks: code for this blog post:
	## https://mathewkiang.com/2020/12/10/my-collaboration-network-for-2010-to-2020-other-plots

	## Imports ----
	library(scholar) # remotes::install_github("jkeirstead/scholar")
	library(tidyverse)
	library(fs)
	library(here)

	## Constants ----