[GRASSLIST:7060] Re: Limitation in r.what

[Daniel Victoria]

Not sure if it's the best way but maybe a bash script will do the trick

for i in <list of raster maps>
do r.what input=$i  east_north=$myeast, $mynorth >> temp
done

This will (hopefully) append the value to the temp file. Just not
entirely sure if the redirect (>>) should be inside the loop or after
the done statement. Can't test it right now

Cheers
Daniel


On 6/7/05, Giacomo Bertoldi <giacomo.bertoldi at ing.unitn.it> wrote:
> I would like to query more than 150 raster layers with r.what, to
> extract the time evolution of a variable in a certan point of a map
> (each layer is for a different time step) but i do receive the error:
> 
> r.what: can only do up to 150 cell files, sorry.
> 
> The call is:
> 
> r.what input=$RASTERMAPS east_north=$myeast, $mynorth
> 
> where $RASTERMAPS is the list of my files.
> 
> How can I overcome this problem?
> 
> Thanks, Giacomo
> 
> ************************************************
>    Giacomo Bertoldi, Ph.D.
>    Dept. of Civil and Environmental Engineering
>    University of Trento - ITALY
>    E-MAIL:    giacomo.bertoldi at ing.unitn.it
> ************************************************
> 
>