PlotPTsmooth | R Documentation |
Project PT axes on the sphere and smooth the image. This function requires function kde2d, from the MASS library.
PlotPTsmooth(paz, pdip, x = 0, y = 0, siz = 1, bcol = "white", border ="black",
IMAGE = TRUE, CONT = TRUE, cont.col = "black",
pal = terrain.colors(100), LABS = FALSE, add = FALSE, NCP=50, NIP=200)
paz |
vector of Axis azimuths, degrees |
pdip |
vector of dip angles, degrees |
x |
x-location of plot center in user coordinates |
y |
y-location of plot center in user coordinates |
siz |
siz of plot in user coordinates |
bcol |
color |
border |
border color |
IMAGE |
logical, TRUE=create an image plot |
CONT |
logical, TRUE=add contour lines |
cont.col |
color of contour lines |
pal |
pallete for image plot |
LABS |
text Label for image |
add |
logical, TRUE=add to plot |
NCP |
integer, Number of points to use for calculating smoothed contours, default=50 |
NIP |
integer, Number of points to use for calculating smoothed image, default=200 |
Program requires MASS libary for 2D smoothing routine kde2d.
For calculating contours the kde2d program creates a smoothed 2D image using NCP points per side. For the images, NIP points are used. To reduce the size of plots, or, if the subplots are very small, reduce NIP to a smaller value for faster plotting.
Graphical Side Effect
Points that fall on the opposite hemisphere are reflected through the origin.
Jonathan M. Lees<jonathan.lees@unc.edu>
kde2d
plot(c(-1,1), c(-1,1), asp=1, type='n')
paz = rnorm(100, mean=297, sd=10)
pdip = rnorm(100, mean=52, sd=8)
PlotPTsmooth(paz, pdip, x=0.5, y=.5, siz=.3, border=NA, bcol='white' ,
LABS=FALSE, add=FALSE, IMAGE=TRUE, CONT=FALSE)
taz = rnorm(100, mean=138, sd=10)
tdip = rnorm(100, mean=12, sd=8)
PlotPTsmooth(taz, tdip, x=-.5, y=.4, siz=.3, border=NA, bcol='white' ,
LABS=FALSE, add=FALSE, IMAGE=TRUE, CONT=TRUE)
########### put them together
plot(c(-1,1), c(-1,1), asp=1, type='n')
PlotPTsmooth(paz, pdip, x=0, y=, siz=1, border=NA, bcol='white' ,
LABS=FALSE, add=FALSE, IMAGE=TRUE, CONT=FALSE)
PlotPTsmooth(taz, tdip, x=0, y=, siz=1, border=NA, bcol='white' ,
LABS=FALSE, add=TRUE, IMAGE=FALSE, CONT=TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.