Google Earth Engine（GEE）——欧空局WorldCover 10 m分辨率土地利用/土地覆盖数据集 2020 V100

此星光明

发布于 2024-02-02 02:18:59

42500

代码可运行

文章被收录于专栏：GEE数据专栏，GEE学习专栏，GEE错误集等专栏GEE数据专栏，GEE学习专栏，GEE错误集等专栏

运行总次数：0

代码可运行

欧空局WorldCover 10 m 2020 V100输入质量¶。 ESA WorldCover 10 m 2020 V100产品是以EPSG:4326投影（地理经纬度CRS）的云优化GeoTIFFs（COG）形式，以3x3度瓦片形式交付。有2651个瓦片，关于访问这个数据集的更多信息可以在这里找到。目前的收集只集中在输入质量层，地图层在谷歌地球引擎中可作为一个图像集。DATA | WORLDCOVER

输入质量层是一个显示地球观测（EO）输入数据质量的每个像素质量指标。该层是一个3波段的GeoTIFF，有

- 波段1：分类工作流程中使用的Sentinel-1 GAMMA0观测值的数量

- 波段2：分类工作流程中使用的哨兵-2 L2A观测值的数量

- 频段3 在分类工作流程中丢弃的无效S2观测值的百分比（0-100）（经过云和云影过滤）。

结合频段2和频段3，可以得出分类工作流程中使用的有效Sentinel-2 L2A观测值的总绝对数。

许可：

The ESA WorldCover product is provided free of charge, without restriction of use. For the full license information see the Creative Commons Attribution 4.0 International License.

Publications, models and data products that make use of these datasets must include proper acknowledgement, including citing the datasets and the journal article as in the following citation.

数据引用

Zanaga, D., Van De Kerchove, R., De Keersmaecker, W., Souverijns, N., Brockmann, C., Quast, R., Wevers, J., Grosu, A.,
Paccini, A., Vergnaud, S., Cartus, O., Santoro, M., Fritz, S., Georgieva, I., Lesiv, M., Carter, S., Herold, M., Li,
Linlin, Tsendbazar, N.E., Ramoino, F., Arino, O., 2021. ESA WorldCover 10 m 2020 v100.
https://doi.org/10.5281/zenodo.5571936

代码：

var iq = ee.ImageCollection("projects/sat-io/open-datasets/ESA_WorldCover_Input_Quality");

/*
ESA 10m World Cover Input Quality Dataset

Per pixel quality indicator showing the quality of the Earth Observation (EO) input data. The layer is a
3 band GeoTIFF with
• Band 1: Number of Sentinel-1 GAMMA0 observations used in the classification workflow
• Band 2: Number of Sentinel-2 L2A observations used in the classification workflow
• Band 3 Percentage (0-100) of invalid S2 observations discarded in the classification workflow
(after cloud and cloud shadow filtering).
Combining Band 2 and Band 3 yields the total absolute number of valid Sentinel-2 L2A observations
used in the classification workflow.
*/

//Count collection objects
print('ESA 10m Input Quality Collection size',iq.size())

//Import palette
var palettes = require('users/gena/packages:palettes')
var vis_b1 = {"bands":["b1"],"min":0,"max":174,palette: palettes.cmocean.Deep[7]};
var vis_b2 = {"bands":["b2"],"min":0,"max":174,palette: palettes.cmocean.Haline[7]};
var vis_b3 = {"bands":["b3"],"min":0,"max":100,palette: palettes.cmocean.Oxy[7]};
var vis_b23 = {"bands":["b23"],"min":0,"max":140,palette: palettes.cmocean.Algae[7]};


var b1_mosaic = iq.select(['b2']).mosaic()
var b3_mosaic = iq.select(['b3']).mosaic()
var combined_mosaic = b1_mosaic.multiply(ee.Image.constant(100).subtract(b3_mosaic)).divide(100).round().rename('b23')

Map.addLayer(iq.mosaic(),vis_b1,'Band 1 Input Quality: Number of Sentinel-1 GAMMA0 observations')
Map.addLayer(iq.mosaic(),vis_b2,'Band 2 Input Quality: Number of Sentinel-2 L2A observations',false)
Map.addLayer(iq.mosaic(),vis_b3,'Band 3 Input Quality: Percentage (0-100) of invalid S2 observations discarded')
Map.addLayer(combined_mosaic,vis_b23,'Band 2 combined with Band 3 Input Quality: total absolute number of valid Sentinel-2 L2A observations',false)