[R-sig-Geo] hist() in raster package

Agustin Lobo alobolistas at gmail.com
Tue Feb 9 14:32:28 CET 2010


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?



More information about the R-sig-Geo mailing list