costrasterGen: Generate a cost Raster

Description

Generate a cost raster from an object of class SpatialPolygons, matrix, or SpatialPointsDataFrame

Usage

costrasterGen(xymat, pols, extent = "polys", projstr)

Arguments

xymat

Matrix of coordinates or a SpatialPointsDataFrame object

pols

SpatialPolygons object

extent

Define extent based on extent of xymat/xyspdf (points) or pols (polys). Default is polys.

projstr

proj4 string defining the inherent projection

Value

RasterLayer

Details

Ensure that the projection of the xymat coordinates and pols match. This can be accomplished by running the projection command on both inputs. If they do not match use the spTransform command.

Examples

Run this code

## Not run: 
# Sr1 <- Polygon(cbind(c(0, 0, 1, 1, 0), c(0, 12, 12, 0, 0)))
# Sr4 <- Polygon(cbind(c(9, 9, 10, 10, 9), c(0, 12, 12, 0, 0)))
# Sr2 <- Polygon(cbind(c(1, 1, 9, 9, 1), c(11, 12, 12, 11, 11)))
# Sr3 <- Polygon(cbind(c(1, 1, 9, 9, 1),c(0, 1, 1, 0, 0)))
# Sr5 <- Polygon(cbind(c(4, 4, 5, 5, 4), c(4, 8, 8, 4, 4)))
# Srs1 <- Polygons(list(Sr1), "s1")
# Srs2 <- Polygons(list(Sr2), "s2")
# Srs3 <- Polygons(list(Sr3), "s3")
# Srs4 <- Polygons(list(Sr4), "s4")
# Srs5 <- Polygons(list(Sr5), "s5")
# 
# pols <- SpatialPolygons(list(Srs1, Srs2, Srs3, Srs4, Srs5), 1:5)
# 
# #using a matrix object
# xymat <- matrix(3, 3, nrow = 1, ncol = 2)
# costras <- costrasterGen(xymat, pols, projstr = NULL)
# 
# #plotting
# plot(costras)
# points(xymat)
# ## End(Not run)