[GRASS-user] v.to.rast shell script

[John C. Tull]

On Sep 19, 2008, at 7:07 PM, John C. Tull wrote:

> Hamish:
>
> On Sep 19, 2008, at 6:45 PM, Hamish wrote:
>
>> John  Tull wrote:
>>> I was wondering if anyone could suggest a simple shell script to
>>> automate the process of running v.to.rast on all layers in
>>> a single vector in GRASS?
>>
>> sure, just a bash for loop + g.mlist will do it:
>>
>> for MAP in `g.mlist vect` ; do
>>  echo v.to.rast in="$MAP" out="`echo $MAP | tr '[:lower:]'  
>> '[:upper:]'`"
>> done
>>
>> in that example the other v.to.rast options would have to be  
>> constant of
>> course.
>>
>> To avoid mistakes I like to capitalize a map name if it is going to  
>> be
>> both a vector and raster map. Probably there is some simple sed  
>> command
>> to just do the first character. otherwise just use "out=$MAP".
>
>
> Thank you for your quick response. This, with the addition of  
> 'type=vect' will give me a list of vectors, but I was wanting a list  
> of layers within a single vector. I am guessing there might be a way  
> to cajole v.info into doing this, but am not sure how at this point.
>
> Thanks,
> John
>


Ok, applying my little grey cells (as Poirot would say) to the problem  
at hand, I believe I have resolved this as follows:

for LAYER in `v.category -g option=report in=vect | awk '{print $1}' |  
uniq` ; do
    v.to.rast in=vect layer="$LAYER" out=vect_$LAYER use=val
done

Hamish: thanks for the basis from which to suss the rest.

John