Changes between Initial Version and Version 1 of Procedures/DeliveryCreation/LiDARElevationDiffs


Ignore:
Timestamp:
Jun 6, 2017 2:47:55 PM (2 years ago)
Author:
dac
Comment:

Initial page

Legend:

Unmodified
Added
Removed
Modified
  • Procedures/DeliveryCreation/LiDARElevationDiffs

    v1 v1  
     1= Calculating Elevation differences for LiDAR data =
     2
     3To calculate the difference between overlapping LiDAR lines a raster is generated from each. For each overlapping pairs the difference in elevation pixels in areas of with low variance (i.e., flat) is calculated.
     4
     5To run for a delivery the following script is used:
     6
     7{{{
     8mkdir -p /tmp/2015_249_lidar_diffs
     9calculate_lidar_elevation_diffs.py -c lidar_genreadme-airborne.cfg \
     10  -p UTM27N
     11  -i delivery/EUFAR15_48-249-lidar-20160223/flightlines/las1.2/
     12  -o /tmp/2015_249_lidar_diffs
     13}}}
     14
     15This will produce an output similar to:
     16
     17```
     18001a 002a 35; 002a 003a 18; 003a 004a 9; 004a 005a 7; 005a 006a 20; 006a 007a 4; 007a 008a 11; 008a 009a 10; 009a 010a 6; 010a 011a 5; 011a 012a 6; 012a 013a 8; 013a 014 Insufficient overlapping pixels; 014 HNI4_cross Insufficient overlapping pixels; HNI4_cross 001b No-overlap; 001b 002b 20; 002b 003b 14; 003b 004b 15; 004b 005b 14; 005b 006b 16; 006b 007b 12; 007b 008b 15; 008b 009b 11; 009b 010b 7; 010b 011b 8; 011b 012b 5; 012b 013b 8; 013b HNI3_cross Insufficient overlapping pixels;
     19```
     20
     21Which can be pasted into the readme config file. The script assumes sequential lines overlap which isn't always the case.
     22
     23To run the script for individual lines use:
     24
     25{{{
     26create_lidar_diff_image.py -d /tmp/2015_249_lidar_diffs/ \
     27   -o /tmp/2015_249_lidar_diffs/001a_002a_diff.tif \
     28   --projection UTM27N \
     29   delivery/EUFAR15_48-249-lidar-20160223/flightlines/las1.2/LDR-EUFAR15_48-2015-249-01.LAS \
     30   delivery/EUFAR15_48-249-lidar-20160223/flightlines/las1.2/LDR-EUFAR15_48-2015-249-02.LAS
     31}}}
     32
     33This will produce the following output:
     34{{{
     35*******************
     36 Min:       0.000
     37 Max:       1.302
     38 Mean:      0.355
     39 StdDev:    0.071
     40 Median:    0.351
     41*******************
     42}}}
     43Copy the median difference into the config file.
     44
     45For areas where there is dense vegetation there might be too few pixels above the variance threshold to calculate a difference. If this is the case a warning will be printed (`Too few valid pixels overlap to be able to output a reliable difference between lines 010a and 011a`). The default variance is 0.1 m, try increasing to 0.2 m using `--max_variance 0.2` and running again to see if this provides sufficient pixels. You can check which areas are being used for the difference and which are masked out (set to NaN) by opening the difference images (e.g., 001a_002a_diff.tif in TuiView or ENVI). If you change the variance make sure it is noted in the ticket.
     46
     47Once the difference table has been populated for all lines the temp directory should be removed.