level.plot: Level plot for gridded data

In GISTools: Some further GIS capabilities for R

Description

Draws a level plot given a SpatialPixelsDataFrame, an index and a shading scheme.

Usage

level.plot(grd, shades, index=1, add=FALSE)

Arguments

grd	A spatialPixelsDataFrame object.
shades	A shading scheme created by shading or auto.shading. If omitted, chosen automatically from grd.
index	Index giving the variable in grd to plot.
add	Whether to add the level plot to an existing plot.

Details

The function returns no value, but draws a level plot on the current graphics device as a side effect.

Value

None (see above).

Author(s)

Chris Brunsdon

Examples

# Data for New Haven to use in example
data(newhaven)
# Do the KDE
breach.dens = kde.points(breach,lims=tracts)
# Plot the result
level.plot(breach.dens)
# Block out the part outside the study area
masker = poly.outer(breach.dens,tracts,extend=100); add.masking(masker)
# Plot census tract boundaries
plot(tracts,add=TRUE)

Example output

Loading required package: maptools
Loading required package: sp
Checking rgeos availability: TRUE
Loading required package: RColorBrewer
Loading required package: MASS
Loading required package: rgeos
rgeos version: 0.3-23, (SVN revision 546)
 GEOS runtime version: 3.4.2-CAPI-1.8.2 r3921 
 Linking to sp version: 1.2-3 
 Polygon checking: TRUE 

NOTE: rgdal::checkCRSArgs: no proj_defs.dat in PROJ.4 shared files

GISTools documentation built on May 6, 2019, 1:09 a.m.