| all <- iconvlist() | |
| names(all) <- all | |
| tries <- llply(all[330:350], function(cur) iconv(song$artist, cur, "UTF-8")) | |
| vapply(tries, function(x) sum(is.na(x)), integer(1)) |
This package contains a set of tools, utilities and stylesheets for working with the Geoscience Australia Australian Bathymetry and Topography Grid.
Here is a sample rendering of this Geoscience Australia Australian Bathymetry and Topography Grid using the process documented in this repository.
| rasterToKML <- function(x, out='Rgrid', png=FALSE, name=out, proj="+proj=longlat") { | |
| require(raster) | |
| require(maptools) | |
| require(rgdal) | |
| require(rgeos) | |
| x <- raster(x) | |
| SG <- as(x, 'SpatialGridDataFrame') | |
| proj4string(SG) <- CRS(proj) | |
| SG <- GE_SpatialGrid(SG) | |
| if(png) png <- sprintf('%s.png', out) else png <- tempfile() |
| # This file contains two functions: | |
| # - combine.samples (used by par.trace.samples) | |
| # - par.trace.samples | |
| # and an example at the bottom. | |
| ## Author: Mikhail Popov (mikhail [at] mpopov.com) | |
| # install.packages("rjags") # JAGS must be installed on system | |
| # install.packages("doMC") # Unix only |
| import os | |
| from osgeo import gdal, ogr, osr | |
| import sys | |
| def initFields(in_lyr, out_lyr): | |
| #Arbitrarily get the first feature | |
| feat = in_lyr[0] | |
| #loop over each field |
| #' Transform raster to data.table | |
| #' | |
| #' @param x Raster* object | |
| #' @param row.names `NULL` or a character vector giving the row names for the data frame. Missing values are not allowed | |
| #' @param optional logical. If `TRUE`, setting row names and converting column names (to syntactic names: see make.names) is optional | |
| #' @param xy logical. If `TRUE`, also return the spatial coordinates | |
| #' @param centroids logical. If TRUE return the centroids instead of all spatial coordinates (only relevant if xy=TRUE) | |
| #' @param sepNA logical. If TRUE the parts of the spatial objects are separated by lines that are NA (only if xy=TRUE and, for polygons, if centroids=FALSE | |
| #' @param ... Additional arguments (none) passed to `raster::as.data.frame` | |
| #' |
| polygonizer <- 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 | |
| # 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) | |
| if (isTRUE(readpoly)) require(rgdal) | |
| if (is.null(pypath)) { |
| sample0110 <- function(size, n, complete.only=FALSE) { | |
| size <- as.integer(size) | |
| n <- as.integer(n) | |
| if(size > 25 || size < 3L) stop( | |
| "Currently size min is 3 and max is 25, though should be possible to allow ", | |
| "smaller and larger with some changes" | |
| ) | |
| # Generate integer pool and weights |
| # http://stackoverflow.com/a/30781090/2725969 | |
| sample0110b <- function(size, n) { | |
| size <- as.integer(size) | |
| n <- as.integer(n) | |
| if(size > 25 || size < 3L) stop("Size out of valid range") | |
| # Generate integer pool and weights |
| library(raster) | |
| library(rasterVis) | |
| library(grid) | |
| r <- raster(system.file("external/test.grd", package="raster")) | |
| ## Main graphic | |
| p1 <- levelplot(r) | |
| ## Inset graphic | |
| p2 <- levelplot(r, | |
| margin = FALSE, |