hist() in raster package
Hi! I'm confused by the following
r <- brick("SDIM0069.tif")
nlayers(r)
[1] 3 Then I do
hist(r)
and get 3 histograms plotted, but I do:
a <- hist(r,plot=F)
and object a seems to hold the histogram of only one layer (actually, the 3rd one):
str(a)
List of 7 $ breaks : num [1:19] 0 5 10 15 20 25 30 35 40 45 ... $ counts : int [1:18] 337 1131 1280 1172 951 815 808 827 741 636 ... $ intensities: num [1:18] 0.00675 0.02264 0.02563 0.02346 0.01904 ... $ density : num [1:18] 0.00675 0.02264 0.02563 0.02346 0.01904 ... $ mids : num [1:18] 2.5 7.5 12.5 17.5 22.5 27.5 32.5 37.5 42.5 47.5 ... $ xname : chr "values" $ equidist : logi TRUE - attr(*, "class")= chr "histogram" Am I doing something wrong? Thanks Agus PS. Is package RemoteSensing being under active development?