How can I reverse the y-axis in a raster plot?

user_80626 0
Viewed 404

I would like to plot a raster with the y-axis increasing from top to bottom, and the raster correspondingly flipped. Is there a way to do this? Specifying ylim in decreasing order gives an error.

For example, I would like this plot with the y-axis running from -100 at the top to 100 at the bottom, with the raster also reversed to put green on top and orange on the bottom:

library(raster)
r <- raster(nrows=10, ncols=10)
r <- setValues(r, 1:ncell(r))
plot(r)



# Fails
plot(r, ylim=c(100, -100))
#> Error in .plotraster2(x, col = col, maxpixels = maxpixels, add = add, : invalid ylim

Created on 2020-10-30 by the reprex package (v0.3.0)

r-raster r
1 Answers

The approach shown below appears to work

Example data

library(raster)
r <- raster(nrows=10, ncols=10, vals=1:100)

Flip the raster vertically, and add the horizontal axis

plot(flip(r, "y"), axes=F)
axis(1)

Get the labels for the vertical axis, and plot them in reverse order

ylabs <- axis(2, labels=FALSE, tick=FALSE)
axis(2, at=ylabs, labels=rev(ylabs))
user_6352450
Related