GEOmap (version 2.1)

getGEOmap: Get Geomap

Description

Get Geomap from ascii files

Usage

getGEOmap(fn)

Arguments

fn

root name

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

Details

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

See Also

plotGEOmapXY, boundGEOmap

Examples

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)