get.shape.state: Retrive SHAPE files

In vetools: Tools for Venezuelan Environmental Data

Description

These functions retreive the necesary SHAPE files to display the Venezuelan political border or any combination of states.

Usage

get.shape.state(abb, shape.file = "venezuelaestados")
get.shape.venezuela(shape.file = "venezuela")

Arguments

abb	a vector of characters containing the two letter abbreviations of the states to load.
shape.file	the base name of the SHAPE file to use.

Details

If the parameter abb is missing, then a data frame is shown and returned containing the states names, abbreviations and SHAPE file IDs.

Value

returns an object of class "SpatialPolygonsDataFrame" that can be plotted using the plot command.

Note

SHAPE files venezuela and "venezuela estados" have a slight size mismatch:

> VE <- get.shape.venezuela()
> VS <- get.shape.state(get.shape.state()$Abb)
> get.shape.range(VE)
      Long.start Long.end Lat.start  Lat.end
SHAPE   -73.3774 -59.7991 0.6498817 12.20108
> get.shape.range(VS)
       Long.start Long.end Lat.start  Lat.end
SHAPE   -73.37749 -59.7991 0.6492503 12.2012

Author(s)

R. Ram<c3><ad>rez. Parque Tecnol<c3><b3>gico Sartenejas, Universidad Sim<c3><b3>n Bol<c3><ad>var. Venezuela.

Wrapped in R by A. M. Sajo-Castelli

References

Maps where constructed and exported from ArcGIS 2.x.

See Also

get.shape.range, get.Grid.size.

Examples

## Get national boudary SHAPE
VE <- get.shape.venezuela()
## Not run: plot(VE, asp=1, axes=T)

## Get list of all available shapes
get.shape.state()

## Get national and statal boudaries SHAPE
VS <- get.shape.state(get.shape.state()$Abb)
## Not run: plot(VS, col="gray80", asp=1, axes=F)

## Retrieve Zone III states
BOAMDA = get.shape.state(c("BO", "AM", "DA"))
## Not run: plot(BOAMDA, add=T, border="darkred", lwd=2, col="pink")

vetools documentation built on May 2, 2019, 10:15 a.m.