Changes between Initial Version and Version 1 of Help/Fullwaveform


Ignore:
Timestamp:
Jun 26, 2012, 4:34:24 PM (12 years ago)
Author:
emca
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Help/Fullwaveform

    v1 v1  
     1 = Processing Full Waveform LiDAR data =
     2
     3Although 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
     13Example commands:
     14
     15{{{
     16python fwf_extract.py  LDR-FW-FW10_02-201010101.LAS
     17}}}
     18 
     19Requests user to input on command line the limits of area to extract and saves full waveform data to current directory
     20
     21{{{
     22python fwf_extract.py  LDR-FW-FW10_02-201010101.LAS -o /lidar/extraction
     23}}}
     24
     25Requests user to input on command line the limits of area to extract and saves full wavefrom data to /lidar/extraction
     26
     27{{{
     28python fwf_extract.py  LDR-FW-FW10_02-201010101.LAS -o /lidar/extraction --area 560002 560000 318002 318000
     29}}}
     30
     31Extracts wavefrom data for area bound by 560002 560000 318002 318000 to /lidar/extraction
     32
     33{{{
     34python fwf_extract.py  LDR-FW-FW10_02-201010101.LAS --header
     35}}}
     36
     37Prints the Header of the LAS file to the screen and exits
     38
     39
     40
     41
     42NOTE: You need to have Python installed on your computer and to add the directroy of your Python distribution to your PATH environment variable. [[BR]]
     43See 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
     51Below are some useful commands for extracting the waveform from the LAS1.3 files.
     52
     53To download the tools and for more information and help see http://www.cs.unc.edu/~isenburg/lastools/
     54
     55
     56{{{
     57las2txt -i LDR-FW-GB12_06-2012-083b-01.LAS -parse MrnxyzWV -o LDR-FW-GB12_06-2012-083b-01.txt
     58}}}
     59
     60converts the LAS file LDR-FW-GB12_06-2012-083b-01.LAS to ASCII text using format MrnxyzWV
     61
     62{{{
     63las2txt -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
     66converts 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{{{
     69laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform
     70}}}
     71
     72compresses 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{{{
     77laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform
     78}}}
     79
     80decomressses 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{{{
     85lasclip -i LDR-FW-GB12_06-2012-083b-01.laz -poly clip.shp -o LDR-FW-GB12_06-2012-083b-01_clip.laz`
     86}}}
     87subsets 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
     90If 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{{{ 
     93laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform
     94laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform
     95las2txt -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
     98If you wish to extract the waveform data for a particular area conatined in clip.shp the workflow would be as follows:
     99
     100{{{
     101laszip -i LDR-FW-GB12_06-2012-083b-01.LAS -waveform
     102laszip -i LDR-FW-GB12_06-2012-083b-01.laz -waveform
     103lasclip -i LDR-FW-GB12_06-2012-083b-01.las -poly clip.shp -o LDR-FW-GB12_06-2012-083b-01_clip.las
     104rename LDR-FW-GB12_06-2012-083b-01.wpd to LDR-FW-GB12_06-2012-083b-01_clip.wpd
     105las2txt -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
     111In 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
     113In !TerraScan, select Manage Trajectories ( /// icon)  from the General tool bar. Under File > Set Directory select the location of the trj files.
     114
     115Alternatively 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
     117Load 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
     119From the main !TerraScan window select Flightline > Deduce Using Time.
     120
     121On 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
     126The 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
     135This 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