World Map centered on Lat-Lon
World Map centered on Lat-Lon with Lambert-Azimuthal projection.
LLsmallcircMap(phicen, lamcen, worldmap, eqlat, eqlon, MAXR = 100, circol = rgb(1, 0.8, 0.8), mapcol = rgb(0, 0, 0), eqcol = rgb(0.8, 0.8, 1) , pch=25, ecex=1 )
phicen
: Center Latitudelamcen
: Center Longitudeworldmap
: GEOmap map structureeqlat
: Latitudes of points, vectoreqlon
: Longitude of points, vectorMAXR
: Maximum radius, degreescircol
: Color for small circlesmapcol
: Color for mapeqcol
: Color for points, single or 2-vectorpch
: Plotting character for pointsecex
: Plotting size for pointsUses a Lamber-Azimuthal projection of the whole globe out to the given radius. If a vector of 2 colors are provided for the eqcol parameter, and the pch is one of (21:25), then a 2-tone points is plotted with ecol[1] on the perimeter, and ecol[2] on the interior.
Graphical side effects
Jonathan M. Leesjonathan.lees@unc.edu
lamaz.eqarea
### earthquake in Noto,Japan data(coastmap) elat = 37.4949 elon = 137.2653 K = list(lon=c(183.3158,188.2173,253.5428,295.3037, 166.4531,110.5354,268.7554,98.9443,212.1384,236.6954), lat=c( 51.8823,-13.9085,34.94591,32.3713,68.0653, -66.2792,38.0557,18.8141,64.8736,44.5855 )) LLsmallcircMap(elat, elon, coastmap, K$lat, K$lon ) LLsmallcircMap(elat, elon, coastmap, K$lat, K$lon, MAXR=80, eqcol=c('blue', 'gold') , mapcol=grey(.8), pch=22, ecex=1.5 )
Useful links