This is a selection of my (open source) tools and code as listed on to my GitHub page. If you find and correct bugs / add functionality please ask to merge this into the master code. If you find my code and tools useful buy me a coffee should we meet.


The SNOTEL R package allows for easy downloading of SNOTEL data directly into your R workspace. It also provides additional functions to summarize the data in terms of snow phenology patterns.


The Google Earth Engine (GEE) code within this project combines MODIS Aqua – Terra snow cover products and derives snow melt – accumulation (snow phenology) dates.


The Ameriflux R toolbox is a collection of functions to facilitate the retrieval and processing and visualization of Ameriflux data.


Functions are provided to (batch) download single pixel or gridded Daymet data directly into your R workspace, or save them as csv/tif files on your computer. Both a batch version as a single download version are provided. Gridded data downloads for a region of interest are specified by a top left / bottom right coordinate pair or a single pixel location.

PhenoCam Installation Tool (PIT)

PhenoCam Installation Tool (PIT) is a set of scripts for Linux/Mac OSX and Windows taking care of the settings as needed by cameras installed by or associated with the PhenoCam network.

NSIDC IMS / CMC snow cover and depth composites

National Snow and Ice Data Center (NSIDC) processing software. Composites raw NSIDC data files into georeferenced (yearly) stacked geotiffs for easy processing. Additional functions are provided to process this data to extract derived data products such as day of first snow melt and the dates spanning the longest snow free period.

ORNL DAAC MODIS Land Product Subsets – Perl Tool

This is a perl client for the Oak Ridge National Laboratories (ORNL) MODIS Land Product Subsets project. It allows for the extraction of time series of several pixels around geographic locations for select MODIS products. The service makes use of the Simple Object Access Protocol (SOAP) web service provided by ORNL (see requirements below).

FOTO toolbox

FOTO (Fourier Transform Textural Ordination) method as used in Couteron et al. 2005, Proisy et al. 2007, Barbier et al. 2010 and Ploton et al. 2012. The FOTO method uses a principal component analysis (PCA) on radially averaged 2D fourier spectra to characterize (greyscale) image texture. Although the techique as presented in these papers is applied on a canopy level, the principle works as well on a smaller scale (at a higher resolution). A proof of concept using a UAV over a grassland can be found in one of my blog posts.

fixed effects Spatial Panel Analysis Method (SPAM)

This is a fixed effects panel analysis function. It uses a method from econometrics to detect trends over time across a window of pixels. This increases the robustness of trend analysis as compared to a pixel by pixel analysis. An added benefit is the speed you gain as you downscale the original data by a factor defined by the panel size used.


Functions to (batch) download single pixel Daymet data directly into your python workspace, or save them as csv files on your computer. Both a batch version as a single download version are provided. Consider downloading gridded data if you download extensive coverage within a single region.