R/linesGEOmapXY.R

In GEOmap: Topographic and Geologic Mapping

linesGEOmapXY<-function(lat=0, lon=0, PROJ=NULL, PMAT=NULL, ... )
  {
    if(missing(PMAT)) { PMAT=NULL }

    pxy = GLOB.XY(lat, lon, PROJ)


    if(!is.null(PMAT))
      {
        tem1 = trans3d(pxy$x, pxy$y, rep(0, length(pxy$y)) , PMAT)
      }
    else
      {
        tem1 =pxy

      }


    lines(tem1, ...)



  }