make.3Dgrid: Methods to prepare 3D prediction locations

Description

Generates a list of objects of type "SpatialPixelsDataFrame" with longitude, latitude and altitude coordinates (these names are used by default for compatibility with the geosamples-class).

Usage

# S4 method for SpatialPixelsDataFrame
make.3Dgrid(obj, 
      proj4s = get("ref_CRS", envir = GSIF.opts), 
      pixsize = get("cellsize", envir = GSIF.opts)[2], 
      resampling_method = "bilinear", 
      NAflag = get("NAflag", envir = GSIF.opts), 
      stdepths = get("stdepths", envir = GSIF.opts), 
      tmp.file = TRUE, show.output.on.console = TRUE, …)
# S4 method for RasterBrick
make.3Dgrid(obj, 
      proj4s = get("ref_CRS", envir = GSIF.opts), 
      pixsize = get("cellsize", envir = GSIF.opts)[2], 
      resampling_method = "bilinear", 
      NAflag = get("NAflag", envir = GSIF.opts), 
      stdepths = get("stdepths", envir = GSIF.opts), 
      tmp.file = TRUE, show.output.on.console = TRUE, …)

Arguments

obj

object of class "SpatialPixelsDataFrame" or "RasterBrick"

proj4s

character; proj4string describing the target coordinate system

pixsize

grid cell size in decimal degrees (set by default at 1/1200 (0.0008333333 or 100 m around equator)

resampling_method

character; resampling method to be passed the reprojection algorithm

NAflag

character; missing value flag

stdepths

numeric; list of standard depths

tmp.file

logical; specifies whether a temporary file name should be generated

show.output.on.console

logical; specifies whether to print out the progress

…

optional arguments that can be passed to the reprojetion algorithm

Value

The output is list of objects of class "SpatialPixelsDataFrame" where the number of elements in the list corresponds to the number of standard depths.

References

Bivand, R.S., Pebesma, E.J., and G?mez-Rubio, V., (2008) Applied Spatial Data Analysis with R. Springer, 378 p.
FWTools (http://fwtools.maptools.org)
gdalUtils package (https://CRAN.R-project.org/package=gdalUtils)
Raster package (https://CRAN.R-project.org/package=raster)

Examples

Run this code

# NOT RUN {
## grids Ebergotzen:
library(plotKML)
library(rgdal)
library(raster)

data(eberg_grid)
gridded(eberg_grid) <- ~x+y
proj4string(eberg_grid) <- CRS("+init=epsg:31467")
##  convert to spatial components:
formulaString <- ~ PRMGEO6+DEMSRT6+TWISRT6+TIRAST6
eberg_spc <- spc(eberg_grid, formulaString)
## create 3D locations in the original coordinate system:
eberg_3Dxy <- sp3D(eberg_spc@predicted)
# }
# NOT RUN {
## wrapper function to create 3D locations in the default WGS84 system:
eberg_3D <- make.3Dgrid(eberg_spc@predicted)
image(eberg_3D[[1]]["PC1"])
## downscale 100 m resolution imagery to 25 m:
data(eberg_grid25)
gridded(eberg_grid25) <- ~x+y
proj4string(eberg_grid25) <- CRS("+init=epsg:31467")
eberg_grid25@data <- cbind(eberg_grid25@data, 
    warp(eberg_grid, pixsize=eberg_grid25@grid@cellsize[1], 
    GridTopology=eberg_grid25@grid, resampling_method="cubicspline")@data)
## this function requires FWTools!
# }