[GRASSLIST:5960] Re: Problem using r.in.ascii

Paul Kelly paul-grass at stjohnspoint.co.uk
Fri Apr 4 03:30:52 EST 2003


Hello

On Thu, 3 Apr 2003, Victor Wren wrote:

> I've got a rather strange problem I can't figure out, using GRASS 5.0.1
>
> I have a large number of sites (on the order of 176,000) that I'm trying to
> use as elevations to create a surface (they were derived from road surveys,
> control points, signposts, and miscellaneous building perimeters, scavenged
> from a DXF file by a script I wrote).  The logical way to do this would be to
> bring them in as sites, and do s.surf.rst, but the surface I get is too
> convoluted (and playing with tension values from 10 to 100 doesn't seem to
> make any noticeable difference).

Did you try using s.surf.idw? I've recently added an improved version to
CVS, which should be able to deal with many hundreds of thousands of sites
quite quickly (provided they aren't located too sparsely) and I would
appreciate some testing for it.

If you don't like using CVS, you can get the changed main.c file from
http://freegis.org/cgi-bin/viewcvs.cgi/*checkout*/grass/src/sites/s.surf.idw/cmd/main.c
Just rename the main.c in src/sites/s.surf.idw/cmd/ and copy this file in,
then recompile with gmake5 -i.
The updated man page is at
http://freegis.org/cgi-bin/viewcvs.cgi/*checkout*/grass/html/html/s.surf.idw.html

Paul




More information about the grass-user mailing list