[GRASS-user] Re: Command to summarize the frequency distribution table

Yasuo SHIMADA sgw00412 at nifty.com
Sun Dec 21 18:36:48 EST 2008


Dear Glynn: 

thank you good advice.

> Reclass maps only work for integers

I see, I understood the cause from your advice.

> Alternatively, if the intervals are equally spaced, you could generate
> an integer version by scaling, e.g.:
> 
> 	r.mapcalc 'newmap = int(oldmap * 10)'
> 

r.mapcalc has a wonderful function, I think.

Okey,it's a nice method. I'll try again. 

Thank you. 

Yasuo SHIMADA
--
> SGW00412 at nifty.com wrote:
> 
> > Thank you useful and kind advice.
> > Sorry, please give advice to me a little more.
> > 
> > I tried the re-classification by using r.reclass method in the raster map of the
> > probability distribution.
> 
> > I want to classify it at 0.1 intervals. So I wrote text files as follows:
> > -----------
> > 0.0 thru 0.1=0 A
> > 0.1 thru 0.2=2 B
> > 0.3 thru 0.4=3 C
> > 0.4 thru 0.5=4 D
> >  ....
> > 0.9 thru 1.0=10 J
> > end
> > -----------
> > 
> > then, I saved as probability_class.txt
> > 
> > r.reclass input=old_rast_map output=new_rast_map < probability_class.txt
> > 
> > then,  following massage printed out on display:
> > -----
> > 0.500000 rounded up to 0
> > 0.600000 rounded up to 0
> 
> > Why ? What's happened ?
> > I do not understand the cause not appropriately classified.
> 
> Reclass maps only work for integers (they are implemented by looking
> up the original cell values in an array, and you can't use
> floating-point values as array indices).
> 
> You can use r.recode to generate a new map using a piecewise-linear
> mapping. The category labels need to be set separately using
> r.category.
> 
> Alternatively, if the intervals are equally spaced, you could generate
> an integer version by scaling, e.g.:
> 
> 	r.mapcalc 'newmap = int(oldmap * 10)'
> 



More information about the grass-user mailing list