[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