|
|
@ -228,7 +228,8 @@ def process(yaml_file):
|
|
|
|
las_data = call_lastools('lasclip', input=input_las, output='-stdout',
|
|
|
|
las_data = call_lastools('lasclip', input=input_las, output='-stdout',
|
|
|
|
args=['-poly', crop_swash_poly], verbose=False)
|
|
|
|
args=['-poly', crop_swash_poly], verbose=False)
|
|
|
|
|
|
|
|
|
|
|
|
las_name = 'crop-1.las'
|
|
|
|
# Export classified, clipped las for delivery to client
|
|
|
|
|
|
|
|
las_name = os.path.join(output_las_dir, las_basename + '.las')
|
|
|
|
with open (las_name, 'wb') as f:
|
|
|
|
with open (las_name, 'wb') as f:
|
|
|
|
f.write(las_data)
|
|
|
|
f.write(las_data)
|
|
|
|
|
|
|
|
|
|
|
@ -238,7 +239,6 @@ def process(yaml_file):
|
|
|
|
args=['-poly', crop_heatmap_poly], verbose=False)
|
|
|
|
args=['-poly', crop_heatmap_poly], verbose=False)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
# Create clipping polygon for heatmap raster
|
|
|
|
# Create clipping polygon for heatmap raster
|
|
|
|
print('Creating heat map cropping polygon...')
|
|
|
|
print('Creating heat map cropping polygon...')
|
|
|
|
shp_name = os.path.join(output_poly_dir, las_basename + '.shp')
|
|
|
|
shp_name = os.path.join(output_poly_dir, las_basename + '.shp')
|
|
|
|