Koen Hufkens khufkens

Founder of BlueGreen Labs | PhD | Addressing climate change through data driven methods

khufkens / backup_snotel.R

Last active March 13, 2025 08:51

Backup USDA SNOTEL

	# This tiny script tries to grab all USDA SNOTEL data
	# for safekeeping. SNOTEL monitors snow and water across
	# the US but mostly in the west.
	# https://www.nrcs.usda.gov/resources/data-and-reports/snow-and-water-interactive-map
	#
	# The data downloaded are the historical reports
	# summarizing daily values of snow pack and additional
	# metrics. Reports lag real-time values with a couple
	# of weeks.
	#

khufkens / clean_gh-pages.sh

Created February 12, 2023 14:31

Create clean gh-pages branch

khufkens / format_rsofun_drivers.R

Created December 5, 2022 09:13

Format rsofun drivers

	#' Format and run p-model drivers on Euler
	#'
	#' Takes site information for a single
	#' or multiple sites and grabs all data
	#' required for a p-model run in rsofun.
	#'
	#' Parameter settings are provided as
	#' arguments, but could be altered after
	#' the fact if desired.
	#'

khufkens / aggregate_netcdf.R

Created November 29, 2022 15:19

NetCDF aggregation using {terra}

	library(terra)
	library(stringr)

	#---- static code ----

	# read in the original data
	r <- rast("LPX-Bern_SH1_gpppft.nc")

	# first subset to only the required years
	# 1978 - 2020

khufkens / batch_unzip.sh

Created November 22, 2022 11:12

Batch unzip in bash

find . -name '*.zip' -exec sh -c 'unzip {};'

khufkens / reference_blue_marble.sh

Created June 8, 2022 18:30

Georeferencing a NASA blue marble world map

gdal_translate -a_ullr -180 90 180 -90 -a_srs EPSG:4326 world.topo.bathy.200408.3x5400x2700.jpg tmp.tif

khufkens / opt_span_example.R

Created April 13, 2022 17:07

Calculate optimal span for MODIS time series

	# Interpolation using a loess smoother
	# and optimal span determination

	library(MODISTools)
	library(phenocamr)
	library(tidyverse)

	VI <- mt_subset(
	product = "MOD13Q1",
	lon = -72.1755,

khufkens / download_plumber2.R

Created November 4, 2021 15:22

Download Plumber 2 data

	library(rvest)

	# base html url
	url <- "http://dapds00.nci.org.au/thredds/fileServer/ks32/CLEX_Data/PLUMBER2/v1-0/Flux/"

	# catalogue
	html_page = read_html("https://dap.nci.org.au/thredds/remoteCatalogService?catalog=http://dapds00.nci.org.au/thredds/catalog/ks32/CLEX_Data/PLUMBER2/v1-0/Flux/catalog.xml")

	# grab files from the catalogue
	files = html_page %>%

khufkens / animation.sh

Last active January 26, 2022 17:57

Twitter movie from PNG sequence

	# convert pngs to mp4
	ffmpeg -r 2 -i tmp/%03d.png -y out.mp4

	# convert mp4 to twitter compatible format
	ffmpeg -i out.mp4 -vf scale=1280:720 -b:v 5000K -minrate 5000K -maxrate 5000K -b:a 128K -r 30 -f mp4 -vcodec libx264 -profile:v high -pix_fmt yuv420p -strict -2 -ac 2 -acodec aac twitter.mp4

khufkens / waffle_plot.R

Last active July 20, 2021 15:43

Complex waffle plot example