draw.brachiopod: Draw Brachiopod

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Draw Brachiopod at a particular location

Usage

1
draw.brachiopod(BK = matrix(), x = 0, y = 0, col = "black", fill = NULL, ...)

Arguments

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

Details

Used for strain illustrations

Value

side effects

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

get.brachiopod

Examples

 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)


  }

geophys documentation built on May 1, 2019, 9:26 p.m.