QGIS » QGIS Application

The HTML file is generated fine, but it's true that it's not as useful as it can be.

I am thinking of this solution: when this is needed, add in the description an output of type OutputFile, with a special name, so SEXTANTE instead of passing it to GRASS, uses it to pipe the result of the corresponding command into that file.

I think it can work

By the way the output of r.quantile is different in windows:

GRASS_INFO_MESSAGE(5216,1): Calculando histograma
GRASS_INFO_END(5216,1)

GRASS_INFO_MESSAGE(5216,2): Computing bins
GRASS_INFO_END(5216,2)

GRASS_INFO_MESSAGE(5216,3): Binning data
GRASS_INFO_END(5216,3)

GRASS_INFO_MESSAGE(5216,4): Sorting bins
GRASS_INFO_END(5216,4)

0:25.000000:1111.000000
1:50.000000:1402.000000
2:75.000000:1735.000000

GRASS_INFO_MESSAGE(5216,5): Calculando quantiles
GRASS_INFO_END(5216,5)

IS this also ok to be parsed by other GRASS commands?

Ok, I added the piping. All those messages are not piped (I gues they are wrriten to err...), so it should work fine

Hi Victor,
Thank you very much! Now is working!

But for some reason, in my case, the html file continues just registering the header info (<h2>r.quantile</h2>).
In the text file everything appears correctly.

Now I just have a doubt. r.quantile returns the upper limit value corresponding to each class, right?
But r.recode, uses r.quantile result as closed interval in the lower value and open interval in the upper value. In other words, assuming that r.quantile returns the following result:

2.000000:6.000000:1
6.000000:8.000000:2
8.000000:12.000000:3
12.000000:20.000000:4
20.000000:872.727295:5

r.recode will use these data as follows:

[2.000000:6.000000[ -> 1
[6.000000:8.000000[ -> 2
[8.000000:12.000000[ -> 3
[12.000000:20.000000[ -> 4
[20.000000:872.727295[ -> 5

But this is not correct. According to quantiles method, the reclassification should be as follows:

]2.000000:6.000000] -> 1
]6.000000:8.000000] -> 2
]8.000000:12.000000] -> 3
]12.000000:20.000000] -> 4
]20.000000:872.727295] -> 5

or am I wrong?

Thanks!

Hi,

As suspected, we can not directly use the output of r.quantile in r.recode. See this thread in the GRASS-user mailing list:

http://osgeo-org.1560.x6.nabble.com/Using-r-quantile-result-with-r-recode-td5046747.html

There are times when this method results, particularly when the boundary values ​​are not values ​​from the input data, but this not always happen, producing erroneous results. For this reason, I increased the priority level for this ticket to blocker.

Any idea about a possible solution?

Thanks!

Pedro Venâncio wrote:

Hi,

As suspected, we can not directly use the output of r.quantile in r.recode. See this thread in the GRASS-user mailing list:

http://osgeo-org.1560.x6.nabble.com/Using-r-quantile-result-with-r-recode-td5046747.html

There are times when this method results, particularly when the boundary values ​​are not values ​​from the input data, but this not always happen, producing erroneous results. For this reason, I increased the priority level for this ticket to blocker.

Any idea about a possible solution?

Thanks!

Hi Pedro,
the SEXTANTE tickets have been merged with the QGIS desktop ones, so I guess we have to lower the priority of this ticket as is not a regression. I understand this an annoying issue, but it is kind of an edge one. Do you agree?

Still empty result on master, see also #12764

well... it seems to work at least on Windows. In my case does not work under Linux, but for others is ok. Before closing I would like to gather more feedback.

On Windows there is the html output, but in the results viewer it is shown the html code instead of the formatted text.
It seems that there are Linux systems where there is the html output (instead of empty) but still it show the raw html in the result viewer.

Fixed in changeset a842d475126d3dd368f7af45278015420149a1ae.