[GRASS-user] Re: Lidar Data automatic processing help

[stephen sefick]

#!/bin/sh
FILES=$(ls /Volumes/biggun/XYZ_DATA_TILES | grep dsm)
for k in $FILES ; do
g.region vect=$k
v.surf.idw input=$k output=dem$k npoints=12 layer=0 column="Point height"
done


will this do what I want?

On Wed, Sep 9, 2009 at 11:51 AM, stephen sefick<ssefick at gmail.com> wrote:
> 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
>



-- 
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