PDAL DSM type set to float, 2.5D dsm res tweak

opendm/dem/pdal.py

@@ -63,7 +63,8 @@ def json_gdal_base(fout, output, radius, resolution=1):
         'resolution': resolution,
         'radius': radius,
         'filename': '{0}.{1}.tif'.format(fout, output[0]),
-        'output_type': output[0]
+        'output_type': output[0],
+        'data_type': 'float'
     })
 
     return json

opendm/mesh.py

@@ -1,5 +1,5 @@
 from __future__ import absolute_import
-import os, shutil, sys, struct, random
+import os, shutil, sys, struct, random, math
 from gippy import GeoImage
 from opendm.dem import commands
 from opendm import system
@@ -19,7 +19,7 @@ def create_25dmesh(inPointCloud, outMesh, dsm_resolution=0.05, depth=8, samples=
     os.mkdir(tmp_directory)
     log.ODM_INFO('Created temporary directory: %s' % tmp_directory)
 
-    radius_steps = [dsm_resolution / 4.0, dsm_resolution / 2.0, dsm_resolution]
+    radius_steps = [dsm_resolution * math.sqrt(2)]
 
     log.ODM_INFO('Creating DSM for 2.5D mesh')
 

scripts/odm_meshing.py

@@ -92,6 +92,7 @@ class ODMeshingCell(ecto.Cell):
                     dsm_resolution=dsm_resolution,
                     depth=self.params.oct_tree,
                     maxVertexCount=self.params.max_vertex,
+                    samples=self.params.samples,
                     verbose=self.params.verbose,
                     max_workers=args.max_concurrency)
           else: