data(GH)
XLIM = c(1226, 1322 )
e = GH$JSTR[[1]][XLIM[1]:XLIM[2]]
n = GH$JSTR[[2]][XLIM[1]:XLIM[2]]
xx = range(e, na.rm =TRUE)
yy = range(n, na.rm =TRUE)
sx = range(c(xx, yy))
x = RPMG::RESCALE(e, 0, 1, sx[1], sx[2])
y = RPMG::RESCALE(n, 0, 1, sx[1], sx[2])
plot(range(x), range(y) , type='n')
lines(x, y, col=grey(0.8) )
DO.PMOT.ARR(x, y)
Run the code above in your browser using DataLab