# Loading the DEM
require(rgdal)
data(EbroDEM1000m)
dem <- EbroDEM1000m
# Plotting the DEM
require(sp)
spplot(dem, scales=list(draw=TRUE, y=list(rot=90)))
# Computing and plotting the hypsometric curve
hypsometric(dem)
Run the code above in your browser using DataLab