@ -406,7 +406,7 @@ def process(yaml_file):
# Interpolate point cloud onto a grid
print('Interpolating to grid...')
xyz_name = os.path.join(tmp_dir, las_basename + '.xyz')
call_lastools('las2dem', input=las_classified_name, output=xyz_name,
call_lastools('blast2dem', input=las_classified_name, output=xyz_name,
args=['-step', 1], verbose=False)
# Make runup clipping mask from gridded point cloud