shp1.filterBounds(roi)(求roi在shp1内所涉及的矢量部分(shp内的n块))
裁剪:image1.clip(roi)(基于roi裁剪image1)
import ee
import geemap
geemap.set_proxy(port=10809)
# geemap.show_youtube('N7rK2aV1R4c')
Map = geemap.Map(center = (38 ,105))
Map
叠加地图
# add earth engine dataset
SRTM = ee.Image('USGS/SRTMGL1_003')
# Set visualization parameters
vis_param = {
'min': 0,
'max': 4000,
'palette': ['006633', 'E5FFCC', '662A00', 'D8D8D8', 'F5F5F5']
}
# add earth engine dem to map
Map.addLayer(SRTM, vis_param, 'SRTM_DEM')
Map
states = ee.FeatureCollection('TIGER/2018/States')
Map.addLayer(states, {}, 'US States')
Map
绘图
Map.draw_features
Map.draw_last_feature
roi = ee.FeatureCollection(Map.draw_features)
select_states = states.filterBounds(roi)
Map.addLayer(select_states, {}, 'Select States')
Map
clipped_image= SRTM.clip(select_states)
Map.addLayer(clipped_image, vis_param, 'Clipped image')
Map