- 微信
- 微博
  
  分享文章到微博
- 复制链接
  
  复制链接到剪贴板

GEE数据集：wissSURFACE3D光栅（DSM）

此星光明发表于 2023/02/24 18:53:02 2023/02/24

【摘要】 swissSURFACE3D光栅（DSM）该数据集目前只对内部人员计划中的人开放。swissSURFACE3D Raster是一个数字表面模型（DSM），它代表了地球的表面，包括可见的和永久的景观元素，如土壤、自然覆盖物和各种建设性工作，但电力线和桅杆除外。 swissSURFACE3D Raster是由swissSURFACE3D的机载LiDAR数据得出的。为了对地表进行建模，使用了以下...

swissSURFACE3D光栅（DSM）
该数据集目前只对内部人员计划中的人开放。

swissSURFACE3D Raster是一个数字表面模型（DSM），它代表了地球的表面，包括可见的和永久的景观元素，如土壤、自然覆盖物和各种建设性工作，但电力线和桅杆除外。 swissSURFACE3D Raster是由swissSURFACE3D的机载LiDAR数据得出的。为了对地表进行建模，使用了以下元素（分类LiDAR点）。

土壤（第2类）
植被(第3类)
建筑物(第6类)
桥梁/人行桥/高架桥（17级）
水(第9类)
为了改善大型河流和湖泊的表面表现，实施了地形景观模型TLM的水道向量。

免责声明：数据集的全部或部分描述是由作者或其作品提供的。

数据预处理
数据集经过了预处理，并以多种格式分发。提供了云端优化的Geotiff文件，并编写了一个程序来获取所有的和当前的31321块瓷砖。

参考文献：

swissSURFACE3D Raster digital surface model (DSM). Last accessed **date** original data from https://www.swisstopo.admin.ch/en/geodata/height/surface3d-raster.html

代码：

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

print('Collection size',swiss3d.size())
Map.setCenter(8.54379, 47.36675,10)

swiss3d = swiss3d.mosaic().setDefaultProjection('EPSG:3857',null,0.5)


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

// Create an "ocean" variable to be used for cartographic purposes
var ocean = swiss3d.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
  swiss3d.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, {}, 'elev palette');

代码链接：https://code.earthengine.google.com/?scriptPath=users/sat-io/awesome-gee-catalog-examples:elevation-bathymetry/swiss3D-RASTER

License & Terms of Use¶

The free geodata and geoservices of swisstopo may be used, distributed and made accessible. Furthermore, they may be enriched and processed and also used commercially. A reference to the source is mandatory. In the case of digital or analogue representations and publications, as well as in the case of dissemination, one of the following source references must be attached in any case: