Retrive SHAPE files

Share:

Description

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

Usage

1
2
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:

1
2
3
4
5
6
7
8
> 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írez. Parque Tecnológico Sartenejas, Universidad Simón Bolí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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## 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")

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.