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

[Jachym Cepicky <jachym.cepicky AT gmail.com>]

ahoj,

On Fri, Jun 19, 2009 at 08:28:37AM +0200, hanoj wrote:
> > 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.

no, vzdyt jo, ne? 

r.mapcalc MASK="if(r#vrstva == 255 && g#vrstva < 34 && g#vrstva > 12 &&
b#vrstva > 0, 1, null())"

pripadne by to slo udelat reclassem (asi podobny postup jako u r.stats)

j


> 
> 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.

> _______________________________________________
> FreeGeoCZ mailing list
> FreeGeoCZ AT fsv.cvut.cz
> http://mailman.fsv.cvut.cz/mailman/listinfo/freegeocz


-- 
Jachym Cepicky
e-mail: jachym.cepicky gmail com
URL: http://les-ejk.cz
GPG: http://www.les-ejk.cz/pgp/JachymCepicky.pgp
Key fingerprint: 0C6D 0EAE 76BD 506C F299  ED8A C8AB 74B8 08D4 E08F

Attachment: signature.asc
Description: Digital signature