Generate a cost Raster

Share:

Description

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

Usage

1
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

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.

Value

RasterLayer

Author(s)

Joseph Stachelek

See Also

spTransform-methods, rasterize

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
## 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)