R: Tape style Lune Plot

TapePlot {RFOC} | R Documentation

## Tape style Lune Plot

### Description

Tape style Lune Plot using Hammer projection

### Usage

    TapePlot(TapeList = list(), add = FALSE, ann = TRUE,
    pcol = c(grey(0), grey(0.85), grey(0.95)))

### Arguments

`TapeList` | List of strokes from TapeBase
`add` | logical, TRUE=add to existing plot
`ann` | logical, TRUE=annotape
`pcol` | 3-vector of colors: inner lines, upper polygon, lower polygon

### Details

Plot an Tape net from the TapeBase function.

### Value

Side effects

### Author(s)

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

### References

Tape, W., and C. Tape (2012), A geometric comparison of source-type plots for moment tensors, Geophys. J. Int., 190, 499-510. https://doi.org/10.1111/j.1365-246X.2012.05490.x

### See Also

TapeBase, HAMMERprojXY

### Examples

    T1 = TapeBase()
    TapePlot(T1)

     data(widdenMoments)
    WM = widdenMoments

             par(mfrow=c(1,1), mai=c(0,0,0,0))
             T1 = TapeBase()
             TapePlot(T1)

             for(i in 1:length(WM$Mxx))
             {
               M1=matrix( c(WM$Mxx[i],WM$Mxy[i],WM$Mxz[i],
          WM$Mxy[i],WM$Myy[i],WM$Myz[i],
           WM$Mxz[i],WM$Myz[i],WM$Mzz[i]), ncol=3, byrow=TRUE)

               E1 = eigen(M1)
               h = SourceType( sort(E1$values, decreasing=TRUE) )
               h$dip = 90-h$phi
               ##  cat(paste(h$dip, h$lam, sep=" "), sep="\n")
               h1 = HAMMERprojXY(h$dip*pi/180, h$lam*pi/180)

               points(h1$x, h1$y,  pch=21, bg="orange" )

             }

[Package RFOC version 3.4-3 Index]