Koen Hufkens khufkens
-
Ghent University
- Antwerp, Belgium
- Sign in to view email
- http://www.khufkens.com
View bad_projection_read.Rmd
--- | |
title: "Faulty projection read" | |
output: html_document | |
--- | |
```{r setup, include=FALSE} | |
knitr::opts_chunk$set(echo = TRUE) | |
``` | |
The default projection data in returned tiles seems to be faulty. Mainly the second parallel North is set to 45 where it should be 60. Below is a reproducible example of the issue. |
View gist:5fe99a9606f544b2d5d6a356a3f93fdd
#!/bin/bash | |
# Insert PDF meta-data into exif header | |
# for pdf files with format author - year - title.pdf. | |
# This allows e-readers to query meta-data for proper | |
# display and sorting of authors and titles | |
# (note: relies on exiftool) | |
for filename in *.pdf; do | |
noext=`echo "${filename%.*}"` |
View el_nino_precip_analysis.r
# map El Nino drought in the Amazon and Africa | |
library(raster) | |
library(tidyverse) | |
library(maps) | |
library(RColorBrewer) | |
# set colour theme | |
colours = brewer.pal(11,"RdBu") | |
# download CHIRP data |
View daymet_server_tests.r
# load required libraries | |
if(!require(devtools)){install.package("devtools")} | |
if(!require(daymetr)){devtools::install_github("khufkens/daymetr")} | |
library("daymetr") | |
# no delay calls | |
no_delay = unlist(lapply(1:100, function(...){ | |
# download data | |
error = try(download_daymet(silent = TRUE)) |
View MCD12Q1_median_class.js
// years to process (from start year t0 to end year t1) | |
var t0 = "2001"; | |
var t1 = "2014"; | |
var LC = ee.ImageCollection('MCD12Q1') | |
.select('Land_Cover_Type_1') | |
.filterDate(t0.concat("-01-01"),t1.concat("-12-31")) | |
.median(); | |
// Create a geometry representing an export region. |
View tropicos_species_distribution.r
#' @param species: genus species or genus | |
#' @param quiet: TRUE / FALSE provides verbose output | |
#' @keywords Tropicos, species distribution | |
#' @examples | |
#' | |
#' # with defaults, outputting a data frame with species distribution | |
#' # for Clematis | |
#' df <- tropicos.species.distribution() | |
#' # returns NA if no data are present | |
#' [requires the rvest package for post-processing] |
View robinson_map.r
library(raster) | |
library(maps) | |
library(maptools) | |
library(sf) | |
library(scales) | |
# set coordinate systems | |
robinson = CRS(" +proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs") | |
latlon = CRS("+init=epsg:4326") |
View 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) |
View extract_MODIS_phenology_layer.r
# calculate adjusted phenology dates | |
# for greenness onset and maximum | |
extract_MODIS_phenology_layer <- function(value="increase",tiles_file="mytiles.txt"){ | |
# load required libraries | |
require(raster) | |
require(MODIS) # to import the SDS layers | |
# list all hdf files, and extract unique years |
View swath2grid.sh
#!/bin/bash | |
# | |
# swath to grid conversion for | |
# MOD04 reflectance data | |
# | |
# get the reprojection information | |
gdal_translate -of VRT HDF4_EOS:EOS_SWATH:"$1":mod04:Mean_Reflectance_Land land.vrt | |
gdal_translate -of VRT HDF4_EOS:EOS_SWATH:"$1":mod04:Mean_Reflectance_Ocean ocean.vrt |
NewerOlder