Hello Morten,
BLAST. The best answer to rasterizing large point clouds is BLAST.
blast2dem -i tiles\*.las -merged -keep_class 2 -o dtm.asc
But you may also try the '-extra_pass' option.
las2dem -i tiles\*.las -merged -keep_class 2 -o dtm.asc -extra_pass
las2dem.exe is a "traditional" in-core algorithm that will first load
all the points into memory, then triangulate them, and then raster the
triangles into a DEM. Triangulating 89,002,420 points results in about
178,004,840 triangles. Since each point needs at least 12 bytes
storage and each triangle at least 28 bytes storage you cannot do this
for that many points/triangles given the 2GB main memory limitation of
a 32 bit windows executable.
blast2dem.exe is a streaming out-of-core algorithm that that will only
keep a fraction of the points in memory at any given time. it is
designed to triangulate points sets that are magnitudes larger than
the available main memory. it creates a seamless TIN that is identical
to that created by las2dem but streams the triangles through memory,
rasterizing the "fionalized" ones, and re-using their memory. If you
are interested this is described in more detail in these two
publications [1,2].
That said, since you are only triangulating the ground points with
classification 2 can also help las2dem.exe by setting the '-
extra_pass' flag. This will cause las2dem.exe to make two passes over
the LAS file. In the first pass it merely counts how many points will
actually be triangulated (i.e. have classification 2) and maybe the
remaining points will fit into memory.
Regards,
Martin @lastools
[1] Martin Isenburg, Yuanxin Liu, Jonathan Shewchuk, Jack Snoeyink,
Streaming Computation of Delaunay Triangulations, Proceedings of
SIGGRAPH'06, pages 1049-1056, July 2006. online at
http://www.cs.unc.edu/~isenburg/papers/ilss-scdt-06.pdf
[2] Martin Isenburg, Yuanxin Liu, Jonathan Shewchuk, Jack Snoeyink,
Tim Thirion, Generating Raster DEM from LiDAR via TIN Streaming,
GIScience'06 Conference Proceedings, pages 186-198, September 2006.
online at
http://www.cs.unc.edu/~isenburg/papers/ilsst-tin2dem-06.pdf
On Mar 13, 2:59 am, MAG <
mor...@mortenaggerholm.dk> wrote:
> I tried to create a DTM from 13 LAS filed using Las2dem.exe. I
> receives the following error:
>
> B:\Test_Area_3\PointCloud>las2dem -i tiles\*.las -merged -keep_class 2
> -o dtm.asc
>
> Please license from '
martin.isenb...@gmail.com' to use LAStools
> commercially.
> WARNING: unlicensed. over 10 million points. inserting black diagonal.
> ERROR: failed malloc for 178004840 TINtriangles.
> ERROR: cannot alloc enough TIN triangles to triangulate 89002420
> points.
> contact
martin.isenb...@gmail.com for help on what to do.
>
> /Morten Aggerholm