DATASET

GHS-BUILT-S R2022A - GHS built-up surface grid, derived from Sentinel2 composite and Landsat, multitemporal (1975-2030)

Collection: GHSL : Global Human Settlement Layer 

Description

The spatial raster dataset depicts the distribution of the built-up (BU) surfaces estimates between 1975 and 2030 in 5 years intervals and two functional use components a) the total BU surface and b) the non-residential (NRES) BU surface. The data is made by spatial-temporal interpolation of five observed collections of multiple-sensor, multiple-platform satellite imageries. Landsat (MSS, TM, ETM sensor) supports the 1975, 1990, 2000, and 2014 epochs. Sentinel2 (S2) composite (GHS-composite-S2 R2020A) supports the 2018 epoch.

The built-up surface fraction (BUFRAC) is estimated at 10m of spatial resolution from the S2 image data, using as learning set a composite of data from GHS-BUILT-S2 R2020A, Facebook, Microsoft, and Open Street Map (OSM) building delineation. The BUFRAC inference is made from the combination of quantized image features (reflectance, derivative of morphological profile DMP) through associative rule learning applied to spatial data analytics, which was introduced as symbolic machine learning (SML). The non-residential (NRES) domain is predicted from S2 image data by observation of radiometric, textural, and morphological features in an object-oriented image processing framework. The multi-temporal dimension is provided by testing by the SML the association between the combination of the quantized radiometric information collected by the Landsat imagery in the past epochs, and the “built-up” (BU) and “non-built-up” (NBU) class abstraction on image segments extracted from S2 images. The spatial-temporal interpolation is solved by rank-optimal spatial allocation using explanatory variables related to the landscape (slope, elevation, distance to water, and distance to vegetation) and related to the observed dynamic of BU surfaces in the past epochs.

Contact

Email
jrc-ghsl-data (at) ec.europa.eu

Contributors

How to cite

Pesaresi, Martino; Politis, Panagiotis (2022): GHS-BUILT-S R2022A - GHS built-up surface grid, derived from Sentinel2 composite and Landsat, multitemporal (1975-2030). European Commission, Joint Research Centre (JRC) [Dataset] doi: 10.2905/D07D81B4-7680-4D28-B896-583745C27085 PID: http://data.europa.eu/89h/d07d81b4-7680-4d28-b896-583745c27085

Keywords

Built-up surface grid GHS BUILT SURFACE GHS-BUILT-S GHSL global map Functional use classification

Data access

GHS-BUILT-S_NRES_GLOBE_R2022A
URL 
  • GHS non-residential (NRES) built-up surface grids derived from joint assessment of Sentinel2 and Landsat sensors, for 1975-2030 (5yrs interval). Values are expressed in unsigned integers and report about the predicted amount of built-up square meters. The data is published at 100m and 1km resolution in World Mollweide (EPSG:54009). For the year 2018 the data is published at 10m as observed from the S2 image data

GHS-BUILT-S_GLOBE_R2022A
URL 
  • GHS built-up surface grids derived from joint assessment of Sentinel2 and Landsat sensors, for 1975-2030 (5yrs interval). Values are expressed in unsigned integers and report about the predicted amount of built-up square meters. The data is published at 100m and 1km resolution in World Mollweide (EPSG:54009). For the year 2018 the data is published at 10m as observed from the S2 image data.

GHSL website
URL 
  • Project Web site

Publications

Publication
Differential area profiles: decomposition properties and efficient computation
Ouzounis G, Pesaresi M, Soille P. Differential area profiles: decomposition properties and efficient computation. IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 34 (8); 2012. p. 1533-1548. JRC59388
  • IEEE COMPUTER SOC, LOS ALAMITOS, USA
  • Abstract

    Differential area profiles (DAPs) are point-based multiscale descriptors used in pattern analysis and image segmentation. They are defined through sets of size-based connected morphological filters that constitute a joint area opening top-hat and area closing bottom-hat scale-space of the input image. The work presented in this paper explores the properties of this image decomposition through sets of area zones. An area zone defines a single plane of the DAP vector field and contains all the peak components of the input image, whose size is between the zone’s attribute extrema. Area zones can be computed efficiently from hierarchical image representation structures, in a way similar to regular attribute filters. Operations on the DAP vector field can then be computed without the need for exporting it first, and an example with the leveling-like convex/concave segmentation scheme is given. This is referred to as the one-pass method and it is demonstrated on the Max-Tree structure. Its computational performance is tested and compared against conventional means for computing differential profiles, relying on iterative application of area openings and closings. Applications making use of the area zone decomposition are demonstrated in problems related to remote sensing and medical image analysis.

Publication
A new compact representation of morphological profiles: report on first massive VHR image processing at the JRC
Pesaresi M, Ouzounis G, Gueguen L. A new compact representation of morphological profiles: report on first massive VHR image processing at the JRC. In Conference Proceedings: Sylvia S. Shen, Paul E. Lewis, editors. Algorithms and Technologies for Multispectral, Hyperspectral, and Ultraspectral Imagery XVIII. Vol. 8390. SPIE; 2012. JRC70542
  • SPIE, U.S.A
  • Abstract

    A new compact representation of differential morphological profile (DMP) vector fields is presented. It is referred to as the CSL model and is conceived to radically reduce the dimensionality of the DMP descriptors. The model maps three characteristic parameters, namely scale, saliency and level, into the RGB space through a HSV transform. The result is a a medium abstraction semantic layer used for visual exploration, image information mining and pattern classification. Fused with the PANTEX built-up presence index, the CSL model converges to an approximate building footprint representation layer in which color represents building class labels. This process is demonstrated on the first high resolution (HR) global human settlement layer (GHSL) computed from multi-modal HR and VHR satellite images. Results of the first massive processing exercise involving several thousands of scenes around the globe are reported along with validation figures.

Publication
A New Approach for the Morphological Segmentation of High-Resolution Satellite Imagery.
Pesaresi M, Benediktsson J. A New Approach for the Morphological Segmentation of High-Resolution Satellite Imagery.. IEEE Transactions on Geoscience and RS 39 (2); 2001. JRC19264
URL 
  • Abstract

    Abstract not available

Publication
Assessment of the Added-Value of Sentinel-2 for Detecting Built-up Areas
Pesaresi M, Corban C, Julea A, Florczyk A, Syrris V, Soille P. Assessment of the Added-Value of Sentinel-2 for Detecting Built-up Areas. REMOTE SENSING 8 (4); 2016. p. 299. JRC99996
  • MDPI AG, BASEL, SWITZERLAND
  • Abstract

    Monitoring of the human-induced changes and the availability of reliable and methodologically consistent urban area maps are essential to support sustainable urban development on a global scale. The Global Human Settlement Layer (GHSL) is a project funded by the European Commission, Joint Research Centre, which aims at providing scientific methods and systems for reliable and automatic mapping of built-up areas from remote sensing data. In the frame of the GHSL, the opportunities offered by the recent availability of Sentinel-2 data are being explored using a novel image classification method, called Symbolic Machine Learning (SML), for detailed urban land cover mapping. In this paper, a preliminary test was implemented with the purpose of: (i) assessing the applicability of the SML classifier on Sentinel-2 imagery; (ii) evaluating the complementarity of Sentinel-1 and Sentinel-2; and (iii) understanding the added-value of Sentinel-2 with respect to Landsat for improving global high-resolution human settlement mapping. The overall objective is to explore areas of improvement, including the possibility of synergistic use of the different sensors. The results showed that noticeable improvement of the quality of the classification could be gained from the increased spatial detail and from the thematic contents of Sentinel-2 compared to the Landsat derived product as well as from the complementarity between Sentinel-1 and Sentinel-2 images.

Publication
A New Method for Earth Observation Data Analytics Based on Symbolic Machine Learning
Pesaresi M; Syrris V; Julea A. A New Method for Earth Observation Data Analytics Based on Symbolic Machine Learning. REMOTE SENSING 8 (5); 2016. p. 399. JRC99747
  • MDPI AG, BASEL, SWITZERLAND
  • Abstract

    This work introduces a new classification method in the remote sensing domain, suitably adapted to dealing with the challenges posed by the big data processing and analytics framework. The method is based on symbolic learning techniques, and it is designed to work in complex and information-abundant environments, where relationships among different data layers are assessed in model-free and computationally-effective modalities. The two main stages of the method are the data reduction-sequencing and the association analysis. The former refers to data representation; the latter searches for systematic relationships between data instances derived from images and spatial information encoded in supervisory signals. Subsequently, a new measure named the evidence-based normalized differential index, inspired by the probability-based family of objective interestingness measures, evaluates these associations. Additional information about the computational complexity of the classification algorithm and some critical remarks are briefly introduced. An application of land cover mapping where the input image features are morphological and radiometric descriptors demonstrates the capacity of the method; in this instructive application, a subset of eight classes from the Corine Land Cover is used as the reference source to guide the training phase.

Publication
GLOBAL FINE-SCALE INFORMATION LAYERS: THE NEED OF A PARADIGM SHIFT
Pesaresi M. GLOBAL FINE-SCALE INFORMATION LAYERS: THE NEED OF A PARADIGM SHIFT. In: Soille P, Marchetti PG, editors. Proceedings of the 2014 conference on Big Data from Space (BiDS`14). Luxembourg (Luxembourg): Publications Office of the European Union; 2014. p. 8-11. JRC92345
  • Publications Office of the European Union, Luxembourg, Luxembourg
  • Abstract

    Global fine-scale information extraction tasks using today’s remote sensing technologies are compared with Big Data Analytics tasks. Issues related to data classification, machine learning and evaluation of the results are discussed. A change of paradigm respect to the classical RS information processing model is proposed in order to cope with the Big RS Data characteristics.

Publication
A Robust Built-up Area Presence Index by Anisotropic Rotation-invariant Textural Measure
Pesaresi M, Gerhardinger A, Kayitakire F. A Robust Built-up Area Presence Index by Anisotropic Rotation-invariant Textural Measure. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing 1 (3); 2008. p. 180-192. JRC37845
  • IEEE, Piscataway, United States of America
URL 
  • Abstract

    A procedure for the calculation of a texture-derived built-up presence index (PanTex) from textural characteristics of panchromatic satellite data is presented. The index is based on fuzzy rule-based composition of anisotropic textural co-occurrence measures derived from the satellite data by the gray-level co-occurrence matrix (GLCM). The strength and weakness of the procedure is analyzed and compared with traditional radiometric and textural approaches with the help of specific examples.

    A variant of the PanTex is proposed including information on vegetation if available (PanTexG).

    The accuracy and robustness of PanTex against seasonal changes, multi-sensor, multi-scene, and data degradation by wavelet-based compression and histogram stretching is discussed with some examples.

Publication
A New Built-Up Presence Index Based On Density of Corners
Gueguen L, Soille P, Pesaresi M. A New Built-Up Presence Index Based On Density of Corners. In Conference Proceedings: Geoscience and Remote Sensing Symposium (IGARSS), 2012 IEEE International. Piscataway (USA): IEEE; 2012. p. 5398-5401. JRC68582
  • IEEE , Piscataway, USA
  • Abstract

    The location and identification of human settlements is a key information in any assessment related to human security and safety decision process, from the preparedness to and mitigation of natural hazards to postdisaster response and reconstruction. The human settlements can be defined as infrastructures build on earth surface for living, trading or any social activity.

    In order to cover and analyze large areas where scarce and heterogeneous data may exist, the automatic information extraction from Earth Observation is the only sustainable strategy. Several studies about the detection of built-up presence from high resolution satellite were performed during the past few years. Among the various methods, a built-up index aggregating anisotropic textural co-occurence measures (PanTex) was demonstrated to be a robust indicator, where 54 high resolution optical images representing cities distributed around the world are processed. The validation of the built-up presence index emphasized its robustness against the type of constructions, the sun angles and the seasons of acquisitions. In this paper, we propose an alternative procedure for the calculation of built-up presence index from panchromatic high resolution satellite images. While radiometric descriptors of human settlements are highly variable across the world and with illumination conditions, the shape of settlements often contains right angles. This property being more stable, we propose a built-up index based on the density of corners. The index is obtained in two steps: - corners are detected by multi scale Harris detector based on differential morphological decompositions; - corners are spatially aggregated to form a density of corners, which is the built-up index. The differential morphological decomposition is a scale-space representation of the image, where image elements are separated by their scales. Then, the Harris corner indicator, which is highly dependent on a scale parameter, can be adapted to a set of scales. The output of the corner detection is a set of points associated to a scale, which represents the right angles in the image. The corners associated to a scale in the range of admissible settlements dimensions are selected and spatially aggregated to derive a density of corners.

    The proposed index is extracted from various high resolution panchromatic images and it is compared to the PanTex. In the experimental section, the high correlation between both indicators proves the suitability of the proposed method for consistently detecting built-up presence. Moreover, it gives a new interpretation of the PanTex which is close to a density of corners. Such an observation is critical for understanding the variablity of the PanTex index in between dense urban and industrial areas (industrial areas being composed of large buildings mechanically contain less corners).

Geographic areas

World

Spatial coverage

Type Value
WKT
POLYGON((-180 90,180 90,180 -90,-180 -90,-180 90))

Temporal coverage

From date To date
1975-01-01 2030-12-31

Additional information

Published by
European Commission, Joint Research Centre
Created date
2022-06-01
Modified date
2022-06-24
Issued date
2022-06-25
Landing page
http://ghsl.jrc.ec.europa.eu/ 
Language(s)
English
Data theme(s)
Regions and cities, Science and technology
Update frequency
irregular
Identifier
http://data.europa.eu/89h/d07d81b4-7680-4d28-b896-583745c27085
Popularity