| library(raster) | |
| library(rhdf5) | |
| library(rgdal) | |
| library(maps) | |
| # manually define the working directory containing the Black Marble VNP46A1 files | |
| wd <- "C:\\myblackmarblefiles" | |
| # List h5 files in the specified directory | |
| setwd(wd) |
| @echo off | |
| set OSGEO4W_ROOT=C:\OSGeo4W64 | |
| set basis_url="https://www.geodatenportal.sachsen-anhalt.de/gfds/ws/wfs/942f5d74-6c2b-263a/GDI-LSA_Schutzgebiete/ows.wfs?REQUEST=GetCapabilities&SERVICE=WFS" | |
| set layer="sg:lau_br_lsa" | |
| call "%OSGEO4W_ROOT%\bin\o4w_env.bat" |
| # Convert raster to polygon using GDAL | |
| ## edited 2019-03-27 to use sf::st_read | |
| polygonise <- function(x, outshape=NULL, gdalformat = 'ESRI Shapefile', | |
| pypath=NULL, readpoly=TRUE, quietish=TRUE) { | |
| ## x: an R Raster layer, or the file path to a raster file recognised by GDAL | |
| ## outshape: the path to the output shapefile (if NULL, a temporary file will be created) | |
| ## gdalformat: the desired OGR vector format. Currently only supports ESRI Shapefile! | |
| ## pypath: the path to gdal_polygonize.py (if NULL, an attempt will be made to determine the location | |
| ## readpoly: should the polygon shapefile be read back into R, and returned by this function? (logical) | |
| ## quietish: should (some) messages be suppressed? (logical) |
| ' This macro creates a "unified inbox" - a search across all inboxes and sent items. | |
| ' Steps to use it: | |
| ' 1. Enable developer tab to see VBA console: | |
| ' Options > Customize Ribbon > | |
| ' Enable the selection field before “Developer” and OK. | |
| ' 2. Add macro: | |
| ' From the Developer tab open Visual Basic dialog | |
| ' Open VBA Add the Sub below to General context | |
| ' 3. Enable macros: | |
| ' Options > Trust Center > Trust Center Settings > Macro Settings |
| # The function that produces the colour matrix | |
| colmat <- function(nbreaks = 3, breakstyle = "quantile", | |
| upperleft = "#0096EB", upperright = "#820050", | |
| bottomleft = "#BEBEBE", bottomright = "#FFE60F", | |
| xlab = "x label", ylab = "y label", plotLeg = TRUE, | |
| saveLeg = FALSE) { | |
| # TODO - replace any tidyr, dplyr etc. functions with data.table # | |
| library(tidyverse) | |
| require(ggplot2) | |
| require(classInt) |
| diverge0 <- function(p, ramp) { | |
| # p: a trellis object resulting from rasterVis::levelplot | |
| # ramp: the name of an RColorBrewer palette (as character), a character | |
| # vector of colour names to interpolate, or a colorRampPalette. | |
| require(RColorBrewer) | |
| require(rasterVis) | |
| if(length(ramp)==1 && is.character(ramp) && ramp %in% | |
| row.names(brewer.pal.info)) { | |
| ramp <- suppressWarnings(colorRampPalette(brewer.pal(11, ramp))) | |
| } else if(length(ramp) > 1 && is.character(ramp) && all(ramp %in% colors())) { |
| stratified <- function(df, group, size, select = NULL, | |
| replace = TRUE, bothSets = FALSE) { | |
| if (is.null(select)) { | |
| df <- df | |
| } else { | |
| if (is.null(names(select))) stop("'select' must be a named list") | |
| if (!all(names(select) %in% names(df))) | |
| stop("Please verify your 'select' argument") | |
| temp <- sapply(names(select), | |
| function(x) df[[x]] %in% select[[x]]) |
Collection of License badges for your Project's README file.
This list includes the most common open source and open data licenses.
Easily copy and paste the code under the badges into your Markdown files.
Translations: (No guarantee that the translations are up-to-date)
| library(foreign) | |
| library(data.table) | |
| library(zoo) | |
| library(Hmisc) | |
| library(gstat) | |
| library(sp) | |
| rm(list=ls()) | |
| # Fill in missing data ---- | |
| sm.thresh <- 0.9 # Min proportion of obs that must be in a station-month |
| #' Source specific lines in an R file | |
| #' | |
| #' @param file character string with the path to the file to source. | |
| #' @param lines numeric vector of lines to source in \code{file}. | |
| source_lines <- function(file, lines){ | |
| source(textConnection(readLines(file)[lines])) | |
| } |
