Abstract Sentinel-2 multi-spectral images collected over periods of several months were used to estimate vegetation height for Gabon and Switzerland. A deep convolutional neural network (CNN) was trained to extract suitable spectral and textural features from reflectance images and to regress per-pixel vegetation height. In Gabon, reference heights for training and validation were derived from airborne LiDAR measurements. In Switzerland, reference heights were taken from an existing canopy height model derived via photogrammetric surface reconstruction. The resulting maps have a mean absolute error (MAE) of 1.7 m in Switzerland and 4.3 m in Gabon (a root mean square error (RMSE) of 3.4 m and 5.6 m, respectively), and correctly estimate vegetation heights up to >50 m. They also show good qualitative agreement with existing vegetation height maps. Our work demonstrates that, given a moderate amount of reference data (i.e., 2000 km2 in Gabon and ≈5800 km2 in Switzerland), high-resolution vegetation height maps with 10 m ground sampling distance (GSD) can be derived at country scale from Sentinel-2 imagery.
[1] Data from spaceborne light detection and ranging (lidar) opens the possibility to map forest vertical structure globally. We present a wall-to-wall, global map of canopy height at 1-km spatial resolution, using 2005 data from the Geoscience Laser Altimeter System (GLAS) aboard ICESat (Ice, Cloud, and land Elevation Satellite). A challenge in the use of GLAS data for global vegetation studies is the sparse coverage of lidar shots (mean = 121 data points degree2 for the L3C campaign). However, GLAS-derived canopy height (RH100) values were highly correlated with other, more spatially dense, ancillary variables available globally, which allowed us to model global RH100 from forest type, tree cover, elevation, and climatology maps. The difference between the model predicted RH100 and footprint level lidar-derived RH100 values showed that error increased in closed broadleaved forests such as the Amazon, underscoring the challenges in mapping tall (>40 m) canopies. The resulting map was validated with field measurements from 66 FLUXNET sites. The modeled RH100 versus in situ canopy height error (RMSE = 6.1 m, R2 = 0.5; or, RMSE = 4.4 m, R2 = 0.7 without 7 outliers) is conservative as it also includes measurement uncertainty and sub pixel variability within the 1-km pixels. Our results were compared against a recently published canopy height map. We found our values to be in general taller and more strongly correlated with FLUXNET data. Our map reveals a global latitudinal gradient in canopy height, increasing towards the equator, as well as coarse forest disturbance patterns.