[GRASS-dev] Re: [GRASS GIS] #359: r.quantile computes incorrect values with default number of bins

GRASS GIS trac at osgeo.org
Sat Feb 20 19:45:47 EST 2010


#359: r.quantile computes incorrect values with default number of bins
-------------------------+--------------------------------------------------
  Reporter:  dylan       |       Owner:  grass-dev at lists.osgeo.org
      Type:  defect      |      Status:  closed                   
  Priority:  major       |   Milestone:  6.5.0                    
 Component:  default     |     Version:  svn-develbranch6         
Resolution:  worksforme  |    Keywords:  r.quantile               
  Platform:  Linux       |         Cpu:  x86-32                   
-------------------------+--------------------------------------------------
Changes (by dylan):

  * milestone:  6.4.0 => 6.5.0

Comment:

 Some updates: Looks like things are working now... with a similar dataset
 I get the following:

 in GRASS:
 {{{
 r.quantile in=rtk_hybrid_beam_0150 quantiles=5
 Computing histogram
  100%
 Computing bins
 Binning data
  100%
 Sorting bins
  100%
 Computing quantiles
 0:20.000000:7096.507812
 1:40.000000:7291.034180
 2:60.000000:7412.145020
 3:80.000000:7514.283691
 }}}

 in R:
 {{{
 library(spgrass6)
 x <- readRAST6('rtk_hybrid_beam_0150')

 quantile(x at data, prob=c(0.2,0.4,0.6,0.8), na.rm=TRUE)
      20%      40%      60%      80%
 7096.507 7291.034 7412.144 7514.284
 }}}


 Seems close enough for me. The test file is attached.

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/359#comment:6>
GRASS GIS <http://grass.osgeo.org>


More information about the grass-dev mailing list