Koen Hufkens khufkens

## equi2polar.sh
#!/bin/bash
# convert equirectangular 2 fisheye

files=`ls *.jpg`

for i in $files;
do
noext=`echo $i | cut -d'.' -f1`

echo $noext

## downsample_extract_hwsd.R
# Downsample HWSD data to the scale of 0.5 degree
# WFDEI / CRU data

library(raster)
library(hwsdr)
library(tidyverse)

# get mu_global at 0.5 degrees from the
# ORNL DAAC (which is easier than reading in
# that monstrous BIL file)

## waffle_plot.R
# data wrangling
library(tidyverse)

# waffle plot
library(waffle)

# combining graphs
library(patchwork)

# text magic

## download_plumber2.R
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 %>%

## animation.sh
# 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

## opt_span_example.R
# Interpolation using a loess smoother
# and optimal span determination

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

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

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

## polar_map.r
#-------- arctic map

library(sp)
library(maps)
library(rgeos)

# function to slice and dice a map and convert it to an sp() object
maps2sp = function(xlim, ylim, l.out = 100, clip = TRUE) {
  stopifnot(require(maps))
  m = map(xlim = xlim, ylim = ylim, plot = FALSE, fill = TRUE)

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

## aggregate_netcdf.R
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
	#!/bin/bash
	# convert equirectangular 2 fisheye

	files=`ls *.jpg`

	for i in $files;
	do
	noext=`echo $i \| cut -d'.' -f1`

	echo $noext
	# Downsample HWSD data to the scale of 0.5 degree
	# WFDEI / CRU data

	library(raster)
	library(hwsdr)
	library(tidyverse)

	# get mu_global at 0.5 degrees from the
	# ORNL DAAC (which is easier than reading in
	# that monstrous BIL file)
	# data wrangling
	library(tidyverse)

	# waffle plot
	library(waffle)

	# combining graphs
	library(patchwork)

	# text magic
	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 %>%
	# 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
	# Interpolation using a loess smoother
	# and optimal span determination

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

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

	library(sp)
	library(maps)
	library(rgeos)

	# function to slice and dice a map and convert it to an sp() object
	maps2sp = function(xlim, ylim, l.out = 100, clip = TRUE) {
	stopifnot(require(maps))
	m = map(xlim = xlim, ylim = ylim, plot = FALSE, fill = TRUE)
	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