GEOmap-package

Topographic and Geologic Mapping

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.

Author(s)

Jonathan M. Lees<jonathan.lees.edu> Maintainer:Jonathan M. Leesjonathan.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.

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

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")
  • Maintainer: Jonathan M. Lees
  • License: GPL (>= 2)
  • Last published: 2024-07-09

Useful links