Geological legend from GEOmap Structure
Create and add Geological legend from GEOmap Structure
geoLEGEND(names, shades, zx, zy, nx, ny, side=1, cex=0.5)
names
: namesof unitsshades
: colorsof unitszx
: width of box, mmzy
: height of box, mmnx
: number of boxes in x-directionny
: number of boxes in y-directionside
: Side of the plot for the legend (1,2,3,4)cex
: Character expansion for text in legendAdds geological legend based on information provided. Legend is placed in margin.
Graphical Side Effects
Jonathan M. Leesjonathan.lees@unc.edu
If plot is resized, should re-run this as the units depend on the screen size information and the transformation of user coordinates.
## Not run: library(RPMG) library(RSEIS) library(GEOmap) library(geomapdata) data(cosogeol) data(cosomap) data(faults) data(hiways) data(owens) proj = cosomap$PROJ XMCOL = setXMCOL() newcol = XMCOL[cosogeol$STROKES$col+1] cosocolnums = cosogeol$STROKES$col cosogeol$STROKES$col = newcol ss = strsplit(cosogeol$STROKES$nam, split="_") geo = unlist(sapply(ss , "[[", 1)) UGEO = unique(geo) mgeo = match( geo, UGEO ) gcol = paste(sep=".", geo, cosogeol$STROKES$col) ucol = unique(gcol) N = length(ucol) spucol = strsplit(ucol,split="\.") names = unlist(sapply(spucol , "[[", 1)) shades = unlist(sapply(spucol , "[[", 2)) ORDN = order(names) ### example: par(mai=c(0.5, 1.5, 0.5, 0.5) ) plotGEOmapXY(cosomap, PROJ=proj, add=FALSE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosogeol, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) geoLEGEND(names[ORDN], shades[ORDN], .28, .14, 4, 16, side=2) #### par(mai=c(0.5, 0.5, 1.0, 0.5) ) plotGEOmapXY(cosomap, PROJ=proj, add=FALSE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosogeol, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) geoLEGEND(names[ORDN], shades[ORDN], .28, .14, 16, 6, side=3) #### par(mai=c(0.5, 0.5, 0.5, 1) ) plotGEOmapXY(cosomap, PROJ=proj, add=FALSE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosogeol, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) geoLEGEND(names[ORDN], shades[ORDN], .28, .14, 3, 16, side=4) #### par(mai=c(1.5, 0.5, 0.5, 0.5) ) plotGEOmapXY(cosomap, PROJ=proj, add=FALSE, ann=FALSE, axes=FALSE) plotGEOmapXY(cosogeol, PROJ=proj, add=TRUE, ann=FALSE, axes=FALSE) geoLEGEND(names[ORDN], shades[ORDN], .28, .14, 16, 3, side=1) ## End(Not run)
Useful links