LIDAR Classification

Purpose

To identify different types of LIDAR return and what they relate to. For instance the basic groups will be: ground, vegetation, buildings and noise. For the purpose of ARSF LIDAR deliveries we just classify noisy points. These typically tend to be due to the atmospheric conditions (cloud and haze), signal noise or single isolated points.

Method

We can use the TerraScan software or the in-house classify_las program for point classification. The TerraScan usage can be found on other pages of the wiki.

classify_las

This is a c++ utility written to do basic classification on LAS files. The current classification algorithms include:

• Isolated points - identify points as noise (classification 7) if there is no other point within x metres of it.
• Absolute elevation - identify points as noise (classification 7) if they are above or below a given elevation.

A fence can be used to only classify points within a rectangular boundary (in any orientation).

Usage

 -lasfilein <LASFILENAME>    The LAS file to be imported and classified
 -lasfileout <LASFILENAME>   The name of the output LAS file 

optional flags

 -fence minx miny maxx maxy  Use a fence whose bounds are given by minx,miny maxx,maxy   
 -fence px py qx qy width    Use a fence defined by centre line pq and width. width is the full width of the rectangle not the distance 
                             from the centre line to the edge.
 -exclude "<macro>"          Exclude points that conform to the given macro from being classified. Macro can be one or more of 
                             the following separated by spaces:
                                i<a  intensity less than a (i.e. do not classify points with intensity less than a)
                                i>b  intensity greater than b (i.e. do not classify points with intensity greater than b)
                                e<c  elevation (z) less than c
                                e>d  elevation (z) greater than d
                                nret=f number of returns = f (i.e. do not want to classify any points which are only single returns nret=1)

method flags

 -isolated x g               Run the isolated points algorithm using a distance of x metres and group size g
 -above x                    Run the absolute elevation algorithm classifying points with elevation above x metres
 -below x                    Run the absolute elevation algorithm classifying points with elevation below x metres   
 -intensity x                Run the intensity based algorithm classifying points below a certain intensity level
 -cloud above|below [a [m [t [i]]]] 
                             Run the cloud search algorithm to classify cloud points either above or below ground. This can be run using the 
                             default values (shown in brackets) or values set on the command line. This method can be run iteratively since 
                             when the mean and std deviation are calculated, noisy points are ignored.
                                a = area over which to calculate the mean (500 metres)
                                m = multiplier for distance away from mean eg mean - m*std deviation (3.0)
                                t = threshold to ignore classification eg areas with std deviations less than t are ignored (4.0 metres)
                                i = intensity value threshold - ignore classifying points with intensity greater than i (10)

The output LAS file will no longer have its points ordered by time, but rather by the quadtree / pointbucket ordering which is used to speed up searching within the software. Any number of methods can be run through on the same command.

Example - to classify isolated points whom have less than 5 points within 4 metres of them

classify_las -lasfilein example.LAS -lasfileout classified.LAS -isolated 4 5

classify_las_files.sh

It is also possible to batch process las files within a folder using classify_las_files.sh.

Uses classify_las to preform the classification. 

Arguments:
    -d:           Directory containing LAS files to be classified or single LAS file to be classified
    -c:            "classification_algorithim(s)_to_be_run" 


Optional:
    -o:           Directory to outout classified LAS files, defaults to input directory  

    -h:           Display usage
    -a:           Details of classification algorithims that are available 
                 (or see http://arsf-dan.nerc.ac.uk/trac/wiki/Processing/LidarClassification)

Example

classify_las_files.sh -d ~airborne/workspace/PROJCODE/leica/proc_laser -c "-isolated 4 4 -below 250 -above 1000"

Full Waveform

If there is full wave form data you will need to copy the classifications from discrete into the fullwaveform files. Classify_las and lag both sort the las files for classification. We need to sort this back to points by time before copying.

You can either create a discrete lidar delivery first (which will sort the las files for you) or run the following command on each line:

lassort -i input.LAS (-o output.LAS)

You can also use lassort.sh which will sort an entire directory for you:

lassort.sh -i inputdir (-o outputdir)

Note that the default is to simply overwrite the files. Specify an output if you don't want to overwrite.

Run the following script after sorting the discrete points:

get_class_las -c classified.LAS -i unclassifiedfw.LAS -o output.LAS

This will copy your classified points into the fw data.

To run this on a whole directory (before discrete lidar delivery creation), run:

batch_fw_class.sh -c classified-directory/ -i unclassified-directory/ -o output-directory/

Note, this script won't sort the LAS files for you.