[GRASS-user] Lidar Data automatic processing help

[stephen sefick]

This is a little script that I have written to import a bunch of tiled
x y z data into GRASS:

#!/bin/sh
FILES=$(ls /Volumes/biggun/XYZ_DATA_TILES | grep dsm)
for k in $FILES ; do
v.in.ascii input=$k output=$k z=3 fs=" "
done

This is how I was doing it before:

#lidar processing
v.in.ascii input=/Volumes/biggun/XYZ_DATA_TILES/dsm16.xyz
output=points_16 fs=" " z=3

#set region to vect
g.region vect=points_16

#interpolate to grid
v.surf.idw input=points_16 at lidar output=dem16 npoints=12 layer=0
column="Point height"

#set color map
r.colors map=dem16 at lidar raster=dem01 at lidar

I would also like to to have the last three commands in the shell
script to pipe the output from 1 to 2 to 3 to 4.  I am just learning
how to write shell scripts and don't know quite how to do this.
thanks for all of your help,


-- 
Stephen Sefick

Let's not spend our time and resources thinking about things that are
so little or so large that all they really do for us is puff us up and
make us feel like gods.  We are mammals, and have not exhausted the
annoying little problems of being mammals.

								-K. Mullis