Click on two points of a plot (map) to obtain an object of class Extent
('bounding box')
drawExtent(show=TRUE, col="red")
logical. If TRUE
, the extent will be drawn on the map
sets the color of the lines of the extent
Extent
# NOT RUN {
r1 <- raster(nrow=10, ncol=10)
r1[] <- runif(ncell(r1))
plot(r1)
# after running the following line, click on the map twice
e <- drawExtent()
# after running the following line, click on the map twice
mean(values(crop(r1, drawExtent())))
# }
Run the code above in your browser using DataLab