getGEOmap: Get Geomap

In GEOmap: Topographic and Geologic Mapping

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)

GEOmap documentation built on Sept. 1, 2023, 5:09 p.m.