Lab 9 - Corridor Analysis & Feature Extraction

GOAL AND BACKGROUND

The goal of this lab was to learn to project and navigate corridor LiDAR data. Feature extraction would also be explored.


METHODS

First, data's metadata was used to determine the projection of the data. The LAS tools in ArcMap were then used to define the projection for the data. The projected data was then opened in LP360, and some measurements were taken such as road widths and locations of power lines that were at risk for vegetation encroachment.

Next, building extraction was explored. Lake County data, classified in previous labs, was used for this. First, building footprints were extracted. Two shapefiles were created from this. One was a perfect outline of building points. The other was outlines with straight lines, which removed the noisy appearance of the footprints by making them more square. Lastly, the buildings heights were analyzed to determine qualification for FEMA's LOMA applications.

RESULTS

Shown below are two instances of vegetation encroachment on power lines found in the corridor data.

Shown below is an example of the difference between the non-squared footprints and the squared footprints, respectively. Notice the jagged lines in the non-squared footprints compared to the straight lines in the squared footprints.

Shown below is a map created of Lake County of buildings in the region that qualify for FEMA's LOMA applications. FEMA recently changed the cutoff from 810 ft ASL to 800 ft ASL. The map shows the buildings that the change affects.

SOURCES

Terrestrial LAS for Algoma, WI, project boundary KMZ, and metadata are from Ayres Associates.

Lab 8 - Vegetation Metrics Modeling

GOAL AND BACKGROUND

The goal of this lab was to learn to learn to calculate vegetation statistics using LiDAR data.


METHODS

First, LP360 was used to create a DTM for the canopy height of the study area. This was done by creating DTM images for both the canopy surface, using first return vegetation points, and ground surface, using last return ground points. A raster calculator was then used to subtract the canopy height from the ground to create a vegetation canopy height raster. There were some negative values in the data from error, and were selected and removed for further calculations.

Next, above ground biomass (AGB) was calculated using model builder in ArcMap. This was done by separating vegetation by species in the study area because species have different parameters for calculations.

This data was then used to calculate the breakdown of the AGB for each species based on stem, branch, and foliage mass. This was done by calculating percentages of each of these based on species and multiplying this by the previously found AGB.

RESULTS

A graph of distributions in vegetation height is shown below.

The canopy height raster is shown below. Negative values are indicated by black symbology.

The AGB map for the entire area is shown below. The species are marked with different colors. The black areas were non-vegetated areas.

Maps for the three parameters: stem, branch, and foliage mass are shown below.

The model created to find the AGB is shown below, followed by a section of the model showing the steps for a single species.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.

Lab 7 - Flood Inundation Modeling

GOAL AND BACKGROUND

The goal of this lab was to learn to use ArcScene to create a flood inundation model for the Eau Claire area.


METHODS

First, using ArcMap, breaklines were created and enforced for water bodies. A DTM image with flattened water was then created. These were brought into ArcScene. An animation was created by raising the water level by 5 feet per frame to simulate the rise in water level due to flooding. The model had to be adjusted using breaklines to ensure that the flooding begins from the river.


RESULTS

The final animation that was created is shown below.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.

Lab 6 - Topo-bathy Applications

GOAL AND BACKGROUND

The goal of this lab was to get experience working with topo-bathy lidar data.


METHODS

First, the data underwent some light QA/QC as some ground points were not well classified. After this was corrected, a shoreline breakline was created in ArcMap. This was then conflated. Next, a breakline was made of the unsubmerged topographical area.

Lastly, a DTM image was exported using the LAS data and the breaklines that were just created.


RESULTS

The image below shows the shoreline breakline that was created and conflated.

The image below shows the breakline that was created for the unsubmerged area.

The next images show the DTM image, followed by the hillshade image that was created.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.

Lab 5 - Breakline Creation, Conflation, and Enforcement

GOAL AND BACKGROUND

The goal of this lab was to learn how to create breaklines within ArcMap using LAS data, to use conflation to constrain downward slopes in rivers, and to enforce breaklines for hydro-flattening of water bodies.


METHODS

The LP360 extension in ArcMap was used for breakline creation. LiDAR data for Eau Claire was used for this section. First, the LAS data was added to ArcMap and displayed using the TIN surface. A ground filter was used to only display ground points. This was necessary so the edges of the river bank were shown instead of vegetation. First, the outline of the river was created, followed by a centerline of the river used to ensure a downward slope.

These shapefiles were brought into LP360. The breakline for the outline of the river was enforced. Then, the river-flattening tool was run, using the river bank shapefile and the river centerline shapefile.

Hydro-flattening was also performed on the Lake County data that the previous labs have used. There was no need to create breaklines for this as they were given, and no rivers in the data meant simpler hydro-flattening could be used. The breaklines were enforced and DTM and contour images were created.



RESULTS

Shown below is the final product in the creation of the centerline and riverbank shapefiles for the Eau Claire data.

Shown below is part of the table showing the constrained Z values along the centerline. The M value is how much the Z was adjusted to keep the river continuously downstream.

Shown below is the final result of the hydro-flattening in the Eau Claire data. Changes in elevation along the river are continuously downward as it flows downstream.

Next, the result from hydro-flattening in the Lake County data is shown. The simple ponds were flattened to be a single elevation.

A contour map of a section of the Lake County data was also created, shown below. It is overlaid on the NAIP imagery of the area.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.

Lab 4 - QA/QC

GOAL AND BACKGROUND

The goal of this lab was to learn how to perform relative and absolute QA/QC, along with doing manual QA/QC on classification.


METHODS

LP360  and the LP360 extension in ArcMap was used for this data processing. First, point cloud density was checked by creating a map that showed low density areas. According to this, point density was acceptable throughout the study area.

Next, relative accuracy was assessed. A map was created that displayed differences in elevation. This was used for swath-to-swath analysis. For this analysis, seamlines were created along flat, unobstructed areas within the overlap lines.

Next, absolute accuracy was assessed. This was done for both vertical and horizontal accuracy. Checkpoint locations for each were imported from a spreadsheet. Vertical checkpoints were referenced for measured values vs. values from the LiDAR data. Horizontal checkpoints were displayed on the map and refenced with the LiDAR data.

Manual QA/QC of classification errors involved looking through the study area and flagging classification errors. Once this was done, these errors were rectified using manual classification tools such as classifying in the profile window.


RESULTS

The map created in relative accuracy assessment showing differences in elevation is shown below. The horizontal lines are overlap of flight lines. The results from swath-to-swath analysis showed few areas above the maximum height value.

The results from the absolute accuracy showed the data was good. Vertical checkpoints were all within 0.10 cm error. Horizontal checkpoint analysis determined the data was in class 0 for both X and Y accuracy.

Manual QA/QC was done through much of the data, an example of this is shown below. The error was first flagged, as shown in the first image, then rectified, as shown in the second. The flag was removed after fixing the error.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.

Lab 3 - Vegetation Classification

GOAL AND BACKGROUND

The goal of this lab was to learn how to classify vegetation in LiDAR data. The vegetation would be classified by low, medium, and high vegetation. High noise would also be classified in this process. QA/QC will be necessary to minimize error.


METHODS

LP360 was used for this data processing. A height filter was used to classify data into either low, medium, or high vegetation. If the points were above the high vegetation threshold, they would be classified as high noise. Only unclassified points were chosen, since vegetation is the last to be classified.


RESULTS

The study area after classification is shown below. The newly classified vegetation is shown as green.

Below are a few examples of closer views of classification within the study area. First is a building that was shown in the previous posts, now with classified vegetation.

Below is a view of some residential housing in the study area. The algorithms had some trouble differentiating between houses and vegetation at some points, especially when it involved overlap. Manual cleanup was used after the algorithm.

SOURCES

Data obtained from Cyril Wilson for use in 358 LiDAR course.