Having problems with blast2dem64

[Malcolm Williamson]

Hello,

I'm running across an error I don't understand on a fairly large set of LiDAR tiles from the USGS. As you can see below, I get a Fatal Error at one point, but then it keeps on running. When it finishes, I end up with a .tif file of only 1 KB, and then a large set of .tmp files. Thoughts?

E:\Photogrammetry\InSAR_more\More_Fayetteville_LiDAR\Denoised>blast2dem64 -i *AR_Benton*.laz -merged -otif -o Fayetteville_AR_LiDAR_DSM_1.tif -nad83_2011 -utm 15N -keep_class 1 2 -first_only
terrain finalization (using 2d quadtree) to level 11 with empty 0.
have bounding box. first pass not needed.
second pass: populating finalization counter grid ...
second pass: done. created 1785807 leaf grid cells and stored 1785807 sprinkle points
spfinalizer: needed 1236.2 seconds
created parent grid cells (new total is 2387329 cells)
third pass: using finalization grid to spatially finalize 2008778376 points...
output precision is float ...
lower left x unspecified. setting it to 357000.
lower left y unspecified. setting it to 3.99551e+06.
step size was unspecified. set to 1.
ncols was unspecified. we set it to 63000.
nrows was unspecified. we set it to 25988.
output file is written as BigTIFF, calculated data size of the raster before compression and without meta data is 6.099209 GB.
FATAL ERROR: vertex not in hash. need pre-order mesh
third pass: done. stored at most -2108908666 chunked points (plus 1785807 sprinkle points)
spfinalizer: needed 2737.89 seconds
duplicate_points 149816
spdelaunay2d: needed 2737.23 seconds
read 2008778376 points. produced TIN with 2008628560 vertices and -277710487 triangles
WARNING: set nverts 2008778376 but v_count 2008628560

E:\Photogrammetry\InSAR_more\More_Fayetteville_LiDAR\Denoised>

Thanks very much,

    -Malcolm

Malcolm D. Williamson
Geospatial Applications and Education Manager
Center for Advanced Spatial Technologies
University of Arkansas
Fayetteville AR 72701
479.501.2734  md...@uark.edu

[Jochen Rapidlasso]

Hi Malcolm,

probably you run into a size issue. Blast is proud to process files up to 2 billion points - you have 2,008,778,376 points, which is quite close to this limit.

Best would be to process without "-merged" and then put the TIFF images together in the next step (using e.g. QGIS).

Or just merge less files, so you stay a bit further from the 2 billion points.

Did you compare with the old blast2dem? Is it able to process there (maybe as PNG or JPG output)?

Cheers,

Jochen @rapidlasso