[GRASSLIST:9280] Re: How to extract elevation of sites?

[Jachym Cepicky]

hi,

for point in `v.out.ascii my_points`; do 
    eval $( echo "$point" | sed -e "s/^/x=/" | sed -e "s/|/; y=/" | sed -e "s/|/; cat=/g" ); 
    r.what in=dem east_north=$x,$y | sed -e "s/||/|$cat|/" >> my_points.txt; 
done

and then run v.in.ascii -zn in=my_poinst.txt out=my_points_z cat=3 z=4

hope it helps

jachym

P.S. If you are using grass 6.1, v.drape could help

On Fri, Dec 02, 2005 at 11:50:49PM +0100, odoardo.zecca at libero.it wrote:
> Hi all.
> 
> I need to extract the elevation of 85 sites from a DEM. How could I do that
> giving a list of coordinates (from a file or standard input) instead of
> interactively use d.what.rast?
> 
> Thank you
> 
> --
> Odoardo Zecca
> via Vallescuropasso, 50		tel.: 328 2883121
> 27043 Broni (PV)		email: odoardo.zecca at libero.it 
> 

-- 
Jachym Cepicky
e-mail: jachym.cepicky at centrum.cz
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/gnupg_public_key/
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
LDF MZLU v Brně
Zemědělská 3
613 00 Brno
e-mail: xcepicky at node.mendelu.cz
URL:    http://mapserver.mendelu.cz
Tel.:   +420 545 134 514