| 1 | = Processing Full Waveform LiDAR data = |
| 2 | |
| 3 | Although there are not many tool capable of processing LAS1.3 full waveform LiDAR data, you may find some of the following useful: |
| 4 | |
| 5 | [wiki:Help/Fullwaveform#ARSFextractiontool ARSF full waveform extraction tool] [[BR]] |
| 6 | [wiki:Help/Fullwaveform#Lastools LAStools] [[BR]] |
| 7 | |
| 8 | [wiki:Help/Fullwaveform#Terrascan TerraScan ] (License required) |
| 9 | |
| 10 | |
| 11 | == ARSF full waveform extraction tool == #ARSFextractiontool |
| 12 | |
| 13 | Example commands: |
| 14 | |
| 15 | {{{ |
| 16 | python fwf_extract.py LDR-FW-FW10_02-201010101.LAS |
| 17 | }}} |
| 18 | |
| 19 | Requests user to input on command line the limits of area to extract and saves full waveform data to current directory |
| 20 | |
| 21 | {{{ |
| 22 | python fwf_extract.py LDR-FW-FW10_02-201010101.LAS -o /lidar/extraction |
| 23 | }}} |
| 24 | |
| 25 | Requests user to input on command line the limits of area to extract and saves full wavefrom data to /lidar/extraction |
| 26 | |
| 27 | {{{ |
| 28 | python fwf_extract.py LDR-FW-FW10_02-201010101.LAS -o /lidar/extraction --area 560002 560000 318002 318000 |
| 29 | }}} |
| 30 | |
| 31 | Extracts wavefrom data for area bound by 560002 560000 318002 318000 to /lidar/extraction |
| 32 | |
| 33 | {{{ |
| 34 | python fwf_extract.py LDR-FW-FW10_02-201010101.LAS --header |
| 35 | }}} |
| 36 | |
| 37 | Prints the Header of the LAS file to the screen and exits |
| 38 | |
| 39 | |
| 40 | |
| 41 | |
| 42 | NOTE: You need to have Python installed on your computer and to add the directroy of your Python distribution to your PATH environment variable. [[BR]] |
| 43 | See the Python user guide for more information: http://docs.python.org/using/windows.html#excursus-setting-environment-variables |
| 44 | |
| 45 | |
| 46 | |
| 47 | == LAStools == #Lastools |
| 48 | |
| 49 | |
| 50 | |
| 51 | Below are some useful commands for extracting the waveform from the LAS1.3 files. |
| 52 | |
| 53 | To download the tools and for more information and help see http://www.cs.unc.edu/~isenburg/lastools/ |
| 54 | |
| 55 | |
| 56 | {{{ |
| 57 | las2txt -i LDR-FW-GB12_06-2012-083b-01.LAS -parse MrnxyzWV -o LDR-FW-GB12_06-2012-083b-01.txt |
| 58 | }}} |
| 59 | |
| 60 | converts the LAS file LDR-FW-GB12_06-2012-083b-01.LAS to ASCII text using format MrnxyzWV |
| 61 | |
| 62 | {{{ |
| 63 | las2txt -i LDR-FW-GB12_06-2012-083b-01.LAS -parse MrnxyzWV -keep_wavepacket 1 -o LDR-FW-GB12_06-2012-083b-01.txt` |
| 64 | }}} |
| 65 | |
| 66 | converts the LAS file LDR-FW-GB12_06-2012-083b-01.LAS to ASCII text using format MrnxyzWV keeping only those points containing waveform information |
| 67 | |
| 68 | {{{ |
| 69 | laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform |
| 70 | }}} |
| 71 | |
| 72 | compresses the LAS file LDR-FW-GB12_06-2012-083b-01.LAS creating two compressed files |
| 73 | * LDR-FW-GB12_06-2012-083b-01.laz - containing the compressed discrete data and wave packet information |
| 74 | * LDR-FW-GB12_06-2012-083b-01.wpz - containing the compressed wave packet data |
| 75 | |
| 76 | {{{ |
| 77 | laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform |
| 78 | }}} |
| 79 | |
| 80 | decomressses the laz file LDR-FW-GB12_06-2012-083b-01.laz creating two file |
| 81 | - LDR-FW-GB12_06-2012-083b-01.las - containing the discrete data and wave packet information |
| 82 | - LDR-FW-GB12_06-2012-083b-01.wpd - containing the wave packet data |
| 83 | |
| 84 | {{{ |
| 85 | lasclip -i LDR-FW-GB12_06-2012-083b-01.laz -poly clip.shp -o LDR-FW-GB12_06-2012-083b-01_clip.laz` |
| 86 | }}} |
| 87 | subsets the laz file LDR-FW-GB12_06-2012-083b-01.laz using the polygon contained in clip.shp and ouputs the result to LDR-FW-GB12_06-2012-083b-01_clip.laz. Note this will only subset the laz file - the wpz file will still contain all points. |
| 88 | |
| 89 | |
| 90 | If you wish to extract the waveform data for a particular area bounded by 407049 200380 407067 200398 the workflow would be as follows: |
| 91 | |
| 92 | {{{ |
| 93 | laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform |
| 94 | laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform |
| 95 | las2txt -i LDR-FW-GB12_06-2012-083b-01.las -clip 407049 200380 407067 200398 -parse MrnxyzWV -o LDR-FW-GB12_06-2012-083b-01_bounding_box.txt |
| 96 | }}} |
| 97 | |
| 98 | If you wish to extract the waveform data for a particular area conatined in clip.shp the workflow would be as follows: |
| 99 | |
| 100 | {{{ |
| 101 | laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform |
| 102 | laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform |
| 103 | lasclip -i LDR-FW-GB12_06-2012-083b-01.las -poly clip.shp -o LDR-FW-GB12_06-2012-083b-01_clip.las |
| 104 | rename LDR-FW-GB12_06-2012-083b-01.wpd to LDR-FW-GB12_06-2012-083b-01_clip.wpd |
| 105 | las2txt -i LDR-FW-GB12_06-2012-083b-01_clip.las -parse MrnxyzWV -o LDR-FW-GB12_06-2012-083b-01_clip.txt |
| 106 | }}} |
| 107 | |
| 108 | |
| 109 | == !TerraScan == #Terrascan |
| 110 | |
| 111 | In order to view the full waveform data you first need to create trajectory files. With the !RawLaser data loaded in ALSPP select Utilities > Generate trj files. This creates an individual trajectory file for each flightline. |
| 112 | |
| 113 | In !TerraScan, select Manage Trajectories ( /// icon) from the General tool bar. Under File > Set Directory select the location of the trj files. |
| 114 | |
| 115 | Alternatively you can use the .sol file. Under Manage Trajectories > File > Import Files > ***.sol. This creates a .trj covering the entire project (which is included with the delivery). |
| 116 | |
| 117 | Load in a full waveform .LAS files. It is advisable to create a small fence and only load points inside the fence as !TerraScan can't seem to cope with entire flightlines. |
| 118 | |
| 119 | From the main !TerraScan window select Flightline > Deduce Using Time. |
| 120 | |
| 121 | On the main !TerraScan toolbox select View Waveform Data to open the waveform viewer (the sideways barchart icon, on the bottom right of the Main toolbar). Under Settings > Display Settings you can adjust how the wave is viewed by varying maximum value (x axis) or sample height (y axis). |
| 122 | |
| 123 | |
| 124 | === !WaveForm Extractions: === |
| 125 | |
| 126 | The PI may request waveform extractions to be included with the delivery. These are ASCII files (one file per return) containing the full waveform information. |
| 127 | |
| 128 | * Draw a fence around the area of interest |
| 129 | * Open the View Waveform Data window |
| 130 | * File > Save As Text |
| 131 | |
| 132 | |
| 133 | === Extract Echoes: === |
| 134 | |
| 135 | This tool can be used to create additional discrete points from the full waveform data. Use the ambient noise value to set the threshold for detecting peaks (on the waveform viewer under Settings > Display Settings). |
| 136 | |
| 137 | * Draw a fence around the area of interest |
| 138 | * Select Extract Echoes from the waveform toolbar |
| 139 | |
| 140 | |
| 141 | |
| 142 | |
| 143 | |
| 144 | |