getGEOmap {GEOmap}

Get Geomap

Description

Get Geomap from ascii files

Usage

getGEOmap(fn)

Arguments

fn

root name

Details

Files are stored as a pair: rootname.strks and rootname.pnts

Value

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

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

plotGEOmapXY, boundGEOmap

Examples

## 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)