Get Geomap
Get Geomap from ascii files
getGEOmap(fn)
fn
: root nameFiles are stored as a pair: rootname.strks and rootname.pnts
STROKES: List of stroke information:
nam: name of stroke
num: number of points
index: index where points start
col: color
style: plotting style: 1=point, 2=line,3=polygon
code: character, geological code
LAT1: bounding box lower left Lat
LAT2: bounding box upper right Lat
LON1: bounding box lower left Lon
LON2: bounding box upper right Lon
POINTS: List of point LL coordinates, list(lat, lon)
PROJ: optional projection parameters
Jonathan M. Leesjonathan.lees@unc.edu
plotGEOmapXY, boundGEOmap
## Not run: library(geomapdata) data(cosomap) data(faults) data(hiways) data(owens) cosogeol = getGEOmap("/home/lees/XMdemo/GEOTHERM/cosogeol") cosogeol = boundGEOmap(cosogeol) proj = cosomap$PROJ plotGEOmapXY(cosomap, PROJ=proj, add=FALSE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosogeol, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosomap, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) plotGEOmapXY(faults, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) ## End(Not run)
Useful links