[GRASS-dev] v.cellstats script

[Martin Wegmann]

On Friday 17 November 2006 16:55, Markus Neteler wrote:
> Martin Wegmann wrote on 11/17/2006 02:34 PM:
> > Hello,
> >
> > Carlos Grohmann and I created a script as substitute for v.cellstats. It
> > would be great if somebody can look through it to give feed-back how to
> > improve it before we post it to the WIKI.
> >
> > regards,
>
> Hi Martin, Carlos,
>
> two quick observations:
>
> This is a bash-ism AFAIK (still needs to be eliminated in some scripts):
>    export LC_NUMERIC=C
>
>
> Before doing that the current region should be saved. In an older mail
> from Glynn it was described how to solve it in an elegant way:
>    #set to current input map region
>    g.region rast=$GIS_OPT_raster
>
> Of course the current region should be restored at the end of the script
> and in case of "trap" (which isn't there, see cleanup() in other scripts).

thanks Markus, I'll will fix it -- in the meantime I found a bug with the 
g.region settings for the coarser raster - only the resolution setting works 
otherwise variance is 0. 

Martin

> Markus
>
> _______________________________________________
> grass-dev mailing list
> grass-dev at grass.itc.it
> http://grass.itc.it/mailman/listinfo/grass-dev