Google Earth Engine（GEE）——全球去除森林和建筑物的哥白尼30米DEM数据集

var fabdem = ee.ImageCollection("projects/sat-io/open-datasets/FABDEM"),
    countries = ee.FeatureCollection("FAO/GAUL/2015/level0");

print('FABDEM Collection size :',fabdem.size())

//Explanation on setting default Projection here https://twitter.com/jstnbraaten/status/1494038930643042309
var elev = fabdem.mosaic().setDefaultProjection('EPSG:3857',null,30)


//you can also use this incase you don't want to specify CRS
//var elev = fabdem.mosaic().setDefaultProjection(glo30.first().projection())

// Add the elevation to the map.  Play with the visualization tools
// to get a better visualization.
Map.addLayer(elev, {}, 'elev', false);

// Use the terrain algorithms to compute a hillshade with 8-bit values.
var shade = ee.Terrain.hillshade(elev);
Map.addLayer(shade, {}, 'hillshade', false);

// Create an "ocean" variable to be used for cartographic purposes
var ocean = elev.lte(0);
Map.addLayer(ocean.mask(ocean), {palette:'000022'}, 'ocean', false);

// Create a custom elevation palette from hex strings.
var elevationPalette = ['006600', '002200', 'fff700', 'ab7634', 'c4d0ff', 'ffffff'];
// Use these visualization parameters, customized by location.
var visParams = {min: 1, max: 3000, palette: elevationPalette};

// Create a mosaic of the ocean and the elevation data
var visualized = ee.ImageCollection([
  // Mask the elevation to get only land
  elev.mask(ocean.not()).visualize(visParams), 
  // Use the ocean mask directly to display ocean.
  ocean.mask(ocean).visualize({palette:'000022'})
]).mosaic();

// Note that the visualization image doesn't require visualization parameters.
Map.addLayer(visualized.clip(countries), {}, 'elev palette');