Nothing
`aGETXprofile` <-
function(jx, jy, jz, LAB="A", myloc=NULL, PLOT=FALSE, asp=1)
{
############### get a cross section through a topographic DEM
############ aGETXprofile : xsec through topo
if(missing(myloc)) { myloc=NULL }
if(missing(LAB)) { LAB="A" }
if(missing(asp)) { asp=NULL }
if(is.null(myloc))
{
myloc = locator(2, type='o')
}
dx = sign(diff(myloc$x))*mean(diff(jx))
dy = sign(diff(myloc$y))*mean(diff(jy))
Lrunvent = lm ( myloc$y ~ myloc$x )
newx = seq(from=myloc$x[1], to=myloc$x[2], by=dx)
JY = Lrunvent$coefficients[1]+Lrunvent$coefficients[2]*newx
newy = seq(from=myloc$y[1], to=myloc$y[2], by=dy)
JX = (newy-Lrunvent$coefficients[1])/Lrunvent$coefficients[2]
allx = c(newx, JX)
ally = c(JY, newy)
ox = order(allx)
allx = allx[ox]
ally = ally[ox]
### points(jx, JY)
boxx = findInterval(allx, jx)
boxy = findInterval(ally, jy, all.inside = FALSE)
flag = boxy>0 & boxy<length(jy)
### points(jx[flag],JY[flag], col=p2)
pts = cbind(boxx[flag], boxy[flag])
LX = allx[flag]
LY = ally[flag]
RX = sqrt((LX-LX[1])^2+ (LY-LY[1])^2)
LZ = jz[pts]
pnt1 = sqrt((myloc$x[1]-LX[1])^2+ (myloc$y[1]-LY[1])^2)
pnt2 = sqrt((myloc$x[2]-LX[1])^2+ (myloc$y[2]-LY[1])^2)
px1 = findInterval(pnt1, RX)
px2 = findInterval(pnt2, RX)
if(PLOT==TRUE)
{
###screens(2)
cdev = dev.cur()
### get(getOption("device"))()
dev.new()
plot(RX, LZ, type='l', xlab="m", ylab="m", ylim=range(jz, na.rm=TRUE), asp=asp)
points( c(pnt1, pnt2), c(LZ[px1], LZ[px2] ), pch=c(6,8), col=c(2,4) )
text(c(pnt1, pnt2),c(LZ[px1], LZ[px2] ) , labels=c(LAB, paste(sep="", LAB, "'")), col=c(2,4), pos=3 )
dev.set(cdev)
}
invisible(list(RX=RX, RZ=LZ, LOC=myloc, LAB=LAB ))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.