Probabilistic change mapping from airborne LiDAR for post-disaster damage assessment
Abstract
When both pre- and post-event LiDAR point clouds are available, change detection can be performed to identify areas that were most affected by a disaster event, and to obtain a map of quantitative changes in terms of height differences. In the case of earthquakes in built-up areas for instance, first responders can use a LiDAR change map to help prioritize search and recovery efforts. The main challenge consists of producing reliable change maps, robust to collection conditions, free of processing artifacts (due for instance to triangulation or gridding), and taking into account the various sources of uncertainty. Indeed, datasets acquired within a few years interval are often of different point density (sometimes an order of magnitude higher for recent data), different acquisition geometries, and very likely suffer from georeferencing errors and geometric discrepancies. All these differences might not be important for producing maps from each dataset separately, but they are crucial when performing change detection. We have developed a novel technique for the estimation of uncertainty maps from the LiDAR point clouds, using Bayesian inference, treating all variables as random. The main principle is to grid all points on a common grid before attempting any comparison, as working directly with point clouds is cumbersome and time consuming. A non-parametric approach based on local linear regression was implemented, assuming a locally linear model for the surface. This enabled us to derive error bars on gridded elevations, and then elevation differences. In this way, a map of statistically significant changes could be computed - whereas a deterministic approach would not allow testing of the significance of differences between the two datasets. This approach allowed us to take into account not only the observation noise (due to ranging, position and attitude errors) but also the intrinsic roughness of the observed surfaces occurring when scanning vegetation. As only elevation differences above a predefined noise level are accounted for (according to a specified confidence interval related to the allowable false alarm rate) the change detection is robust to all these sources of noise. To first validate the approach, we built small-scale models and scanned them using a terrestrial laser scanner to establish 'ground truth'. Changes were manually applied to the models then new scans were performed and analyzed. Additionally, two airborne datasets of the Monterey Peninsula, California, were processed and analyzed. The first one was acquired during 2010 (with relatively low point density, 1-3 pts/m2), and the second one was acquired during 2012 (with up to 30 pts/m2). To perform the comparison, a new point cloud registration technique was developed and the data were registered to a common 1 m grid. The goal was to correct systematic shifts due to GPS and INS errors, and focus on the actual height differences regardless of the absolute planimetric accuracy of the datasets. Though no major disaster event occurred between the two acquisition dates, sparse changes were detected and interpreted mostly as construction and natural landscape evolution.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2013
- Bibcode:
- 2013AGUFM.G33A0982J
- Keywords:
-
- 0520 COMPUTATIONAL GEOPHYSICS Data analysis: algorithms and implementation;
- 1986 INFORMATICS Statistical methods: Inferential;
- 4324 NATURAL HAZARDS Spatial decision support systems;
- 4346 NATURAL HAZARDS Emergency response and evacuations