Title: A Hidden Markov Forest Model for Terrain-Aware Flood Inundation Mapping from Earth Imagery
Flood inundation mapping from Earth imagery plays a vital role in rapid disaster response and national water forecasting. However, the problem is non-trivial due to significant imagery noise and obstacles, complex spatial dependency on 3D terrains, spatial non-stationarity, and high computational cost. Existing machine learning approaches are mostly terrain-unaware and are prone to produce spurious results due to imagery noise and obstacles, requiring significant efforts in post-processing. Recently, several terrain- aware methods were proposed that incorporate complex spatial dependency (e.g., water flow directions on 3D terrains) but they assume that the inferred flood surface level is spatially stationary, making them insufficient for a large heterogeneous geographic area. To address these limitations, this paper proposes a novel spatial learning framework called hidden Markov forest, which decomposes a large heterogeneous area into local stationary zones, represents spatial dependency on 3D terrains via zonal trees (forest), and jointly infers the class map in different zonal trees with spatial regularization. We design efficient inference algorithms based on dynamic programming and multi-resolution filtering. Evaluations on real-world datasets show that our method outperforms baselines and our proposed computational refinement significantly reduces the time cost.  more » « less
Award ID(s):
2106461
PAR ID:
10437003
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the 2023 SIAM International Conference on Data Mining (SDM)
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Given earth imagery with spectral features on a terrain surface, this paper studies surface segmentation based on both explanatory features and surface topology. The problem is important in many spatial and spatiotemporal applications such as flood extent mapping in hydrology. The problem is uniquely challenging for several reasons: first, the size of earth imagery on a terrain surface is often much larger than the input of popular deep convolutional neural networks; second, there exists topological structure dependency between pixel classes on the surface, and such dependency can follow an unknown and non-linear distribution; third, there are often limited training labels. Existing methods for earth imagery segmentation often divide the imagery into patches and consider the elevation as an additional feature channel. These methods do not fully incorporate the spatial topological structural constraint within and across surface patches and thus often show poor results, especially when training labels are limited. Existing methods on semi-supervised and unsupervised learning for earth imagery often focus on learning representation without explicitly incorporating surface topology. In contrast, we propose a novel framework that explicitly models the topological skeleton of a terrain surface with a contour tree from computational topology, which is guided by the physical constraint (e.g., water flow direction on terrains). Our framework consists of two neural networks: a convolutional neural network (CNN) to learn spatial contextual features on a 2D image grid, and a graph neural network (GNN) to learn the statistical distribution of physics-guided spatial topological dependency on the contour tree. The two models are co-trained via variational EM. Evaluations on the real-world flood mapping datasets show that the proposed models outperform baseline methods in classification accuracy, especially when training labels are limited. 
    more » « less
  2. The accurate and prompt mapping of flood-affected regions is important for effective disaster management, including damage assessment and relief efforts. While high-resolution optical imagery from satellites during disasters presents an opportunity for automated flood inundation mapping, existing segmentation models face challenges due to noises such as cloud cover and tree canopies. Thanks to the digital elevation model (DEM) data readily available from sources such as United States Geological Survey (USGS), terrain guidance was utilized by recent graphical models such as hidden Markov trees (HMTs) to improve segmentation quality. Unfortunately, these methods either can only handle a small area where water levels at different locations are assumed to be consistent or require restricted assumptions such as there is only one river channel. This article presents an algorithm for flood extent mapping on large-scale Earth imagery, applicable to a large geographic area with multiple river channels. Since water level can vary a lot from upstream to downstream, we propose to detect river pixels to partition the remaining pixels into localized zones, each with a unique water level. In each zone, water at all locations flows to the same river entry point. Pixels in each zone are organized by an HMT to capture water flow directions guided by elevations. Moreover, a novel regularization scheme is designed to enforce inter-zone consistency by penalizing pixel-pairs of adjacent zones that violate terrain guidance. Efficient parallelization is made possible by coloring the zone adjacency graph to identify zones and zone-pairs that have no dependency and hence can be processed in parallel, and incremental one-pass terrain-guided scanning is conducted wherever applicable to reuse computations. Experiments demonstrate that our solution is more accurate than existing solutions and can efficiently and accurately map out flooding pixels in a giant area of size 24,805 × 40,129. Despite the imbalanced workloads caused by a few large zonal HMTs dominating the serial computing time, our parallelization approach is effective and manages to achieve up to 14.3× speedup on a machine with Intel Xeon Gold 6126 CPU @ 2.60 GHz (24 cores, 48 threads) using 32 threads. 
    more » « less
  3. Improving high-resolution (meter-scale) mapping of snow-covered areas in complex and forested terrains is critical to understanding the responses of species and water systems to climate change. Commercial high-resolution imagery from Planet Labs, Inc. (Planet, San Francisco, CA, USA) can be used in environmental science, as it has both high spatial (0.7–3.0 m) and temporal (1–2 day) resolution. Deriving snow-covered areas from Planet imagery using traditional radiometric techniques have limitations due to the lack of a shortwave infrared band that is needed to fully exploit the difference in reflectance to discriminate between snow and clouds. However, recent work demonstrated that snow cover area (SCA) can be successfully mapped using only the PlanetScope 4-band (Red, Green, Blue and NIR) reflectance products and a machine learning (ML) approach based on convolutional neural networks (CNN). To evaluate how additional features improve the existing model performance, we: (1) build on previous work to augment a CNN model with additional input data including vegetation metrics (Normalized Difference Vegetation Index) and DEM-derived metrics (elevation, slope and aspect) to improve SCA mapping in forested and open terrain, (2) evaluate the model performance at two geographically diverse sites (Gunnison, Colorado, USA and Engadin, Switzerland), and (3) evaluate the model performance over different land-cover types. The best augmented model used the Normalized Difference Vegetation Index (NDVI) along with visible (red, green, and blue) and NIR bands, with an F-score of 0.89 (Gunnison) and 0.93 (Engadin) and was found to be 4% and 2% better than when using canopy height- and terrain-derived measures at Gunnison, respectively. The NDVI-based model improves not only upon the original band-only model’s ability to detect snow in forests, but also across other various land-cover types (gaps and canopy edges). We examined the model’s performance in forested areas using three forest canopy quantification metrics and found that augmented models can better identify snow in canopy edges and open areas but still underpredict snow cover under forest canopies. While the new features improve model performance over band-only options, the models still have challenges identifying the snow under trees in dense forests, with performance varying as a function of the geographic area. The improved high-resolution snow maps in forested environments can support studies involving climate change effects on mountain ecosystems and evaluations of hydrological impacts in snow-dominated river basins. 
    more » « less
  4. The large-scale deployment of wind energy underscores the critical need for accurate resource characterization to reduce uncertainty in power estimates and to enable the installation of wind farms in increasingly complex terrains. Accurate wind resource assessment in peri-urban and moderately complex terrains remains a significant challenge due to spatial heterogeneity in surface terrain features and atmospheric thermal stability. This study investigates the influence of surface complexity and atmospheric stratification on vertical wind profiles at a utility-scale wind turbine site in Cedar Rapids, Iowa. One year of multi-level wind data from a 106-meter-tall meteorological tower were analyzed to quantify variations in the wind shear exponent α, wind direction veer, and horizontal turbulence intensity (TI) across open-field and complex-surface wind sectors and four thermal stability classes, defined by the bulk Richardson number Rib. The results show that the wind shear exponent α increases systematically with atmospheric stability. Over the open-field terrain, α ranges from 0.11 in unstable conditions to 0.45 in strongly stable conditions, compared to 0.17 and 0.40 over the complex surface. A pronounced diurnal variation in α was observed, particularly during the summer months. Wind veer was greatest and exceeded 30° under strongly stable conditions over open terrain. Elevated TI values peaked at 32 m in height due to flow separation and wake turbulence from nearby vegetation and sloping terrain. These findings highlight the importance of incorporating terrain-induced and thermally driven variability into wind resource assessments to improve power prediction and turbine siting in complex heterogeneous terrain environments. 
    more » « less
  5. The remote-sensing based Flood Crop Loss Assessment Service System (RF-CLASS) is a web service based system developed and managed by the Center for Spatial Information Science and Systems (CSISS). The system uses Moderate Resolution Imaging Spectroradiometer (MODIS)-based flood data, which was implemented by the Dartmouth Flood Observatory (DFO), to provide an estimation of crop loss from floods. However, due to the spectral similarity between water and shadow, a noticeable amount of false classification of shadow can be found in the DFO flood products. Traditional methods can be utilized to remove cloud shadow and part of mountain shadow. This paper aims to develop an algorithm to filter out noise from permanent mountain shadow in the flood layer. The result indicates that mountain shadow was significantly removed by using the proposed approach. In addition, the gold standard test indicated a small number of actual water surfaces were misidentified by the proposed algorithm. Furthermore, experiments also suggest that increasing the spatial resolution of the slope helped reduce more noise in mountains. The proposed algorithm achieved acceptable overall accuracy (>80%) in all different filters and higher overall accuracies were observed when using lower slope filters. This research is one of the very first discussions on identifying false flood classification from terrain shadow by using the highly efficient method. 
    more » « less