<br><font size=2 face="sans-serif">Hi Folks,</font>
<br>
<br><font size=3>I aggregated all of the bare earth lidar points for the
state of North Carolina into a single file &nbsp;Imported all 8.5 billion
points for NC into one point layer with no topology. &nbsp;I was sort of
curious to see if I could generate a seamless 20 ft elevation grid for
the State of North Carolina from a single data layerusing bspline and&nbsp;
tried the following command:<br>
<br>
GRASS 6.5.svn (ncstpft_nad83):~ &gt; v.surf.bspline input=all_nc_be_pts2
raster=all_nc_be_20ft_bspline sie=40 sin=40 layer=0<br>
WARNING: Coor files of vector map &lt;all_nc_be_pts2@statewide&gt; is larger
than<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; it should be (158913789952
bytes excess)<br>
Cells for raster map &lt;all_nc_be_20ft_bspline&gt; will be interpolated<br>
ERROR: G_calloc: unable to allocate 18446744071812941729 * 8 bytes of<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; memory at dalloc.c:66</font>
<br>
<br><font size=3>This was for a region with the following settings:</font>
<br>
<br><font size=3>projection: 99 (Lambert Conformal Conic) &nbsp;- North
Carolina State Plane Feet NAD83</font>
<br><font size=3>zone: &nbsp; &nbsp; &nbsp; 0</font>
<br><font size=3>datum: &nbsp; &nbsp; &nbsp;nad83</font>
<br><font size=3>ellipsoid: &nbsp;grs80</font>
<br><font size=3>north: &nbsp; &nbsp; &nbsp;1050000</font>
<br><font size=3>south: &nbsp; &nbsp; &nbsp;30000</font>
<br><font size=3>west: &nbsp; &nbsp; &nbsp; 400000</font>
<br><font size=3>east: &nbsp; &nbsp; &nbsp; 3060000</font>
<br><font size=3>nsres: &nbsp; &nbsp; &nbsp;20</font>
<br><font size=3>ewres: &nbsp; &nbsp; &nbsp;20</font>
<br><font size=3>rows: &nbsp; &nbsp; &nbsp; 51000</font>
<br><font size=3>cols: &nbsp; &nbsp; &nbsp; 133000</font>
<br><font size=3>cells: &nbsp; &nbsp; &nbsp;6783000000</font>
<br>
<br><font size=3><br>
<br>
I then zoomed into the northeast corner of the state and set the region
to :<br>
north:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1047363<br>
south:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 795102<br>
west:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2505067<br>
east:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3007693<br>
nsres:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20.00007928<br>
ewres:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 20.00023875<br>
rows:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12613<br>
cols:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 25131<br>
cells:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 316977303<br>
<br>
and tried again for this area<br>
<br>
GRASS 6.5.svn (ncstpft_nad83):~ &gt; v.surf.bspline input=all_nc_be_pts2
raster=ne_nc_be_20ft_bspline sie=40 sin=40 layer=0<br>
WARNING: Coor files of vector map &lt;all_nc_be_pts2@statewide&gt; is larger
than<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; it should be (158913789952
bytes excess)<br>
Cells for raster map &lt;ne_nc_be_20ft_bspline&gt; will be interpolated<br>
subregion 1 of 5459<br>
<br>
I killed that one and set the sin and sie to 400 rather than 40 .&nbsp;
This dropped the regions down to 66.&nbsp; <br>
<br>
GRASS 6.5.svn (ncstpft_nad83):~ &gt; v.surf.bspline input=all_nc_be_pts2
raster=ne_nc_be_20ft_bspline sie=400 sin=400 layer=0<br>
WARNING: Coor files of vector map &lt;all_nc_be_pts2@statewide&gt; is larger
than<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; it should be (158913789952
bytes excess)<br>
Cells for raster map &lt;ne_nc_be_20ft_bspline&gt; will be interpolated<br>
subregion 1 of 66<br>
</font>
<br><font size=3>This is about halfway done, so I should see the results
in another 1.5 days or so. &nbsp;</font>
<br>
<br><font size=3>I have two questions. &nbsp;First, is there a problem
with importing large point datasets that is being highlighted by the warning?
&nbsp;I did notice that the size of the coords and hist files did not match
for the point data set.</font>
<br>
<br><font size=3>&nbsp; &nbsp; &nbsp; &nbsp; Second, &nbsp;is about
v.surf.bspline and sin and sie. &nbsp;I have seen in the documentation
that the sie an sin needs to be twice as large as the spacing between points
as a good starting point. &nbsp; The density of the lidar data sets that
make up this ground point layer varies between 1m and 5 m . &nbsp; Obviously,
making sin and sie larger reduces the number of subregions, but how does
that affect the accuracy of the surface generated? &nbsp;If I'm trying
to generate a 20ft grid, is using sin and sie options of 400 absurd? &nbsp;(
I.E., &nbsp;It works faster but the results are less accurate)</font>
<br><font size=3><br>
</font>
<br><font size=3>The version of GRASS used was the weekly snapshot from
12/25/2010. &nbsp;The computer has a 1.6 GHz Xeon quad-core cpu with 16
GB RAM.</font>
<br>
<br>
<br><font size=3>Doug</font>
<br>
<br>
<br>
<br>
<br><font size=2 face="sans-serif">Doug Newcomb &nbsp; &nbsp; &nbsp; &nbsp;
&nbsp; &nbsp; <br>
USFWS<br>
Raleigh, NC<br>
919-856-4520 ext. 14 doug_newcomb@fws.gov<br>
---------------------------------------------------------------------------------------------------------<br>
The opinions I express are my own and are not representative of the official
policy of the U.S.Fish and Wildlife Service or Dept. of the Interior. &nbsp;
Life is too short for undocumented, proprietary data formats.</font>