GEE 夜间灯光数据裁剪教程：以无锡市为例

本教程将演示如何使用 Google Earth Engine (GEE) 和 Python 代码裁剪夜间灯光数据，以无锡市为例。

1. 准备工作

在开始之前，请确保您已完成以下准备工作：

• 已注册 Google Earth Engine 账号并成功登录。* 已安装 geemap Python 库。您可以使用 pip install geemap 命令进行安装。

2. 代码实现

以下是使用 GEE 裁剪无锡市夜间灯光数据的 Python 代码：pythonimport geemapimport ee

初始化 Earth Engineee.Initialize()

加载无锡市的边界数据wuxi = ee.FeatureCollection('users/wuxi_boundary')

加载 DMSP-OLS 夜间灯光数据集ntl = ee.ImageCollection('NOAA/DMSP-OLS/NIGHTTIME_LIGHTS')

定义裁剪函数def clip_image(image): return image.clip(wuxi)

裁剪夜间灯光数据集ntl_wuxi = ntl.map(clip_image)

可视化裁剪后的夜间灯光数据集vis_params = { 'min': 0, 'max': 63, 'palette': ['000000', 'FFFF00', 'FF0000']}

创建地图Map = geemap.Map()Map.centerObject(wuxi, 9)Map.addLayer(ntl_wuxi.median(), vis_params, '夜间灯光')Map.addLayer(wuxi, {}, '无锡市边界')Map.add_legend(legend_title='夜间灯光强度')Map.addLayerControl()Map

3. 代码说明

1. 首先，我们导入 geemap 和 ee 库。2. 然后，初始化 Earth Engine 并加载无锡市边界数据和 DMSP-OLS 夜间灯光数据集。3. 接着，我们定义一个名为 clip_image 的函数，该函数使用无锡市边界裁剪输入的影像。4. 然后，我们使用 map 函数将 clip_image 函数应用于整个夜间灯光数据集，得到裁剪后的数据集 ntl_wuxi。5. 最后，我们使用 geemap 库可视化裁剪后的夜间灯光数据，并在地图上添加无锡市边界、图例和图层控制器。

4. 总结

通过本教程，您学习了如何使用 GEE 和 Python 代码裁剪夜间灯光数据。您可以将此代码应用于其他地区或其他类型的地理空间数据。