GEOmap-package {GEOmap} | R Documentation |
Topographic and Geologic Mapping
Description
Topographic and Geologic Mapping
Details
Set of routines for making Map Projections (forward and inverse), Topographic Maps, Perspective plots, geological databases, interactive plotting and selection of focus regions.
Note
- High level plotting:
-
BASICTOPOMAP DOTOPOMAPI geoLEGEND GEOsymbols locworld plotGEOmap plotGEOmapXY linesGEOmapXY rectGEOmapXY textGEOmapXY pointsGEOmapXY insideGEOmapXY plotUTM plotworldmap XSECDEM
- PLOTTING:
-
circle addLLXY addTIX antipolygon zebra demcmap setXMCOL shade.col
- Geological Map Symbols:
-
bcars faultdip faultperp horseshoe normalfault OverTurned perpen teeth thrust SynAnticline SSfault
- Data manipulation:
-
getGEOmap boundGEOmap SELGEOmap geoarea GEOTOPO getGEOperim GETXprofile Lintersect LOCPOLIMAP pline selectPOLImap setplotmat SETPOLIMAP settopocol subsetTOPO
- Misc:
-
getgreatarc ccw difflon DUMPLOC getsplineG inpoly inside PointsAlong polyintern
- Projections:
-
setPROJ projtype GLOB.XY XY.GLOB MAPconstants GCLCFR lambert.cc.ll lambert.cc.xy lambert.ea.ll lambert.ea.xy lcgc merc.sphr.ll merc.sphr.xy utmbox utm.elps.ll utm.elps.xy utm.sphr.ll utm.sphr.xy stereo.sphr.ll stereo.sphr.xy equid.cyl.ll equid.cyl.xy
Author(s)
Jonathan M. Lees<jonathan.lees.edu> Maintainer:Jonathan M. Lees<jonathan.lees@unc.edu>
References
Snyder, John P., Map Projections- a working manual, USGS, Professional Paper, 1987.
Lees, J. M., Geotouch: Software for Three and Four Dimensional GIS in the Earth Sciences, Computers & Geosciences, 26, 7, 751-761, 2000.
See Also
RSEIS
Examples
################ projections
proj = setPROJ(type = 2, LAT0 =23, LON0 = 35)
### get lat-lon
LL = XY.GLOB(200, 300, proj)
## find x-y again, should be the same
XY = GLOB.XY(LL$lat, LL$lon, proj)
XY
################
library(geomapdata)
data(worldmap)
KAMlat = c(48.5, 65)
KAMlon = c(150, 171)
PLOC=list(LON=KAMlon,LAT=KAMlat)
PLON = seq(from=KAMlon[1], to=KAMlon[2], by=2)
PLAT = seq(from=KAMlat[1], to=KAMlat[2], by=2)
proj = setPROJ(2, LON0=mean(KAMlon), LAT0=mean(KAMlat))
xy = GLOB.XY(KAMlat, KAMlon , proj)
kbox=list(x=range(xy$x, na.rm=TRUE), y=range(xy$y, na.rm=TRUE))
plot(kbox$x,kbox$y, type='n', axes=FALSE, xlab="", ylab="", asp=1)
plotGEOmapXY(worldmap, LIM=c(KAMlon[1], KAMlat[1], KAMlon[2],
KAMlat[2]), add=TRUE, PROJ=proj, axes=FALSE, xlab="", ylab="" )
sqrTICXY(kbox , proj, side=c(1,2,3,4), LLgrid=TRUE, col=grey(.7) )
title("Crude Map of Kamchatka")