使用ICESat-2(Ice, Cloud, and land Elevation Satellite-2)的ALT06数据来计算冰川高度是一个常见的应用。在GEE上实现这个应用需要以下步骤:
- 打开GEE平台,载入所需数据
var icesat2 = ee.ImageCollection("NASA/ICEBridge/EXPLORER"); // 载入ICESat-2数据
- 选择数据源和区域
// 选择数据源和日期范围
var icesat2_filter = icesat2.filterDate('2018-10-14', '2018-10-20')
.filterBounds(geometry)
// 选择一个区域
var geometry = ee.Geometry.Point([83.7212, 28.3975]).buffer(5000); // 选择以该点为中心的半径为5000米的区域
- 定义函数来创建高度图像
// 创建高度图像
function createHeightMap(image) {
// 选取高度数据和可用返回波束(ATL06)的像素
var height = image.select('h_mean');
var quality = image.select('atl06_quality_summary').eq(0);
var mask = height.mask().and(quality);
return height.updateMask(mask);
}
- 编写代码来获取冰川高度
// 获取冰川高度
var heightMap = icesat2_filter.map(createHeightMap).mosaic();
- 可选步骤:应用除地球曲率外的其他因素来调整高度
// 使用外部因素调整高度
var iceThickness = heightMap.expression(
'(density * g * thickness + height * earth_radius) / earth_radius - height',
{density: 917, g: 9.81, thickness: 100}
);
可以使用以下链接来了解更多有关ICESat-2 ALT06数据和在GEE上计算冰川高度的信息:
- https://nsidc.org/data/icesat-2/data-sets
- https://developers.google.com/earth-engine/datasets/catalog/NASA_ICEBridge_ALT06_300m
- https://ieeexplore.ieee.org/abstract/document/9341058
- https://www.mdpi.com/2072-4292/11/22/2633/htm