Hi!
I'm confused by the following
> r <- brick("SDIM0069.tif")
> nlayers(r)
[1] 3Then 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? _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
