Where we live – A summary of the achievements and planned evolution of the Global Urban Footprint

The TerraSAR-X (TSX) mission provides a distinguished collection of high resolution satellite images that shows great promise for a global monitoring of human settlements. Hence, the German Aerospace Center (DLR) has developed the Urban Footprint Processor (UFP) that represents an operational framework for the mapping of built-up areas based on a mass processing and analysis of TSX imagery. The UFP includes functionalities for data management, feature extraction, unsupervised classification, mosaicking, and post-editing. Based on >180.000 TSX StripMap scenes, the UFP was used in 2016 to derive a global map of human presence on Earth in a so far unique spatial resolution of 12 m per grid cell: the Global Urban Footprint (GUF). This work provides a comprehensive summary of the major achievements related to the Global Urban Footprint initiative, with dedicated sections focusing on aspects such as UFP methodology, basic product characteristics (specification, accuracy, global figures on urbanization derived from GUF), the user community, and the already initiated future roadmap of follow-on activities and products. The active community of >250 institutions already working with the GUF data documents the relevance and suitability of the GUF initiative and the underlying high-resolution SAR imagery with respect to the provision of key information on the human presence on earth and the global human settlements properties and patterns, respectively

Digital elevation model correction in urban areas using extreme gradient boosting, land cover and terrain parameters

The accuracy of digital elevation models (DEMs) in urban areas is influenced by numerous factors including land cover and terrain irregularities. Moreover, building artifacts in global DEMs cause artificial blocking of surface flow pathways. This compromises their quality and adequacy for hydrological and environmental modelling in urban landscapes where precise and accurate terrain information is needed. In this study, the extreme gradient boosting (XGBoost) ensemble algorithm is adopted for enhancing the accuracy of two medium-resolution 30m DEMs over Cape Town, South Africa: Copernicus GLO-30 and ALOS World 3D (AW3D). XGBoost is a scalable, portable and versatile gradient boosting library that can solve many environmental modelling problems. The training datasets are comprised of eleven predictor variables including elevation, urban footprints, slope, aspect, surface roughness, topographic position index, terrain ruggedness index, terrain surface texture, vector roughness measure, forest cover and bare ground cover. The target variable (elevation error) was calculated with respect to highly accurate airborne LiDAR. After training and testing, the model was applied for correcting the DEMs at two implementation sites. The correction achieved significant accuracy gains which are competitive with other proposed methods. The root mean square error (RMSE) of Copernicus DEM improved by 46 to 53% while the RMSE of AW3D DEM improved by 72 to 73%. These results showcase the potential of gradient boosted trees for enhancing the quality of DEMs, and for improved hydrological modelling in urban catchments.Comment: 8 page

Rapid Damage Mapping for the 2015 M_w 7.8 Gorkha Earthquake Using Synthetic Aperture Radar Data from COSMO–SkyMed and ALOS-2 Satellites

The 25 April 2015 M_w 7.8 Gorkha earthquake caused more than 8000 fatalities and widespread building damage in central Nepal. The Italian Space Agency’s COSMO–SkyMed Synthetic Aperture Radar (SAR) satellite acquired data over Kathmandu area four days after the earthquake and the Japan Aerospace Exploration Agency’s Advanced Land Observing Satellite-2 SAR satellite for larger area nine days after the mainshock. We used these radar observations and rapidly produced damage proxy maps (DPMs) derived from temporal changes in Interferometric SAR coherence. Our DPMs were qualitatively validated through comparison with independent damage analyses by the National Geospatial-Intelligence Agency and the United Nations Institute for Training and Research’s United Nations Operational Satellite Applications Programme, and based on our own visual inspection of DigitalGlobe’s WorldView optical pre- versus postevent imagery. Our maps were quickly released to responding agencies and the public, and used for damage assessment, determining inspection/imaging priorities, and reconnaissance fieldwork

Urban morphology parameters from global digital elevation models: implications for aerodynamic roughness for wind-speed estimation

Urban morphology and aerodynamic roughness parameters are derived from three global digital elevation models (GDEM): Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), Shuttle Radar Topography Mission (SRTM), and TanDEM-X. Initially, each is compared to benchmark elevation data in London (UK). A moving window extracts ground heights from the GDEMs, generating terrain models with root-mean-square accuracy of up to 3 m. Subtraction of extracted ground heights provides roughness-element heights only, allowing for calculation of morphology parameters. The parameters are calculated for eight directional sectors of 1 km grid-squares. Apparent merging of roughness elements in all GDEMs causes height-based parameter underestimation, whilst plan and frontal areas are over- and under-estimated, respectively. Combined, these lead to an underestimation of morphometrically-derived aerodynamic roughness parameters. Parameter errors are least for the TanDEM-X data. Further comparison in five cities (Auckland, Greater London, New York, Sao Paulo, Tokyo) provides basis for empirical corrections to TanDEM-X-derived geometric parameters. These reduce the error in parameters across the cities and for a separate location. Meteorological observations in central London give insight to wind-speed estimation accuracy using roughness parameters from the different elevation databases. The proposed corrections to TanDEM-X parameters lead to improved wind-speed estimates, which combined with the improved spatial representation of parameters across cities demonstrates their potential for use in future studies

Mapping Local Climate Zones for a Worldwide Database of the Form and Function of Cities

Progress in urban climate science is severely restricted by the lack of useful information that describes aspects of the form and function of cities at a detailed spatial resolution. To overcome this shortcoming we are initiating an international effort to develop the World Urban Database and Access Portal Tools (WUDAPT) to gather and disseminate this information in a consistent manner for urban areas worldwide. The first step in developing WUDAPT is a description of cities based on the Local Climate Zone (LCZ) scheme, which classifies natural and urban landscapes into categories based on climate-relevant surface properties. This methodology provides a culturally-neutral framework for collecting information about the internal physical structure of cities. Moreover, studies have shown that remote sensing data can be used for supervised LCZ mapping. Mapping of LCZs is complicated because similar LCZs in different regions have dissimilar spectral properties due to differences in vegetation, building materials and other variations in cultural and physical environmental factors. The WUDAPT protocol developed here provides an easy to understand workflow; uses freely available data and software; and can be applied by someone without specialist knowledge in spatial analysis or urban climate science. The paper also provides an example use of the WUDAPT project results

Large-scale 3D Modelling of the Built Environment - Joint Analysis of TanDEM-X, Sentinel-2 and Open Street Map Data

Continental to global scale mapping of the human settlement extent based on earth observation satellite data has made considerable progress. Nevertheless, the current approaches only provide a two-dimensional representation of the built environment. Therewith, a full characterization is restricted in terms of the urban morphology and built-up density, which can only be gained by a detailed examination of the vertical settlement extent. This paper introduces a methodology for the extraction of three-dimensional (3D) information on human settlements by analyzing the digital elevation and radar intensity data collected by the German TanDEM-X satellite mission in combination with multispectral Sentinel-2 imagery and data from the Open Street Map initiative and the Global Urban Footprint human settlement mask. The first module of the underlying processor generates a normalized digital surface model from the TanDEM-X digital elevation model for all regions marked as a built-up area by the Global Urban Footprint. The second module generates a building mask based on a joint processing of Open Street Map, TanDEM-X/TerraSAR-X radar images, the calculated normalized digital surface model and Sentinel-2 imagery. Finally, a third module allocates the local relative heights of the normalized digital surface model to the building structures provided by the building mask. The outcome of the procedure is a 3D map of the built environment showing the estimated local height for all identified vertical building structures at 12 m spatial resolution. The results of a first validation campaign based on reference data collected for the seven cities of Amsterdam (NL), Indianapolis (US), Kigali (RW), Munich (DE), New York (US), Vienna (AT), and Washington (US) indicate the potential of the proposed methodology to accurately estimate the distribution of building heights within the built-up area

Towards a 20m global building map from Sentinel-1 SAR Data

This study introduces a technique for automatically mapping built-up areas using synthetic aperture radar (SAR) backscattering intensity and interferometric multi-temporal coherence generated from Sentinel-1 data in the framework of the Copernicus program. The underlying hypothesis is that, in SAR images, built-up areas exhibit very high backscattering values that are coherent in time. Several particular characteristics of the Sentinel-1 satellite mission are put to good use, such as its high revisit time, the availability of dual-polarized data, and its small orbital tube. The newly developed algorithm is based on an adaptive parametric thresholding that first identifies pixels with high backscattering values in both VV and VH polarimetric channels. The interferometric SAR coherence is then used to reduce false alarms. These are caused by land cover classes (other than buildings) that are characterized by high backscattering values that are not coherent in time (e.g., certain types of vegetated areas). The algorithm was tested on Sentinel-1 Interferometric Wide Swath data from five different test sites located in semiarid and arid regions in the Mediterranean region and Northern Africa. The resulting building maps were compared with the Global Urban Footprint (GUF) derived from the TerraSAR-X mission data and, on average, a 92% agreement was obtained.Peer ReviewedPostprint (published version