GEE案例:Landsat 5、7、8影像构建1985-2023年rsei生态遥感指数详细代码

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是使用Google Earth Engine (GEE) 对Landsat-8多幅影像进行批量地形校正的代码示例: ```javascript // Load Landsat 8 surface reflectance data. var collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') .filterDate('2019-01-01', '2019-12-31') .filterBounds(geometry) .select(['B4', 'B3', 'B2', 'B5', 'B6', 'B7']); // Function to apply terrain correction to a single image. var applyTerrainCorrection = function(image) { // Load elevation data. var elevation = ee.Image('USGS/SRTMGL1_003'); // Compute terrain correction parameters. var tc = ee.Terrain.hillShadow(elevation, image.metadata('sunAzimuth'), image.metadata('sunElevation')); // Apply terrain correction to the image. var corrected = image.addBands(tc.select(['hillshade'])); return corrected; }; // Apply terrain correction to the image collection. var correctedCollection = collection.map(applyTerrainCorrection); // Export the corrected images to Google Drive. Export.image.toDrive({ image: correctedCollection, description: 'Landsat8_terrain_corrected', folder: 'GEE_exports', scale: 30, region: geometry }); ``` 在上述代码中,我们首先加载了Landsat 8表面反射率数据,并使用`filterDate()`和`filterBounds()`函数对数据进行筛选。然后,我们定义了一个名为`applyTerrainCorrection`的函数,该函数接受一个影像作为输入,并返回已进行地形校正的影像。在函数中,我们加载了SRTM高程数据,并使用`ee.Terrain.hillShadow()`函数计算了地形校正参数。最后,我们将地形校正应用到了影像中,并使用`addBands()`函数将地形校正参数作为新波段添加到了影像中。 接下来,我们使用`map()`函数将`applyTerrainCorrection`函数应用到整个影像集合中,并将结果存储在名为`correctedCollection`的变量中。最后,我们使用`Export.image.toDrive()`函数将地形校正后的影像导出到Google Drive中。 需要注意的是,上述代码中的`geometry`变量需要根据实际情况进行设置,以指定影像的空间范围。此外,如果需要,还可以调整`scale`参数以控制导出影像的空间分辨率。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

此星光明

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值