Module 2.3: Surfaces - Accuracy in DEMs

The purpose of this week's lab was to determine the vertical accuracy of DEMs, including determining bias. I was given a high resolution DEM for a section of North Carolina created from LIDAR data and a table of field data of elevation at the ground surface collected using high accuracy survey methods. The table contained 5 land cover types (a-e). The table had coordinates for each point which were converted to a point shapefile using the XY Table to Point tool. The points within the DEM were selected and saved as a separate shapefile for the analysis.

To get the values of the raster beneath the points, I used the Extract Multi Values to Points. This tools grabs the elevation value of the pixel in the DEM directly beneath each sample point and then adds a new field to the point shapefile with that value. Because these values were in feet, I added a new field to the attribute table and converted the feet to meters. I then calculated accuracy for each land cover type as well as the entire data set. The calculations included the root mean square (RMSE), accuracy at the 68th percentile, accuracy at the 95th percentile, and Bias (or mean error). The summary of the results of these calculations is below:


RMSE is the most common metric used to evaluate accuracy. The lower the RMSE, the more accurate the result. RMSE was calculated using the following formula from Fisher and Tate (2006)*:

I calculated bias using the mean error (ME) using the following equation from Fisher and Tate (2006)*:

Using this calculation, bias or ME can either be positive or negative. This represents the overall under- or overestimation of the elevations in the DEM. Positive values indicate where values were overestimated whereas negative values indicate areas where values were underestimated. Both types of bias exist in this data set and vary by land cover type.

*Fisher, P.F., and N.J. Tate. 2006. Causes and consequences of error in digital elevation models. Progress in Physical Geography, 30(4), 467-489.


Comments

Post a Comment

Popular posts from this blog

Module 6: Proportional Symbol and Bivariate Choropleth Mapping

Image
The objectives of this week's lab was to explore proportional symbol and bivariate choropleth mapping techniques, as well as customizing legends for each technique. One of the data sets for this lab included positive and negative values. The final map utilizes proportional symbols to show both positive and negative values. This map is essentially two proportional symbol maps overlaid on each other, one for positive job growth and another for negative job growth, with the symbols for both variables sized the same. Utilizing proportional symbols to show both positive and negative values is an effective way to visualize areas that show the greatest and lowest values. These areas are easy to see quickly on the map and give the map reader a simple way to visualize the variable being mapped. Fig 1. Proportional symbol map using positive and negative values.  Another data set for this lab examined contained two data variables: percent of the US population in each state that is obe
Read more

Module 2: Coordinate Systems

Image
In this week's lab I explore different coordinate systems and projections in order to better understand the nature of distortions introduced by the systems. For the last part of the lab, I was tasked with selecting a state in the contiguous US and selecting an appropriate projection for the area. I selected the state of Wyoming as my area of interest. Wyoming has 4 State Plane zones and two UTM zones, so neither a single State Plane nor UTM system would work for this area without distorting the size or shape of the state. Therefore, the best coordinate system for this area is Custom Coordinate system. I chose NAD 1983 WyLam (Meters), a custom system specifically for the state of Wyoming, as the projected coordinate system. This coordinate system is a Lambert secant-case normal-aspect conic type projection which preserves shape and is appropriate for a general reference map. The standard parallels are placed at 41N and 45N, the boundaries of the state of Wyoming, which minimizes
Read more

Module 4: Color Concepts & Choropleth Mapping

Image
In this week's lab, I explored using color and selecting color ramps in ColorBrewer, as well as applying these concepts to create meaningful choropleth maps. I created a map of population change from 2010 to 2014 in Georgia counties. I normalized the data and produced a choropleth map of population change using a diverging color scheme with 6 classes divided around a midpoint of 0. Fig 1. Choropleth map of population change in Georgia counties using a diverging color scheme. For the Georgia population data set, I chose to divide the data set into 6 classes around a midpoint of 0: 3 showing a decrease in population and 3 showing an increase in the population. The classes were grouped based on a manual classification closely based on the Natural Break classification, where the midpoint was adjusted to 0 with the other groups unchanged. Natural Breaks classification creates classes based on similarities in values in the data, thereby revealing patterns in the data. By adjust
Read more