skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Accuracy of Structure-from-Motion/Multiview Stereo Terrain Models: A Practical Assessment for Applications in Field Geology
We assess the accuracy of Structure-from-Motion/Multiview stereo (SM) terrain models acquired ad hoc or without high-resolution ground control to analyze their usage as a base for inexpensive 3D bedrock geologic mapping. Our focus is on techniques that can be utilized in field projects without the use of heavy and/or expensive equipment or the placement of ground control in logistically challenging sites (e.g., steep cliff faces or remote settings). We use a Terrestrial Light Detection and Ranging (LiDAR) survey as a basis for the comparison of two types of SM models: (1) models developed from images acquired in a chartered airplane flight with ground control referenced by natural objects located on Google Earth scenes; and (2) drone flights with a georeference established solely from camera positions located by conventional, differentially corrected Global Navigation Satellite systems (GNSS). We find that all our SM models are indistinguishable in scale from the LiDAR reference model. The SM models do, however, show rigid body translations and rotations, with translations generally within the 1–5 m size of the natural objects used for ground control, the resolution of the GNSS receivers, or both. The rigid body rotations can be attributed to a poor imaging plan, which can be avoided with survey planning. Analyses of point densities in various models show a limitation of Terrestrial LiDAR point clouds as a mapping base due to the rapid falloff of resolution with distance. In contrast, SM models are characterized by relatively uniform point densities controlled by camera optics, the numbers of images, and the distance from the target. This uniform density is the product of the Multiview stereo step in SM processing that fills areas between key points and is important for bedrock geologic mapping because it affords direct interpretation on a point cloud at a relatively uniform scale throughout a model. Our results indicate that these simple methods allow SM model construction to be accurate to the range of conventional GNSS with resolutions to the submeter, even cm, scale depending on data acquisition parameters. Thus, SM models can, and should, serve as a base for high-resolution geologic mapping, particularly in a steep terrain where conventional techniques fail. Our SM models appear to provide accurate visualizations of geologic features over km scales that allow detailed geologic mapping in 3D with a relative accuracy to the decimeter or centimeter level and absolute positioning in the 2–5 m precision of GNSS; a geometric precision that will allow unprecedented new studies of any geologic system where geometry is the fundamental data.  more » « less
Award ID(s):
2049603
PAR ID:
10502768
Author(s) / Creator(s):
;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Geosciences
Volume:
13
Issue:
7
ISSN:
2076-3263
Page Range / eLocation ID:
217
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We use structure from motion–multiview stereo (SM) terrain models developed from ground-based images and images acquired from uncrewed aircraft (aka drones) as a base map for three-dimensional (3-D) mapping on the walls of a deep canyon in the Panamint Range of eastern California, USA. The ability to manipulate the 3-D model with views from arbitrary look directions and broad scale range revealed structures that were invisible to conventional two-dimensional (2-D) mapping because of both the scale of the structures and their exposure on vertical to near-vertical cliff faces. The analysis supports field evidence for four phases of ductile deformation, with only one of the younger phases documented on early geologic maps of the area. The oldest deformational event (D1) produced the main metamorphic fabric and pre-dates Late Cretaceous plutons. This deformation produced a 200–250-m-thick high-strain zone localized along marbles at the top of the Kingston Peak Formation and lower Noonday Formation. Geometric analysis from the model suggests strongly that large sheath folds at scales of 100–300 m are developed within these marbles. Large measured finite strains indicate displacement across this apparent shear zone of at least 4–5 km and displacements of tens of kilometers are allowable, yet the structure is invisible to conventional mapping because the high-strain zone is stratabound. The main fabric shows two clear overprints and a third that is likely an even younger deformation. D2 and D3 generated tight to close, recumbent folds and open to tight, upright folds, respectively, both folding the main foliation with localized development of crenulation cleavages axial planar to the folds. An additional overprint shows no clear cross-cutting relationship with D2 or D3 fabrics and could be a manifestation of either of those events, although the deformation is spatially limited to a narrow shear zone beneath a brittle, dextral-normal fault with the same kinematics as a mylonitic fabric in a Cretaceous granite in the footwall. This observation suggests an extensional, core complex–style deformation to produce this structure. We suggest that 3-D mapping has the potential to revolutionize geologic mapping studies, particularly where steep topography provides 3-D views that are virtually invisible on conventional 2-D maps. Previously bewildering geologic puzzles can be solved by the ability to visualize large cliff exposures from arbitrary angles and map the features in true 3-D at resolutions to the centimeter level. Although this study emphasized intermediate scales imaged by a drone, our methods here are easily extended to larger scales using a crewed aircraft for imaging. We suggest these methods should be used routinely in frontier areas with steep terrain where aviation is already in use for access, but the methods can be employed anywhere steep terrain “hides” major rock exposures on conventional 2-D maps. 
    more » « less
  2. The University of Alaska Fairbanks T Field is a legacy farm field that is part of the National Science Foundation (NSF) Funded Permafrost Grown project. We are studying the long-term effects of permafrost thaw following initial clearing for cultivation purposes. In this regard, we have acquired very high resolution light detection and ranging (LiDAR) data and digital photography from a DJI M300 drone using a Zenmuse L1 and a MicaSense RedEdge-P camera. The Zenmuse L1 integrates a Livox Lidar module, a high-accuracy inertial measurement units (IMU), and a camera with a 1-inch CMOS on a 3-axis stabilized gimbal. The MicaSense RedEdge-P camera has five multispectral bands and a high-resolution panchromatic band. The drone was configured to fly in real-time kinematic (RTK) mode at an altitude of 60 meters above ground level using the DJI D-RTK 2 base station. Data was acquired using a 50% sidelap and a 70% frontlap for the Zenmuse L1 and an 80% sidelap and a 75% frontlap for the MicaSense. Additional ground control was established with a Leica GS18 global navigation satellite system (GNSS) and all data have been post-processed to World Geodetic System 1984 (WGS84) universal transverse mercator (UTM) Zone 6 North using ellipsoid heights. Data outputs include a two-class-classified LiDAR point cloud, digital surface model, digital terrain model, an orthophoto mosaic, and a multispectral orthoimage consisting of five bands. Image acquisition occurred on 18 August 2023. 
    more » « less
  3. The Midnight Sun Golf Course in Fairbanks, Alaska is a legacy farm field that is part of the National Science Foundation (NSF) Funded Permafrost Grown project. This 65 hectare (ha) parcel was initially cleared for agriculture purposes but changed land-use practices to a golf course around 25 years ago. The land-use conversion was in part due to ice-rich permafrost thaw following clearing. We are studying the long-term effects of permafrost thaw following initial clearing for cultivation purposes. We are working with the current landowners to provide information regarding ongoing thermokarst development on the property and to conduct studies in reforested portions of the land area to understand land clearing and reforestation on permafrost-affected soils. In this regard, we have acquired very high resolution light detection and ranging (LiDAR) data and digital photography from a DJI M300 drone using a Zenmuse L1. The Zenmuse L1 integrates a Livox Lidar module, a high-accuracy inertial measurement units (IMU), and a camera with a 1-inch CMOS on a 3-axis stabilized gimbal. The drone was configured to fly in real-time kinematic (RTK) mode at an altitude of 60 meters above ground level using the DJI D-RTK 2 base station. Data was acquired using a 50% sidelap and a 70% frontlap. Additional ground control was established with a Leica GS18 global navigation satellite system (GNSS) and all data have been post-processed to World Geodetic System 1984 (WGS84) universal transverse mercator (UTM) Zone 6 North using ellipsoid heights. Data outputs include a two-class classified LiDAR point cloud, digital surface model, digital terrain model, and an orthophoto mosaic. Image acquisition occurred on 10 September 2023. The input images are available for download at http://arcticdata.io/data/10.18739/A2PC2TB1T. 
    more » « less
  4. Abstract. Here we describe the curriculum and outcomes from a data-intensivegeomorphic analysis course, “Geoscience Field Issues Using High-ResolutionTopography to Understand Earth Surface Processes”, which pivoted to virtualin 2020 due to the COVID-19 pandemic. The curriculum covers technologies formanual and remotely sensed topographic data methods, including (1) GlobalPositioning Systems and Global Navigation Satellite System (GPS/GNSS)surveys, (2) Structure from Motion (SfM) photogrammetry, and (3) ground-based(terrestrial laser scanning, TLS) and airborne lidar. Course content focuseson Earth-surface process applications but could be adapted for othergeoscience disciplines. Many other field courses were canceled in summer2020, so this course served a broad range of undergraduate and graduatestudents in need of a field course as part of degree or researchrequirements. Resulting curricular materials are available freely within theNational Association of Geoscience Teachers' (NAGT's) “Teaching with Online Field Experiences” collection. Theauthors pre-collected GNSS data, uncrewed-aerial-system-derived (UAS-derived) photographs, and ground-based lidar, which students then used in courseassignments. The course was run over a 2-week period and had synchronousand asynchronous components. Students created SfM models that incorporatedpost-processed GNSS ground control points and created derivative SfM and TLSproducts, including classified point clouds and digital elevation models(DEMs). Students were successfully able to (1) evaluate the appropriatenessof a given survey/data approach given site conditions, (2) assess pros andcons of different data collection and post-processing methods in light offield and time constraints and limitations of each, (3) conduct error andgeomorphic change analysis, and (4) propose or implement a protocol to answera geomorphic question. Overall, our analysis indicates the course had asuccessful implementation that met student needs as well as course-specificand NAGT learning outcomes, with 91 % of students receiving an A, B, or Cgrade. Unexpected outcomes of the course included student self-reflectionand redirection and classmate support through a daily reflection anddiscussion post. Challenges included long hours in front of a computer,computing limitations, and burnout because of the condensed nature of thecourse. Recommended implementation improvements include spreading the courseout over a longer period of time or adopting only part of the course andproviding appropriate computers and technical assistance. This paperand published curricular materials should serve as an implementation andassessment guide for the geoscience community to use in virtual or in-personhigh-resolution topographic data courses that can be adapted for individuallabs or for an entire field or data course. 
    more » « less
  5. Abstract Seismic hazard assessment, such as the U.S. Geological Survey (USGS) National Seismic Hazard Model (NSHM), relies on estimates of fault slip rate based on geology and/or geodetic observations such as the Global Navigation Satellite System (GNSS), including the Global Positioning System. Geodetic fault slip rates may be estimated within a 3D spherical block model, in which the crust is divided into microplates bounded by mapped faults; fault slip rates are determined by the relative rotations of adjacent microplates. Uncertainty in selecting appropriate block-bounding faults and in forming closed microplates has limited the interpretability of block models for seismic hazard modeling. By introducing an automated block closure algorithm and regularizing the resulting densely spaced block model with total variation regularization, I develop the densest and most complete block model of the western continental United States to date. The model includes 853 blocks bounded by 1017 geologically identified fault sections from the USGS NSHM Fault Sections database. Microplate rotations and fault slip rates are constrained by 4979 GNSS velocities and 1243 geologic slip rates. I identify a regularized solution that fits the GNSS velocity field with a root mean square misfit of 1.9 mm/yr and reproduces 57% of geologic slip rates within reported geologic uncertainty and model sensitivity, consistent with other geodetic-based models in this Focus Section. This block model includes slip on faults that are not included in the USGS NSHM Fault sections database (but are required to form closed blocks) for an estimate of “off-fault” deformation of 3.62×1019  N·m/yr, 56% of the total calculated moment accumulation rate in the model. 
    more » « less