Description Usage Arguments Details Value Author(s) See Also Examples
Draw Brachiopod at a particular location
1 | draw.brachiopod(BK = matrix(), x = 0, y = 0, col = "black", fill = NULL, ...)
|
BK |
brachiopod matrix |
x |
x-coordinate location on plot |
y |
y-coordinate location on plot |
col |
col for outline |
fill |
color for fill |
... |
other par parameters for plotting |
Used for strain illustrations
side effects
Jonathan M. Lees<jonathan.lees@unc.edu>
get.brachiopod
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 | gb = get.brachiopod()
draw.brachiopod(gb)
BB = get.brachiopod()
phi0=0
phi = rep(phi0, 6)
shears = seq(from=0, to=1.5, length=6)
shearmats = list()
for(i in 1:length(shears))
{
shr = shears[i]
shearmats[[i]] = matrix(c(1, shr, 0, 1), ncol=2)
}
phi0=0
phi = rep(phi0, 6)
phi=runif(6, 0, 90)
shears = rep(0, 6)
for(i in 1:length(phi))
{
shr = shears[i]
shearmats[[i]] = matrix(c(1, shr, 0, 1), ncol=2)
}
x = rep(c(25,75), times=3)
y = as.vector( t(matrix( rep(c( 75, 50,25 ), times=2) , ncol=2)) )
siz = rep(30, 6)
plot(c(0, 100), c(0, 100), asp=1, type='n', ann=FALSE, axes=FALSE)
for(i in 1:length(x))
{
cosp=cos(phi[i]*pi/180)
sinp=sin(phi[i]*pi/180)
####### make a rotation matrix
#### rmat=matrix(c(cosp,-sinp, x, sinp, cosp, y), ncol=2)
####
rmat1=matrix(c(cosp,-sinp, sinp, cosp), ncol=2)
total = rmat1
rmat= rbind( total, c(x[i],y[i]) )
n = length(BB[,1])
RACK1 = cbind(siz[i]*BB[,1:2], rep(1,n))
draw.brachiopod(RACK1)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.