GEOmap (version 2.5-11)

getGEOmap: Get Geomap

Description

Get Geomap from ascii files

Usage

getGEOmap(fn)

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

Arguments

fn

root name

Author

Jonathan M. Lees<[e-mail]>

Details

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

See Also

plotGEOmapXY, boundGEOmap

Examples

if (FALSE) {
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)
  
}