GEE Python版/geemap 学习笔记 04 使用绘图工具与GEE数据交互

创建交互式地图

import os
import geemap
import ee
os.environ["HTTP_PROXY"] = "http://127.0.0.1:1234"
os.environ["HTTPS_PROXY"] = "http://127.0.0.1:1234"
Map = geemap.Map()
Map

 添加GEE数据集与底图

# Add Earth Engine dataset
dem = ee.Image('USGS/SRTMGL1_003')
landcover = ee.Image("ESA/GLOBCOVER_L4_200901_200912_V2_3").select('landcover')
landsat7 = ee.Image('LANDSAT/LE7_TOA_5YEAR/1999_2003').select(
    ['B1', 'B2', 'B3', 'B4', 'B5', 'B7']
)
states = ee.FeatureCollection("TIGER/2018/States")

# Set visualization parameters.
vis_params = {
    'min': 0,
    'max': 4000,
    'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5'],
}

# Add Earth Engine layers to Map
Map.addLayer(dem, vis_params, 'SRTM DEM', True, 0.5)
Map.addLayer(landcover, {}, 'Land cover')
Map.addLayer(
    landsat7,
    {'bands': ['B4', 'B3', 'B2'], 'min': 20, 'max': 200, 'gamma': 2.0},
    'Landsat 7',
) 
Map.addLayer(states, {}, "US States")
Map

 绘制形状

 检索刚刚绘制的所有形状的信息

Map.draw_features

检索最后一个形状的信息

Map.draw_last_feature

 

 过滤你所不需要的影像

roi = ee.FeatureCollection(Map.draw_features)
selected_states = states.filterBounds(roi)
Map.addLayer(selected_states, {}, "Selected_states")

 

 

clipped_image = dem.clip(selected_states)
Map.addLayer(clipped_image, vis_params, "Clipped_image")