Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I'm trying to retrieve the following pixel statistics:

  1. Proportion of zero value pixels:

    [nº of zero value pixels / (total nº of pixels - nº of NA pixels)]

  2. Proportion of pixels smaller than 3rd quartile (= 0.75):

    [nº of pixels < 3rd quartile / (total nº of pixels - nº of NA pixels)]

Here is a reproducible example:

f <- system.file("external/test.grd", package="raster")  
r <- raster(f) #r is the object  
share|improve this question
up vote 4 down vote accepted

One way of doing this is to work with the matrix:

#Raster to matrix
#Proportion of zero value pixels
#Proportion of values < 3rd Quantile     
value=as.numeric(summary(as.vector(r_matrix))[5])#position of the 3rd quantile value
share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.