GEOTOPO {GEOmap} | R Documentation |
Extract subset of a topographic database, interpolate and plot using the persp program.
GEOTOPO(TOPO, PLOC, PROJ, calcol=NULL, npoints=500)
TOPO |
list of x,y,z for a DEM |
PLOC |
Location list, includes vectors LON and Lat |
PROJ |
projection |
calcol |
color table for coloring elevations above sea level |
npoints |
number of points in x-y grid, default=500 |
The return matrix PMAT is a rotation matrix used for adding geographic (projected) data onto the perspective plot.
PMAT |
Matrix from persp, used for adding other geographic information |
Jonathan M. Lees<jonathan.lees.edu>
subsetTOPO, settopocol, subsetTOPO, persp, DOTOPOMAPI
library(geomapdata) data(ETOPO5) PLOC=list(LON=c(137.008, 141.000),LAT=c(34.000, 36.992), x=c(137.008, 141.000), y=c(34.000, 36.992) ) PROJ = setPROJ(type=2, LAT0=mean(PLOC$y) , LON0=mean(PLOC$x) ) COLS = settopocol() JMAT = GEOTOPO(ETOPO5, PLOC, PROJ, COLS$calcol)