36 GIS Applications, LiDAR, and Digital Landscapes (Posters)

Tuesday, 17 April 2018
8:00 AM - 5:30 PM
First Floor Lobby

SUPERCOMPUTING EXPERIMENTS TOWARDS MULTISCALE MACHINE-ASSISTED LANDFORM IDENTIFICATION FROM LIDAR DATA

Paper number: 36-4
Program: Theme Sessions
Time: 8:00 AM – 5:30 PM

Authors

Abstract

We have undertaken a project to detect landforms for soil mapping from large lidar point cloud elevation datasets. The goals of the project are to (a) from the lidar point cloud data, create bare earth digital elevation models (DEM) at any desired project scale and with any desired footprint, and make the process tractable at large scale; (b) detect ridges and gullies, the latter of which are critical for erosion analysis but, until lidar data became available, mappable only in the field; (c) extend the analysis to the range of landforms used in soil survey. That landforms occur across a range of scales and may be superimposed is a particular challenge. For example, an area in Illinois may include dunes, lake plains, and megascale glacial lineations superimposed upon a moraine. DEM construction is based on last returns from lidar point cloud data sets using IDW interpolation. Initial experiments using the open-source PointstoGrid algorithm implemented in GRASS processed the 6440 km2 Embarrass watershed of east central Illinois (976 GB of LAS data) in ~8 hr. Instrument and other noise was best removed with a series of Median and Lee filter operations. Implementation of the algorithm in C across 200 parallel processors reduced the entire processing time to 42 minutes.

For landform identification, we are experimenting with image analysis techniques first developed to differentiate ridge and valley morphology of human brains from Magnetic Resonance Imagery. The Difference of Gaussians (DoG) method generates a Gaussian Scale Space representation of the DEM by resampling it at progressively coarser grids, in our case from 3.95 ft to 101.12 ft. Each iteration is saved and subtracted from its previous (finer) scale to generate a new set of images, which represent the DoG at a given scale. Curvature across multiple scales can then be extracted by performing the processing on several of these DoG layers. The process generates ~ 1TB of data per county, including interim products. On the supercomputer, this takes 10 minutes. The next step is to interpret that curvature as multiscalar landforms.

Geological Society of America Abstracts with Programs. Vol. 50, No. 4

doi: 10.1130/abs/2018NC-313039

© Copyright 2018 The Geological Society of America (GSA), all rights reserved. Permission is hereby granted to the author(s) of this abstract to reproduce and distribute it freely, for noncommercial purposes. Permission is hereby granted to any individual scientist to download a single copy of this electronic file and reproduce up to 20 paper copies for noncommercial purposes advancing science and education, including classroom use, providing all reproductions include the complete content shown here, including the author information. All other forms of reproduction and/or transmittal are prohibited without written permission from GSA Copyright Permissions.