[GRASS-user] r.reclass using pipe from r.stats?

[Glynn Clements]

Milton Cezar Ribeiro wrote:

> I have a raster map with all categories (~50)
> and I have about 500 spatial subsets of the same
> original map. As is expected, each of spatial subset
> map have less categories than the original map.
> 
> I can check the list of categories for each subset
> map using r.stats. But I need to generate a new
> set of maps with the spatial range of original map,
> but keeping only those classes of my original map.
> 
> Just supose that "r.stats mysubset1" result something
> like:
> 
> 2
> 4
> 9
> 11
> 43
> 
> So, I would like that my new map with my full spatial
> range keep only these 5 categories, and all other
> categories come to be "no data".
> 
> I know if I generate a rule file with
> 2=2
> 4=9
> 11=11
> 43=43
> I can get the desired results, but as the number
> of spatial subsets is very large I would like turn it
> more operational.
> 
> Is there a suggestion of how can I do this using pipe like
> r.stats mysubset1 | r.reclass mysubset1 mysubset1_fullrange rules= -??

	r.stats mysubset1 | awk '{print $1,"=",$1} ; END {print "* = NULL"}' | r.reclass ...

-- 
Glynn Clements <glynn at gclements.plus.com>