danh

Python scripts for processing lidar surveys from the Central Coast.

You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

Go to file

Dan Howe fcc6b810c8 Use temp files instead of stdout piping, to support older lastools versions		6 years ago
param-files	Update inputs for corrected input las files	6 years ago
README.md	Update README	6 years ago
generate_heatmaps.py	Update docstring to explain relative paths in arcpy	6 years ago
las_manipulation.py	Use temp files instead of stdout piping, to support older lastools versions	6 years ago
las_outputs.py	Export classified, clipped las for delivery	6 years ago
nielsen_volumes.py	Suppress numpy floating point arithmetic warnings	6 years ago
plot_heatmaps.py	Restore export to jpg (removed in previous commit)	6 years ago
survey_tools.py	Correctly handle profiles that are outside survey area	6 years ago

README.md

Central Coast aerial lidar processing

Python scripts for cleaning, classifying, and extracting survey profiles from las point clouds.

1. Enter input parameters

Input parameters are entered into excel spreadsheets, e.g. param-files/survey-1.xlsx. The python scripts take their input in yaml format. The spreadsheets can be converted to yaml by running:

python excel_to_yaml.py

2. Manipulate las data

las_manipulation.py does the following to the input las file:

clip to beach extent
classify ground points (so vegetation and buildings can be removed later)
draw a polygon to remove swash zone interference

Example usage:

# Process single survey at specific beach
python las_manipulation.py param-files/survey-1-avoca.yaml

# Process single survey at multiple beaches
python las_manipulation.py param-files/survey-1-avoca.yaml param-files/survey-1-pearl.yaml

# Process all surveys at specific beach
python las_manipulation.py param-files/*avoca.yaml

# Process all beaches for specific survey date
python las_manipulation.py param-files/survey-1*.yaml

3. Check swash clipping mask

The swash clipping mask should be checked visually to ensure it is correct.

4. Export results

las_outputs.py does the following, based on the outputs of las_manipulation.py:

remove swash zone
extract elevations from point cloud at profile locations
complete volume analysis based on Nielsen et al. (1992)
save figures of profile cross sections
save raster of ground surface

Example usage:

python las_outputs.py param-files/survey-1-avoca.yaml

5. Generate erosion heatmaps

generate_heatmaps.py creates a heatmap raster by calculating the elevation difference between the current and previous surveys.

generate_heatmaps.py relies on arcpy, so this script must be run using the version of python that was installed with ArcGIS.

Example usage:

C:/Python27/ArcGIS10.4/python generate_heatmaps.py param-files/survey-2-avoca.yaml

6. Export erosion heatmaps

plot_heatmaps.py exports a jpg erosion heatmap at all of the beaches covered by the specified heatmap tif.

plot_heatmaps.py relies on arcpy, so this script must be run using the version of python that was installed with ArcGIS.

Example usage:

C:/Python27/ArcGIS10.4/python plot_heatmaps.py param-files/survey-2-avoca.yaml