ggplot_projection_shapefile: Plot a raster (either a raster layer or a inla projection...

In timcdlucas/INLAutils: Utility Functions for 'INLA'

Plot a raster (either a raster layer or a inla projection matrix) with a shape file

Description

Plot a raster (either a raster layer or a inla projection matrix) with a shape file

Usage

ggplot_projection_shapefile(
  raster = NULL,
  projector = NULL,
  spatialpolygons = NULL,
  mesh = NULL,
  shapecol = "white"
)

Arguments

raster	Either a RasterLayer or a matrix (e.g. from inla.mesh.project)
projector	If raster is a matrix, a projector object (from inla.mesh.projector) is also needed.
spatialpolygons	A SpatialPolygonsDataFrame object.
mesh	An inla.mesh object to be plotted. Still rough.
shapecol	Colour for the shape file outlines.

Examples

## Not run: 
set.seed(2)
library(INLA)

# Create inla projector
n <- 20
loc <- matrix(runif(n*2), n, 2)
mesh <- inla.mesh.create(loc, refine=list(max.edge=0.05))
projector <- inla.mesh.projector(mesh)

field <- cos(mesh$loc[,1]*2*pi*3)*sin(mesh$loc[,2]*2*pi*7)
projection <- inla.mesh.project(projector, field)

# And the shape file
crds <- loc[chull(loc), ]
SPls <- SpatialPolygons(list(Polygons(list(Polygon(crds)), ID = 'a')))

ggplot_projection_shapefile(projection, projector, SPls) + theme_minimal()


# Alternatively plot a raster
library(raster)
raster <- raster(projection)
extent(raster) <- c(range(projector$x), range(projector$y))
ggplot_projection_shapefile(raster, spatialpolygons = SPls)

## End(Not run)

timcdlucas/INLAutils documentation built on Nov. 29, 2022, 5:41 a.m.