Re: [FreeGeoCZ] rast in Grass (r.mask a r.in.gdal)

[hanoj <ehanoj AT gmail.com>]

> ahoj,
*** Cau!

> On Thu, Jun 18, 2009 at 12:46:21PM +0200, hanoj wrote:
> > 2) r.mask
> > Tento modul umi, nasadit masku na rastr podle indexovych hodnot. Ja
> > bych ale potreboval masku stanovit podle RGB hodnot, protoze moje
> > hledana RGB (0,255,0) ma ruzne indexy napric rastry. Je na to nejaky
> > figl?
> > Pozadi: jedna se o stovky rastru puvodne vektorovych (rozliseni 2m/px
> > na uzemi cele republiky), ktere ve for-cyklu postupne vektorizuji.
> >
> > Jeste me ted napada obejit r.mask pomoci r.patch a cyklus provest jen
> > nad jednim rastrem za pomoci g.region. Vytvareni tak velkeho rastru
> > ale zkousel az naposled.
> >
>
> to by melo jit mapcalcem
>
> r.mapcalc MASK=r#vrstva
>
> vezme cervenou
*** Vezme cervenou, ale jako kanal ze vsech barev (tzn cervenou i z
sede skaly barev atd.). Ja hledam jen konkretni RGB hodnotu.

Nakonec jsem vyuzil vlastnosti cele sady, ze zelena RGB (0,255,0) je
vzdycky posledni ve vyctu  indexu barev. Ne moc hezke, ale pro tento
pripad pouzitelne.

# find green color (latest row is green)
mapmask=`r.stats -c -p -l input=$map | wc -l`

# first num is 0 than (( result - 1 ))
(( mapmask-- ))

#mast set to green color with data
r.mask -o --quiet input=$map maskcats=$mapmask


PS: pokus o vektorizaci http://googluj.cz/1432 pro OSM

diky

Dan B.