R/varsquig.R

In RSEIS: Seismic Time Series Analysis Tools

varsquig<-function(x,y, L=locator(2) , FLIP=FALSE, filcol="blue", tracecol="red" , var=0, xpd=TRUE)
{
  if(missing(FLIP) ) { FLIP=FALSE  }
  if(missing(L) ) { L=locator(2)  }
  if(missing(filcol) ) { filcol="red"  }
  if(missing(tracecol) ) { tracecol="red"  }
  
  if(missing(var) ) {  var=TRUE }
  ##  plot(x, y , type='n')

  if(is.numeric(var))
    {
      zee = var*(max(y, na.rm=TRUE)-mean(y, na.rm=TRUE))
    }
  else
    {
      zee = 0
    }

  if(FLIP==TRUE)
    {

      rx = RPMG::RESCALE(y, min(L$x), max(L$x), min(y), max(y))
      ry = RPMG::RESCALE(x, max(L$y), min(L$y), min(x), max(x))
      zer =RPMG::RESCALE(zee, min(L$x), max(L$x), min(y), max(y))

      yup = rx>zer
      g = rx
      g[!yup] = zer 
      g[1] = zer
      g[length(rx)] = zer
      if(var) polygon(g, ry, col=filcol, border=NA)

    }
  else
    {
      rx = RPMG::RESCALE(x, L$x[1], L$x[2], min(x), max(x))
      ry = RPMG::RESCALE(y, L$y[1], L$y[2], min(y), max(y))
      zer =RPMG::RESCALE(zee, L$y[1], L$y[2], min(y), max(y))

      yup = ry>zer
      g = ry
      g[!yup] = zer 
      g[1] = zer
      g[length(ry)] = zer
      if(var)  polygon(rx, g , col=filcol, border=NA, xpd=xpd)

    }


  lines(rx, ry, col=tracecol, xpd=xpd)

### lines(g, ry, col=col)


}