profilerraster
How can I add/display scale to two marginal plot?
RAD1998.all <- stack(list.ras)
list.files <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = F, pattern = ".asc.bil"))
#list.files
outlist <- substr(list.files, 31, 43)
#outlist
png(paste("1998",outlist[[1]],".png",sep=""))
levelplot(RAD1998.all[[1]])
dev.off()
There is a code example for adding a spatial scalebar in the FAQ section of rasterVis:
https://oscarperpinan.github.io/rastervis/FAQ.html#scalebar
f <- system.file("external/test.grd", package="raster")
r <- raster(f)
levelplot(r) +
layer({
xs <- seq(181000, 181400, by=100)
grid.rect(x=xs, y=330500,
width=100, height=30,
gp=gpar(fill=rep(c('transparent', 'black'), 2)),
default.units='native')
grid.text(x= xs - 50, y=330560, seq(0, 400, by=100),
gp=gpar(cex=0.5), rot=30,
default.units='native')
})
Use lines3d from the rgl package but make sure to plot the raster with adjust=FALSE or the axes won't line up.
Where r is a raster and xyz is a data frame with x, y, z (height) locations, I do:
plot3D(r, col = rainbow, adjust=FALSE)
lines3d(xyz)
and get
You can get the heights of points from the raster using extract if you don't have them already, and you might want to add a little offset to the heights to keep them just clear of the terrain, so they don't get clipped by any interpolation. You can see my line gets a bit clipped because it's a rough transect with only ten points.
You can use any rgl package functions to augment that plot, such as axis3d("x") to show the x-axis.
Best Answer
I have committed changes in the repository that implement the axis.
Use this code to install the latest development version:
Then, use the new argument
axis.marginto enable the axis.