Mapping 10-m global impervious surface area (GISA-10m) using multi-source geospatial data
- 1School of Remote Sensing and Information Engineering, Wuhan University, Wuhan, 430079, P.R. China
- 2State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan, 430079, P.R. China
- 1School of Remote Sensing and Information Engineering, Wuhan University, Wuhan, 430079, P.R. China
- 2State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan, 430079, P.R. China
Abstract. Artificial impervious surface area (ISA) documents human footprints. Accurate, timely, and detailed ISA datasets are therefore essential for global climate change and urban planning. However, due to the lack of sufficient training samples and operational mapping methods, global ISA mapping at 10-m resolution is still lacking. To this end, we proposed a global ISA mapping method leveraging multi-source geospatial data. Based on the existing satellite-derived ISA maps and the crowdsourcing OpenStreetMap (OSM), 58 million training samples were extracted via a series of temporal, spatial, spectral, and geometric rules. Combined with over 2.7 million Sentinel optical and radar images on the Google Earth Engine, we produced the 10 m global ISA dataset (GISA-10m). Based on the test samples that are independent to the training set, GISA-10m embraced an overall accuracy greater than 86 %. In addition, the GISA-10m was comprehensively compared with the existing global ISA datasets, and the superiority of GISA-10m was demonstrated. It was found that China and the United States embraced the largest ISA and road area. The global rural ISA was 2.2 times that of urban while rural road area was 1.5 times larger than that of urban region. The global road area accounted for 14.2 % of the global ISA, 57.9 % of which was located in the top ten countries. Generally, the produced GISA-10m dataset and the proposed sampling and mapping method are able to achieve rapid and efficient global mapping, and have potential for detecting other land covers. It was also indicated that global ISA mapping can be improved by incorporating refined OSM data. GISA-10m can be used as a fundamental parameter for Earth system science, and provide valuable support for of urban planning and water cycle study. The GSIA-10m can be freely downloaded from http://doi.org/10.5281/zenodo.5791855 (Huang et al, 2021).
- Preprint
(4677 KB) -
Supplement
(705 KB) - BibTeX
- EndNote
Xin Huang et al.
Status: open (until 07 Jun 2022)
-
CC1: 'Comment on essd-2021-458', Chong Liu, 17 Apr 2022
reply
Accurate mapping of artifical impervious surface using remote sensing is challenging, especially at continental and global scales. The authors here provides an really exciting ISA map dataset which makes the above-mentioned challenge partially addressed. Given a 10 m spatial resolution, the GISA-10m is able to detect some subtle patterns that cannot be extracted by previously 30~300m products. As a public user, I only have one concern as follows.
The city group often include multiple cities with different scales. So, I believe it is important for potential users to know the accuracy of the GISA-10m for the cities with different scales, i.e., small, middle, and big city. Is it possible to compare the overall accuracy of the GISA-10m across different city sizes?
-
RC1: 'Comment on essd-2021-458', Anonymous Referee #1, 03 May 2022
reply
This manuscript proposed an efficient method to produce the 10-m global impervious surface areas (GISA-10m) based on the existing ISA maps, Sentinel-1/2 images, and OSM data. Compared to existing global GISA products, GISA-10m can provide higher spatial resolution while keeping higher accuracy. The inter-comparison with existing datasets demonstrated the superiority of GISA-10m. Analysis of ISA on rural and urban areas further revealed the urbanization level and landscape of different countries in more details. In particular, an interesting point of GISA-10m is that it is able to delineate the area of roads across the world, making GISA-10m valuable for relevant urban studies. In general, this manuscript is well presented and makes novel contributions. However, some issues should be clarified to improve this manuscript. Specific comments include the following aspects:
1) In Section 3.1.3, the authors used 200 trees for training the random forest classifier, while the effect of the number of trees is not analyzed. Besides, the key parameter, e.g., the number of features used for training each tree, is not clarified. Please provide this information for better understanding.
2) Line 90: do you mean by "operating by"?
3) Line 103: relevant reference should be provided to support “the terrain distortion caused by the combination of two orbits”.
4) L120: it should be "Landsat 8".
5) L121: I found both "GLCFCS" and "GLC_FCS" in the manuscript. Please explain.
6) Line 201: the original OSM data are provided in vector form. When this data was converted to 10-m raster, whether the majority rule was applied? The majority rule refers to “a pixel (10m × 10m) was labelled as ISA if more than half of its area was cover by ISA, otherwise it was identified as NISA”. Please clarified this issue.
7) L295: why the total number of visually interpreted samples was 10800 when 200 samples were selected in 59 grids? Please check.
8) Section 3.3: it is better to move this section to Section 5.1, since the detailed discussion has been presented in Section 5.1.
9) Figure 6: is it possible to compare the continental accuracies of other datasets presented in Table 3? The comparison at continental level may give a clear difference of different datasets.
10) Figure 8: not much information in it.
11) Figure 11: it is better to use (a), (b), (c) to distinguish each subgraph.
12) Figure 17: this figure is not clear enough for presenting 30 grids. It is suggested to add legend and put this figure to the supplementary materials.
13) Line 372: “extracted” or “detected”?
14) Table 9: whether test samples used in Table 9 are from visually interpreted samples? Please clarify this.
-
RC2: 'Comment on essd-2021-458', Xian Guo, 12 May 2022
reply
The mapping of 10-m impervious surfaces at the global scale using multiple geodata sources is interesting. The authors applied temporal-spatial-spectral-geometrical rules to generate samples, and validation of the results is comprehensive and adequate. They also attempted to delineate the spatial distribution of impervious surface in urban and non-urban areas. The manuscript fits the journal's scope and the dataset is valuable, which is suitable for publication in ESSD. However, the paper still has some flaws (see my comments below) which should be further clarified or discussed before acceptance.
My major concern lies in the completeness and correctness of the OSM data. How about the effect of the geographic bias in spatial distribution of OSM data? More analysis is needed to discuss this issue.
Line 10: “global ISA mapping” should be “global ISA datasets”
Line 21: “refined OSM data” -> “OSM data”.
Line 80-85: The GISA-10m dataset attempted to further delineate road regions from the ISA. This should be mentioned in the introduction and abstract.
Line 152: “multiple sources” is not clear, and can be modified as “multi-source datasets”.
Figure 1. It would be better to label each step, e.g., "Step 1. Training sample generation".
Line 157. The authors selected the GlobeLand30 in 2010 but chosed other data (e.g., GISA and FROM-GLC) in 2016. Would the temporal gap between these data impact the quality of training data?
Line 171. How did you define edge pixels? I think the edge pixels are different between 30-m and 10-m images, as a non-edge pixel in a 30m image may be edge pixels in a 10m image. Could you clarify this issue?
Line 197. Why buildings with area less than 100 m2 were excluded?
Line 210. Why did the authors remove the OSM samples intersected with those from other global datasets?
Line 235. Please explain why these features were chosen.
Line 286. How many RF models were built?
Line 293. How did the authors select the ISA test points? If the points were mostly located in urban areas, it might bias the assessment result. Could you provide the ISA density around these ISA points?
Figure 9. It's interesting to see the accuracy in rural and arid areas. How about urban areas?
Line 379. How did you divide the rural and urban areas?
Line 380. What do you mean by Global ISA?
Figure 14. The title of subgraph seems incorrect.
Figures 16 and 17 may be moved to the supplements.
Line 523. "difference" or " differences"
Line 500: “distinguish well ISA from NISA” -> “distinguish ISA from NISA effectively”
Xin Huang et al.
Data sets
Mapping 10-m global impervious surface area (GISA-10m) using multi-source geospatial data Xin Huang; Jie Yang; Wenrui Wang; Zhengrong Liu https://zenodo.org/record/5791855
Xin Huang et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
345 | 102 | 13 | 460 | 30 | 5 | 8 |
- HTML: 345
- PDF: 102
- XML: 13
- Total: 460
- Supplement: 30
- BibTeX: 5
- EndNote: 8
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1