Thanks for the quick reply. I seem to have gotten it working, however
verrryyyyy sloowwwlllly on a 2.4 GHZ machine, 1 GB RAM. Any chance the rst
code has been optimized for 64 bit?
The below syntax took about 30 minutes to run....one tile of LIDAR
data...Andy, I think you and I will be working on the same dataset (entire
neuse) at one run...read the paper on RST and the various parameters, but
can't seem to find a groove that doesn't bring my workstation to it's knees.
Read in archives of people getting millions of points done in minutes!
GRASS 6.1.cvs (counties):~ > g.region vect=grid_0793 res=10
GRASS 6.1.cvs (counties):~ > g.region -pm
projection: 99 (Lambert Conformal Conic)
zone: 0
datum: nad83
ellipsoid: a=6378137 es=0.00669438
north: 740004.9853608
south: 729995.4277972
west: 2089999.43912146
east: 2100009.54270455
nsres: 3.04787138
ewres: 3.04803764
rows: 1001
cols: 1001
v.surf.rst input=grid_0793 layer=0 zcolumn=flt1 dmax=24.998895
dmin=4.999779 elev=elev_0973 zmult=1.0 tension=40. smooth=0.1 segmax=10
npmin=100
Authors: original version - H.Mitasova, L.Mitas, I. Kosinovsky, D.P. Gerdes
See manual pages for reference and publications T3
Percent complete: Reading lines from vector map ...
Reading nodes from vector map ...
Warning: ignoring 610 points -- too dense
(Why are we dropping 610 points? Could these all be in one place?)
The number of points from vector file is 308552 The number of points outside
of region 0 The number of points being used is 307942
Processing all selected output files
will require 4949732 bytes of disk space for temp files Percent complete:
A[100]=-0.100000
A[100]=-0.100000
A[100]=-0.100000
A[100]=-0.100000
A[100]=-0.100000
A[100]=-0.100000
history initiated
Thomas Colson
North Carolina State University
Department of Forestry and Environmental Resources
(919) 673 8023
tom_colson@ncsu.edu
Calendar:
www4.ncsu.edu/~tpcolson
-----Original Message-----
From: Andrew Danner [mailto:adanner@cs.duke.edu]
Sent: Wednesday, June 15, 2005 3:42 PM
To: tom_colson@ncsu.edu
Cc: grasslist@baylor.edu
Subject: Re: [GRASSLIST:7173] Gettting off the ground with v.surf.rst
Tom,
Some of your parameters seem a bit odd, mainly dmin, segmax, and npmin seem
way too high. What is your grid resolution from g.region -p?
Typically dmin is set to be 1/2 of the grid resolution so if you want a 20
meter resolution surface, try a dmin of 10. Decreasing dmin typically
increases the number of points used for the interpolation. Since your dmin
is high and your point count is low, I would look into adjusting dmin first.
Once you fix that though, you may find that the interpolation step takes
forever, mostly because segmax and npmin are set quite high as well. For
LIDAR data from NCfloodmaps or LDART I typically use npmin=100 and
segmax=20. The running time of the interpolation is roughly proportional to
(n/segmax)*npmin^3, so cutting both npmin and segmax in half makes the
program run about 4 times faster.
See if that helps.
-Andy
On Wed, 2005-06-15 at 14:52 -0400, Thomas Colson wrote:
this may appear to be a stupid question, but I'm having a hard time
getting the "desired" output from v.surf.rst, and I don't want to run
downstairs and bother Dr. Mitasova with such a menial question....
I have a x,y,z vector of LIDAR elevation data, and I'm trying to get a
raster elevation surface based on the z col. I've tried many different
variations of the output below, with the same errors "The number of
points being used is 3 " and so forth. I'm sure this has been done
many times by someone in the know, any chance of an example command
syntax so I can get this running?
Thanks.
v.surf.rst input=lidar3 layer=0 zcolumn=flt1 dmax=100000 dmin=10000
elev=surf7 zmult=1.0 tension=40. smooth=0.1 segmax=400 npmin=2000
Authors: original version - H.Mitasova, L.Mitas, I. Kosinovsky, D.P.
Gerdes See manual pages for reference and publications T3
Percent complete: Reading lines from vector map ...
Reading nodes from vector map ...
Warning: strip exists with insufficient data
Warning: ignoring 348352 points -- too dense
WARNING: 3 points given for interpolation (after thinning) is less
than given NPMIN=2000 Warning : there is less than 2050 points for
interpolation, no segmentation is necessary, to run the program
faster, set segmax=2050 (see
manual)
The number of points from vector file is 348355 The number of points
outside of region 0 The number of points being used is 3
Processing all selected output files
will require 1612 bytes of disk space for temp files Percent complete:
history initiated